Search

Polyploidy is an evolutionary innovation for many animals and all flowering plants, but its impact on selection and domestication remains elusive. Here we analyze genome evolution and diversification for all five allopolyploid cotton species, including economically important Upland and Pima cottons. Although these polyploid genomes are conserved in gene content and synteny, they have diversified by subgenomic transposon exchanges that equilibrate genome size, evolutionary rate heterogeneities and positive selection between homoeologs within and among lineages. These differential evolutionary trajectories are accompanied by gene-family diversification and homoeolog expression divergence among polyploid lineages. Selection and domestication drive parallel gene expression similarities in fibers of two cultivated cottons, involving coexpression networks and N6-methyladenosine RNA modifications. Furthermore, polyploidy induces recombination suppression, which correlates with altered epigenetic landscapes and can be overcome by wild introgression. These genomic insights will empower efforts to manipulate genetic recombination and modify epigenetic landscapes and target genes for crop improvement., Sequencing and genomic diversification of five allopolyploid cotton species provide insights into polyploid genome evolution and epigenetic landscapes for cotton improvement.

Sequence assembly of large and repeat-rich plant genomes has been challenging, requiring substantial computational resources and often several complementary sequence assembly and genome mapping approaches. The recent development of fast and accurate long-read sequencing by circular consensus sequencing (CCS) on the PacBio platform may greatly increase the scope of plant pan-genome projects. Here, we compare current long-read sequencing platforms regarding their ability to rapidly generate contiguous sequence assemblies in pan-genome studies of barley (Hordeum vulgare). Most long-read assemblies are clearly superior to the current barley reference sequence based on short-reads. Assemblies derived from accurate long reads excel in most metrics, but the CCS approach was the most cost-effective strategy for assembling tens of barley genomes. A downsampling analysis indicated that 20-fold CCS coverage can yield very good sequence assemblies, while even five-fold CCS data may capture the complete sequence of most genes. We present an updated reference genome assembly for barley with near-complete representation of the repeat-rich intergenic space. Long-read assembly can underpin the construction of accurate and complete sequences of multiple genomes of a species to build pan-genome infrastructures in Triticeae crops and their wild relatives.

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., Pecan is an important specialty crop that has experienced extensive interspecific hybridization and nearly-obligate outcrossing. Here, the authors assemble diploid genomes of four outbred genotypes, identify interspecific introgressions through comparative genomics analyses, and map QTLs associated with pest resistance.

Significance Understanding how individual genetic loci contribute to trait variation across geographic space is of fundamental importance for understanding evolutionary adaptations. Our study demonstrates that most loci underlying locally adaptive trait variation have beneficial effects in some geographic regions while conferring little or no detectable cost in other parts of the geographic range of switchgrass over two field seasons of study. Thus, loci that contribute to local adaptation vary in the degree to which they are costly in alternative environments but typically confer greater benefits than costs. Further, our study suggests that breeding locally adapted varieties of switchgrass will be a boon to the biofuel industry, as locally adaptive loci could be combined to increase local yields in switchgrass., Local adaptation is the process by which natural selection drives adaptive phenotypic divergence across environmental gradients. Theory suggests that local adaptation results from genetic trade-offs at individual genetic loci, where adaptation to one set of environmental conditions results in a cost to fitness in alternative environments. However, the degree to which there are costs associated with local adaptation is poorly understood because most of these experiments rely on two-site reciprocal transplant experiments. Here, we quantify the benefits and costs of locally adaptive loci across 17° of latitude in a four-grandparent outbred mapping population in outcrossing switchgrass (Panicum virgatum L.), an emerging biofuel crop and dominant tallgrass species. We conducted quantitative trait locus (QTL) mapping across 10 sites, ranging from Texas to South Dakota. This analysis revealed that beneficial biomass (fitness) QTL generally incur minimal costs when transplanted to other field sites distributed over a large climatic gradient over the 2 y of our study. Therefore, locally advantageous alleles could potentially be combined across multiple loci through breeding to create high-yielding regionally adapted cultivars.

Centro de Investigación Forestal (CIFOR), Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable ele ments, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants., This work was supported by NSF DEB-1541005 and 1542609 and start-up funds from UF to S.F.M.; microMORPH Cross-Disciplinary Training Grant, Sigma-Xi Grant-In-Aid of Research, and Society for the Study of Evolution Rosemary Grant Award to S.B.C.; NSF DEB-1239992 to N.J.W.; the Emil Aaltonen Foundation and the University of Turku to S.O.; and NSF DEB-1541506 to J.G.B. and S.F.M. The work conducted by the U.S. Department of Energy Joint Genome Institute was supported by the Office of Science of the U.S. Department of Energy under contract no. DE-AC02-05CH11231., 12 Pág. Supplementary material for this article is available at http://advances.sciencemag.org/cgi/ content/full/7/27/eabh2488/DC1

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.

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent., Healey and colleagues presented a reference sequence assembly of Corymbia citriodora (spotted gum), a tree which is crucial for timber, pulp, and paper, as well as carbon sequestration and essential oil production.

Background: Seagrasses (Alismatales) are the only fully marine angiosperms. Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings. Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding. Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: As an important marine angiosperm, the improved Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

16 Pags.- 8 Figs. This article is licensed under a Creative Commons Attribution 4.0 International License., Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution., The work conducted by the US DOE Joint Genome Institute is supported by the Office of Science of the US Department of Energy under Contract no. DE-AC02-05CH11231 and User Agreement No. FP00006675 between the University of Zaragoza (User) and the DOE (Community Science Program 503504 proposal); S.J. was supported by the DOE Office of Science Graduate Student Research Program; P.C., B.C.-M., A.D.P., and R.S. were supported by a Spanish Ministry of Economy and Competitiveness grant CGL2016-79790-P; P.C., A.D.P., and R.S. were supported by a Bioflora A01-17 grant co-funded by the Spanish Aragon Government and the European Social Fund. R.H., E.W., and J.L. were supported by the National Science Centre Poland (grants DEC-2012/04/A/NZ3/00572 and DEC-2014/14/M/NZ2/00519). J.H.D. and C.N. were supported by Leverhulme Trust (grant no. 10754). A.C.E. and A.D. were supported by the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013)/ERC grant agreement no [GA335691 “Effectomics”], the Austrian Science Fund (FWF): [P27429-B22, P27818-B22, I 3033-B22], and the Austrian Academy of Science (OEAW).

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.

The ages and sizes of a sex-determination region (SDR) are difficult to determine in non-model species. Due to the lack of recombination and enrichment of repetitive elements in SDRs, the quality of assembly with short sequencing reads is universally low. Unique features present in the SDRs help provide clues about how SDRs are established and how they evolve in the absence of recombination. Several Populus species have been reported with a male heterogametic configuration of sex (XX/XY system) mapped on chromosome 19, but the exact location of the SDR has been inconsistent among species, and thus far, none of these SDRs has been fully assembled in a genomic context. Here we identify the Y-SDR from a Y-linked contig directly from a long-read PacBio assembly of a Populus trichocarpa male individual. We also identified homologous gene sequences in the SDR of P. trichocarpa and the SDR of the W chromosome in Salix purpurea. We show that inverted repeats (IRs) found in the Y-SDR and the W-SDR are lineage-specific. We hypothesize that, although the two IRs are derived from the same orthologous gene within each species, they likely have independent evolutionary histories. Furthermore, the truncated inverted repeats in P. trichocarpa may code for small RNAs that target the homologous gene for RNA-directed DNA methylation. These findings support the hypothesis that diverse sex-determining systems may be achieved through similar evolutionary pathways, thereby providing a possible mechanism to explain the lability of sex-determination systems in plants in general.

Pecan (Carya illinoinensis) is the most economically important member of the Carya genus and has been collected and evaluated across its broad geographic range in the process of crop improvement. In this study we obtained complete chloroplast genome sequences from two pecan genotypes, 87MX3-2.11 and the ‘Lakota’ cultivar (160,545 and 160,819 bp in length, respectively). The chloroplast genome of C. illinoinensis maintains the conserved structure typical of Juglandaceae and other land plants and is a circular molecule that includes a large single-copy (LSC) and a small single-copy (SSC) region, separated by a pair of inverted repeats (IRa and IRb). There were 124 genes found on the 87MX3-2.11 chloroplast genome and 123 on ‘Lakota’ (including multiple copies of the same gene), with 108 and 107 unique genes, respectively (counting only one copy of each gene). Different genes are found among C. illinoinensis, C. sinensis, and Juglans chloroplast genomes. C. illinoinensis is missing rps16 gene and has fewer copies of some tRNA genes, with ‘Lakota’ lacking a start codon of rps12 gene, compared with other related species. The nucleotide divergence between the two pecan chloroplast genomes reflects the genetic diversity of geographically separated populations of the species. Genomic divergence was also confirmed by the phylogenetic relationship of 19 whole chloroplast genome sequences representing Juglandaceae taxa. The complete chloroplast genome sequences in this study provide a foundation for understanding the influences of geographical adaptation, gene flow, and horticultural trait inheritance, in order to develop functional genomic tools for regional selection and pecan breeding.

BackgroundSex chromosomes in a wide variety of species share common characteristics, including the presence of suppressed recombination surrounding sex determination loci. They have arisen independently in numerous lineages, providing a conclusive example of convergent evolution. Mammalian sex chromosomes contain multiple palindromic repeats across the non-recombining region that facilitate sequence conservation through gene conversion, and contain genes that are crucial for sexual reproduction. Plant sex chromosomes are less well understood, and in particular it is not clear how coding sequence conservation is maintained in the absence of homologous recombination.ResultsHere we present the first evidence of large palindromic structures in a plant sex chromosome, based on a highly contiguous assembly of the W chromosome of the dioecious shrubSalix purpurea. Two consecutive palindromes span over a region of 200 kb, with conspicuous 20 kb stretches of highly conserved sequences among the four arms. The closely-related speciesS. suchowensisalso has two copies of a portion of the palindrome arm and provides strong evidence for gene conversion. Four genes in the palindrome are homologous to genes in the SDR of the closely-related genusPopulus, which is located on a different chromosome. These genes show distinct, floral-biased expression patterns compared to paralogous copies on autosomes.ConclusionThe presence of palindromic structures in sex chromosomes of mammals and plants highlights the intrinsic importance of these features in adaptive evolution in the absence of recombination. Convergent evolution is driving both the independent establishment of sex chromosomes as well as their fine-scale sequence structure.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the ‘pan-genome’1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley—comprising landraces, cultivars and a wild barley—that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

BackgroundPlants 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 eco-evolutionary context, particularly in long-lived perennials.ResultsHere, we generated a new high-quality reference genome from the oldest branch of a wildPopulus trichocarpatree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we used 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.ConclusionTaken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

Diverse wild and weedy crop relatives hold genetic variants underlying key evolutionary innovations of crops under domestication. Here, we provide genome resources and probe the genetic basis of domestication traits in green millet (Setaria viridis), a close wild relative of foxtail millet (S. italica). Specifically, we develop and exploit a platinum-quality genome assembly andde novoassemblies for 598 wild accessions to identify loci underlying a) response to climate, b) a key ‘loss of shattering’ trait that permits mechanical harvest, and c) leaf angle, a major predictor of yield in many grass crops. With CRISPR-Cas9 genome editing, we validatedLess Shattering1(SvLES1) as a novel gene for seed shattering, which is rendered non-functional via a retrotransposon insertion inSiLes1, the domesticated loss-of-shattering allele ofS. italica. Together these results and resources projectS. viridisas a key model species for complex trait dissection and biotechnological improvement of panicoid crops.

Allohexaploid (2n = 6x = 42) intermediate wheatgrass (Thinopyrum intermedium), abbreviated IWG, is an outcrossing perennial grass belonging to the tertiary gene pool of wheat. Perenniality would be valuable option for grain production, but attempts to introgress this complex trait from wheat-Thinopyrum hybrids have not been commercially successful. Efforts to breed IWG itself as a dual-purpose forage and grain crop have demonstrated useful progress and applications, but grain yields are significantly less than wheat. Therefore, genetic and physical maps have been developed to accelerate domestication of IWG. Herein, these maps were used to identify quantitative trait loci (QTLs) and candidate genes associated with IWG grain production traits in a family of 266 full-sib progenies derived from two heterozygous parents, M26 and M35. Transgressive segregation was observed for 17 traits related to seed size, shattering, threshing, inflorescence capacity, fertility, stem size, and flowering time. A total of 111 QTLs were detected in 36 different regions using 3826 genotype-by-sequence markers in 21 linkage groups. The most prominent QTL had a LOD score of 15 with synergistic effects of 29% and 22% over the family means for seed retention and percentage of naked seeds, respectively. Many QTLs aligned with one or more IWG gene models corresponding to 42 possible domestication orthogenes including the wheat Q and RHT genes. A cluster of seed-size and fertility QTLs showed possible alignment to a putative Z self-incompatibility gene, which could have detrimental grain-yield effects when genetic variability is low. These findings elucidate pathways and possible hurdles in the domestication of IWG.

Background: The availability of thousands of complete rice genome sequences from diverse varieties and accessions has laid the foundation for in-depth exploration of the rice genome. One drawback to these collections is that most of these rice varieties have long life cycles, and/or low transformation efficiencies, which limits their usefulness as model organisms for functional genomics studies. In contrast, the rice variety Kitaake has a rapid life cycle (9 weeks seed to seed) and is easy to transform and propagate. For these reasons, Kitaake has emerged as a model for studies of diverse monocotyledonous species. Results: Here, we report the de novo genome sequencing and analysis of Oryza sativa ssp. japonica variety KitaakeX, a Kitaake plant carrying the rice XA21 immune receptor. Our KitaakeX sequence assembly contains 377.6 Mb, consisting of 33 scaffolds (476 contigs) with a contig N50 of 1.4 Mb. Complementing the assembly are detailed gene annotations of 35,594 protein coding genes. We identified 331,335 genomic variations between KitaakeX and Nipponbare (ssp. japonica), and 2,785,991 variations between KitaakeX and Zhenshan97 (ssp. indica). We also compared Kitaake resequencing reads to the KitaakeX assembly and identified 219 small variations. The high-quality genome of the model rice plant KitaakeX will accelerate rice functional genomics. Conclusions: The high quality, de novo assembly of the KitaakeX genome will serve as a useful reference genome for rice and will accelerate functional genomics studies of rice and other species., BMC Genomics, 20, ISSN:1471-2164

Background The process of crop domestication often consists of two stages: initial domestication, where the wild species is first cultivated by humans, followed by diversification, when the domesticated species are subsequently adapted to more environments and specialized uses. Selective pressure to increase sugar accumulation in certain varieties of the cereal crop Sorghum bicolor is an excellent example of the latter; this has resulted in pronounced phenotypic divergence between sweet and grain-type sorghums, but the genetic mechanisms underlying these differences remain poorly understood. Results Here we present a new reference genome based on an archetypal sweet sorghum line and compare it to the current grain sorghum reference, revealing a high rate of nonsynonymous and potential loss of function mutations, but few changes in gene content or overall genome structure. We also use comparative transcriptomics to highlight changes in gene expression correlated with high stalk sugar content and show that changes in the activity and possibly localization of transporters, along with the timing of sugar metabolism play a critical role in the sweet phenotype. Conclusions The high level of genomic similarity between sweet and grain sorghum reflects their historical relatedness, rather than their current phenotypic differences, but we find key changes in signaling molecules and transcriptional regulators that represent new candidates for understanding and improving sugar metabolism in this important crop. Electronic supplementary material The online version of this article (10.1186/s12864-019-5734-x) contains supplementary material, which is available to authorized users.

American chestnut was once a foundation species of eastern North American forests, but was rendered functionally extinct in the early 20th century by an exotic fungal blight (Cryphonectria parasitica). Over the past 30 years, the American Chestnut Foundation (TACF) has pursued backcross breeding to generate hybrids that combine the timber‐type form of American chestnut with the blight resistance of Chinese chestnut based on a hypothesis of major gene resistance. To accelerate selection within two backcross populations that descended from two Chinese chestnuts, we developed genomic prediction models for five presence/absence blight phenotypes of 1,230 BC3F2 selection candidates and average canker severity of their BC3F3 progeny. We also genotyped pure Chinese and American chestnut reference panels to estimate the proportion of BC3F2 genomes inherited from parent species. We found that genomic prediction from a method that assumes an infinitesimal model of inheritance (HBLUP) has similar accuracy to a method that tends to perform well for traits controlled by major genes (Bayes C). Furthermore, the proportion of BC3F2 trees' genomes inherited from American chestnut was negatively correlated with the blight resistance of these trees and their progeny. On average, selected BC3F2 trees inherited 83% of their genome from American chestnut and have blight resistance that is intermediate between F1 hybrids and American chestnut. Results suggest polygenic inheritance of blight resistance. The blight resistance of restoration populations will be enhanced through recurrent selection, by advancing additional sources of resistance through fewer backcross generations, and by potentially by breeding with transgenic blight‐tolerant trees.

Life in high salinity environments poses challenges to cells in a variety of ways: maintenance of ion homeostasis and nutrient acquisition, often while concomitantly enduring saturating irradiances. Dunaliella salina has an exceptional ability to thrive even in saturated brine solutions. This ability has made it a model organism for studying responses to abiotic stress factors. Here we describe the occurrence of unique gene families, expansion of gene families, or gene losses that might be linked to osmoadaptive strategies. We discovered multiple unique genes coding for several of the homologous superfamily of the Ser-Thr-rich glycosyl-phosphatidyl-inositol-anchored membrane family and of the glycolipid 2-alpha-mannosyltransferase family, suggesting that such components on the cell surface are essential to life in high salt. Gene expansion was found in families that participate in sensing of abiotic stress and signal transduction in plants. One example is the patched family of the Sonic Hedgehog receptor proteins, supporting a previous hypothesis that plasma membrane sterols are important for sensing changes in salinities in D. salina. We also investigated genome-based capabilities regarding glycerol metabolism and present an extensive map for core carbon metabolism. We postulate that a second broader glycerol cycle exists that also connects to photorespiration, thus extending the previously described glycerol cycle. Further genome-based analysis of isoprenoid and carotenoid metabolism revealed duplications of genes for 1-deoxy-D-xylulose-5-phosphate synthase (DXS) and phytoene synthase (PSY), with the second gene copy of each enzyme being clustered together. Moreover, we identified two genes predicted to code for a prokaryotic-type phytoene desaturase (CRTI), indicating that D. salina may have eukaryotic and prokaryotic elements comprising its carotenoid biosynthesis pathways. In brief, our genomic data provide the basis for further gene discoveries regarding sensing abiotic stress, the metabolism of this halophilic alga, and its potential in biotechnological applications.

We report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G.max. The G.max assemblies provided are for widely used US cultivars: the northern line Williams82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G.soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7snpsperkb between these lines and G.soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G.soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G.soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

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., Drought is a major factor limiting crop productivity. Here, via eQTL analysis and comparative genomics, the authors show compensatory evolution between trans-regulatory loci and transcription factor binding sites that shape the drought response networks in the model C4 grass Panicum hallii.

Environmental heterogeneity can drive patterns of functional trait variation and lead to the formation of locally adapted ecotypes. Plant ecotypes are often differentiated by suites of correlated root and shoot traits that share common genetic, developmental, and physiological relationships. For instance, although plant water loss is largely governed by shoot systems, root systems determine water access and constrain shoot water status. To evaluate the genetic basis of root and shoot trait divergence, we developed a recombinant inbred population derived from mesic and xeric ecotypes of the perennial grass Panicum hallii. Our study sheds light on the genetic architecture underlying the relationships between root and shoot traits. We identified several genomic “hotspots” which control suites of correlated root and shoot traits, thus indicating genetic coordination between plant organ systems in the process of ecotypic divergence. Genomic regions of colocalized quantitative trait locus (QTL) for the majority of shoot and root growth related traits were independent of colocalized QTL for shoot and root resource acquisition traits. The allelic effects of individual QTL underscore ecological specialization for drought adaptation between ecotypes and reveal possible hybrid breakdown through epistatic interactions. These results have implications for understanding the factors constraining or facilitating local adaptation in plants.

Recent climate trends are driving rapid shifts in global precipitation patterns, leading to changes in soil water availability that can impact plant performance and distribution. Soil water availability is an especially important driver of contemporary evolution and ecotype formation in plant populations. In the process of ecotype formation, populations can diverge across many traits and exhibit different niche characteristics, which requires coordination between plant organ systems. For instance, although plant water loss is largely governed by shoot systems, root systems determine water access and constrain shoot water status. Understanding the genetic architecture of root traits and their relationship with shoot traits helps to develops a more complete picture of the adaptive differences that arise between ecotypes in response to changes in water availability. Panicum hallii is an emerging model system for C4 perennial grasses, including the important biofuel crop switchgrass ( Panicum virgatum ). Here, we produced an intercross between individuals of the xeric and mesic ecotypes of P. hallii , utilized a single seed decent method to generate a population of recombinant inbred lines (RIL) at the F7 generation and subsequently constructed a new genetic map based on whole genome re-sequencing. Utilizing extensive phenotyping of root and shoot traits and a quantitative genetic approach, we identified several genomic 9hotspots9 which control suites of correlated root and shoot traits, thus indicating genetic coordination between plant organ systems in the process of ecotypic divergence. In addition, we found that genomic regions of colocalized QTL for the majority of shoot and root growth related traits were independent of colocalized QTL for shoot and root resource acquisition traits. Finally, we confirmed major finding from the greenhouse in a field setting.

Author(s): Filiault, D; Ballerini, E; Mandakova, T; Akoz, G; Derieg, N; Schmutz, J; Jenkins, J; Grimwood, J; Shu, S; Hayes, R; Hellsten, U; Barry, K; Yan, J; Mihaltcheva, S; Karafiatova, M; Nizhynska, V; Lysak, M; Hodges, S; Nordborg, M | Abstract: The columbine genus Aquilegia is a classic example of an adaptive radiation, involving a wide variety of pollinators and habitats. Here we present the genome assembly of A. coerulea "Goldsmith", complemented by high-coverage sequencing data from 10 wild species covering the world-wide distribution. Our analysis reveals extensive allele sharing among species, and sheds light on the complex process of radiation. We also present the remarkable discovery that the evolutionary history of an entire chromosome differed from that of the rest of the genome --- a phenomenon which we do not fully understand, but which highlights the need to consider chromosomes in an evolutionary context.

Background Biofuels derived from lignocellulosic plant material are an important component of current renewable energy strategies. Improvement efforts in biofuel feedstock crops have been primarily focused on increasing biomass yield with less consideration for tissue quality or composition. Four primary components found in the plant cell wall contribute to the overall quality of plant tissue and conversion characteristics, cellulose and hemicellulose polysaccharides are the primary targets for fuel conversion, while lignin and ash provide structure and defense. We explore the genetic architecture of tissue characteristics using a quantitative trait loci (QTL) mapping approach in Panicum hallii, a model lignocellulosic grass system. Diversity in the mapping population was generated by crossing xeric and mesic varietals, comparative to northern upland and southern lowland ecotypes in switchgrass. We use near-infrared spectroscopy with a primary analytical method to create a P. hallii specific calibration model to quickly quantify cell wall components. Results Ash, lignin, glucan, and xylan comprise 68% of total dry biomass in P. hallii: comparable to other feedstocks. We identified 14 QTL and one epistatic interaction across these four cell wall traits and found almost half of the QTL to localize to a single linkage group. Conclusions Panicum hallii serves as the genomic model for its close relative and emerging biofuel crop, switchgrass (P. virgatum). We used high throughput phenotyping to map genomic regions that impact natural variation in leaf tissue composition. Understanding the genetic architecture of tissue traits in a tractable model grass system will lead to a better understanding of cell wall structure as well as provide genomic resources for bioenergy crop breeding programs. Electronic supplementary material The online version of this article (10.1186/s13068-018-1033-z) contains supplementary material, which is available to authorized users.

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

To identify switchgrass homologs of rice genes, known/predicted to control biomass and stress response-related traits, we screened 96,000 clones from two switchgrass bacterial artificial chromosome (BAC) libraries. Full-length sequencing of 311 BAC clones revealed sequence for ∼3.2 % (51.7 Mb) of the switchgrass genome, coding for 3948 genes. A comparison with Arabidopsis and five grass genomes revealed that switchgrass genes share the highest number of homologs with rice (95.5 %) followed by foxtail millet (91.7 %) and Sorghum (91.5 %). One hundred eighteen of the annotated genes are unique to switchgrass. Gene annotation and ontology analysis revealed 695 genes belonging to gene families targeted in the screening. These include 350 kinase, 203 glycosyltransferase (GT), 109 glycoside hydrolase (GH), and 33 ethylene responsive transcription factor (ERF) family genes. Rice homologs of 65 genes, identified here, have demonstrated roles in bioenergy-relevant traits. These include 14 GT2 family genes involved in the synthesis of cellulose and hemicelluloses. Comparative expression analysis in six switchgrass organs revealed a conserved expression pattern for three cellulose synthase (CesA1, CesA2, and CesA9) and five cellulose-synthase-like genes (CslA2, CslA11, CslC1, CslD4, and CslE6). CslF genes that encode mixed linkage glucans are expressed in wider range of tissues in switchgrass compared with rice.

The halotolerant alga Dunaliella salina is a model for stress tolerance and is used commercially for production of beta-carotene (=pro-vitamin A). The presented draft genome of the genuine strain CCAP19/18 will allow investigations into metabolic processes involved in regulation of stress responses, including carotenogenesis and adaptations to life in high-salinity environments.

Botryococcus braunii has long been known as a prodigious producer of liquid hydrocarbon oils that can be converted into combustion engine fuels. This draft genome for the B race of B. braunii will allow researchers to unravel important hydrocarbon biosynthetic pathways and identify possible regulatory networks controlling this unusual metabolism.

Setaria viridis is a rapid-life-cycle model panicoid grass. To identify genes that may contribute to inflorescence architecture and thus have the potential to influence grain yield in related crops such as maize, we conducted an N-nitroso-N-methylurea (NMU) mutagenesis of S. viridis and screened for visible inflorescence mutant phenotypes. Of the approximately 2,700 M2 families screened, we identified four recessive sparse panicle mutants (spp1-spp4) characterized by reduced and uneven branching of the inflorescence. To identify the gene underlying the sparse panicle1 (spp1) phenotype, we performed bulked segregant analysis and deep sequencing to fine map it to an approximately 1 Mb interval. Within this interval, we identified disruptive mutations in two genes. Complementation tests between spp1 and spp3 revealed they were allelic, and deep sequencing of spp3 identified an independent disruptive mutation in SvAUX1 (AUXIN1), one of the two genes in the ∼1 Mb interval and the only gene disruption shared between spp1 and spp3. SvAUX1 was found to affect both inflorescence development and root gravitropism in S. viridis. A search for orthologous mutant alleles in maize confirmed a very similar role of ZmAUX1 in maize, which highlights the utility of S. viridis in accelerating functional genomic studies in maize.

The anatomically simple plants that first colonized land must have acquired molecular and biochemical adaptations to drought stress. Abscisic acid (ABA) coordinates responses leading to desiccation tolerance in all land plants. We identified ABA nonresponsive mutants in the model bryophyte Physcomitrella patens and genotyped a segregating population to map and identify the ABA NON-RESPONSIVE (ANR) gene encoding a modular protein kinase comprising an N-terminal PAS domain, a central EDR domain, and a C-terminal MAPKKK-like domain. anr mutants fail to accumulate dehydration tolerance-associated gene products in response to drought, ABA, or osmotic stress and do not acquire ABA-dependent desiccation tolerance. The crystal structure of the PAS domain, determined to 1.7-A resolution, shows a conserved PAS-fold that dimerizes through a weak dimerization interface. Targeted mutagenesis of a conserved tryptophan residue within the PAS domain generates plants with ABA nonresponsive growth and strongly attenuated ABA-responsive gene expression, whereas deleting this domain retains a fully ABA-responsive phenotype. ANR orthologs are found in early-diverging land plant lineages and aquatic algae but are absent from more recently diverged vascular plants. We propose that ANR genes represent an ancestral adaptation that enabled drought stress survival of the first terrestrial colonizers but were lost during land plant evolution.

Background A landmark in soybean research, Glyma1.01, the first whole genome sequence of variety Williams 82 (Glycine max L. Merr.) was completed in 2010 and is widely used. However, because the assembly was primarily built based on the linkage maps constructed with a limited number of markers and recombinant inbred lines (RILs), the assembled sequence, especially in some genomic regions with sparse numbers of anchoring markers, needs to be improved. Molecular markers are being used by researchers in the soybean community, however, with the updating of the Glyma1.01 build based on the high-resolution linkage maps resulting from this research, the genome positions of these markers need to be mapped. Results Two high density genetic linkage maps were constructed based on 21,478 single nucleotide polymorphism loci mapped in the Williams 82 x G. soja (Sieb. & Zucc.) PI479752 population with 1083 RILs and 11,922 loci mapped in the Essex x Williams 82 population with 922 RILs. There were 37 regions or single markers where marker order in the two populations was in agreement but was not consistent with the physical position in the Glyma1.01 build. In addition, 28 previously unanchored scaffolds were positioned. Map data were used to identify false joins in the Glyma1.01 assembly and the corresponding scaffolds were broken and reassembled to the new assembly, Wm82.a2.v1. Based upon the plots of the genetic on physical distance of the loci, the euchromatic and heterochromatic regions along each chromosome in the new assembly were delimited. Genomic positions of the commonly used markers contained in BARCSOYSSR_1.0 database and the SoySNP50K BeadChip were updated based upon the Wm82.a2.v1 assembly. Conclusions The information will facilitate the study of recombination hot spots in the soybean genome, identification of genes or quantitative trait loci controlling yield, seed quality and resistance to biotic or abiotic stresses as well as other genetic or genomic research. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2344-0) contains supplementary material, which is available to authorized users.

Polyploid species have long been thought to be recalcitrant to whole-genome assembly. By combining high-throughput sequencing, recent developments in parallel computing, and genetic mapping, we derive, de novo, a sequence assembly representing 9.1 Gbp of the highly repetitive 16 Gbp genome of hexaploid wheat, Triticum aestivum, and assign 7.1 Gb of this assembly to chromosomal locations. The genome representation and accuracy of our assembly is comparable or even exceeds that of a chromosome-by-chromosome shotgun assembly. Our assembly and mapping strategy uses only short read sequencing technology and is applicable to any species where it is possible to construct a mapping population. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0582-8) contains supplementary material, which is available to authorized users.

International audience; Genetic maps are key tools in genetic research as they constitute the framework for many applications, such as quantitative trait locus analysis, and support the assembly of genome sequences. The resequencing of the two parents of a cross between Eucalyptus urophylla and Eucalyptus grandis was used to design a single nucleotide polymorphism (SNP) array of 6000 markers evenly distributed along the E. grandis genome. The genotyping of 1025 offspring enabled the construction of two high-resolution genetic maps containing 1832 and 1773 markers with an average marker interval of 0.45 and 0.5 cM for E. grandis and E. urophylla, respectively. The comparison between genetic maps and the reference genome highlighted 85% of collinear regions. A total of 43 noncollinear regions and 13 nonsynthetic regions were detected and corrected in the new genome assembly. This improved version contains 4943 scaffolds totalling 691.3 Mb of which 88.6% were captured by the 11 chromosomes. The mapping data were also used to investigate the effect of population size and number of markers on linkage mapping accuracy. This study provides the most reliable linkage maps for Eucalyptus and version 2.0 of the E. grandis genome.

Fixed chromosomal inversions can reduce gene flow and promote speciation in two ways: by suppressing recombination and by carrying locally favored alleles at multiple loci. However, it is unknown whether favored mutations slowly accumulate on older inversions or if young inversions spread because they capture preexisting adaptive Quantitative Trait Loci (QTLs). By genetic mapping, chromosome painting and genome sequencing we have identified a major inversion controlling ecologically important traits in Boechera stricta. The inversion arose since the last glaciation and subsequently reached local high frequency in a hybrid speciation zone. Furthermore, the inversion shows signs of positive directional selection. To test whether the inversion could have captured existing, linked QTLs, we crossed standard, collinear haplotypes from the hybrid zone and found multiple linked phenology QTLs within the inversion region. These findings provide the first direct evidence that linked, locally adapted QTLs may be captured by young inversions during incipient speciation.

Scott Jackson, Jeremy Schmutz, Phillip McClean and colleagues report the genome sequence of the common bean (Phaseolus vulgaris) and resequenced wild individuals and landraces from Mesoamerican and Andean gene pools, showing that common bean underwent two independent domestications. Supplementary information The online version of this article (doi:10.1038/ng.3008) contains supplementary material, which is available to authorized users., Common bean (Phaseolus vulgaris L.) is the most important grain legume for human consumption and has a role in sustainable agriculture owing to its ability to fix atmospheric nitrogen. We assembled 473 Mb of the 587-Mb genome and genetically anchored 98% of this sequence in 11 chromosome-scale pseudomolecules. We compared the genome for the common bean against the soybean genome to find changes in soybean resulting from polyploidy. Using resequencing of 60 wild individuals and 100 landraces from the genetically differentiated Mesoamerican and Andean gene pools, we confirmed 2 independent domestications from genetic pools that diverged before human colonization. Less than 10% of the 74 Mb of sequence putatively involved in domestication was shared by the two domestication events. We identified a set of genes linked with increased leaf and seed size and combined these results with quantitative trait locus data from Mesoamerican cultivars. Genes affected by domestication may be useful for genomics-enabled crop improvement. Supplementary information The online version of this article (doi:10.1038/ng.3008) contains supplementary material, which is available to authorized users.

Background Theobroma cacao L. cultivar Matina 1-6 belongs to the most cultivated cacao type. The availability of its genome sequence and methods for identifying genes responsible for important cacao traits will aid cacao researchers and breeders. Results We describe the sequencing and assembly of the genome of Theobroma cacao L. cultivar Matina 1-6. The genome of the Matina 1-6 cultivar is 445 Mbp, which is significantly larger than a sequenced Criollo cultivar, and more typical of other cultivars. The chromosome-scale assembly, version 1.1, contains 711 scaffolds covering 346.0 Mbp, with a contig N50 of 84.4 kbp, a scaffold N50 of 34.4 Mbp, and an evidence-based gene set of 29,408 loci. Version 1.1 has 10x the scaffold N50 and 4x the contig N50 as Criollo, and includes 111 Mb more anchored sequence. The version 1.1 assembly has 4.4% gap sequence, while Criollo has 10.9%. Through a combination of haplotype, association mapping and gene expression analyses, we leverage this robust reference genome to identify a promising candidate gene responsible for pod color variation. We demonstrate that green/red pod color in cacao is likely regulated by the R2R3 MYB transcription factor TcMYB113, homologs of which determine pigmentation in Rosaceae, Solanaceae, and Brassicaceae. One SNP within the target site for a highly conserved trans-acting siRNA in dicots, found within TcMYB113, seems to affect transcript levels of this gene and therefore pod color variation. Conclusions We report a high-quality sequence and annotation of Theobroma cacao L. and demonstrate its utility in identifying candidate genes regulating traits.

International audience; The shift from outcrossing to selfing is common in flowering plants, but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible

The oomycete vegetable pathogen Phytophthora capsici has shown remarkable adaptation to fungicides and new hosts. Like other members of this destructive genus, P. capsici has an explosive epidemiology, rapidly producing massive numbers of asexual spores on infected hosts. In addition, P. capsici can remain dormant for years as sexually-recombined oospores, making it difficult to produce crops at infested sites, and allowing outcrossing populations to maintain significant genetic variation. Genome sequencing, development of a high-density genetic map, and integrative genomic/genetic characterization of P. capsici field isolates and intercross progeny revealed significant mitotic loss of heterozygosity (LOH) and higher levels of SNVs than those reported for humans, plants, and P. infestans. LOH was detected in clonally propagated field isolates and sexual progeny, cumulatively affecting >30% of the genome. LOH altered genotypes for more than 11,000 single nucleotide variant (SNV) sites and showed a strong association with changes in mating type and pathogenicity. Overall, it appears that LOH may provide a rapid mechanism for fixing alleles and may be an important component of adaptability for P. capsici.

SummaryGenomic comparisons of chordates, hemichordates, and echinoderms can inform hypotheses for the evolution of these strikingly different phyla from the last common deuterostome ancestor [1–5]. Because hox genes play pivotal developmental roles in bilaterian animals [6–8], we analyzed the Hox complexes of two hemichordate genomes. We find that Saccoglossus kowalevskii and Ptychodera flava both possess 12-gene clusters, with mir10 between hox4 and hox5, in 550 kb and 452 kb intervals, respectively. Genes hox1–hox9/10 of the clusters are in the same genomic order and transcriptional orientation as their orthologs in chordates, with hox1 at the 3′ end of the cluster. At the 5′ end, each cluster contains three posterior genes specific to Ambulacraria (the hemichordate-echinoderm clade), two forming an inverted terminal pair. In contrast, the echinoderm Strongylocentrotus purpuratus contains a 588 kb cluster [9] of 11 orthologs of the hemichordate genes, ordered differently, plausibly reflecting rearrangements of an ancestral hemichordate-like ambulacrarian cluster. Hox clusters of vertebrates and the basal chordate amphioxus [10] have similar organization to the hemichordate cluster, but with different posterior genes. These results provide genomic evidence for a well-ordered complex in the deuterostome ancestor for the hox1–hox9/10 region, with the number and kind of posterior genes still to be elucidated.

We generated a high-quality reference genome sequence for foxtail millet (Setaria italica). The ~400-Mb assembly covers ~80% of the genome and >95% of the gene space. The assembly was anchored to a 992-locus genetic map and was annotated by comparison with >1.3 million expressed sequence tag reads. We produced more than 580 million RNA-Seq reads to facilitate expression analyses. We also sequenced Setaria viridis, the ancestral wild relative of S. italica, and identified regions of differential single-nucleotide polymorphism density, distribution of transposable elements, small RNA content, chromosomal rearrangement and segregation distortion. The genus Setaria includes natural and cultivated species that demonstrate a wide capacity for adaptation. The genetic basis of this adaptation was investigated by comparing five sequenced grass genomes. We also used the diploid Setaria genome to evaluate the ongoing genome assembly of a related polyploid, switchgrass (Panicum virgatum).