Chromosome-scale comparative sequence analysis unravels molecular mechanisms of genome dynamics between two wheat cultivars

Fine-tuning of heading time by earliness per se effect due to multi-allelic variants in VRN-B3 locus of hexaploid wheat

MAIN CONCLUSION

Wheat VRN-B3 contains multi-allelic variants conferring earliness per se effects and can generate a new allele by recombining multiple alleles, highlighting its importance for the fine-tuning of heading time. Fine-tuning of heading time is required for breeding well-adapted varieties of wheat in regional environments. The VRN-B3 locus, which encodes the FT-B1 gene, is known as a vernalization gene. In this study, we analyzed two alleles in the VRN-B3 region: QHt-7B_Zen of a Japanese variety (Zenkouji-komugi) and QHt-7B_spelt of a spelt wheat strain (st. Rumania). Phenotypic evaluation of near-isogeneic lines (NILs) of 'Chinese Spring' (CS) under long-day (16 h) conditions showed that QHt-7B_Zen and QHt-7B_spelt conferred approximately 3.9 d earlier and 3.0 d later heading time compared with CS, respectively. Differences in heading times among NILs were also observed for fully vernalized plants under long-day conditions, indicating their earliness per se effect. Both QTLs behaved as single genes with incomplete dominant effects, and fine-mapping showed that FT-B1 was responsible for heading time. Droplet digital PCR analysis revealed that QHt-7B_Zen contained three copies of FT-B1, similar to the CS. QHt-7B_spelt had one FT-B1 copy with 14 substitutions, a 15 bp insertion in the 4.8 kb promoter region, and one amino acid substitution in the third exon, which could be designated as a novel allele, Vrn-B3f. Furthermore, a new allele with two FT-B1 copies conferring an intermediate heading time between the parents was created by the recombination of FT-B1 copies of CS and NIL for QHt-7B_spelt. Our findings indicate that fine-tuning heading time is possible through the versatility of the VRN-B3 locus, which can generate multi-allelic variants that have both vernalization and earliness per se effects in wheat.

Evolution of Plant Genome Size and Composition

The rapid development of sequencing technology has led to an explosion of plant genome data, opening up more opportunities for research in the field of comparative evolutionary analysis of plant genomes. In this review, we focus on changes in plant genome size and composition, examining the effects of polyploidy, whole-genome duplication, and alternations in transposable elements on plant genome architecture and evolution, respectively. In addition, to address gaps in the available information, we also collected and analyzed 234 representative plant genome data as a supplement. We aim to provide a comprehensive, up-to-date summary of information on plant genome architecture and evolution in this review.

NRC Immune receptor networks show diversified hierarchical genetic architecture across plant lineages

Plants' complex immune systems include nucleotide-binding domain and leucine-rich repeat-containing (NLR) proteins, which help recognize invading pathogens. In solanaceous plants, the NRC (NLR required for cell death) family includes helper NLRs that form a complex genetic network with multiple sensor NLRs to provide resistance against pathogens. However, the evolution and function of NRC networks outside solanaceous plants are currently unclear. Here, we conducted phylogenomic and macroevolutionary analyses comparing NLRs identified from different asterid lineages and found that NRC networks expanded significantly in most lamiids but not in Ericales and campanulids. Using transient expression assays in Nicotiana benthamiana, we showed that NRC networks are simple in Ericales and campanulids, but have high complexity in lamiids. Phylogenetic analyses grouped the NRC helper NLRs into three NRC0 subclades that are conserved, and several family-specific NRC subclades of lamiids that show signatures of diversifying selection. Functional analyses revealed that members of the NRC0 subclades are partially interchangeable, whereas family-specific NRC members in lamiids lack interchangeability. Our findings highlight the distinctive evolutionary patterns of the NRC networks in asterids and provide potential insights into transferring disease resistance across plant lineages.

Finding needles in a haystack: identification of inter-specific introgressions in wheat genebank collections using low-coverage sequencing data

Recently, entire genebank collections of wheat have been extensively characterized with sequencing data. We have identified introgressions using these genotyping-by-sequencing and whole-genome sequencing data. On the basis of our results, we provide information about predicted introgressions at 1-Mb resolution for 9,172 wheat samples as a resource for breeders and scientists. We recommend that all plant genetic resources, including genebank collections, be characterized using a combination of variant calling and introgression prediction. This is necessary to identify potential duplicates in collections efficiently and reliably, and to select promising germplasms with potentially beneficial introgressions for further characterization and prospective breeding application.

A survey of lineage-specific genes in Triticeae reveals de novo gene evolution from genomic raw material

Diploid plant genomes typically contain ~35,000 genes, almost all belonging to highly conserved gene families. Only a small fraction are lineage-specific, which are found in only one or few closely related species. Little is known about how genes arise de novo in plant genomes and how often this occurs; however, they are believed to be important for plants diversification and adaptation. We developed a pipeline to identify lineage-specific genes in Triticeae, using newly available genome assemblies of wheat, barley, and rye. Applying a set of stringent criteria, we identified 5942 candidate Triticeae-specific genes (TSGs), of which 2337 were validated as protein-coding genes in wheat. Differential gene expression analyses revealed that stress-induced wheat TSGs are strongly enriched in putative secreted proteins. Some were previously described to be involved in Triticeae non-host resistance and cold response. Additionally, we show that 1079 TSGs have sequence homology to transposable elements (TEs), ~68% of them deriving from regulatory non-coding regions of Gypsy retrotransposons. Most importantly, we demonstrate that these TSGs are enriched in transmembrane domains and are among the most highly expressed wheat genes overall. To summarize, we conclude that de novo gene formation is relatively rare and that Triticeae probably possess ~779 lineage-specific genes per haploid genome. TSGs, which respond to pathogen and environmental stresses, may be interesting candidates for future targeted resistance breeding in Triticeae. Finally, we propose that non-coding regions of TEs might provide important genetic raw material for the functional innovation of TM domains and the evolution of novel secreted proteins.

Resistance of QYm.nau-2D to wheat yellow mosaic virus was derived from an alien introgression into common wheat

The QYm.nau-2D locus conferring wheat yellow mosaic virus resistance is an exotic introgression and we developed 11 diagnostic markers tightly linked to QYm.nau-2D. Wheat yellow mosaic virus (WYMV) is a serious disease of winter wheat in China. Breeding resistant varieties is the most effective strategy for WYMV control. A WYMV resistant locus QYm.nau-2D on the chromosome arm 2DL has been repeatedly reported but the mapped region is large. In the present study, we screened recombinants using a biparental population and mapped QYm.nau-2D into an 18.8 Mb physical interval. By genome-wide association studies of 372 wheat varieties for WYMV resistance in four environments, we narrowed down QYm.nau-2D into a 16.4 Mb interval. Haplotype analysis indicated QYm.nau-2D were present as six different states due to recombination during hybridization breeding. QYm.nau-2D was finally mapped into a linkage block of 11.2 Mb. Chromosome painting using 2D specific probes and collinearity analysis among the published sequences corresponding to QYm.nau-2D region indicated the block was an exotic introgression. The Illumina-sequenced reads of four diploid Aegilops species were mapped to the sequence of Fielder, a variety having the introgression. The mapping reads were significantly increased at the putative introgression regions of Fielder. Ae. uniaristata (NN) had the highest mapping reads, suggesting that QYm.nau-2D was possibly an introgression from genome N. We investigated the agronomic performances of different haplotypes and observed no linkage drag of the alien introgression for the 15 tested traits. For marker-assisted selection of QYm.nau-2D, we developed 11 diagnostic markers tightly linked to the locus. This research provided a case study of an exotic introgression, which has been utilized in wheat improvement for WYMV resistance.

Patterns of genetic variation and QTLs controlling grain traits in a collection of global wheat germplasm revealed by high-quality SNP markers

BACKGROUND

Establish a molecular breeding program involved assembling a diverse germplasm collection and generating accurate genotypes to characterize their genetic potential and associate them with agronomic traits. In this study, we acquired over eight hundred wheat accessions from international gene banks and assessed their genetic relatedness using high-quality SNP genotypes. Understanding the scope of genomic variation in this collection allows the breeders to utilize the genetic resources efficiently while improving wheat yield and quality.

RESULTS

A wheat diversity panel comprising 39 durum wheat, 60 spelt wheat, and 765 bread wheat accessions was genotyped on iSelect 90 K wheat SNP arrays. A total of 57,398 SNP markers were mapped to IWGSC RefSeq v2.1 assembly, over 30,000 polymorphic SNPs in the A, B, D genomes were used to analyze population structure and diversity, the results revealed the separation of the three species and the differentiation of CIMMYT improved breeding lines and landraces or widely grown cultivars. In addition, several chromosomal regions under selection were detected. A subset of 280 bread wheat accessions was evaluated for grain traits, including grain length, width, surface area, and color. Genome-wide association studies (GWAS) revealed that several chromosomal regions were significantly linked to known quantitative trait loci (QTL) controlling grain-related traits. One of the SNP peaks at the end of chromosome 7A was in strong linkage disequilibrium (LD) with WAPO-A1, a gene that governs yield components.

CONCLUSIONS

Here, the most updated and accurate physical positions of SNPs on 90 K genotyping array are provided for the first time. The diverse germplasm collection and associated genotypes are available for the wheat researchers to use in their molecular breeding program. We expect these resources to broaden the genetic basis of original breeding and pre-breeding materials and ultimately identify molecular markers associated with important agronomic traits which are evaluated in diverse environmental conditions.

ggComp enables dissection of germplasm resources and construction of a multiscale germplasm network in wheat

Accurate germplasm characterization is a vital step for accelerating crop genetic improvement, which remains largely infeasible for crops such as bread wheat (Triticum aestivum L.), which has a complex genome that undergoes frequent introgression and contains many structural variations. Here, we propose a genomic strategy called ggComp, which integrates resequencing data with copy number variations and stratified single-nucleotide polymorphism densities to enable unsupervised identification of pairwise germplasm resource-based Identity-By-Descent (gIBD) blocks. The reliability of ggComp was verified in wheat cultivar Nongda5181 by dissecting parental-descent patterns represented by inherited genomic blocks. With gIBD blocks identified among 212 wheat accessions, we constructed a multi-scale genomic-based germplasm network. At the whole-genome level, the network helps to clarify pedigree relationship, demonstrate genetic flow, and identify key founder lines. At the chromosome level, we were able to trace the utilization of 1RS introgression in modern wheat breeding by hitchhiked segments. At the single block scale, the dissected germplasm-based haplotypes nicely matched with previously identified alleles of "Green Revolution" genes and can guide allele mining and dissect the trajectory of beneficial alleles in wheat breeding. Our work presents a model-based framework for precisely evaluating germplasm resources with genomic data. A database, WheatCompDB (http://wheat.cau.edu.cn/WheatCompDB/), is available for researchers to exploit the identified gIBDs with a multi-scale network.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems

Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.

Long-range assembly of sequences helps to unravel the genome structure and small variation of the wheat-Haynaldia villosa translocated chromosome 6VS.6AL

Genomics studies in wild species of wheat have been limited due to the lack of references; however, new technologies and bioinformatics tools have much potential to promote genomic research. The wheat-Haynaldia villosa translocation line T6VS·6AL has been widely used as a backbone parent of wheat breeding in China. Therefore, revealing the genome structure of translocation chromosome 6VS·6AL will clarify how this chromosome formed and will help to determine how it affects agronomic traits. In this study, chromosome flow sorting, NGS sequencing and Chicago long-range linkage assembly were innovatively used to produce the assembled sequences of 6VS·6AL, and gene prediction and genome structure characterization at the molecular level were effectively performed. The analysis discovered that the short arm of 6VS·6AL was actually composed of a large distal segment of 6VS, a small proximal segment of 6AS and the centromere of 6A, while the collinear region in 6VS corresponding to 230-260 Mb of 6AS-Ta was deleted when the recombination between 6VS and 6AS occurred. In addition to the molecular mechanism of the increased grain weight and enhanced spike length produced by the translocation chromosome, it may be correlated with missing GW2-V and an evolved NRT-V cluster. Moreover, a fine physical bin map of 6VS was constructed by the high-throughput developed 6VS-specific InDel markers and a series of newly identified small fragment translocation lines involving 6VS. This study will provide essential information for mining of new alien genes carried by the 6VS·6AL translocation chromosome.

Long-read sequence assembly: a technical evaluation in barley

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.

Introgressing the Aegilops tauschii genome into wheat as a basis for cereal improvement

Increasing crop production is necessary to feed the world's expanding population, and crop breeders often utilize genetic variations to improve crop yield and quality. However, the narrow diversity of the wheat D genome seriously restricts its selective breeding. A practical solution is to exploit the genomic variations of Aegilops tauschii via introgression. Here, we established a rapid introgression platform for transferring the overall genetic variations of A. tauschii to elite wheats, thereby enriching the wheat germplasm pool. To accelerate the process, we assembled four new reference genomes, resequenced 278 accessions of A. tauschii and constructed the variation landscape of this wheat progenitor species. Genome comparisons highlighted diverse functional genes or novel haplotypes with potential applications in wheat improvement. We constructed the core germplasm of A. tauschii, including 85 accessions covering more than 99% of the species' overall genetic variations. This was crossed with elite wheat cultivars to generate an A. tauschii-wheat synthetic octoploid wheat (A-WSOW) pool. Laboratory and field analysis with two examples of the introgression lines confirmed its great potential for wheat breeding. Our high-quality reference genomes, genomic variation landscape of A. tauschii and the A-WSOW pool provide valuable resources to facilitate gene discovery and breeding in wheat.

A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum

Identifying the molecular basis of resistance to pathogens is critical to promote a chemical-free cropping system. In plants, nucleotide-binding leucine-rich repeat constitute the largest family of disease resistance (R) genes, but this resistance can be rapidly overcome by the pathogen, prompting research into alternative sources of resistance. Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the most important diseases of common bean. This study aimed to identify the molecular basis of Co-x, an anthracnose R gene conferring total resistance to the extremely virulent C. lindemuthianum strain 100. To that end, we sequenced the Co-x 58 kb target region in the resistant JaloEEP558 (Co-x) common bean and identified KTR2/3, an additional gene encoding a truncated and chimeric CRINKLY4 kinase, located within a CRINKLY4 kinase cluster. The presence of KTR2/3 is strictly correlated with resistance to strain 100 in a diversity panel of common beans. Furthermore, KTR2/3 expression is up-regulated 24 hours post-inoculation and its transient expression in a susceptible genotype increases resistance to strain 100. Our results provide evidence that Co-x encodes a truncated and chimeric CRINKLY4 kinase probably resulting from an unequal recombination event that occurred recently in the Andean domesticated gene pool. This atypical R gene may act as a decoy involved in indirect recognition of a fungal effector.

DNA Modification Patterns within the Transposable Elements of the Fig (Ficus carica L.) Genome

Transposable element activity can be harmful to the host’s genome integrity, but it can also provide selective advantages. One strategy to cope with transposons is epigenetic control through DNA base modifications. We report the non-canonic DNA modification dynamics of fig (Ficus carica L.) by exploiting high-quality genome reference and related N4-methylcytosine (4mC) and N6-methyladenine (6mA) data. Overall, 1.49% of transposon nucleotides showed either 4mC or 6mA modifications: the 4mC/6mA ratio was similar in Class I and Class II transposons, with a prevalence of 4mC, which is comparable to coding genes. Different percentages of 4mC or 6mA were observed among LTR-retrotransposon lineages and sub-lineages. Furthermore, both the Copia and Gypsy retroelements showed higher modification rates in the LTR and coding regions compared with their neighbour regions. Finally, the unconventional methylation of retrotransposons is unrelated to the number of close genes, suggesting that the 4mC and 6mA frequency in LTR-retrotransposons should not be related to transcriptional repression in the adjacency of the element. In conclusion, this study highlighted unconventional DNA modification patterns in fig transposable elements. Further investigations will focus on functional implications, in regards to how modified retroelements affect the expression of neighbouring genes, and whether these epigenetic markers can spread from repeats to genes, shaping the plant phenotype.

Chromosome genomics uncovers plant genome organization and function

The identification of causal genomic loci and their interactions underlying various traits in plants has been greatly aided by progress in understanding the organization of the nuclear genome. This provides clues to the responses of plants to environmental stimuli at the molecular level. Apart from other uses, these insights are needed to fully explore the potential of new breeding techniques that rely on genome editing. However, genome analysis and sequencing is not straightforward in the many agricultural crops and their wild relatives that possess large and complex genomes. Chromosome genomics streamlines this task by dissecting the genome to single chromosomes whose DNA is then used instead of nuclear DNA. This results in a massive and lossless reduction in DNA sample complexity, reduces the time and cost of the experiment, and simplifies data interpretation. Flow cytometric sorting of condensed mitotic chromosomes makes it possible to purify single chromosomes in large quantities, and as the DNA remains intact this process can be coupled successfully with many techniques in molecular biology and genomics. Since the first experiments with flow cytometric sorting in the late 1980s, numerous applications have been developed, and chromosome genomics has been having a significant impact in many areas of research, including the sequencing of complex genomes of important crops and gene cloning. This review discusses these applications, describes their contribution to advancements in plant genome analysis and gene cloning, and outlines future directions.

Evidence for the Accumulation of Nonsynonymous Mutations and Favorable Pleiotropic Alleles During Wheat Breeding

Plant breeding leads to the genetic improvement of target traits by selecting a small number of genotypes from among typically large numbers of candidate genotypes after careful evaluation. In this study, we first investigated how mutations at conserved nucleotide sites normally viewed as deleterious, such as nonsynonymous sites, accumulated in a wheat, Triticum aestivum, breeding lineage. By comparing a 150 year old ancestral and modern cultivar, we found recent nucleotide polymorphisms altered amino acids and occurred within conserved genes at frequencies expected in the absence of purifying selection. Mutations that are deleterious in other contexts likely had very small or no effects on target traits within the breeding lineage. Second, we investigated if breeders selected alleles with favorable effects on some traits and unfavorable effects on others and used different alleles to compensate for the latter. An analysis of a segregating population derived from the ancestral and modern parents provided one example of this phenomenon. The recent cultivar contains the Rht-B1b green revolution semi-dwarfing allele and compensatory alleles that reduce its negative effects. However, improvements in traits other than plant height were due to pleiotropic loci with favorable effects on traits and to favorable loci with no detectable pleiotropic effects. Wheat breeding appears to tolerate mutations at conserved nucleotide sites and to only select for alleles with both favorable and unfavorable effects on traits in exceptional situations.

Impact of Transposable Elements on Methylation and Gene Expression across Natural Accessions of Brachypodium distachyon

Transposable elements (TEs) constitute a large fraction of plant genomes and are mostly present in a transcriptionally silent state through repressive epigenetic modifications, such as DNA methylation. TE silencing is believed to influence the regulation of adjacent genes, possibly as DNA methylation spreads away from the TE. Whether this is a general principle or a context-dependent phenomenon is still under debate, pressing for studying the relationship between TEs, DNA methylation, and nearby gene expression in additional plant species. Here, we used the grass Brachypodium distachyon as a model and produced DNA methylation and transcriptome profiles for 11 natural accessions. In contrast to what is observed in Arabidopsis thaliana, we found that TEs have a limited impact on methylation spreading and that only few TE families are associated with a low expression of their adjacent genes. Interestingly, we found that a subset of TE insertion polymorphisms is associated with differential gene expression across accessions. Thus, although not having a global impact on gene expression, distinct TE insertions may contribute to specific gene expression patterns in B. distachyon.

GenoTypeMapper: graphical genotyping on genetic and sequence-based maps

BACKGROUND

The rising availability of assemblies of large genomes (e.g. bread and durum wheat, barley) and their annotations deliver the basis to graphically present genome organization of parents and progenies on a physical scale. Genetic maps are a very important tool for breeders but often represent distorted models of the actual chromosomes, e.g., in centromeric and telomeric regions. This biased picture might lead to imprecise assumptions and estimations about the size and complexity of genetic regions and the selection of suitable molecular markers for the incorporation of traits in breeding populations or near-isogenic lines (NILs). Some software packages allow the graphical illustration of genotypic data, but to the best of our knowledge, suitable software packages that allow the comparison of genotypic data on the physical and genetic scale are currently unavailable.

RESULTS

We developed a simple Java-based-software called GenoTypeMapper (GTM) for comparing genotypic data on genetic and physical maps and tested it for effectiveness on data of two NILs that carry QTL-regions for drought stress tolerance from wild emmer on chromosome 2BS and 7AS. Both NILs were more tolerant to drought stress than their recurrent parents but exhibited additional undesirable traits such as delayed heading time.

CONCLUSIONS

In this article, we illustrate that the software easily allows users to display and identify additional chromosomal introgressions in both NILs originating from the wild emmer parent. The ability to detect and diminish linkage drag can be of particular interest for pre-breeding purposes and the developed software is a well-suited tool in this respect. The software is based on a simple allele-matching algorithm between the offspring and parents of a crossing scheme. Despite this simple approach, GTM seems to be the only software that allows us to analyse, illustrate and compare genotypic data of offspring of different crossing schemes with up to four parents in two different maps. So far, up to 500 individuals with a maximum number of 50,000 markers can be examined with the software. The main limitation that hampers the performance of the software is the number of markers that are examined in parallel. Since each individual must be analysed separately, a maximum of ten individuals can currently be displayed in a single run. On a computer with an Intel five processor of the 8th generation, GTM can reliably either analyse a single individual with up to 12,000 markers or ten individuals with up to 3,600 markers in less than five seconds. Future work aims to improve the performance of the software so that more complex crossing schemes with more parents and more markers can be analysed.

Structural Variations Affecting Genes and Transposable Elements of Chromosome 3B in Wheats

Structural variations (SVs) such as copy number and presence–absence variations are polymorphisms that are known to impact genome composition at the species level and are associated with phenotypic variations. In the absence of a reference genome sequence, their study has long been hampered in wheat. The recent production of new wheat genomic resources has led to a paradigm shift, making possible to investigate the extent of SVs among cultivated and wild accessions. We assessed SVs affecting genes and transposable elements (TEs) in a Triticeae diversity panel of 45 accessions from seven tetraploid and hexaploid species using high-coverage shotgun sequencing of sorted chromosome 3B DNA and dedicated bioinformatics approaches. We showed that 23% of the genes are variable within this panel, and we also identified 330 genes absent from the reference accession Chinese Spring. In addition, 60% of the TE-derived reference markers were absent in at least one accession, revealing a high level of intraspecific and interspecific variability affecting the TE space. Chromosome extremities are the regions where we observed most of the variability, confirming previous hypotheses made when comparing wheat with the other grasses. This study provides deeper insights into the genomic variability affecting the complex Triticeae genomes at the intraspecific and interspecific levels and suggests a phylogeny with independent hybridization events leading to different hexaploid species.

Longin CF, Ganal M, Himmelbach A, Reif JC, Schulthess AW. Exome association analysis sheds light onto leaf rust (Puccinia triticina) resistance genes currently used in wheat breeding (Triticum aestivum L.)

Resistance breeding is crucial for a sustainable control of leaf rust (Puccinia triticina) in wheat (Triticum aestivum L.) while directly targeting functional variants is the Holy Grail for efficient marker-assisted selection and map-based cloning. We assessed the limits and prospects of exome association analysis for severity of leaf rust in a large hybrid wheat population of 1574 single-crosses plus their 133 parents. After imputation and quality control, exome sequencing revealed 202 875 single-nucleotide polymorphisms (SNPs) covering 19.7% of the high-confidence annotated gene space. We performed intensive data mining and found significant associations for 2171 SNPs corresponding to 50 different loci. Some of these associations mapped in the proximity of the already known resistance genes Lr21, Lr34-B, Lr1 and Lr10, while other associated genomic regions, such as those on chromosomes 1A and 3D, harboured several annotated genes putatively involved in resistance. Validation with an independent population helped to narrow down the list of putative resistance genes that should be targeted by fine-mapping. We expect that the proposed strategy of intensive data mining coupled with validation will significantly influence research in plant genetics and breeding.

Allele mining of TaGRF-2D gene 5'-UTR in Triticum aestivum and Aegilops tauschii genotypes

The low diversity of the D-subgenome of bread wheat requires the involvement of new alleles for breeding. In grasses, the allelic state of Growth Regulating Factor (GRF) gene is correlated with nitrogen uptake. In this study, we characterized the sequence of TaGRF-2D and assessed its diversity in bread wheat and goatgrass Aegilops tauschii (genome DD). In silico analysis was performed for reference sequence searching, primer pairs design and sequence assembly. The gene sequence was obtained using Illumina and Sanger sequencing. The complete sequences of TaGRF-2D were obtained for 18 varieties of wheat. The polymorphism in the presence/absence of two GCAGCC repeats in 5' UTR was revealed and the GRF-2D-SSR marker was developed. Our results showed that the alleles 5' UTR-250 and 5' UTR-238 were present in wheat varieties, 5' UTR-250 was presented in the majority of wheat varieties. In Ae. tauschii ssp. strangulata (likely donor of the D-subgenome of polyploid wheat), most accessions carried the 5' UTR-250 allele, whilst most Ae. tauschii ssp. tauschii have 5' UTR-244. The developed GRF-2D-SSR marker can be used to study the genetic diversity of wheat and Ae. tauschii.

Identification and characterization of large-scale genomic rearrangements during wheat evolution

Following allopolyploidization, nascent polyploid wheat species react with massive genomic rearrangements, including deletion of transposable element-containing sequences. While such massive rearrangements are considered to be a prominent process in wheat genome evolution and speciation, their structure, extent, and underlying mechanisms remain poorly understood. In this study, we retrieved ~3500 insertions of a specific variant of Fatima, one of the most dynamic gypsy long-terminal repeat retrotransposons in wheat from the recently available high-quality genome drafts of Triticum aestivum (bread wheat) and Triticum turgidum ssp. dicoccoides or wild emmer, the allotetraploid mother of all modern wheats. The dynamic nature of Fatima facilitated the identification of large (i.e., up to ~ 1 million bases) Fatima-containing insertions/deletions (indels) upon comparison of bread wheat and wild emmer genomes. We characterized 11 such indels using computer-assisted analysis followed by PCR validation, and found that they might have occurred via unequal intra-strand recombination or double-strand break (DSB) events. Additionally, we observed one case of introgression of novel DNA fragments from an unknown source into the wheat genome. Our data thus indicate that massive large-scale DNA rearrangements might play a prominent role in wheat speciation.

Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis

The diploid wild cotton species Gossypium australe possesses excellent traits including resistance to disease and delayed gland morphogenesis, and has been successfully used for distant breeding programmes to incorporate disease resistance traits into domesticated cotton. Here, we sequenced the G. australe genome by integrating PacBio, Illumina short read, BioNano (DLS) and Hi-C technologies, and acquired a high-quality reference genome with a contig N50 of 1.83 Mb and a scaffold N50 of 143.60 Mb. We found that 73.5% of the G. australe genome is composed of various repeat sequences, differing from those of G. arboreum (85.39%), G. hirsutum (69.86%) and G. barbadense (69.83%). The G. australe genome showed closer collinear relationships with the genome of G. arboreum than G. raimondii and has undergone less extensive genome reorganization than the G. arboreum genome. Selection signature and transcriptomics analyses implicated multiple genes in disease resistance responses, including GauCCD7 and GauCBP1, and experiments revealed induction of both genes by Verticillium dahliae and by the plant hormones strigolactone (GR24), salicylic acid (SA) and methyl jasmonate (MeJA). Experiments using a Verticillium-resistant domesticated G. barbadense cultivar confirmed that knockdown of the homologues of these genes caused a significant reduction in resistance against Verticillium dahliae. Moreover, knockdown of a newly identified gland-associated gene GauGRAS1 caused a glandless phenotype in partial tissues using G. australe. The G. australe genome represents a valuable resource for cotton research and distant relative breeding as well as for understanding the evolutionary history of crop genomes.

Variation in phosphorus and sulfur content shapes the genetic architecture and phenotypic associations within the wheat grain ionome

Dissection of the genetic basis of wheat ionome is crucial for understanding the physiological and biochemical processes underlying mineral accumulation in seeds, as well as for efficient crop breeding. Most of the elements essential for plants are metals stored in seeds as chelate complexes with phytic acid or sulfur-containing compounds. We assume that the involvement of phosphorus and sulfur in metal chelation is the reason for strong phenotypic correlations within ionome. Adjustment of element concentrations for the effect of variation in phosphorus and sulfur seed content resulted in drastic change of phenotypic correlations between the elements. The genetic architecture of wheat grain ionome was characterized by quantitative trait loci (QTL) analysis using a cross between durum and wild emmer wheat. QTL analysis of the adjusted traits and two-trait analysis of the initial traits paired with either P or S considerably improved QTL detection power and accuracy, resulting in the identification of 105 QTLs and 617 QTL effects for 11 elements. Candidate gene search revealed some potential functional associations between QTLs and corresponding genes within their intervals. Thus, we have shown that accounting for variation in P and S is crucial for understanding of the physiological and genetic regulation of mineral composition of wheat grain ionome and can be implemented for other plants.

Genome of an iconic Australian bird: High-quality assembly and linkage map of the superb fairy-wren (Malurus cyaneus)

The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.

Where the Wild Things Are: Transposable Elements as Drivers of Structural and Functional Variations in the Wheat Genome

Transposable elements (TEs) are major contributors to genome plasticity and thus are likely to have a dramatic impact on genetic diversity and speciation. Recent technological developments facilitated the sequencing and assembly of the wheat genome, opening the gate for whole genome analysis of TEs in wheat, which occupy over 80% of the genome. Questions that have been long unanswered regarding TE dynamics throughout the evolution of wheat, are now being addressed more easily, while new questions are rising. In this review, we discuss recent advances in the field of TE dynamics in wheat and possible future directions.

Tools and Strategies for Long-Read Sequencing and De Novo Assembly of Plant Genomes

The commercial release of third-generation sequencing technologies (TGSTs), giving long and ultra-long sequencing reads, has stimulated the development of new tools for assembling highly contiguous genome sequences with unprecedented accuracy across complex repeat regions. We survey here a wide range of emerging sequencing platforms and analytical tools for de novo assembly, provide background information for each of their steps, and discuss the spectrum of available options. Our decision tree recommends workflows for the generation of a high-quality genome assembly when used in combination with the specific needs and resources of a project.

Frequent intra- and inter-species introgression shapes the landscape of genetic variation in bread wheat

BACKGROUND

Bread wheat is one of the most important and broadly studied crops. However, due to the complexity of its genome and incomplete genome collection of wild populations, the bread wheat genome landscape and domestication history remain elusive.

RESULTS

By investigating the whole-genome resequencing data of 93 accessions from worldwide populations of bread wheat and its diploid and tetraploid progenitors, together with 90 published exome-capture data, we find that the B subgenome has more variations than A and D subgenomes, including SNPs and deletions. Population genetics analyses support a monophyletic origin of domesticated wheat from wild emmer in northern Levant, with substantial introgressed genomic fragments from southern Levant. Southern Levant contributes more than 676 Mb in AB subgenomes and enriched in the pericentromeric regions. The AB subgenome introgression happens at the early stage of wheat speciation and partially contributes to their greater genetic diversity. Furthermore, we detect massive alien introgressions that originated from distant species through natural and artificial hybridizations, resulting in the reintroduction of ~ 709 Mb and ~ 1577 Mb sequences into bread wheat landraces and varieties, respectively. A large fraction of these intra- and inter-introgression fragments are associated with quantitative trait loci of important traits, and selection events are also identified.

CONCLUSION

We reveal the significance of multiple introgressions from distant wild populations and alien species in shaping the genetic components of bread wheat, and provide important resources and new perspectives for future wheat breeding.

The Coiled-Coil NLR Rph1, Confers Leaf Rust Resistance in Barley Cultivar Sudan

Unraveling and exploiting mechanisms of disease resistance in cereal crops is currently limited by their large repeat-rich genomes and the lack of genetic recombination or cultivar (cv)-specific sequence information. We cloned the first leaf rust resistance gene Rph1 (Rph1 a) from cultivated barley (Hordeum vulgare) using "MutChromSeq," a recently developed molecular genomics tool for the rapid cloning of genes in plants. Marker-trait association in the CI 9214/Stirling doubled haploid population mapped Rph1 to the short arm of chromosome 2H in a physical region of 1.3 megabases relative to the barley cv Morex reference assembly. A sodium azide mutant population in cv Sudan was generated and 10 mutants were confirmed by progeny-testing. Flow-sorted 2H chromosomes from Sudan (wild type) and six of the mutants were sequenced and compared to identify candidate genes for the Rph1 locus. MutChromSeq identified a single gene candidate encoding a coiled-coil nucleotide binding site Leucine-rich repeat (NLR) receptor protein that was altered in three different mutants. Further Sanger sequencing confirmed all three mutations and identified an additional two independent mutations within the same candidate gene. Phylogenetic analysis determined that Rph1 clustered separately from all previously cloned NLRs from the Triticeae and displayed highest sequence similarity (89%) with a homolog of the Arabidopsis (Arabidopsis thaliana) disease resistance protein 1 protein in Triticum urartu In this study we determined the molecular basis for Rph1-mediated resistance in cultivated barley enabling varietal improvement through diagnostic marker design, gene editing, and gene stacking technologies.

Contribution of recent technological advances to future resistance breeding

The development of durable host resistance strategies to control crop diseases is a primary need for sustainable agricultural production in the future. This article highlights the potential of recent progress in the understanding of host resistance for future cereal breeding. Much of the novel work is based on advancements in large-scale sequencing and genomics, rapid gene isolation techniques and high-throughput molecular marker technologies. Moreover, emerging applications on the pathogen side like effector identification or field pathogenomics are discussed. The combination of knowledge from both sides of cereal pathosystems will result in new approaches for resistance breeding. We describe future applications and innovative strategies to implement effective and durable strategies to combat diseases of major cereal crops while reducing pesticide dependency.

Detecting Large Chromosomal Modifications Using Short Read Data From Genotyping-by-Sequencing

Markers linked to agronomic traits are of the prerequisite for molecular breeding. Genotyping-by-sequencing (GBS) data enables to detect small polymorphisms including single nucleotide polymorphisms (SNPs) and short insertions or deletions (InDels) that can be used, for instance, for marker-assisted selection, population genetics, and genome-wide association studies (GWAS). Here, we aim at detecting large chromosomal modifications in barley and wheat based on GBS data. These modifications could be duplications, deletions, substitutions including introgressions as well as alterations of DNA methylation. We demonstrate that GBS coverage analysis is capable to detect Hordeum vulgare/Hordeum bulbosum introgression lines. Furthermore, we identify large chromosomal modifications in barley and wheat collections. Hence, large chromosomal modifications, including introgressions and copy number variations (CNV), can be detected easily and can be used as markers in research and breeding without additional wet-lab experiments.