Identification of loci associated with water use efficiency and symbiotic nitrogen fixation in soybean.

Soybean (Glycine max) production is greatly affected by persistent and/or intermittent droughts in rainfed soybean-growing regions worldwide. Symbiotic N2 fixation (SNF) in soybean can also be significantly hampered even under moderate drought stress. The objective of this study was to identify genomic regions associated with shoot carbon isotope ratio (δ13C) as a surrogate measure for water use efficiency (WUE), nitrogen isotope ratio (δ15N) to assess relative SNF, N concentration ([N]), and carbon/nitrogen ratio (C/N). Genome-wide association mapping was performed with 105 genotypes and approximately 4 million single-nucleotide polymorphism markers derived from whole-genome resequencing information. A total of 11, 21, 22, and 22 genomic loci associated with δ13C, δ15N, [N], and C/N, respectively, were identified in two environments. Nine of these 76 loci were stable across environments, as they were detected in both environments. In addition to the 62 novel loci identified, 14 loci aligned with previously reported quantitative trait loci for different C and N traits related to drought, WUE, and N2 fixation in soybean. A total of 58 Glyma gene models encoding for different genes related to the four traits were identified in the vicinity of the genomic loci.

Parallel tuning of semi-dwarfism via differential splicing of Brachytic1 in commercial maize and smallholder sorghum.

In the current genomic era, the search and deployment of new semi-dwarf alleles have continued to develop better plant types in all cereals. We characterized an agronomically optimal semi-dwarf mutation in Zea mays L. and a parallel polymorphism in Sorghum bicolor L. We cloned the maize brachytic1 (br1-Mu) allele by a modified PCR-based Sequence Amplified Insertion Flanking Fragment (SAIFF) approach. Histology and RNA-Seq elucidated the mechanism of semi-dwarfism. GWAS linked a sorghum plant height QTL with the Br1 homolog by resequencing a West African sorghum landraces panel. The semi-dwarf br1-Mu allele encodes an MYB transcription factor78 that positively regulates stalk cell elongation by interacting with the polar auxin pathway. Semi-dwarfism is due to differential splicing and low functional Br1 wild-type transcript expression. The sorghum ortholog, SbBr1, co-segregates with the major plant height QTL qHT7.1 and is alternatively spliced. The high frequency of the Sbbr1 allele in African landraces suggests that African smallholder farmers used the semi-dwarf allele to improve plant height in sorghum long before efforts to introduce Green Revolution-style varieties in the 1960s. Surprisingly, variants for differential splicing of Brachytic1 were found in both commercial maize and smallholder sorghum, suggesting parallel tuning of plant architecture across these systems.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches.

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina).

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics.

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

Genomic prediction of switchgrass winter survivorship across diverse lowland populations.

In the North-Central United States, lowland ecotype switchgrass can increase yield by up to 50% compared with locally adapted but early flowering cultivars. However, lowland ecotypes are not winter tolerant. The mechanism for winter damage is unknown but previously has been associated with late flowering time. This study investigated heading date (measured for two years) and winter survivorship (measured for three years) in a multi-generation population generated from two winter-hardy lowland individuals and diverse southern lowland populations. Sequencing data (311,776 markers) from 1,306 individuals were used to evaluate genome-wide trait prediction through cross-validation and progeny prediction (n = 52). Genetic variance for heading date and winter survivorship was additive with high narrow-sense heritability (0.64 and 0.71, respectively) and reliability (0.68 and 0.76, respectively). The initial negative correlation between winter survivorship and heading date degraded across generations (F1r = -0.43, pseudo-F2r = -0.28, pseudo-F2 progeny r = -0.15). Within-family predictive ability was moderately high for heading date and winter survivorship (0.53 and 0.52, respectively). A multi-trait model did not improve predictive ability for either trait. Progeny predictive ability was 0.71 for winter survivorship and 0.53 for heading date. These results suggest that lowland ecotype populations can obtain sufficient survival rates in the northern United States with two or three cycles of effective selection. Despite accurate genomic prediction, naturally occurring winter mortality successfully isolated winter tolerant genotypes and appears to be an efficient method to develop high-yielding, cold-tolerant switchgrass cultivars.

A Genome-Wide Association Study Identifying Single-Nucleotide Polymorphisms for Iron and Zinc Biofortification in a Worldwide Barley Collection.

Micronutrient deficiency affects half of the world's population, mostly in developing countries. Severe health issues such as anemia and inadequate growth in children below five years of age and pregnant women have been linked to mineral deficiencies (mostly zinc and iron). Improving the mineral content in staple crops, also known as mineral biofortification, remains the best approach to address mineral malnutrition. Barley is a staple crop in some parts of the world and is a healthy choice since it contains ß-glucan, a high dietary protein. Barley mineral biofortification, especially with zinc and iron, can be beneficial since barley easily adapts to marginalized areas and requires less input than other frequently consumed cereals. In this study, we analyzed zinc and iron content in 496 barley samples. The samples were genotyped with an Illumina 50 K SNP chip. Genome-wide association studies (GWAS) identified 62 SNPs and 68 SNPs (p < 0.001) associated with iron and zinc content in grains, respectively. After a Bonferroni correction (p < 0.005), there were 12 SNPs (single-nucleotide polymorphism) associated with Zn and 6 for iron. SNP annotations revealed proteins involved in membrane transport, Zn and Fe binding, linked to nutrient remobilization in grains. These results can be used to develop biofortified barley via marker-assisted selection (MAS), which could alleviate mineral malnutrition.

A generalist-specialist trade-off between switchgrass cytotypes impacts climate adaptation and geographic range.

Polyploidy results from whole-genome duplication and is a unique form of heritable variation with pronounced evolutionary implications. Different ploidy levels, or cytotypes, can exist within a single species, and such systems provide an opportunity to assess how ploidy variation alters phenotypic novelty, adaptability, and fitness, which can, in turn, drive the development of unique ecological niches that promote the coexistence of multiple cytotypes. Switchgrass, Panicum virgatum, is a widespread, perennial C4 grass in North America with multiple naturally occurring cytotypes, primarily tetraploids (4×) and octoploids (8×). Using a combination of genomic, quantitative genetic, landscape, and niche modeling approaches, we detect divergent levels of genetic admixture, evidence of niche differentiation, and differential environmental sensitivity between switchgrass cytotypes. Taken together, these findings support a generalist (8×)­specialist (4×) trade-off. Our results indicate that the 8× represent a unique combination of genetic variation that has allowed the expansion of switchgrass' ecological niche and thus putatively represents a valuable breeding resource.

The recent evolutionary rescue of a staple crop depended on over half a century of global germplasm exchange.

Rapid environmental change can lead to population extinction or evolutionary rescue. The global staple crop sorghum (Sorghum bicolor) has recently been threatened by a global outbreak of an aggressive new biotype of sugarcane aphid (SCA; Melanaphis sacchari). We characterized genomic signatures of adaptation in a Haitian breeding population that had rapidly adapted to SCA infestation, conducting evolutionary population genomics analyses on 296 Haitian lines versus 767 global accessions. Genome scans and geographic analyses suggest that SCA adaptation has been conferred by a globally rare East African allele of RMES1, which spread to breeding programs in Africa, Asia, and the Americas. De novo genome sequencing revealed potential causative variants at RMES1. Markers developed from the RMES1 sweep predicted resistance in eight independent commercial and public breeding programs. These findings demonstrate the value of evolutionary genomics to develop adaptive trait technology and highlight the benefits of global germplasm exchange to facilitate evolutionary rescue.

Four chromosome scale genomes and a pan-genome annotation to accelerate pecan tree breeding.

Genome-enabled biotechnologies have the potential to accelerate breeding efforts in long-lived perennial crop species. Despite the transformative potential of molecular tools in pecan and other outcrossing tree species, highly heterozygous genomes, significant presence-absence gene content variation, and histories of interspecific hybridization have constrained breeding efforts. To overcome these challenges, here, we present diploid genome assemblies and annotations of four outbred pecan genotypes, including a PacBio HiFi chromosome-scale assembly of both haplotypes of the 'Pawnee' cultivar. Comparative analysis and pan-genome integration reveal substantial and likely adaptive interspecific genomic introgressions, including an over-retained haplotype introgressed from bitternut hickory into pecan breeding pedigrees. Further, by leveraging our pan-genome presence-absence and functional annotation database among genomes and within the two outbred haplotypes of the 'Lakota' genome, we identify candidate genes for pest and pathogen resistance. Combined, these analyses and resources highlight significant progress towards functional and quantitative genomics in highly diverse and outbred crops.

Genetic dissection of natural variation in oilseed traits of camelina by whole-genome resequencing and QTL mapping.

Camelina [Camelina sativa (L.) Crantz] is an oilseed crop in the Brassicaceae family that is currently being developed as a source of bioenergy and healthy fatty acids. To facilitate modern breeding efforts through marker-assisted selection and biotechnology, we evaluated genetic variation among a worldwide collection of 222 camelina accessions. We performed whole-genome resequencing to obtain single nucleotide polymorphism (SNP) markers and to analyze genomic diversity. We also conducted phenotypic field evaluations in two consecutive seasons for variations in key agronomic traits related to oilseed production such as seed size, oil content (OC), fatty acid composition, and flowering time. We determined the population structure of the camelina accessions using 161,301 SNPs. Further, we identified quantitative trait loci (QTL) and candidate genes controlling the above field-evaluated traits by genome-wide association studies (GWAS) complemented with linkage mapping using a recombinant inbred line (RIL) population. Characterization of the natural variation at the genome and phenotypic levels provides valuable resources to camelina genetic studies and crop improvement. The QTL and candidate genes should assist in breeding of advanced camelina varieties that can be integrated into the cropping systems for the production of high yield of oils of desired fatty acid composition.

The tepary bean genome provides insight into evolution and domestication under heat stress.

Tepary bean (Phaseolus acutifolis A. Gray), native to the Sonoran Desert, is highly adapted to heat and drought. It is a sister species of common bean (Phaseolus vulgaris L.), the most important legume protein source for direct human consumption, and whose production is threatened by climate change. Here, we report on the tepary genome including exploration of possible mechanisms for resilience to moderate heat stress and a reduced disease resistance gene repertoire, consistent with adaptation to arid and hot environments. Extensive collinearity and shared gene content among these Phaseolus species will facilitate engineering climate adaptation in common bean, a key food security crop, and accelerate tepary bean improvement.

Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

A genome resource for green millet Setaria viridis enables discovery of agronomically valuable loci.

Wild and weedy relatives of domesticated crops harbor genetic variants that can advance agricultural biotechnology. Here we provide a genome resource for the wild plant green millet (Setaria viridis), a model species for studies of C4 grasses, and use the resource to probe domestication genes in the close crop relative foxtail millet (Setaria italica). We produced a platinum-quality genome assembly of S. viridis and de novo assemblies for 598 wild accessions and exploited these assemblies to identify loci underlying three traits: response to climate, a 'loss of shattering' trait that permits mechanical harvest and leaf angle, a predictor of yield in many grass crops. With CRISPR-Cas9 genome editing, we validated Less Shattering1 (SvLes1) as a gene whose product controls seed shattering. In S. italica, this gene was rendered nonfunctional by a retrotransposon insertion in the domesticated loss-of-shattering allele SiLes1-TE (transposable element). This resource will enhance the utility of S. viridis for dissection of complex traits and biotechnological improvement of panicoid crops.

A genome assembly and the somatic genetic and epigenetic mutation rate in a wild long-lived perennial Populus trichocarpa.

BACKGROUND: Plants can transmit somatic mutations and epimutations to offspring, which in turn can affect fitness. Knowledge of the rate at which these variations arise is necessary to understand how plant development contributes to local adaption in an ecoevolutionary context, particularly in long-lived perennials. RESULTS: Here, we generate a new high-quality reference genome from the oldest branch of a wild Populus trichocarpa tree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we use a multi-omics approach to quantify age-related somatic changes at the genetic, epigenetic, and transcriptional level. We show that the per-year somatic mutation and epimutation rates are lower than in annuals and that transcriptional variation is mainly independent of age divergence and cytosine methylation. Furthermore, a detailed analysis of the somatic epimutation spectrum indicates that transgenerationally heritable epimutations originate mainly from DNA methylation maintenance errors during mitotic rather than during meiotic cell divisions. CONCLUSION: Taken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

Genome-wide analysis of the transcriptional response to drought stress in root and leaf of common bean.

Genes related to the response to drought stress in leaf and root tissue of drought-susceptible (DS) and tolerant (DT) genotypes were characterized by RNA-Seq. In total, 54,750 transcripts, representative of 28,590 genes, were identified; of these, 1,648 were of high-fidelity (merge of 12 libraries) and described for the first time in the Andean germplasm. From the 1,239 differentially expressed genes (DEGs), 458 were identified in DT, with a predominance of genes in categories of oxidative stress, response to stimulus and kinase activity. Most genes related to oxidation-reduction terms in roots were early triggered in DT (T75) compared to DS (T150) suggestive of a mechanism of tolerance by reducing the damage from ROS. Among the KEGG enriched by DEGs up-regulated in DT leaves, two related to the formation of Sulfur-containing compounds, which are known for their involvement in tolerance to abiotic stresses, were common to all treatments. Through qPCR, 88.64% of the DEGs were validated. A total of 151,283 variants were identified and functional effects estimated for 85,780. The raw data files were submitted to the NCBI database. A transcriptome map revealed new genes and isoforms under drought. These results supports a better understanding of the drought tolerance mechanisms in beans.

The genomic landscape of molecular responses to natural drought stress in Panicum hallii.

Environmental stress is a major driver of ecological community dynamics and agricultural productivity. This is especially true for soil water availability, because drought is the greatest abiotic inhibitor of worldwide crop yields. Here, we test the genetic basis of drought responses in the genetic model for C4 perennial grasses, Panicum hallii, through population genomics, field-scale gene-expression (eQTL) analysis, and comparison of two complete genomes. While gene expression networks are dominated by local cis-regulatory elements, we observe three genomic hotspots of unlinked trans-regulatory loci. These regulatory hubs are four times more drought responsive than the genome-wide average. Additionally, cis- and trans-regulatory networks are more likely to have opposing effects than expected under neutral evolution, supporting a strong influence of compensatory evolution and stabilizing selection. These results implicate trans-regulatory evolution as a driver of drought responses and demonstrate the potential for crop improvement in drought-prone regions through modification of gene regulatory networks.

Population genomics and climate adaptation of a C4 perennial grass, Panicum hallii (Poaceae).

BACKGROUND: Understanding how and why genetic variation is partitioned across geographic space is of fundamental importance to understanding the nature of biological species. How geographical isolation and local adaptation contribute to the formation of ecotypically differentiated groups of plants is just beginning to be understood through population genomic studies. We used whole genome sequencing combined with association study of climate to discover the drivers of differentiation in the perennial C4 grass Panicum hallii. RESULTS: Sequencing of 89 natural accessions of P.hallii revealed complex population structure across the species range. Major population genomic separation was found between subspecies P.hallii var. hallii and var. filipes as well as between at least four major unrecognized subgroups within var. hallii. At least 139 genomic SNPs were significantly associated with temperature or precipitation across the range and these SNPs were enriched for non-synonymous substitutions. SNPs associated with temperature and aridity were more often found in or near genes than expected by chance and enriched for putative involvement in dormancy processes, seed maturation, response to hyperosmosis and salinity, abscisic acid metabolism, hormone metabolism, and drought recovery. CONCLUSIONS: Both geography and climate adaptation contribute significantly to patterns of genome-wide variation in P.hallii. Population subgroups within P.hallii may represent early stages in the formation of ecotypes. Climate associated loci identified here represent promising targets for future research in this and other perennial grasses.

White seed color in common bean (Phaseolus vulgaris) results from convergent evolution in the P (pigment) gene.

The presence of seed color in common bean (Phaseolus vulgaris) requires the dominant-acting P (pigment) gene, and white seed is a recessive phenotype in all domesticated races of the species. P was classically associated with seed size, thus describing it as the first genetic marker for a quantitative trait. The molecular structure of P was characterized to understand the selection of white seeds during bean diversification and the relationship of P to seed weight. P was identified by homology searches, a genome-wide association study (GWAS) and gene remodeling, and confirmed by gene silencing. Allelic variation was assessed by a combination of resequencing and marker development, and the relationship between P and seed weight was assessed by a GWAS study. P is a member of clade B of subclass IIIf of plant basic helix-loop-helix (bHLH) proteins. Ten race-specific P alleles conditioned the white seed phenotype, and each causative mutation affected at least one bHLH domain required for color expression. GWAS analysis confirmed the classic association of P with seed weight. In common bean, white seeds are the result of convergent evolution and, among plant species, orthologous convergence on a single transcription factor gene was observed.