Plant genetic resources are considered a treasure trove of valuable, untapped diversity that holds the key to breeding the crops of the future. However, the use of these resources in breeding is often limited due to the lack of comprehensive phenotypic characterization. The present study provides extensive historical phenotypic data from nine genebanks as a MIAPPE compliant data set. We compiled and curated phenotypic data from 43,293 wheat accessions, encompassing 460,399 data points across 52 traits, including the three core traits of plant height, heading time, and thousand kernel weight from seven decades. The exceptional quality of the presented dataset was highlighted by predominantly high heritabilities. Phenotypic data of such quantity and quality is a crucial resource for unlocking the valuable diversity of plant genetic resources for agricultural advancement.
Wheat genetic resources hold the diversity required to mitigate agricultural challenges from climate change and reduced inputs. Using DArTseq, we genotype 461 wheat landraces and cultivars and evaluate them for powdery mildew resistance. By developing a k-mer-based GWAS approach with fully assembled genomes of Triticum aestivum and its progenitors, we uncover 25% more resistance-associated k-mers than single-reference methods, outperforming SNP-based GWAS in both loci detection and mapping precision. In total, we detect 34 powdery mildew resistance loci, including 27 potentially novel regions. Our approach underscores the importance of integrating multiple reference genomes to unlock the potential of wheat germplasm.
Background: Wheat landraces and cultivars stored in gene banks worldwide represent a valuable source of genetic diversity for discovering genes critical for agriculture, which is increasingly constrained by climate change and inputs reduction. We assembled and genotyped, using DArTseq technology, a panel of 461 accessions representative of the genetic diversity of Swiss wheat material. The collection was evaluated for powdery mildew resistance under field conditions for two consecutive years and at the seedling stage with 10 different wheat powdery mildew isolates. Results: To identify the genetic basis of mildew resistance in wheat, we developed a k-mer-based GWAS approach using multiple fully-assembled genomes including Triticum aestivum as well as four progenitor genomes. Compared to approaches based on single reference genomes, we unambiguously mapped an additional 25% resistance-associated k-mers. Our approach outperformed SNP-based GWAS in terms of number of loci identified and precision of mapping. In total, we detected 34 (Pm) powdery mildew resistance loci, including seven previously-described and more importantly 27 novel loci active at the seedling stage. Furthermore, we identified a region associated with adult plant resistance, which was not detected with SNP-based approaches. Conclusions: The described non-reference-based approach highlights the potential of integrating multiple wheat reference genomes with k-mer GWAS to harness the untapped genetic diversity present in germplasm collections. ### Competing Interest Statement The authors have declared no competing interest.
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.
The YABBY gene family is one of the plant transcription factors present in all seed plants. The family members were extensively studied in various plants and shown to play important roles in plant growth and development, such as the polarity establishment in lateral organs, the formation and development of leaves and flowers, and the response to internal plant hormone and external environmental stress signals. In this study, a total of 364 YABBY genes were identified from 37 Brassicaceae genomes, of which 15 were incomplete due to sequence gaps, and nine were imperfect (missing C2C2 zinc-finger or YABBY domain) due to sequence mutations. Phylogenetic analyses resolved these YABBY genes into six compact clades except for a YAB3-like gene identified in Aethionema arabicum. Seventeen Brassicaceae species each contained a complete set of six basic YABBY genes (i.e., 1 FIL, 1 YAB2, 1 YAB3, 1 YAB5, 1 INO and 1 CRC), while 20 others each contained a variable number of YABBY genes (5–25) caused mainly by whole-genome duplication/triplication followed by gene losses, and occasionally by tandem duplications. The fate of duplicate YABBY genes changed considerably according to plant species, as well as to YABBY gene type. These YABBY genes were shown to be syntenically conserved across most of the Brassicaceae species, but their functions might be considerably diverged between species, as well as between paralogous copies, as demonstrated by the promoter and expression analysis of YABBY genes in two Brassica species (B. rapa and B. oleracea). Our study provides valuable insights for understanding the evolutionary story of YABBY genes in Brassicaceae and for further functional characterization of each YABBY gene across the Brassicaceae species.
Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.
High-quality nanopore genome assemblies were generated for two Brassica nigra genotypes (Ni100 and CN115125); a member of the agronomically important Brassica species. The N50 contig length for the two assemblies were 17.1 Mb (58 contigs) and 0.29 Mb (963 contigs), respectively, reflecting recent improvements in the technology. Comparison with a de novo short read assembly for Ni100 corroborated genome integrity and quantified sequence related error rates (0.002%). The contiguity and coverage allowed unprecedented access to low complexity regions of the genome. Pericentromeric regions and coincidence of hypo-methylation enabled localization of active centromeres and identified a novel centromere-associated ALE class I element which appears to have proliferated through relatively recent nested transposition events (<1 million years ago). Computational abstraction was used to define a post-triplication Brassica specific ancestral genome and to calculate the extensive rearrangements that define the genomic distance separating B. nigra from its diploid relatives.
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.
It is only recently, with the advent of long-read sequencing technologies, that we are beginning to uncover previously uncharted regions of complex and inherently recursive plant genomes. To comprehensively study and exploit the genome of the neglected oilseed Brassica nigra , we generated two high-quality nanopore de novo genome assemblies. The N50 contig lengths for the two assemblies were 17.1 Mb (12 contigs), one of the best among 324 sequenced plant genomes, and 0.29 Mb (424 contigs), respectively, reflecting recent improvements in the technology. Comparison with a de novo short-read assembly corroborated genome integrity and quantified sequence-related error rates (0.2%). The contiguity and coverage allowed unprecedented access to low-complexity regions of the genome. Pericentromeric regions and coincidence of hypomethylation enabled localization of active centromeres and identified centromere-associated ALE family retro-elements that appear to have proliferated through relatively recent nested transposition events (<1 Ma). Genomic distances calculated based on synteny relationships were used to define a post-triplication Brassica- specific ancestral genome, and to calculate the extensive rearrangements that define the evolutionary distance separating B. nigra from its diploid relatives.
Rapeseed or canola (B. napus) is the second most important oilseed crop of the world. It is also a favourite plant for basic and breeding research. Due to its origin and evolution, rapeseed has a complex polyploid genome. Recent sequencing of the corresponding genomes provides the basis for a better understanding and exploitation of the genetic diversity involved in major rapeseed traits. However, directed selection for major quality characteristics, i.e. minimal erucic acid content and low glucosinolate level, has caused genetic bottlenecks limiting genetic variation in the current gene pools of cultivated oilseed rape (OSR). Therefore, broadening genetic diversity is an important aim of research and a necessary prerequisite for further progress by OSR breeding. In agricultural production, rapeseed is nowadays an indispensable component of crop rotations in major growing areas such as Australia, Western Canada, Central China and many countries of the European Union. In many cases, OSR is the only leaf crop among dominating cereal species. Therefore, OSR as a component of crop rotations helps to maintain soil fertility and contributes to sustainable production therefore. As a major cash crop OSR substantially contributes to farmers’ incomes and therefore helps to stabilize rural populations. Beyond that, as major globally traded agricultural commodities rapeseed/ canola and rapeseed/canola oil and meal significantly input the national products of a number of countries e.g. Canada. Rapeseed/ canola is a raw material for vegetable oil and extraction meal as feed, food and fuel. The oil is mainly used as a high-value salad oil for dressings etc. due to its high contents of oleic acid (ca. 60%) and poly-unsaturated linolenic acid (omega-3, ca. 10%). Nevertheless, a large part is also used as a mobility fuel for diesel cars and tractors, particularly in Germany and Europe. The extraction meal (and protein) from oil processing is now recognized as a highly valuable animal feed, particularly for ruminants (cattle) but also for monogastric farm animals (pigs, poultry). Furthermore, the interest in rapeseed protein for the purpose of human nutrition is increasing. Optimal contents of the major compounds mentioned before represent the main requirements for rapeseed/canola varieties today. Consequently, quality characteristics are major criteria for variety testing and registration therefore. Other major requirements for modern rapeseed varieties are yield W. Friedt (&) Department of Plant Breeding, University of Giessen, 35392 Giessen, Germany e-mail: Wolfgang.Friedt@agrar.uni-giessen.de J. Tu T. Fu College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, China © Springer Nature Switzerland AG 2018 S. Liu et al. (eds.), The Brassica napus Genome, Compendium of Plant Genomes, https://doi.org/10.1007/978-3-319-43694-4_1 1 and various agronomic traits securing seed yield, i.e. resistance against fungal diseases and insect pests. Because of the environmental concerns and the rejection of agrochemicals such as insecticides (e.g. ban of neonicotinoids in the EU), genetic approaches of establishing resistant crop cultivars constantly gain importance. OSR breeding has long been a relatively ordinary process of repeated selection for resistance, quality and yield, the propagation of improved populations and their release as new open pollinated (OP) varieties. Since this approach is not very effective, breeders have been interested in breeding hybrids instead. Today, F1 hybrids represent the major variety type. They are produced with genetic male sterility systems, most of which are based on cytoplasmic mutations causing male sterility (cms). The higher performance of hybrids is caused by “heterotic effect” which largely depends on the genetic distance between the parents. Therefore, distant genetic pools need to be generated for the development of female and male parents. The future potential of variety design is consequently determined by the usefulness of the genetic pools and the performance of hybrid parents extracted from them. Modern breeding tools based on biotechnology and genomics can substantially contribute to a better exploitation of useful genetic diversity, i.e. specific genes and genetic networks. Better varieties are a precondition for further crop improvement. Future quality OSR cultivars deserve high yield potential, combined with good stability due to disease and pest resistance. This will be the basis for exploiting the great agronomical and industrial advantages of the rapeseed plant. 1.1 General Relevance of Brassica napus Scientific research on Brassica napus is particularly relevant because of the following reasons: (i) It is a major crop plant and member of a large plant family (Brassicaceae) including many crops and weeds with many unique features. (ii) It also represents one of the few dicot plants which are particularly important today to break plant production systems (rotations) often dominated by cereals and
The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.
Summary Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.
Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.
Rapeseed (Brassica napus), an important oilseed crop, has adapted to diverse climate zones and latitudes by forming three main ecotype groups, namely winter, semi-winter, and spring types. However, genetic variations underlying the divergence of these ecotypes are largely unknown. Here, we report the global pattern of genetic polymorphisms in rapeseed determined by resequencing a worldwide collection of 991 germplasm accessions. A total of 5.56 and 5.53 million single-nucleotide polymorphisms (SNPs) as well as 1.86 and 1.92 million InDels were identified by mapping reads to the reference genomes of "Darmor-bzh'' and "Tapidor,'' respectively. We generated a map of allelic drift paths that shows splits and mixtures of the main populations, and revealed an asymmetric evolution of the two subgenomes of B. napus by calculating the genetic diversity and linkage disequilibrium parameters. Selective-sweep analysis revealed genetic changes in genes orthologous to those regulating various aspects of plant development and response to stresses. A genome-wide association study identified SNPs in the promoter regions of FLOWERING LOCUS T and FLOWERING LOCUS C orthologs that corresponded to the different rapeseed ecotype groups. Our study provides important insights into the genomic footprints of rapeseed evolution and flowering-time divergence among three ecotype groups, and will facilitate screening of molecular markers for accelerating rapeseed breeding.
Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (Ar ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with Ar , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.
More and more RING finger genes were found to be implicated in various important biological processes. In the present study, a total of 731 RING domains in 715 predicted proteins were identified in Brassica rapa genome (AA, 2n = 20), which were further divided into eight types: RING-H2 (371), RING-HCa (215), RING-HCb (47), RING-v (44), RING-C2 (38), RING-D (10), RING-S/T (5) and RING-G (1). The 715 RING finger proteins were further classified into 51 groups according to the presence of additional domains. 700 RING finger protein genes were mapped to the 10 chromosomes of B. rapa with a range of 47 to 111 genes for each chromosome. 667 RING finger protein genes were expressed in at least one of the six tissues examined, indicating their involvement in various physiological and developmental processes in B. rapa . Hierarchical clustering analysis of RNA-seq data divided them into seven major groups, one of which includes 231 members preferentially expressed in leaf, and constitutes then a panel of gene candidates for studying the genetic and molecular mechanisms of leafy head traits in Brassica crops. Our results lay the foundation for further studies on the classification, evolution and putative functions of RING finger protein genes in Brassica species.
A large panel of hexaploid wheat synthetics was developed. Their tetraploid parents consisted of either four extracted wheat tetraploids (ETWs) or four natural present-day tetraploids, and their diploid parents consisted of twenty accessions of Aegilops tauschii . Analysis of meiotic behaviour of the synthetics showed that chromosome pairing is highly variable and depends on the progenitor. The meiotic behaviour in the four ETWs was compared to that of the natural tetraploid wheats. It appears there was no evolution at the hexaploid level of the meiotic genes carried by the A and B genomes. We also reach the conclusion that the neo-allohexaploids at the origin of present-day wheat had a meiotic behaviour close to that of the present-day hexaploid wheat. It is likely that other neo-hexaploids with an impaired meiosis were formed, but they had no future due to their more or less rapid disappearance due to increasing aneuploidy level and structural changes, mainly Robertsonian translocations.