Search NASA⌕ Search

SEARCH · Search NASA

Results for “genome structure and evolution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Genomic analysis of Klebsiella aerogenes circulating in New Mexico

Klebsiella aerogenes is an opportunistic pathogen and a growing cause of healthcare-associated infections, characterized by multidrug resistance and the emergence of global high-risk clones. However, regional genomic surveillance data remain limited. Here, we sought to characterize the population structure, transmission dynamics and resistance mechanisms of clinical K. aerogenes in Albuquerque, New Mexico. We sequenced 177 clinical isolates collected between 2021 and 2023. We also developed a novel, species-specific PopPUNK database to facilitate rapid, high-resolution typing. The New Mexico K. aerogenes population was diverse but dominated by two global pandemic lineages, ST93 (47.5%) and ST4 (7.9%), which were significantly enriched for the virulence factors yersiniabactin and colibactin. Genomic evidence for recent local transmission was rare, with only four putative transmission pairs identified. The resistome was characterized by intrinsic and adaptive mutations. Nearly all isolates possessed gyrA mutations associated with decreased fluoroquinolone susceptibility. Mutations in the AmpC regulator AmpD and the outer membrane porin Omp36 were common, particularly within the dominant ST93 lineage. These mutations have been associated with increased AmpC-mediated carbapenem resistance. Our findings underscore the critical importance of genomic surveillance to monitor the transmission and evolution of adaptive resistance.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning enables identification of an alternative yeast galactose utilization pathway

How genomic differences contribute to phenotypic differences is a major question in biology. The recently characterized genomes, isolation environments, and qualitative patterns of growth on 122 sources and conditions of 1,154 strains from 1,049 fungal species (nearly all known) in the yeast subphylum Saccharomycotina provide a powerful, yet complex, dataset for addressing this question. We used a random forest algorithm trained on these genomic, metabolic, and environmental data to predict growth on several carbon sources with high accuracy. Known structural genes involved in assimilation of these sources and presence/absence patterns of growth in other sources were important features contributing to prediction accuracy. By further examining growth on galactose, we found that it can be predicted with high accuracy from either genomic (92.2%) or growth data (82.6%) but not from isolation environment data (65.6%). Prediction accuracy was even higher (93.3%) when we combined genomic and growth data. After the GALactose utilization genes, the most important feature for predicting growth on galactose was growth on galactitol, raising the hypothesis that several species in two orders, Serinales and Pichiales (containing the emerging pathogen Candida auris and the genus Ogataea, respectively), have an alternative galactose utilization pathway because they lack the GAL genes. Growth and biochemical assays confirmed that several of these species utilize galactose through an alternative oxidoreductive D-galactose pathway, rather than the canonical GAL pathway. Machine learning approaches are powerful for investigating the evolution of the yeast genotype–phenotype map, and their application will uncover novel biology, even in well-studied traits.

59 BASIC BIOLOGICAL SCIENCES↗

Probing interspecies metabolic interactions within a synthetic binary microbiome using genome-scale modeling

Metabolic interactions within a microbial community play a key role in determining the structure, function, and composition of the community. However, due to the complexity and intractability of natural microbiomes, limited knowledge is available on interspecies interactions within a community. In this work, using a binary synthetic microbiome, a methanotroph-photoautotroph (M-P) coculture, as the model system, we examined different genome-scale metabolic modeling (GEM) approaches to gain a better understanding of the metabolic interactions within the coculture, how they contribute to the enhanced growth observed in the coculture, and how they evolve over time. Using batch growth data of the model M-P coculture, we compared three GEM approaches for microbial communities. Two of the methods are existing approaches: SteadyCom, a steady state GEM, and dynamic flux balance analysis (DFBA) Lab, a dynamic GEM. We also proposed an improved dynamic GEM approach, DynamiCom, for the M-P coculture. SteadyCom can predict the metabolic interactions within the coculture but not their dynamic evolutions; DFBA Lab can predict the dynamics of the coculture but cannot identify interspecies interactions. DynamiCom was able to identify the cross-fed metabolite within the coculture, as well as predict the evolution of the interspecies interactions over time. A new dynamic GEM approach, DynamiCom, was developed for a model M-P coculture. Constrained by the predictions from a validated kinetic model, DynamiCom consistently predicted the top metabolites being exchanged in the M-P coculture, as well as the establishment of the mutualistic N-exchange between the methanotroph and cyanobacteria. The interspecies interactions and their dynamic evolution predicted by DynamiCom are supported by ample evidence in the literature on methanotroph, cyanobacteria, and other cyanobacteria-heterotroph cocultures.

59 BASIC BIOLOGICAL SCIENCES↗

Computer Modeling of Protocellular Functions: Peptide Insertion in Membranes

Lipid vesicles became the precursors to protocells by acquiring the capabilities needed to survive and reproduce. These include transport of ions, nutrients and waste products across cell walls and capture of energy and its conversion into a chemically usable form. In modem organisms these functions are carried out by membrane-bound proteins (about 30% of the genome codes for this kind of proteins). A number of properties of alpha-helical peptides suggest that their associations are excellent candidates for protobiological precursors of proteins. In particular, some simple a-helical peptides can aggregate spontaneously and form functional channels. This process can be described conceptually by a three-step thermodynamic cycle: 1 - folding of helices at the water-membrane interface, 2 - helix insertion into the lipid bilayer and 3 - specific interactions of these helices that result in functional tertiary structures. Although a crucial step, helix insertion has not been adequately studied because of the insolubility and aggregation of hydrophobic peptides. In this work, we use computer simulation methods (Molecular Dynamics) to characterize the energetics of helix insertion and we discuss its importance in an evolutionary context. Specifically, helices could self-assemble only if their interactions were sufficiently strong to compensate the unfavorable Free Energy of insertion of individual helices into membranes, providing a selection mechanism for protobiological evolution.

Rodriquez-Gomez, D.↗

Discovery of an ultraspecific triuret hydrolase (TrtA) establishes the triuret biodegradation pathway

Triuret (carbonyldiurea) is an impurity found in industrial urea fertilizer (<0.1% w/w) that is applied, worldwide, around 300 million pounds each year on agricultural lands. In addition to anthropogenic sources, endogenous triuret has been identified in amoeba and human urine, the latter being diagnostic for hypokalemia. The present study is the first to describe the metabolic breakdown of triuret, which funnels into biuret metabolism. We identified the gene responsible for triuret decomposition ( trtA ) in bacterial genomes, clustered with biuH , which encodes biuret hydrolase and has close protein sequence homology. TrtA is a member of the isochorismatase-like hydrolase (IHL) protein family, similarly to BiuH, and has a catalytic efficiency (k cat/ K M ) of 6 x 10 5 M –1 s –1 , a K M for triuret of 20 μM, and exquisite substrate specificity. Indeed, TrtA has four orders of magnitude less activity with biuret. Crystal structures of TrtA in apo and holo form were solved and compared with the BiuH structure. The high substrate selectivity was found to be conveyed by second shell residues around each active site. Mutagenesis of residues conserved in TrtA to the alternate consensus found in BiuHs revealed residues critical to triuret hydrolase activity but no single mutant evolved more biuret activity, and likely a combination of mutations is required to interconvert between TrtA, BiuH functions. TrtA-mediated triuret metabolism is relatively rare in recorded genomes (1–2%), but is largely found in plant-associated, nodulating, and endophytic bacteria. This study suggests functions for triuret hydrolase in certain eukaryotic intermediary processes and prokaryotic intermediary or biodegradative metabolism.

36 MATERIALS SCIENCE↗

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

35S rDNA evolution↗

Monolignol radical-radical coupling networks in western red cedar and Arabidopsis and their evolutionary implications

The discovery of a nine-member multigene dirigent family involved in control of monolignol radical-radical coupling in the ancient gymnosperm, western red cedar, suggested that a complex multidimensional network had evolved to regulate such processes in vascular plants. Accordingly, in this study, the corresponding promoter regions for each dirigent multigene member were obtained by genome-walking, with Arabidopsis being subsequently transformed to express each promoter fused to the beta-glucuronidase (GUS) reporter gene. It was found that each component gene of the proposed network is apparently differentially expressed in individual tissues, organs and cells at all stages of plant growth and development. The data so obtained thus further support the hypothesis that a sophisticated monolignol radical-radical coupling network exists in plants which has been highly conserved throughout vascular plant evolution.

NASA Discipline Plant Biology↗

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta)

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. Here, we report that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum . The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f 4 ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.

59 BASIC BIOLOGICAL SCIENCES↗

Population Genomics of Pseudocercospora griseola Reveals New Groups in the Middle American Clade and the Presence of the Endophytic Bacterium Achromobacter xylosoxidans

Angular leaf spot (ALS), caused by Pseudocercospora griseola is an important disease of common beans. P. griseola, is highly variable and has co-evolved with its host. In this study, 48 isolates of P. griseola from Puerto Rico, Guatemala, Honduras and Tanzania were sequenced (3RADseq), resulting in the de novo assembly of 42,214 contigs. Phylogenomic, population genetic structure and principal component analyses using 1,260 SNPs divided these isolates into two populations, Andean and Middle American, while the Middle American population was further divided into three sub-populations. There were moderate to high levels of differentiation between P. griseola populations, with pairwise Fst values ranging from 0.11 to 0.95. The Andean population was composed of isolates from Tanzania, and was separated from the Middle American population (Fst = 0.95). The Middle American population was separated into 3 subpopulations including isolates from: 1. Guatemala and Honduras, 2. Tanzania, and 3. Puerto Rico. Pathogenicity testing of 27 isolates from Puerto Rico, using 12 common bean differential lines, identified ten races, but these races were not associated with SNPs found in virulence genes. DNA of an endophytic bacterium (Achromobacter xylosoxidans) was found in seven mildly virulent isolates suggesting a possible role of the bacterium in the observed virulence patterns. To understand the evolution and diversity of P. griseola, further study of the virulence genes and the interactions among the endophytic bacterium, the fungus, and the host plant is required. Such information is critical to inform breeding strategies for the development of resistant germplasm and cultivars.

Serrato-Diaz, Luz M. [U.S. Department of Agricultu↗

Genomic Features of Antimicrobial Resistance in Staphylococcus pseudintermedius Isolated from Dogs with Pyoderma in Argentina and the United States: A Comparative Study

Staphylococcus pseudintermedius is the most common opportunistic pathogen in dogs and methicillin resistance (MRSP) has been identified as an emerging problem in canine pyoderma. Here, we evaluated the antimicrobial resistance (AMR) features and phylogeny of S. pseudintermedius isolated from canine pyoderma cases in Argentina (n = 29) and the United States (n = 29). 62% of isolates showed multi-drug resistance. The AMR genes found: mecA, blaZ, ermB, dfrG, catA, tetM, aac(6')-aph(2"), in addition to tetK and lnuA (only found in U.S. isolates). Two point mutations were detected: grlA(S80I)-gyrA(S84L), and grlA(D84N)-gyrA(S84L) in one U.S. isolate. A mutation in rpoB (H481N) was found in two isolates from Argentina. SCCmec type III, SCCmec type V, ΨSCCmec 57395 were identified in the Argentinian isolates; and SCCmec type III, SCCmec type IVg, SCCmec type V, and SCCmec type VII variant in the U.S. cohort. Sequence type (ST) ST71 belonging to a dominant clone was found in isolates from both countries, and ST45 only in Argentinian isolates. This is the first study to comparatively analyze the population structure of canine pyoderma-associated S. pseudintermedius isolates in Argentina and in the U.S. It is important to maintain surveillance on S. pseudintermedius populations to monitor AMR and gain further understanding of its evolution and dissemination.

60 APPLIED LIFE SCIENCES↗

Small Ion Channel Linking Molecular Simulations and Electrophysiology

Ion channels are pore-forming protein assemblies that mediate the transport of small ions across cell membranes. Otherwise, membrane bilayers would be almost impermeable to ions incapable to traverse the low dielectric constant, hydrophobic membrane core. Ion channels are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses, cardiac functions, muscle contraction and apoptosis. On the other extreme of biological complexity, viral ion channels (viroporins) influence many stages of the virus infection cycle either through regulating virus replication, such as entry, assembly and release or modulating the electrochemical balance in the subcellular compartments of host cells. Ion channels were crucial components of protocells. Their emergence facilitated adaptation of nascent life to different environmental conditions. The earliest ion channels must have been much simpler than most of their modern ancestors. Viral channels are among only a few naturally occurring models to study the structure, function and evolution of primordial channels. Experimental studies of these properties are difficult and often unreliable. In principle, computational methods, and molecular dynamics (MD) simulations in particular, can aid in providing information about both the structure and the function of ion channels. However, MD suffers from its own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. It is, therefore, essential to determine the reliability of MD simulations. We propose to do so on the basis of two criteria. One is channel stability on time scales that extend for several microseconds or longer. The other is the ability to reproduce the measured ionic conductance as a function of applied voltage. If both the stability and the calculated ionic conductance are satisfactory it will greatly increase our confidence that the structure and the function of a channel are described sufficiently accurately. To our knowledge, long time scale stability (approx.10 micro-sec) and the correct electrophysiology have been shown so far for only one channel - the synthetic LS3 hexamer). In this presentation, this approach will be discussed in application to two viral channels - Vpu, encoded by the HIV-1 genome and p7 of hepatitis C.

Pohorille, Andrzej↗

Re-evaluating the Systematics of Dendrolycopodium Using Restriction-Site Associated DNA-Sequencing

The genus Dendrolycopodium (Lycopodiaceae) includes four to five species across North America and East Asia. Species identification in Dendrolycopodium is difficult due to limited or inconsistent characters. In addition, plants with intermediate morphologies regularly occur, potentially indicative of interspecific hybridization. To determine the species relationships in Dendrolycopodium and investigate the existence of hybrids, we generated a draft genome assembly for D. obscurum and carried out double-digest restriction-site associated DNA sequencing (RADSeq) on 86 Dendrolycopodium specimens. Our sampling includes all the described species and 11 individuals with intermediate morphology. We find that the genus can be divided into four clades that largely correspond to the described taxa, as well as evidence of interspecific hybridization. Within these clades, our STRUCTURE analysis suggests that there are multiple finer subgroups, with evidence of hybridization and introgression between these subgroups. Given the limited availability of specimens collected from Asia, the status of the various Asian species remains uncertain and will require further study. In summary, our study confirms several hybrid relationships in Dendrolycopodium and provides a clear phylogenetic framework for future taxonomic revision.

59 BASIC BIOLOGICAL SCIENCES↗

Nasa GeneLab Computomics Reveal Horizontal Gene Transfer on International Space Station Environmental Metagenomes

Prokaryotic lifeforms can be observed to demonstrate many keen adaptive advantages, perhaps facilitated by a nature simplistic relative to divergent domains of life. In particular, decompartmentalized gene expression facilitates adaptation by allowing free exchange of genetic material, albeit at the cost of increased susceptibility to genetic damage. Thus, these lifeforms must compensate by embracing diverse investment strategies in an attempt to “brute force” the evolvability equation through precipitous genesis, lean metabolic efficiency, and sheer population. This prokaryotic archetype also enables symbiotic relationships with secondary mobile genetic elements known as plasmids, which have been shown to drive evolution on rapid temporal scales through processes such as conjugation and transformation. This study attempts to decipher whether these mechanisms of horizontal gene transfer (HGT) are major factors in determining prokaryote fitness within a unique isolated environment, the International Space Station (ISS). The ISS Microbial Tracking (MT) project has generated a wealth of data concerning the successive reigns of microbial genera that appear to thrive amidst harsh conditions for life. Despite relatively higher doses of ionizing radiation as compared to Earth, complications associated with microgravity, and the anti-microbial mélange deployed, microbial life still persists in this environment. The NASA GeneLab serves as a data repository and analysis platform to enable researchers to access space flight factor related data. With the use of GeneLab’s modern computational suites (computomics), phylogenetic and functional genomic investigations of HGT events were conducted on the data generated from the MT-1 project. The putative data concerning the plasmid population (plasmidome) of the ISS was algorithmically derived and compared to those of habitats with similar environmental dynamics- such as living quarters and hospitals- to investigate whether these HGT elements may play crucial role(s) in shaping the microbiome of this closed habitat that serves as the only inhabited structure in space.

Bense, Nicholas↗

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

60 APPLIED LIFE SCIENCES↗

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium

The “genomic shock” hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession “Bhyb26.” We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species’ similarity in transposable element load may account for the subtlety of the observed genome dominance.

59 BASIC BIOLOGICAL SCIENCES↗

Expanding the genomic encyclopedia of Actinobacteria with 824 isolate reference genomes

The phylum Actinobacteria includes important human pathogens like Mycobacterium tuberculosis and Corynebacterium diphtheriae and renowned producers of secondary metabolites of commercial interest, yet only a small part of its diversity is represented by sequenced genomes. Here, we present 824 actinobacterial isolate genomes in the context of a phylum-wide analysis of 6,700 genomes including public isolates and metagenome-assembled genomes (MAGs). We estimate that only 30%–50% of projected actinobacterial phylogenetic diversity possesses genomic representation via isolates and MAGs. A comparison of gene functions reveals novel determinants of host-microbe interaction as well as environment-specific adaptations such as potential antimicrobial peptides. We identify plasmids and prophages across isolates and uncover extensive prophage diversity structured mainly by host taxonomy. Analysis of >80,000 biosynthetic gene clusters reveals that horizontal gene transfer and gene loss shape secondary metabolite repertoire across taxa. Our observations illustrate the essential role of and need for high-quality isolate genome sequences.

59 BASIC BIOLOGICAL SCIENCES↗

Community‐Level Metabolic Shifts Following Land Use Change in the Amazon Rainforest Identified by a Supervised Machine Leaning Approach

ABSTRACT The Amazon rainforest has been subjected to high rates of deforestation, mostly for pasturelands, over the last few decades. This change in plant cover is known to alter the soil microbiome and the functions it mediates, but the genomic changes underlying this response are still unresolved. In this study, we used a combination of deep shotgun metagenomics complemented by a supervised machine learning approach to compare the metabolic strategies of tropical soil microbial communities in pristine forests and long‐term established pastures in the Amazon. Machine learning‐derived metagenome analysis indicated that microbial community structures (bacteria, archaea and viruses) and the composition of protein‐coding genes were distinct in each plant cover type environment. Forest and pasture soils had different genomic diversities for the above three taxonomic groups, characterised by their protein‐coding genes. These differences in metagenome profiles in soils under forests and pastures suggest that metabolic strategies related to carbohydrate and energy metabolisms were altered at community level. Changes were also consistent with known modifications to the C and N cycles caused by long‐term shifts in aboveground vegetation and were also associated with several soil physicochemical properties known to change with land use, such as the C/N ratio, soil temperature and exchangeable acidity. In addition, our analysis reveals that these alterations in land use can also result in changes to the composition and diversity of the soil DNA virome. Collectively, our study indicates that soil microbial communities shift their overall metabolic strategies, driven by genomic alterations observed in pristine forests and long‐term established pastures with implications for the C and N cycles.

carbon and nitrogen cycles↗

Genome analyses reveal population structure and a purple stigma color gene candidate in finger millet

Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a ‘poor man’s crop’, finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21 st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

59 BASIC BIOLOGICAL SCIENCES↗