Megabase level sequencing reveals contrasted organization and evolution patterns of the wheat gene and transposable element spaces

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

An integrated strategy involving high-throughput sequencing to characterize an unknown GM wheat event in Canada

Glyphosate-resistant wheat plants were discovered in southern Alberta in 2017, representing an unauthorized GM release in Canada. The Canadian Food Inspection Agency undertook a series of experiments to characterize and identify this unknown GM wheat, as well as to develop and validate construct-specific and event-specific qPCR assays. Results of PCR-based assays and Sanger sequencing indicated the presence of CaMV 35S promoter (p35S), Rice Actin 1 intron (RactInt1), CP4-EPSPS gene and nopaline synthase terminator (tNOS) elements in the unknown GM wheat. Genome walking and bead capture strategies, combined with high-throughput sequencing, were used to identify the 5' and 3' wheat junctions and the subsequent mapping of the insert to chromosome 3B of the wheat genome. A probable transformation vector, pMON25497, was recognized, and further testing identified the unknown GM wheat as MON71200 event, one of two events obtained with the pMON25497 vector. The two construct-specific assays targeted the junctions of the RactInt1 and the CP4-EPSPS elements and the CP4-EPSPS and tNOS elements, while the event-specific assay was located at the 3' junction into the wheat genome. Both construct-specific and event-specific assays had limits of detection of 0.10% of MON71200 in a seed pool. As expected, the two construct-specific assays cross-reacted with other wheat and corn events containing the same elements in the same order. No cross-reactivity was observed for the event-specific assay. The integrated strategy employed in this study can serve as a model for other cases when facing similar challenges involving unknown GM events.

Identification and Expressional Analysis of siRNAs Responsive to Fusarium graminearum Infection in Wheat

The outbreak of Fusarium head blight (FHB) poses a serious threat to wheat production as it leads to both significant yield losses and accumulation of several mycotoxins including deoxynivalenol (DON) in the grains, which are harmful to human and livestock. To date, hundreds of FHB-resistance-related quantitative trait loci (QTLs) have been reported, but only a few of them have been cloned and used for breeding. Small interfering RNAs (siRNA) have been reported in plants to mediate host defense against pathogens, but they have rarely been reported in wheat-FHB interaction. In order to identify the key siRNAs that can potentially be used in the improvement of resistance to FHB, siRNAs from the spikes of an FHB-resistant variety Sumai 3 and an FHB-susceptible variety of Chinese Spring (CS) were sequenced after F. graminearum infection and mock inoculation, respectively. The expression patterns of the siRNAs of interest were analyzed. A total of 4019 siRNAs of high-confidence were identified, with 131 being CS-specific, 309 Sumai 3-specific and 3071 being common in both varieties. More than 87% of these siRNAs were 24 nt in length. An overall down-regulation trend was found for siRNAs in the spikes of both varieties after being infected with F. graminearum. The expression patterns for Triticum aestivum Dicer-like 3 (TaDCL3) that synthesizes 24 nt siRNAs were validated by qRT-PCR, which were positively correlated with those of the siRNAs. A total of 85% of the differentially expressed genes putatively targeted by the siRNAs were significantly up-regulated after infection, showing a negative correlation with the overall down-regulated expression of siRNAs. Interestingly, the majority of the up-regulated genes are annotated as disease resistance. These results suggested that the inhibition of siRNA by F. graminearum up-regulated the disease resistance genes, which were putatively suppressed by siRNAs through RNA-directed DNA methylation (RdDM). Consequently, the resistant capability to F. graminearum infection was enhanced. This study provides novel clues for investigating the function of siRNA in wheat-F. graminearum interaction.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15 Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T. aestivum cultivars (6x = AABBDD) and a single genome for Triticum durum (4x = AABB), Triticum dicoccoides (4x = AABB), Triticum urartu (2x = AA), and Aegilops tauschii (2x = DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

A Reappraisal of Polyploidy Events in Grasses (Poaceae) in a Rapidly Changing World

Around 80% of megaflora species became extinct at the Cretaceous–Paleogene (K–Pg) boundary. Subsequent polyploidy events drove the survival of thousands of plant species and played a significant historical role in the development of the most successful modern cereal crops. However, current and rapid global temperature change poses an urgent threat to food crops worldwide, including the world’s big three cereals: rice, wheat, and maize, which are members of the grass family, Poaceae. Some minor cereals from the same family (such as teff) have grown in popularity in recent years, but there are important knowledge gaps regarding the similarities and differences between major and minor crops, including how polyploidy affects their biological processes under natural and (a)biotic stress conditions and thus the potential to harness polyploidization attributes for improving crop climate resilience. This review focuses on the impact of polyploidy events on the Poaceae family, which includes the world’s most important food sources, and discusses the past, present, and future of polyploidy research for major and minor crops. The increasing accessibility to genomes of grasses and their wild progenitors together with new tools and interdisciplinary research on polyploidy can support crop improvement for global food security in the face of climate change.

tRNA-derived fragments from wheat are potentially involved in susceptibility to Fusarium head blight

BACKGROUND

Fusarium head blight (FHB) caused by Fusarium graminearum is a devastating fungal disease of wheat. The mechanism underlying F. graminearum-wheat interaction remains largely unknown. tRNA-derived fragments (tRFs) are RNase-dependent small RNAs derived from tRNAs, and they have not been reported in wheat yet, and whether tRFs are involved in wheat-F. graminearum interactions remains unknown.

RESULTS

Herein, small RNAs from the spikelets inoculated with F. graminearum and mock from an FHB-susceptible variety Chinese Spring (CS) and an FHB-resistant variety Sumai3 (SM) were sequenced respectively. A total of 1249 putative tRFs were identified, in which 15 tRFs was CS-specific and 12 SM-specific. Compared with mock inoculation, 39 tRFs were significantly up-regulated across both wheat varieties after F. graminearum challenge and only nine tRFs were significantly down-regulated. tRFGlu, tRFLys and tRFThr were dramatically induced by F. graminearum infection, with significantly higher fold changes in CS than those in SM. The expression patterns of the three highly induced tRFs were further validated by stem-loop qRT-PCR. The accumulation of tRFs were closely related to ribonucleases T2 family members, which were induced by F. graminearum challenge. The tRFs' targets in host were predicted and were validated by RNA sequencing.

CONCLUSION

Integrative analysis of the differentially expressed tRFs and their candidate targets indicated that tRFGlu, tRFLys and tRFThr might negatively regulate wheat resistance to FHB. Our results unvealed the potential roles of tRFs in wheat-F. graminearum interactions.

Characterization and Use in Wheat Breeding of Leaf Rust Resistance Genes from Durable Varieties

Simple Summary

Wheat leaf rust is one of the most significant diseases worldwide, incited by a parasitic fungus which infects leaves, affecting grain yield. This pathogen is spread by the wind over large areas through microscopic spores. This huge number of spores favors the selection of virulent forms; therefore, there is a continuous need for new resistance genes to control this disease without fungicides. These resistant genes are naturally found in resistant wheat varieties and can be introduced by standard crosses. In this work, seven resistant genes were introduced into several commercial susceptible varieties. The selection of resistance genes was assisted by DNA markers that are close to these genes on the chromosome. Additionally, the selection of desirable traits from the commercial variety was also assisted by DNA markers to accelerate the process. In field testing, the varieties developed here were resistant to leaf rust, and suitable for commercial use.

Abstract

Leaf rust is one of the most significant diseases of wheat worldwide. In Argentina, it is one of the main reasons for variety replacement that becomes susceptible after large-scale use. Some varieties showed durable resistance to this disease, including Buck Manantial and Sinvalocho MA. RILs (Recombinant Inbred Lines) were developed for each of these varieties and used in genetics studies to identify components of resistance, both in greenhouse inoculations using leaf rust races, and in field evaluations under natural population infections. In Buck Manantial, the APR gene LrBMP1 was associated with resistance in field tests. In crosses involving Sinvalocho MA, four genes were previously identified and associated with resistance in field testing: APR (Adult Plant Resistance) gene LrSV1, the APR genetic system LrSV2 + LrcSV2 and the ASR (All Stage Resistance) gene LrG6. Using backcrosses, LrBMP1 was introgressed in four commercial susceptible varieties and LrSV1, LrSV2 + LrcSV2 and LrG6 were simultaneously introgressed in three susceptible commercial varieties. The use of molecular markers for recurrent parent background selection allowed us to select resistant lines with more than 80% similarity to commercial varieties. Additionally, progress towards positional cloning of the genetic system LrSV2 + LrcSV2 for leaf rust APR is reported.

TaNAC032 transcription factor regulates lignin-biosynthetic genes to combat Fusarium head blight in wheat

Fusarium head blight (FHB) is a destructive disease affecting cereal crops globally due to mycotoxin contamination of grains that reduce yield and quality. Among hundreds of QTLs identified for resistance, the QTL-Fhb1 is of significant interest even today, for its major contribution to FHB resistance. Previously, QTL-Fhb1 dissection based on a combined metabolo-genomics approach, identified a few potential resistance genes, including a NAC like transcription factor for FHB resistance. Sequencing and phylogenetic analysis confirmed NAC to be the wheat TaNAC032. Also, the quantitative RT-PCR studies revealed a greater induced expression of TaNAC032 in resistant NIL in comparison to susceptible NIL upon Fusarium graminearum (Fg) infection. The virus-induced gene silencing (VIGS) based functional validation of TaNAC032 in resistant NIL confirmed increased disease severity and fungal biomass. Metabolic profiling revealed low abundances of resistance-related (RR) metabolites in TaNAC032 silenced NIL-R compared to non-silenced. Silenced plants showed decreased transcript abundances of RR metabolite biosynthetic genes associated with a reduction in total lignin content in rachis, confirming the regulatory role of TaNAC032 in wheat in response to Fg infection. If TaNA032 is mutated in an FHB susceptible cultivar, it can be edited to enhance FHB resistance.

Rhizospheric bacterial community structure of Triticum and Aegilops revealed by pyrosequencing analysis of the 16S rRNA gene: dominance of the A genome over the B and D genomes

This study examined the relationship between host plant and rhizospheric bacterial community structure, including composition and diversity, in Triticum and Aegilops species (12 and two accessions, respectively) as well as three closely related species, barley, rye and oat (four accessions), to explore the possibility that wheat root and rhizosphere interaction can be utilized for wheat breeding and biotechnology in the future. For this purpose, DNA was isolated from rhizospheric soil samples and one control non-rhizospheric soil sample, and the 16S rRNA gene region was amplified and subjected to DNA pyrosequencing. A total of 132,888 amplicons were analyzed. Bacterial composition at the phylum level was similar among the 18 rhizospheric samples; however, the proportion of Acidobacteria was much lower in these samples than in the control non-rhizospheric soil sample, indicating that rhizospheres influenced the bacterial composition even at the higher taxonomic level. Across host plant genome types (three levels of ploidy and three major genomes, A, B and D), there was no detectable difference in phylum composition or species diversity. Estimated bacterial species diversity was higher in the control soil sample than in plant rhizospheric soils, implying that bacterial species diversity was reduced in rhizospheres. A PCoA plot and UPGMA dendrogram based on the bacterial species composition showed that control soil was distantly located from the plant rhizospheric samples and that Triticum, Aegilops and related species were well separated. PERMANOVA analysis detected statistically significant differentiation among these four groups. Clustering of Triticum species suggested that the A genome was dominant over the B and D genomes, with respect to the influence on rhizospheric bacterial species composition. Although the cause was not investigated in this study, these results clearly indicated that the genetic constitution of the plant host exerted a strong influence on rhizospheric bacterial community structure.

Differential chromatin accessibility landscape reveals structural and functional features of the allopolyploid wheat chromosomes

Background

Our understanding of how the complexity of the wheat genome influences the distribution of chromatin states along the homoeologous chromosomes is limited. Using a differential nuclease sensitivity assay, we investigate the chromatin states of the coding and repetitive regions of the allopolyploid wheat genome.

Results

Although open chromatin is found to be significantly enriched around genes, the majority of MNase-sensitive regions are located within transposable elements (TEs). Chromatin of the smaller D genome is more accessible than that of the larger A and B genomes. Chromatin states of different TEs vary among families and are influenced by the TEs’ chromosomal position and proximity to genes. While the chromatin accessibility of genes is influenced by proximity to TEs, and not by their position on the chromosomes, we observe a negative chromatin accessibility gradient along the telomere-centromere axis in the intergenic regions, positively correlated with the distance between genes. Both gene expression levels and homoeologous gene expression bias are correlated with chromatin accessibility in promoter regions. The differential nuclease sensitivity assay accurately predicts previously detected centromere locations. SNPs located within more accessible chromatin explain a higher proportion of genetic variance for a number of agronomic traits than SNPs located within more closed chromatin.

Conclusions

Chromatin states in the wheat genome are shaped by the interplay of repetitive and gene-encoding regions that are predictive of the functional and structural organization of chromosomes, providing a powerful framework for detecting genomic features involved in gene regulation and prioritizing genomic variation to explain phenotypes.

The mechanisms and significance of the positional control of centromeres and telomeres in plants

The centromere and telomere are universal heterochromatic domains; however, the proper positioning of those domains in nuclear space during the mitotic interphase differs among eukaryotes. Consequently, the question arises how and why this difference occurs. Studies over the past 2 decades have identified several nuclear membrane proteins, nucleolar proteins, and the structural maintenance of a chromosome complex as factors involved in the positional control of centromeres and/or telomeres during the mitotic interphase in yeasts, animals, and plants. In this review, with a primary focus on plants, the roles of those factors are summarized, and the biological significance of proper centromere and telomere positionings during the mitotic interphase is discussed in an effort to provide guidance for this question.

Fine mapping and gene cloning in the post-NGS era: advances and prospects

Improvement in traits of agronomic importance is the top breeding priority of crop improvement programs. Majority of these agronomic traits show complex quantitative inheritance. Identification of quantitative trait loci (QTLs) followed by fine mapping QTLs and cloning of candidate genes/QTLs is central to trait analysis. Advances in genomic technologies revolutionized our understanding of genetics of complex traits, and genomic regions associated with traits were employed in marker-assisted breeding or cloning of QTLs/genes. Next-generation sequencing (NGS) technologies have enabled genome-wide methodologies for the development of ultra-high-density genetic linkage maps in different crops, thus allowing placement of candidate loci within few kbs in genomes. In this review, we compare the marker systems used for fine mapping and QTL cloning in the pre- and post-NGS era. We then discuss how different NGS platforms in combination with advanced experimental designs have improved trait analysis and fine mapping. We opine that efficient genotyping/sequencing assays may circumvent the need for cumbersome procedures that were earlier used for fine mapping. A deeper understanding of the trait architectures of agricultural significance will be crucial to accelerate crop improvement.

What Can Long Terminal Repeats Tell Us About the Age of LTR Retrotransposons, Gene Conversion and Ectopic Recombination?

LTR retrotransposons constitute a significant part of plant genomes and their evolutionary dynamics play an important role in genome size changes. Current methods of LTR retrotransposon age estimation are based only on LTR (long terminal repeat) divergence. This has prompted us to analyze sequence similarity of LTRs in 25,144 LTR retrotransposons from fifteen plant species as well as formation of solo LTRs. We found that approximately one fourth of nested retrotransposons showed a higher LTR divergence than the pre-existing retrotransposons into which they had been inserted. Moreover, LTR similarity was correlated with LTR length. We propose that gene conversion can contribute to this phenomenon. Gene conversion prediction in LTRs showed potential converted regions in 25% of LTR pairs. Gene conversion was higher in species with smaller genomes while the proportion of solo LTRs did not change with genome size in analyzed species. The negative correlation between the extent of gene conversion and the abundance of solo LTRs suggests interference between gene conversion and ectopic recombination. Since such phenomena limit the traditional methods of LTR retrotransposon age estimation, we recommend an improved approach based on the exclusion of regions affected by gene conversion.

Combining a New Exome Capture Panel With an Effective varBScore Algorithm Accelerates BSA-Based Gene Cloning in Wheat

The discovery of functional genes underlying agronomic traits is of great importance for wheat improvement. Here we designed a new wheat exome capture probe panel based on IWGSC RefSeq v1.0 genome sequence information and developed an effective algorithm, varBScore, that can sufficiently reduce the background noise in gene mapping and identification. An effective method, termed bulked segregant exome capture sequencing (BSE-Seq) for identifying causal mutations or candidate genes was established by combining the use of a newly designed wheat exome capture panel, sequencing of bulked segregant pools from segregating populations, and the robust algorithm varBScore. We evaluated the effectiveness of varBScore on SNP calling using the published dataset for mapping and cloning the yellow rust resistance gene Yr7 in wheat. Furthermore, using BSE-Seq, we rapidly identified a wheat yellow leaf mutant gene, ygl1, in an ethyl methanesulfonate (EMS) mutant population and found that a single mutation of G to A at 921 position in the wild type YGL1 gene encoding magnesium-chelatase subunit chlI caused the leaf yellowing phenotype. We further showed that mutation of YGL1 through CRISPR/Cas9 gene editing led to a yellow phenotype on the leaves of transgenic wheat, indicating that ygl1 is the correct causal gene responsible for the mutant phenotype. In summary, our approach is highly efficient for discovering causal mutations and gene cloning in wheat.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (Jb), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (JrJvsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

In-silico prediction of novel genes responsive to drought and salinity stress tolerance in bread wheat (Triticum aestivum)

Common wheat (Triticum aestivum) is the most widely grown cereal crop and is cultivated extensively in dry regions. Water shortage, resulting from either drought or salinity, leads to slow growth and loss of wheat yield. In order to predict new genes responsive to the drought and salt stresses in wheat, 6,717 expressed sequence tags (ESTs), expressed in drought and salinity stress conditions were collected from the National Center for Biotechnology Information (NCBI). The downloaded ESTs were clustered and assembled into 354 contigs; 14 transcription factor families in 29 contigs were identified. In addition, 119 contigs were organized in five enzyme classes. Biological functions were obtained for only 324 of the 354 contigs using gene ontology. In addition, using Kyoto Encyclopedia of Genes and Genomes database, 191 metabolic pathways were identified. The remaining contigs were used for further analysis and the search for new genes responsive to drought and salt stresses. These contigs were mapped on the International Wheat Genome Sequencing Consortium RefSeq v1.0 assembly, the most complete version of the reference sequence of the bread wheat variety Chinese Spring. They were found to have from one to three locations on the subgenomes A, B, and D. Full-length gene sequences were designed for these contigs, which were further validated using promoter analysis. These predicted genes may have applications in molecular breeding programs and wheat drought and salinity research.

Plasticity in Triticeae centromere DNA sequences: a wheat × tall wheatgrass (decaploid) model

Centromeres mediate chromosome attachment to microtubules and maintain the integrity of chromosomes for proper segregation of the sister chromatids during cell division. Advances in the assembly of Triticeae genome sequences combined with the capacity to recover hybrid species derived from very distantly related species provides potential experimental systems for linking retrotransposon amplification and repositioning of centromeres via non-mendelian inheritance in partial amphiploid breeds. The decaploid tall wheatgrass (Thinopyrum ponticum) is one of the most successfully used perennial species in wheat breeding for generating translocation lines with valuable agronomic traits. We found that wheat centromere retrotransposons CRW and Quinta widely occur within the tall wheatgrass genome. In addition, one of the genome donors to Th. ponticum, Pseudoroegneria stipifolia (StSt), has been shown to have Abigail and a satellite repeat, CentSt. We also found two other centromeric retrotransposons, Abia and CL135 in Th. ponticum by ChIP-seq. Examination of partial amphiploid lines that were generated in the 1970s demonstrated extensive modification in centromere sequences using CentSt, Abigail and Abia as probes. We also detected that St-genome chromosomes were more enriched with Abigail and CentSt, whereas E-genome chromosomes were enriched with CRW and Quinta in tall wheatgrass and its closer relatives. It can be concluded that bursts of transposition of retrotransposons and repositioning of centromeres via non-mendelian segregation are common in partial amphiploids derived from interspecific hybrids. Practically speaking, our study reveals that the existence of homologous centromere functional sequences in both a donor and its receptor can substantially contribute to the successful transfer of alien genes into crop species. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/sra/SRR9089557; https://www.ncbi.nlm.nih.gov/sra/SRR9089558; https://www.ncbi.nlm.nih.gov/sra/SRR9089559; https://www.ncbi.nlm.nih.gov/sra/SRR9089560; https://www.ncbi.nlm.nih.gov/sra/SRR9089561; https://www.ncbi.nlm.nih.gov/sra/SRR9089562; https://www.ncbi.nlm.nih.gov/sra/SRR9089563; https://www.ncbi.nlm.nih.gov/sra/SRR9089564; https://www.ncbi.nlm.nih.gov/nuccore/MK999394; https://www.ncbi.nlm.nih.gov/nuccore/MK999395; https://www.ncbi.nlm.nih.gov/nuccore/MK999396.

Serpins: Genome-Wide Characterisation and Expression Analysis of the Serine Protease Inhibitor Family in Triticum aestivum

The serine protease inhibitor (serpin) gene family is the largest family of protease inhibitors. Serine protease inhibitors have an active, but under-characterized, role in grain development and defense against pathogen attack in cereal crops. By exploiting publicly available genomic, transcriptomic and proteomic data for wheat (Triticum aestivum), we have identified and annotated the entire ’serpinome’ of wheat and constructed a high-quality and robust phylogenetic tree of the gene family, identifying paralogous and homeologous clades from the hexaploid wheat genome, including the Serpin-Z group that have been well characterized in barley. Using publicly available RNAseq data (http://www.wheat-expression.com/), expression profiles of the wheat serpins were explored across a variety of tissues from the developing grain, spikelet and spike. We show that the SERPIN-Z clade, among others, are highly expressed during grain development, and that there is homeologous and paralogous functional redundancy in this gene family. Further to their role in grain development, serpins play an important but under-explored role in response to fungal pathogens. Using 13 RNAseq datasets of wheat tissues infected by fungal pathogens, we identified 37 serpins with a significant disease response. The majority of the disease-responsive serpins were upregulated by Fusarium graminearum, a destructive fungal pathogen that attacks the spike and developing grain of wheat. As serpins are ubiquitous in wheat grain, the genes encoding serpins may be linked to grain development, with their disease response a result of pleiotropy.

A deletion mutation in TaHRC confers Fhb1 resistance to Fusarium head blight in wheat

Fusarium head blight (FHB), which is mainly caused by Fusarium graminearum, is a destructive wheat disease that threatens global wheat production. Fhb1, a quantitative trait locus discovered in Chinese germplasm, provides the most stable and the largest effect on FHB resistance in wheat. Here we show that TaHRC, a gene that encodes a putative histidine-rich calcium-binding protein, is the key determinant of Fhb1-mediated resistance to FHB. We demonstrate that TaHRC encodes a nuclear protein conferring FHB susceptibility and that a deletion spanning the start codon of this gene results in FHB resistance. Identical sequences of the TaHRC-R allele in diverse accessions indicate that Fhb1 had a single origin, and phylogenetic and haplotype analyses suggest that the TaHRC-R allele most likely originated from a line carrying the Dahongpao haplotype. This discovery opens a new avenue to improve FHB resistance in wheat, and possibly in other cereal crops, by manipulating TaHRC sequence through bioengineering approaches.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm

BACKGROUND

The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity.

RESULTS

We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs.

CONCLUSIONS

Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

RepetDB: a unified resource for transposable element references

Background

Thanks to their ability to move around and replicate within genomes, transposable elements (TEs) are perhaps the most important contributors to genome plasticity and evolution. Their detection and annotation are considered essential in any genome sequencing project. The number of fully sequenced genomes is rapidly increasing with improvements in high-throughput sequencing technologies. A fully automated de novo annotation process for TEs is therefore required to cope with the deluge of sequence data.

However, all automated procedures are error-prone, and an automated procedure for TE identification and classification would be no exception. It is therefore crucial to provide not only the TE reference sequences, but also evidence justifying their classification, at the scale of the whole genome. A few TE databases already exist, but none provides evidence to justify TE classification. Moreover, biological information about the sequences remains globally poor.

Results

We present here the RepetDB database developed in the framework of GnpIS, a genetic and genomic information system. RepetDB is designed to store and retrieve detected, classified and annotated TEs in a standardized manner. RepetDB is an implementation with extensions of InterMine, an open-source data warehouse framework used here to store, search, browse, analyze and compare all the data recorded for each TE reference sequence. InterMine can display diverse information for each sequence and allows simple to very complex queries. Finally, TE data are displayed via a worldwide data discovery portal. RepetDB is accessible at urgi.versailles.inra.fr/repetdb.

Conclusions

RepetDB is designed to be a TE knowledge base populated with full de novo TE annotations of complete (or near-complete) genome sequences. Indeed, the description and classification of TEs facilitates the exploration of specific TE families, superfamilies or orders across a large range of species. It also makes possible cross-species searches and comparisons of TE family content between genomes.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

Mapping a gene on wheat chromosome 4BL involved in a complementary interaction with adult plant leaf rust resistance gene LrSV2

A complementary gene to LrSV2 for specific adult plant leaf rust resistance in wheat was mapped on chromosome 4BL, tightly linked to Lr12 / 31. LrSV2 is a race-specific adult plant leaf rust (Puccinia triticina) resistance gene on subdistal chromosome 3BS detected in the cross of the traditional Argentinean wheat (Triticum aestivum) variety Sinvalocho MA and the experimental line Gama6. The analysis of the cross of R46 [recombinant inbred line (RIL) derived from Sinvalocho MA carrying LrSV2 gene and the complementary gene Lrc-SV2 identified in the current paper] and the commercial variety Relmo Siriri (not carrying neither of these two genes) allowed the detection of the unlinked complementary gene Lrc-SV2 because the presence of one dominant allele of both is necessary to express the LrSV2-specific adult plant resistance. Lrc-SV2 was mapped within a 1-cM interval on chromosome 4BL using 100 RILs from the cross Sinvalocho MA × Purple Straw. This genetic system resembles the Lr27+31 seedling resistance reported in the Australian varieties Gatcher and Timgalen where interacting genes map at similar chromosomal positions. However, in high-resolution maps, Lr27 and LrSV2 were already mapped to adjacent intervals on 3BS and Lrc-SV2 map position on 4BL is distal to the reported Lr12/31-flanking microsatellites.

Genome-wide analyses of miniature inverted-repeat transposable elements reveals new insights into the evolution of the Triticum-Aegilops group

The sequence drafts of wild emmer and bread wheat facilitated high resolution, genome-wide analysis of transposable elements (TEs), which account for up to 90% of the wheat genome. Despite extensive studies, the role of TEs in reshaping nascent polyploid genomes remains to be fully understood. In this study, we retrieved miniature inverted-repeat transposable elements (MITEs) from the recently published genome drafts of Triticum aestivum, Triticum turgidum ssp. dicoccoides, Aegilops tauschii and the available genome draft of Triticum urartu. Overall, 239,126 MITE insertions were retrieved, including 3,874 insertions of a newly identified, wheat-unique MITE family that we named "Inbar". The Stowaway superfamily accounts for ~80% of the retrieved MITE insertions, while Thalos is the most abundant family. MITE insertions are distributed in the seven homologous chromosomes of the wild emmer and bread wheat genomes. The remarkably high level of insertions in the B sub-genome (~59% of total retrieved MITE insertions in the wild emmer genome draft, and ~41% in the bread wheat genome draft), emphasize its highly repetitive nature. Nearly 52% of all MITE insertions were found within or close (less than 100bp) to coding genes, and ~400 MITE sequences were found in the bread wheat transcriptome, indicating that MITEs might have a strong impact on wheat genome expression. In addition, ~40% of MITE insertions were found within TE sequences, and remarkably, ~90% of Inbar insertions were located in retrotransposon sequences. Our data thus shed new light on the role of MITEs in the diversification of allopolyploid wheat species.

Birth and Death of LTR-Retrotransposons in Aegilops tauschii

Long terminal repeat-retrotransposons (LTR-RTs) are a major component of all flowering plant genomes. To analyze the time dynamics of LTR-RTs, we modeled the insertion rates of the 35 most abundant LTR-RT families in the genome of Aegilops tauschii, one of the progenitors of wheat. Our model of insertion rate (birth) takes into account random variation in LTR divergence and the deletion rate (death) of LTR-RTs. Modeling the death rate is crucial because ignoring it would underestimate insertion rates in the distant past. We rejected the hypothesis of constancy of insertion rates for all 35 families and showed by simulations that our hypothesis test controlled the false-positive rate. LTR-RT insertions peaked from 0.064 to 2.39 MYA across the 35 families. Among other effects, the average age of elements within a family was negatively associated with recombination rate along a chromosome, with proximity to the closest gene, and weakly associated with the proximity to its 5' end. Elements within a family that were near genes colinear with genes in the genome of tetraploid emmer wheat tended to be younger than those near noncolinear genes. We discuss these associations in the context of genome evolution and stability of genome sizes in the tribe Triticeae. We demonstrate the general utility of our models by analyzing the two most abundant LTR-RT families in Arabidopsis lyrata, and show that these families differed in their insertion dynamics. Our estimation methods are available in the R package TE on CRAN.

Impact of transposable elements on genome structure and evolution in bread wheat

BACKGROUND

Transposable elements (TEs) are major components of large plant genomes and main drivers of genome evolution. The most recent assembly of hexaploid bread wheat recovered the highly repetitive TE space in an almost complete chromosomal context and enabled a detailed view into the dynamics of TEs in the A, B, and D subgenomes.

RESULTS

The overall TE content is very similar between the A, B, and D subgenomes, although we find no evidence for bursts of TE amplification after the polyploidization events. Despite the near-complete turnover of TEs since the subgenome lineages diverged from a common ancestor, 76% of TE families are still present in similar proportions in each subgenome. Moreover, spacing between syntenic genes is also conserved, even though syntenic TEs have been replaced by new insertions over time, suggesting that distances between genes, but not sequences, are under evolutionary constraints. The TE composition of the immediate gene vicinity differs from the core intergenic regions. We find the same TE families to be enriched or depleted near genes in all three subgenomes. Evaluations at the subfamily level of timed long terminal repeat-retrotransposon insertions highlight the independent evolution of the diploid A, B, and D lineages before polyploidization and cases of concerted proliferation in the AB tetraploid.

CONCLUSIONS

Even though the intergenic space is changed by the TE turnover, an unexpected preservation is observed between the A, B, and D subgenomes for features like TE family proportions, gene spacing, and TE enrichment near genes.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

A Quantitative Proteomics View on the Function of Qfhb1, a Major QTL for Fusarium Head Blight Resistance in Wheat

Fusarium head blight (FHB) is a highly detrimental disease of wheat. A quantitative trait locus for FHB resistance, Qfhb1, is the most utilized source of resistance in wheat-breeding programs, but very little is known about its resistance mechanism. In this study, we elucidated a prospective FHB resistance mechanism by investigating the proteomic signatures of Qfhb1 in a pair of contrasting wheat near-isogenic lines (NIL) after 24 h of inoculation of wheat florets by Fusarium graminearum. Statistical comparisons of the abundances of protein spots on the 2D-DIGE gels of contrasting NILs (fhb1+ NIL = Qfhb1 present; fhb1- NIL = Qfhb1 absent) enabled us to select 80 high-ranking differentially accumulated protein (DAP) spots. An additional evaluation confirmed that the DAP spots were specific to the spikelet from fhb1- NIL (50 spots), and fhb1+ NIL (seven spots). The proteomic data also suggest that the absence of Qfhb1 makes the fhb1- NIL vulnerable to Fusarium attack by constitutively impairing several mechanisms including sucrose homeostasis by enhancing starch synthesis from sucrose. In the absence of Qfhb1, Fusarium inoculations severely damaged photosynthetic machinery; altered the metabolism of carbohydrates, nitrogen and phenylpropanoids; disrupted the balance of proton gradients across relevant membranes; disturbed the homeostasis of many important signaling molecules induced the mobility of cellular repair; and reduced translational activities. These changes in the fhb1- NIL led to strong defense responses centered on the hypersensitive response (HSR), resulting in infected cells suicide and the consequent initiation of FHB development. Therefore, the results of this study suggest that Qfhb1 largely functions to either alleviate HSR or to manipulate the host cells to not respond to Fusarium infection.

Molecular Characterization and Expression of PFT, an FHB Resistance Gene at the Fhb1 QTL in Wheat

Fusarium head blight (FHB) is a destructive fungal disease in wheat worldwide. Efforts have been carried out to combat this disease, and the pore-forming toxin-like (PFT) gene at the quantitative trait locus (QTL) Fhb1 was isolated and found to confer resistance to FHB in Sumai 3. In this study, we characterized PFT in 348 wheat accessions. Four haplotypes of PFT were identified. The wild haplotype of PFT had higher resistance than other haplotypes and explained 13.8% of phenotypic variation in FHB resistance by association analysis. PFT was highly expressed during early flowering and increased after Fusarium graminearum treatment in Sumai 3. Analysis of the 5' flanking sequence of PFT predicted that the cis elements of the PFT promoter were related to hormones and biological defense responses. However, PFT existed not only in the FHB-resistant accessions but also in some susceptible accessions. These results suggested that FHB resistance in a diverse range of wheat genotypes is partially conditioned by PFT. The profiling of FHB resistance and the PFT locus in this large collection of wheat germplasm may prove helpful for incorporating FHB resistance into wheat breeding programs, although more work is needed to reveal the exact role of the QTL Fhb1 in conferring resistance to fungal spread.

Association genetics studies on frost tolerance in wheat (Triticum aestivum L.) reveal new highly conserved amino acid substitutions in CBF-A3, CBF-A15, VRN3 and PPD1 genes

BACKGROUND

Understanding the genetic basis of frost tolerance (FT) in wheat (Triticum aestivum L.) is essential for preventing yield losses caused by frost due to cellular damage, dehydration and reduced metabolism. FT is a complex trait regulated by a number of genes and several gene families. Availability of the wheat genomic sequence opens new opportunities for exploring candidate genes diversity for FT. Therefore, the objectives of this study were to identity SNPs and insertion-deletion (indels) in genes known to be involved in frost tolerance and to perform association genetics analysis of respective SNPs and indels on FT.

RESULTS

Here we report on the sequence analysis of 19 candidate genes for FT in wheat assembled using the Chinese Spring IWGSC RefSeq v1.0. Out of these, the tandem duplicated C-repeat binding factors (CBF), i.e. CBF-A3, CBF-A5, CBF-A10, CBF-A13, CBF-A14, CBF-A15, CBF-A18, the vernalisation response gene VRN-A1, VRN-B3, the photoperiod response genes PPD-B1 and PPD-D1 revealed association to FT in 235 wheat cultivars. Within six genes (CBF-A3, CBF-A15, VRN-A1, VRN-B3, PPD-B1 and PPD-D1) amino acid (AA) substitutions in important protein domains were identified. The amino acid substitution effect in VRN-A1 on FT was confirmed and new AA substitutions in CBF-A3, CBF-A15, VRN-B3, PPD-B1 and PPD-D1 located at highly conserved sites were detected. Since these results rely on phenotypic data obtained at five locations in 2 years, detection of significant associations of FT to AA changes in CBF-A3, CBF-A15, VRN-A1, VRN-B3, PPD-B1 and PPD-D1 may be exploited in marker assisted breeding for frost tolerance in winter wheat.

CONCLUSIONS

A set of 65 primer pairs for the genes mentioned above from a previous study was BLASTed against the IWGSC RefSeq resulting in the identification of 39 primer combinations covering the full length of 19 genes. This work demonstrates the usefulness of the IWGSC RefSeq in specific primer development for highly conserved gene families in hexaploid wheat and, that a candidate gene association genetics approach based on the sequence data is an efficient tool to identify new alleles of genes important for the response to abiotic stress in wheat.

Ten steps to get started in Genome Assembly and Annotation

As a part of the ELIXIR-EXCELERATE efforts in capacity building, we present here 10 steps to facilitate researchers getting started in genome assembly and genome annotation. The guidelines given are broadly applicable, intended to be stable over time, and cover all aspects from start to finish of a general assembly and annotation project.

Intrinsic properties of genomes are discussed, as is the importance of using high quality DNA. Different sequencing technologies and generally applicable workflows for genome assembly are also detailed. We cover structural and functional annotation and encourage readers to also annotate transposable elements, something that is often omitted from annotation workflows. The importance of data management is stressed, and we give advice on where to submit data and how to make your results Findable, Accessible, Interoperable, and Reusable (FAIR).

High throughput SNP discovery and genotyping in hexaploid wheat

Because of their abundance and their amenability to high-throughput genotyping techniques, Single Nucleotide Polymorphisms (SNPs) are powerful tools for efficient genetics and genomics studies, including characterization of genetic resources, genome-wide association studies and genomic selection. In wheat, most of the previous SNP discovery initiatives targeted the coding fraction, leaving almost 98% of the wheat genome largely unexploited. Here we report on the use of whole-genome resequencing data from eight wheat lines to mine for SNPs in the genic, the repetitive and non-repetitive intergenic fractions of the wheat genome. Eventually, we identified 3.3 million SNPs, 49% being located on the B-genome, 41% on the A-genome and 10% on the D-genome. We also describe the development of the TaBW280K high-throughput genotyping array containing 280,226 SNPs. Performance of this chip was examined by genotyping a set of 96 wheat accessions representing the worldwide diversity. Sixty-nine percent of the SNPs can be efficiently scored, half of them showing a diploid-like clustering. The TaBW280K was proven to be a very efficient tool for diversity analyses, as well as for breeding as it can discriminate between closely related elite varieties. Finally, the TaBW280K array was used to genotype a population derived from a cross between Chinese Spring and Renan, leading to the construction a dense genetic map comprising 83,721 markers. The results described here will provide the wheat community with powerful tools for both basic and applied research.

ImmunoFISH: Simultaneous Visualisation of Proteins and DNA Sequences Gives Insight Into Meiotic Processes in Nuclei of Grasses

ImmunoFISH is a method combining immunolabelling (IL) with fluorescent in situ hybridisation (FISH) to simultaneously detect the nuclear distribution of proteins and specific DNA sequences within chromosomes. This approach is particularly important when analysing meiotic cell division where morphogenesis of individual proteins follows stage-specific changes and is accompanied by a noticeable chromatin dynamism. The method presented here is simple and provides reliable results of high quality signal, low background staining and can be completed within 2 days following preparation. Conventional widefield epifluorescent or laser scanning microscopy can be used for high resolution and three-dimensional analysis. Fixation and preparation techniques were optimised to best preserve nuclear morphology and protein epitopes without the need for any antigen retrieval. Preparation of plant material involved short cross-linking fixation of meiotic tissues with paraformaldehyde (PFA) followed by enzyme digestion and slide-mounting. In order to avoid rapid sample degradation typical of shortly fixed plant materials, and to be able to perform IL later, slides were snap-frozen and stored at -80°C. Ultra-freezing produced a remarkable degree of structural preservation for up to 12 months, whereby sample quality was similar to that of fresh material. Harsh chemicals and sample dehydration were avoided throughout the procedure and permeability was ensured by a 0.1-0.3% detergent treatment. The ImmunoFISH method was developed specifically for studying meiosis in Triticeae, but should also be applicable to other grass and plant species.

Current Strategies of Polyploid Plant Genome Sequence Assembly

Polyploidy or duplication of an entire genome occurs in the majority of angiosperms. The understanding of polyploid genomes is important for the improvement of those crops, which humans rely on for sustenance and basic nutrition. As climate change continues to pose a potential threat to agricultural production, there will increasingly be a demand for plant cultivars that can resist biotic and abiotic stresses and also provide needed and improved nutrition. In the past decade, Next Generation Sequencing (NGS) has fundamentally changed the genomics landscape by providing tools for the exploration of polyploid genomes. Here, we review the challenges of the assembly of polyploid plant genomes, and also present recent advances in genomic resources and functional tools in molecular genetics and breeding. As genomes of diploid and less heterozygous progenitor species are increasingly available, we discuss the lack of complexity of these currently available reference genomes as they relate to polyploid crops. Finally, we review recent approaches of haplotyping by phasing and the impact of third generation technologies on polyploid plant genome assembly.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

Functional and DNA-protein binding studies of WRKY transcription factors and their expression analysis in response to biotic and abiotic stress in wheat (Triticum aestivum L.)

WRKY, a plant-specific transcription factor family, plays vital roles in pathogen defense, abiotic stress, and phytohormone signalling. Little is known about the roles and function of WRKY transcription factors in response to rust diseases in wheat. In the present study, three TaWRKY genes encoding complete protein sequences were cloned. They belonged to class II and III WRKY based on the number of WRKY domains and the pattern of zinc finger structures. Twenty-two DNA-protein binding docking complexes predicted stable interactions of WRKY domain with W-box. Quantitative real-time-PCR using wheat near-isogenic lines with or without Lr28 gene revealed differential up- or down-regulation in response to biotic and abiotic stress treatments which could be responsible for their functional divergence in wheat. TaWRKY62 was found to be induced upon treatment with JA, MJ, and SA and reduced after ABA treatments. Maximum induction of six out of seven genes occurred at 48 h post inoculation due to pathogen inoculation. Hence, TaWRKY (49, 50, 52, 55, 57, and 62) can be considered as potential candidate genes for further functional validation as well as for crop improvement programs for stress resistance. The results of the present study will enhance knowledge towards understanding the molecular basis of mode of action of WRKY transcription factor genes in wheat and their role during leaf rust pathogenesis in particular.

Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure

While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.

The role of differential gene content in the evolution and function of eukaryotic genomes remains poorly explored. Here the authors assemble and annotate the Brachypodium distachyon pan-genome consisting of 54 diverse lines and reveal the differential present genes as a major driver of phenotypic variation.

Chromosome Evolution in Connection with Repetitive Sequences and Epigenetics in Plants

Chromosome evolution is a fundamental aspect of evolutionary biology. The evolution of chromosome size, structure and shape, number, and the change in DNA composition suggest the high plasticity of nuclear genomes at the chromosomal level. Repetitive DNA sequences, which represent a conspicuous fraction of every eukaryotic genome, particularly in plants, are found to be tightly linked with plant chromosome evolution. Different classes of repetitive sequences have distinct distribution patterns on the chromosomes. Mounting evidence shows that repetitive sequences may play multiple generative roles in shaping the chromosome karyotypes in plants. Furthermore, recent development in our understanding of the repetitive sequences and plant chromosome evolution has elucidated the involvement of a spectrum of epigenetic modification. In this review, we focused on the recent evidence relating to the distribution pattern of repetitive sequences in plant chromosomes and highlighted their potential relevance to chromosome evolution in plants. We also discussed the possible connections between evolution and epigenetic alterations in chromosome structure and repatterning, such as heterochromatin formation, centromere function, and epigenetic-associated transposable element inactivation.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Resolving fine-grained dynamics of retrotransposons: comparative analysis of inferential methods and genomic resources

Transposable elements support genome diversification, but comparison of their proliferation and genomic distribution within and among species is necessary to characterize their role in evolution. Such inferences are challenging because of potential bias with incomplete sampling of repetitive genome regions. Here, using the assembled genome as well as genome skimming datasets in Arabis alpina, we assessed the limits of current approaches inferring the biology of transposable elements. Long terminal repeat retrotransposons (LTR-RTs) identified in the assembled genome were classified into monophyletic lineages (here called tribes), including families of similar copies in Arabis along with elements from related Brassicaceae. Inference of their dynamics using divergence of LTRs in full-length copies and mismatch distribution of genetic variation among all copies congruently highlighted recent transposition bursts, although ancient proliferation events were apparent only with mismatch distribution. Similar inferences of LTR-RT dynamics based on random sequences from genome skimming were highly correlated with assembly-based estimates, supporting accurate analyses from shallow sequencing. Proportions of LTR-RT copies next to genes from both assembled genomes and genome skimming were congruent, pointing to tribes being over- or under-represented in the vicinity of genes. Finally, genome skimming at low coverage revealed accurate inferences of LTR-RT dynamics and distribution, although only the most abundant families appeared robustly analysed at 0.1X. Examining the pitfalls and benefits of approaches relying on different genomic resources, we highlight that random sequencing reads represent adequate data suitably complementing biased samples of LTR-RT copies retrieved from assembled genomes towards comprehensive surveys of the biology of transposable elements.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

Detection of alien genetic introgressions in bread wheat using dot-blot genomic hybridisation

Simple, reliable methods for the identification of alien genetic introgressions are required in plant breeding programmes. The use of genomic dot-blot hybridisation allows the detection of small Hordeum chilense genomic introgressions in the descendants of genetic crosses between wheat and H. chilense addition or substitution lines in wheat when molecular markers are difficult to use. Based on genomic in situ hybridisation, DNA samples from wheat lines carrying putatively H. chilense introgressions were immobilised on a membrane, blocked with wheat genomic DNA and hybridised with biotin-labelled H. chilense genomic DNA as a probe. This dot-blot screening reduced the number of plants necessary to be analysed by molecular markers or in situ hybridisation, saving time and money. The technique was sensitive enough to detect a minimum of 5 ng of total genomic DNA immobilised on the membrane or about 1/420 dilution of H. chilense genomic DNA in the wheat background. The robustness of the technique was verified by in situ hybridisation. In addition, the detection of other wheat relative species such as Hordeum vulgare, Secale cereale and Agropyron cristatum in the wheat background was also reported.

The Diversity of Sequence and Chromosomal Distribution of New Transposable Element-Related Segments in the Rye Genome Revealed by FISH and Lineage Annotation

Transposable elements (TEs) in plant genomes exhibit a great variety of structure, sequence content and copy number, making them important drivers for species diversity and genome evolution. Even though a genome-wide statistic summary of TEs in rye has been obtained using high-throughput DNA sequencing technology, the accurate diversity of TEs in rye, as well as their chromosomal distribution and evolution, remains elusive due to the repetitive sequence assembling problems and the high dynamic and nested nature of TEs. In this study, using genomic plasmid library construction combined with dot-blot hybridization and fluorescence in situ hybridization (FISH) analysis, we successfully isolated 70 unique FISH-positive TE-related sequences including 47 rye genome specific ones: 30 showed homology or partial homology with previously FISH characterized sequences and 40 have not been characterized. Among the 70 sequences, 48 sequences carried Ty3/gypsy-derived segments, 7 sequences carried Ty1/copia-derived segments and 15 sequences carried segments homologous with multiple TE families. 26 TE lineages were found in the 70 sequences, and among these lineages, Wilma was found in sequences dispersed in all chromosome regions except telomeric positions; Abiba was found in sequences predominantly located at pericentromeric and centromeric positions; Wis, Carmilla, and Inga were found in sequences displaying signals dispersed from distal regions toward pericentromeric positions; except DNA transposon lineages, all the other lineages were found in sequences displaying signals dispersed from proximal regions toward distal regions. A high percentage (21.4%) of chimeric sequences were identified in this study and their high abundance in rye genome suggested that new TEs might form through recombination and nested transposition. Our results also gave proofs that diverse TE lineages were arranged at centromeric and pericentromeric positions in rye, and lineages like Abiba might play a role in their structural organization and function. All these results might help in understanding the diversity and evolution of TEs in rye, as well as their driving forces in rye genome organization and evolution.

Vacuum and Co-cultivation Agroinfiltration of (Germinated) Seeds Results in Tobacco Rattle Virus (TRV) Mediated Whole-Plant Virus-Induced Gene Silencing (VIGS) in Wheat and Maize

Tobacco rattle virus (TRV)-mediated virus-induced gene silencing (VIGS) has been frequently used in dicots. Here we show that it can also be used in monocots, by presenting a system involving use of a novel infiltration solution (containing acetosyringone, cysteine, and Tween 20) that enables whole-plant level VIGS of (germinated) seeds in wheat and maize. Using the established system, phytoene desaturase (PDS) genes were successfully silenced, resulting in typical photo-bleaching symptoms in the leaves of treated wheat and maize. In addition, three wheat homoeoalleles of MLO, a key gene repressing defense responses to powdery mildew in wheat, were simultaneously silenced in susceptible wheat with this system, resulting in it becoming resistant to powdery mildew. The system has the advantages generally associated with TRV-mediated VIGS systems (e.g., high-efficiency, mild virus infection symptoms, and effectiveness in different organs). However, it also has the following further advantages: (germinated) seed-stage agroinfiltration; greater rapidity and convenience; whole-plant level gene silencing; adequately stable transformation; and suitability for studying functions of genes involved in seed germination and early plant development stages.

Chromosome-specific sequencing reveals an extensive dispensable genome component in wheat

The hexaploid wheat genotype Chinese Spring (CS) has been used worldwide as the reference base for wheat genetics and genomics, and significant resources have been used by the international community to generate a reference wheat genome based on this genotype. By sequencing flow-sorted 3B chromosome from a hexaploid wheat genotype CRNIL1A and comparing the obtained sequences with those available for CS, we detected that a large number of sequences in the former were missing in the latter. If the distribution of such sequences in the hexaploid wheat genome is random, CRNILA sequences missing in CS could be as much as 159.3 Mb even if only fragments of 50 bp or longer were considered. Analysing RNA sequences available in the public domains also revealed that dispensable genes are common in hexaploid wheat. Together with those extensive intra- and interchromosomal rearrangements in CS, the existence of such dispensable genes is another factor highlighting potential issues with the use of reference genomes in various studies. Strong deviation in distributions of these dispensable sequences among genotypes with different geographical origins provided the first evidence indicating that they could be associated with adaptation in wheat.