Reconciling the evolutionary origin of bread wheat (Triticum aestivum)

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

Genomic content of wheat has a higher influence than plant domestication status on the ability to interact with Pseudomonas plant growth-promoting rhizobacteria

Plant evolutionary history has had profound effects on belowground traits, which is likely to have impacted the ability to interact with microorganisms, but consequences on root colonization and gene expression by plant growth-promoting rhizobacteria (PGPR) remain poorly understood. Here, we tested the hypothesis that wheat genomic content and domestication are key factors determining the capacity for PGPR interaction. Thus, 331 wheat representatives from eight Triticum or Aegilops species were inoculated under standardized conditions with the generalist PGPR Pseudomonas ogarae F113, using an autofluorescent reporter system for monitoring F113 colonization and expression of phl genes coding for the auxinic inducing signal 2,4-diacetylphloroglucinol. The interaction with P. ogarae F113 was influenced by ploidy level, presence of genomes AA, BB, DD, and domestication. While root colonization was higher for hexaploid and tetraploid species, and phl expression level higher for hexaploid wheat, the diploid Ae. tauschii displayed higher phl induction rate (i.e., expression:colonisation ratio) on roots. However, a better potential of interaction with F113 (i.e., under non-stress gnotobiotic conditions) did not translate, after seed inoculation, into better performance of wheat landraces in non-sterile soil under drought. Overall, results showed that domestication and especially plant genomic content modulate the PGPR interaction potential of wheats.

Genome shock in a synthetic allotetraploid wheat invokes subgenome-partitioned gene regulation, meiotic instability, and karyotype variation

It is becoming increasingly evident that interspecific hybridization at the homoploid level or coupled with whole-genome duplication (i.e. allopolyploidization) has played a major role in biological evolution. However, the direct impacts of hybridization and allopolyploidization on genome structure and function, phenotype, and fitness remains to be fully understood. Synthetic hybrids and allopolyploids are trackable experimental systems that can be used to address this issue. In this study, we resynthesized a pair of reciprocal F1 hybrids and corresponding reciprocal allotetraploids using the two diploid progenitor species of bread wheat (Triticum aestivum, BBAADD), namely T. urartu (AA) and Aegilops tauschii (DD). By comparing phenotypes related to growth, development, and fitness, and by analysing genome expression in both hybrids and allotetraploids in relation to the parents, we found that the types and trends of karyotype variation in the immediately formed allotetraploids were correlated with both instability of meiosis and chromosome- and subgenome-biased expression. We determined clear advantages of allotetraploids over diploid F1 hybrids in several morphological traits including fitness that mirrored the tissue- and developmental stage-dependent subgenome-partitioning of the allotetraploids. The allotetraploids were meiotically unstable primarily due to homoeologous pairing that varied dramatically among the chromosomes. Nonetheless, the manifestation of organismal karyotype variation and the occurrence of meiotic irregularity were not concordant, suggesting a role of functional constraints probably imposed by subgenome- and chromosome-biased gene expression. Our results provide new insights into the direct impacts and consequences of hybridization and allopolyploidization that are relevant to evolution and likely to be informative for future crop improvement approaches using synthetic polyploids.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

Genetic diversity analysis in wheat cultivars using SCoT and ISSR markers, chloroplast DNA barcoding and grain SEM

BACKGROUND

Wheat is a major cereal that can narrow the gap between the increasing human population and food production. In this connection, assessing genetic diversity and conserving wheat genetic resources for future exploitation is very important for breeding new cultivars that may withstand the expected climate change. The current study evaluates the genetic diversity in selected wheat cultivars using ISSR and SCoT markers, the rbcL and matK chloroplast DNA barcoding, and grain surface sculpture characteristics. We anticipate that these objectives may prioritize using the selected cultivars to improve wheat production. The selected collection of cultivars may lead to the identification of cultivars adapted to a broad spectrum of climatic environments.

RESULTS

Multivariate clustering analyses of the ISSR and SCoT DNA fingerprinting polymorphism grouped three Egyptian cultivars with cultivar El-Nielain from Sudan, cultivar Aguilal from Morocco, and cultivar Attila from Mexico. In the other group, cultivar Cook from Australia and cultivar Chinese-166 were differentiated from four other cultivars: cultivar Cham-10 from Syria, cultivar Seri-82 from Mexico, cultivar Inqalab-91 from Pakistan, and cultivar Sonalika from India. In the PCA analysis, the Egyptian cultivars were distinct from the other studied cultivars. The rbcL and matK sequence variation analysis indicated similarities between Egyptian cultivars and cultivar Cham-10 from Syria and cultivar Inqalab-91 from Pakistan, whereas cultivar Attila from Mexico was distinguished from all other cultivars. Combining the data of ISSR and SCoT with the rbcL and matK results retained the close resemblance among the two Egyptian cultivars EGY1: Gemmeiza-9 and EGY3: Sakha-93, and the Moroccan cultivar Aguilal, and the Sudanese cultivar El-Nielain and between Seri-82, Inqalab-91, and Sonalika cultivars. The analysis of all data distinguished cultivar Cham-10 from Syria from all other cultivars, and the analysis of grain traits indicated a close resemblance between cv. Cham-10 from and the two Egyptian cultivars Gemmeiza-9 and Sakha-93.

CONCLUSIONS

The analysis of rbcL and matK chloroplast DNA barcoding agrees with the ISSR and the SCoT markers in supporting the close resemblance between the Egyptian cultivars, particularly Gemmeiza-9 and Sakha-93. The ISSR and SCoT data analyses significantly expressed high differentiation levels among the examined cultivars. Cultivars with closer resemblance may be recommended for breeding new wheat cultivars adapted to various climatic environments.

The TaTCP4/10-B1 cascade regulates awn elongation in wheat (Triticum aestivum L.)

Awns are important morphological markers for wheat and exert a strong physiological effect on wheat yield. The awn elongation suppressor B1 has recently been cloned through association and linkage analysis in wheat. However, the mechanism of awn inhibition centered around B1 remains to be clarified. Here, we identified an allelic variant in the coding region of B1 through analysis of re-sequencing data; this variant causes an amino acid substitution and premature termination, resulting in a long-awn phenotype. Transcriptome analysis indicated that B1 inhibited awn elongation by impeding cytokinin- and auxin-promoted cell division. Moreover, B1 directly repressed the expression of TaRAE2 and TaLks2, whose orthologs have been reported to promote awn development in rice or barley. More importantly, we found that TaTCP4 and TaTCP10 synergistically inhibited the expression of B1, and a G-to-A mutation in the B1 promoter attenuated its inhibition by TaTCP4/10. Taken together, our results reveal novel mechanisms of awn development and provide genetic resources for trait improvement in wheat.

Wheat genome architecture influences interactions with phytobeneficial microbial functional groups in the rhizosphere

Wheat has undergone a complex evolutionary history, which led to allopolyploidization and the hexaploid bread wheat Triticum aestivum. However, the significance of wheat genomic architecture for beneficial plant-microbe interactions is poorly understood, especially from a functional standpoint. In this study, we tested the hypothesis that wheat genomic architecture was an overriding factor determining root recruitment of microorganisms with particular plant-beneficial traits. We chose five wheat species representing genomic profiles AA (Triticum urartu), BB {SS} (Aegilops speltoides), DD (Aegilops tauschii), AABB (Triticum dicoccon) and AABBDD (Triticum aestivum) and assessed by quantitative polymerase chain reaction their ability to interact with free-nitrogen fixers, 1-aminocyclopropane-1-carboxylate deaminase producers, 2,4-diacetylphloroglucinol producers and auxin producers via the phenylpyruvate decarboxylase pathway, in combination with Illumina MiSeq metabarcoding analysis of N fixers (and of the total bacterial community). We found that the abundance of the microbial functional groups could fluctuate according to wheat genomic profile, as did the total bacterial abundance. N fixer diversity and total bacterial diversity were also influenced significantly by wheat genomic profile. Often, rather similar results were obtained for genomes DD (Ae. tauschii) and AABBDD (T. aestivum), pointing for the first time that the D genome could be particularly important for wheat-bacteria interactions.

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

Does one subgenome become dominant in the formation and evolution of a polyploid?

BACKGROUND

Polyploids are common in flowering plants and they tend to have more expanded ranges of distributions than their diploid progenitors. Possible mechanisms underlying polyploid success have been intensively investigated. Previous studies showed that polyploidy generates novel changes and that subgenomes in allopolyploid species often differ in gene number, gene expression levels and levels of epigenetic alteration. It is widely believed that such differences are the results of conflicts among the subgenomes. These differences have been treated by some as subgenome dominance, and it is claimed that the magnitude of subgenome dominance increases in polyploid evolution.

SCOPE

In addition to changes which occurred during evolution, differences between subgenomes of a polyploid species may also be affected by differences between the diploid donors and changes which occurred during polyploidization. The variable genome components in many plant species are extensive, which would result in exaggerated differences between a subgenome and its progenitor when a single genotype or a small number of genotypes are used to represent a polyploid or its donors. When artificially resynthesized polyploids are used as surrogates for newly formed genotypes which have not been exposed to evolutionary selection, differences between diploid genotypes available today and those involved in the formation of the natural polyploid genotypes must also be considered.

CONCLUSIONS

Contrary to the now widely held views that subgenome biases in polyploids are the results of conflicts among the subgenomes and that one of the parental subgenomes generally retains more genes which are more highly expressed, available results show that subgenome biases mainly reflect legacy from the progenitors and that they can be detected before the completion of polyploidization events. Further, there is no convincing evidence that the magnitudes of subgenome biases have significantly changed during evolution for any of the allopolyploid species assessed.

Dynamic chromatin regulatory programs during embryogenesis of hexaploid wheat

BACKGROUND

Plant and animal embryogenesis have conserved and distinct features. Cell fate transitions occur during embryogenesis in both plants and animals. The epigenomic processes regulating plant embryogenesis remain largely elusive.

RESULTS

Here, we elucidate chromatin and transcriptomic dynamics during embryogenesis of the most cultivated crop, hexaploid wheat. Time-series analysis reveals stage-specific and proximal-distal distinct chromatin accessibility and dynamics concordant with transcriptome changes. Following fertilization, the remodeling kinetics of H3K4me3, H3K27ac, and H3K27me3 differ from that in mammals, highlighting considerable species-specific epigenomic dynamics during zygotic genome activation. Polycomb repressive complex 2 (PRC2)-mediated H3K27me3 deposition is important for embryo establishment. Later H3K27ac, H3K27me3, and chromatin accessibility undergo dramatic remodeling to establish a permissive chromatin environment facilitating the access of transcription factors to cis-elements for fate patterning. Embryonic maturation is characterized by increasing H3K27me3 and decreasing chromatin accessibility, which likely participates in restricting totipotency while preventing extensive organogenesis. Finally, epigenomic signatures are correlated with biased expression among homeolog triads and divergent expression after polyploidization, revealing an epigenomic contributor to subgenome diversification in an allohexaploid genome.

CONCLUSIONS

Collectively, we present an invaluable resource for comparative and mechanistic analysis of the epigenomic regulation of crop embryogenesis.

Phylogenomics as an effective approach to untangle cross-species hybridization event: A case study in the family Nymphaeaceae

Hybridization is common and considered as an important evolutionary force to increase intraspecific genetic diversity. Detecting hybridization events is crucial for understanding the evolutionary history of species and further improving molecular breeding. The studies on identifying hybridization events through the phylogenomic approach are still limited. We proposed the conception and method of identifying allopolyploidy events by phylogenomics. The reconciliation and summary of nuclear multi-labeled gene family trees were adopted to untangle hybridization events from next-generation data in our novel phylogenomic approach. Given horticulturalists’ relatively clear cultivated crossbreeding history, the water lily family is a suitable case for examining recent allopolyploidy events. Here, we reconstructed and confirmed the well-resolved nuclear phylogeny for the Nymphaeales family in the context of geological time as a framework for identifying hybridization signals. We successfully identified two possible allopolyploidy events with the parental lineages for the hybrids in the family Nymphaeaceae based on summarization from multi-labeled gene family trees of Nymphaeales. The lineages where species Nymphaea colorata and Nymphaea caerulea are located may be the progenitors of horticultural cultivated species Nymphaea ‘midnight’ and Nymphaea ‘Woods blue goddess’. The proposed hybridization hypothesis is also supported by horticultural breeding records. Our methodology can be widely applied to identify hybridization events and theoretically facilitate the genome breeding design of hybrid plants.

Exploration of the genomic atlas of Dof transcription factor family across genus Oryza provides novel insights on rice breeding in changing climate

DNA-binding with one finger (Dof) transcription factors have been demonstrated to regulate various stresses and developmental processes in plants. Their identification and comparative evolutionary analyses in cultivated and wild species of genus oryza were yet to be explored. In this context, we report a comprehensive genomics atlas of DNA-binding with one finger (Dof) family genes in 13 diverse rice genomes (five cultivated and eight rice wild-relatives) through a genome-wide scanning approach. A galore of 238 Dof genes, identified across the genus Oryza, are categorized into seven distinct subgroups by comparative phylogenetic analysis with the model plant Arabidopsis. Conserved motifs and gene structure analyses unveiled the prevalence of species- and subgroups-specific structural and functional diversity that is expediating with the evolutionary period. Our results indicate that Dof genes might have undergone strong purifying selections and segmental duplications to expand their gene family members in corresponding Oryza genomes. We speculate that miR2927 potentially targets the Dof domain to regulate gene expression under different climatic conditions, which are supported by in-silico and wet-lab experiments-based expression profiles. In a nutshell, we report several superior haplotypes significantly associated with early flowering in a treasure trove of 3,010 sequenced rice accessions and have validated these haplotypes with two years of field evaluation-based flowering data of a representative subpanel. Finally, we have provided some insights on the resolution of Oryza species phylogeny discordance and divergence highlighting the mosaic evolutionary history of the genus Oryza. Overall, this study reports a complete genomic landscape of the Dof family in cultivated and wild Oryza species that could greatly facilitate in fast-track development of early maturing and climate-resilient rice cultivars through modern haplotype-led breeding.

The grain quality of wheat wild relatives in the evolutionary context

We evaluated the potential of wheat wild relatives for the improvement in grain quality characteristics including micronutrients (Fe, Zn) and gluten and identified diploid wheats and the timopheevii lineage as the most promising resources. Domestication enabled the advancement of civilization through modification of plants according to human requirements. Continuous selection and cultivation of domesticated plants induced genetic bottlenecks. However, ancient diversity has been conserved in crop wild relatives. Wheat (Triticum aestivum L.; Triticum durum Desf.) is one of the most important staple foods and was among the first domesticated crop species. Its evolutionary diversity includes diploid, tetraploid and hexaploid species from the Triticum and Aegilops taxa and different genomes, generating an AA, BBAA/GGAA and BBAADD/GGAAA^mA^m genepool, respectively. Breeding and improvement in wheat altered its grain quality. In this review, we identified evolutionary patterns and the potential of wheat wild relatives for quality improvement regarding the micronutrients Iron (Fe) and Zinc (Zn), the gluten storage proteins α-gliadins and high molecular weight glutenin subunits (HMW-GS), and the secondary metabolite phenolics. Generally, the timopheevii lineage has been neglected to date regarding grain quality studies. Thus, the timopheevii lineage should be subject to grain quality research to explore the full diversity of the wheat gene pool.

Hybridization, missing wild ancestors and the domestication of cultivated diploid bananas

Hybridization and introgressions are important evolutionary forces in plants. They contribute to the domestication of many species, including understudied clonal crops. Here, we examine their role in the domestication of a clonal crop of outmost importance, banana (Musa ssp.). We used genome-wide SNPs generated for 154 diploid banana cultivars and 68 samples of the wild M. acuminata to estimate and geo-localize the contribution of the different subspecies of M. acuminata to cultivated banana. We further investigated the wild to domesticate transition in New Guinea, an important domestication center. We found high levels of admixture in many cultivars and confirmed the existence of unknown wild ancestors with unequal contributions to cultivated diploid. In New Guinea, cultivated accessions exhibited higher diversity than their direct wild ancestor, the latter recovering from a bottleneck. Introgressions, balancing selection and positive selection were identified as important mechanisms for banana domestication. Our results shed new lights on the radiation of M. acuminata subspecies and on how they shaped banana domestication. They point candidate regions of origin for two unknown ancestors and suggest another contributor in New Guinea. This work feed research on the evolution of clonal crops and has direct implications for conservation, collection, and breeding.

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.

Genome-wide identification of the jumonji C domain- containing histone demethylase gene family in wheat and their expression analysis under drought stress

Methylation and demethylation of histone play a crucial role in regulating chromatin formation and gene expression. The jumonji C (JmjC) domain-containing proteins are demethylases that are involved in regulating epigenetic modification in plants. In our study, the JmjC genes in Triticum aestivum L., Triticum turgidum L., Triticum dicoccoides L., Triticum urartu L., and Aegilops tauschii L. were identified. Phylogenetic relationship and colinearity analysis revealed that the wheat JmjC genes were conserved in A, B, and D subgenomes during evolution. Cis-acting elements analysis showed that elements related to stress response, hormone response, and light response were found in wheat JmjC genes. The expression of JmjC genes was affected by tissue types and developmental stages, and members of the same subfamily tended to have similar expression patterns in wheat. They also showed a unique expression pattern in root during PEG (Polyethylene glycol) treatment. In conclusion, comprehensive analysis indicated that three members (Tr-1A-JMJ2, Tr-1B-JMJ2, and Tr-1D-JMJ2) might be regulated by several hormones and function in the early stages of drought stress, while eight members (Tr-1B-JMJ3, Tr-4B-JMJ1, Tr-7A-JMJ1, etc.) displayed a significantly high expression after 24 h of PEG treatment, indicating a role in the later stages of drought stress. This research presents the first genome-wide study of the JmjC family in wheat, and lays the foundation for promoting the study of their functional characterization in wheat drought resistance.

Dispersed emergence and protracted domestication of polyploid wheat uncovered by mosaic ancestral haploblock inference

Major crops are all survivors of domestication bottlenecks. Studies have focused on the genetic loci related to the domestication syndrome, while the contribution of ancient haplotypes remains largely unknown. Here, an ancestral genomic haploblock dissection method is developed and applied to a resequencing dataset of 386 tetraploid/hexaploid wheat accessions, generating a pan-ancestry haploblock map. Together with cytoplastic evidences, we reveal that domesticated polyploid wheat emerged from the admixture of six founder wild emmer lineages, which contributed the foundation of ancestral mosaics. The key domestication-related loci, originated over a wide geographical range, were gradually pyramided through a protracted process. Diverse stable-inheritance ancestral haplotype groups of the chromosome central zone are identified, revealing the expanding routes of wheat and the trends of modern wheat breeding. Finally, an evolution model of polyploid wheat is proposed, highlighting the key role of wild-to-crop and interploidy introgression, that increased genomic diversity following bottlenecks introduced by domestication and polyploidization.

The contribution of ancient haplotypes to domestication is largely unknown. Here, the authors develop an ancestral genomic haploblock dissection method to generate a mosaic pan-ancestry genomic map and reveal that the domesticated polyploidy wheat emerged from the admixture of six founder wild emmer linages.

Development and Characterization of Triticum aestivum-Aegilops longissima 6Sl Recombinants Harboring a Novel Powdery Mildew Resistance Gene Pm6Sl

Powdery mildew of wheat is a foliar disease that is spread worldwide. Cultivation of resistant varieties is the most effective, economical, and environmentally friendly strategy to curb this disease. Powdery mildew resistance genes (Pm) are the primary resources for resistance breeding, and new Pm genes are in constant demand. Previously, we identified Aegilops longissima chromosome 6S^l#3 as a carrier of powdery mildew resistance and designated the resistance gene as Pm6Sl. Here, we reported the design of 24 markers specific to 6S^l#3 on the basis of the full-length cDNA sequences of 6S^l#3 donor Ae. longissma accession TA1910, and the development of wheat-Ae. longissima 6S^l#3 introgression stocks by ph1b-induced homoeologous recombination. Further, 6S^l#3 introgression lines were identified and characterized by integration analysis of powdery mildew responses, in situ hybridization, and molecular markers and Pm6Sl was mapped to a distal interval of 42.80 Mb between markers Ael58410 and Ael57699 in the long arm of 6S^l#3. Two resistant recombinants, R43 (T6BS.6BL-6S^l#3L) and T27 (Ti6AS.6AL-6S^l#3L-6AL), contained segments harboring Pm6Sl with less than 8% of 6S^l#3 genomic length, and two markers were diagnostic for Pm6Sl. This study broadened powdery mildew resistance gene resources for wheat improvement and provided a fundamental basis for fine mapping and cloning of Pm6Sl to further understand its molecular mechanism of disease resistance.

Whole-genome resequencing of the wheat A subgenome progenitor Triticum urartu provides insights into its demographic history and geographic adaptation

Triticum urartu is the progenitor of the A subgenome in tetraploid and hexaploid wheat. Uncovering the landscape of genetic variations in T. urartu will help us understand the evolutionary and polyploid characteristics of wheat. Here, we investigated the population genomics of T. urartu by genome-wide sequencing of 59 representative accessions collected around the world. A total of 42.2 million high-quality single-nucleotide polymorphisms and 3 million insertions and deletions were obtained by mapping reads to the reference genome. The ancient T. urartu population experienced a significant reduction in effective population size (Ne) from ∼3 000 000 to ∼140 000 and subsequently split into eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations during the Younger Dryas period. A map of allelic drift paths displayed splits and mixtures between different geographic groups, and a strong genetic drift towards hexaploid wheat was also observed, indicating that the direct donor of the A subgenome originated from northwestern Syria. Genetic changes were revealed between the eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations in genes orthologous to those regulating plant development and stress responses. A genome-wide association study identified two single-nucleotide polymorphisms in the exonic regions of the SEMI-DWARF 37 ortholog that corresponded to the different T. urartu ecotype groups. Our study provides novel insights into the origin and genetic legacy of the A subgenome in polyploid wheat and contributes a gene repertoire for genomics-enabled improvements in wheat breeding.

Dissection of Structural Reorganization of Wheat 5B Chromosome Associated With Interspecies Recombination Suppression

Chromosomal rearrangements that lead to recombination suppression can have a significant impact on speciation, and they are also important for breeding. The regions of recombination suppression in wheat chromosome 5B were identified based on comparisons of the 5B map of a cross between the Chinese Spring (CS) variety of hexaploid wheat and CS-5Bdic (genotype CS with 5B substituted with its homologue from tetraploid Triticum dicoccoides) with several 5B maps of tetraploid and hexaploid wheat. In total, two regions were selected in which recombination suppression occurred in cross CS × CS-5Bdic when compared with other maps: one on the short arm, 5BS_RS, limited by markers BS00009810/BS00022336, and the second on the long arm, 5BL_RS, between markers Ra_c10633_2155 and BS00087043. The regions marked as 5BS_RS and 5BL_RS, with lengths of 5 Mb and 3.6 Mb, respectively, were mined from the 5B pseudomolecule of CS and compared to the homoeologous regions (7.6 and 3.8 Mb, respectively) of the 5B pseudomolecule of Zavitan (T. dicoccoides). It was shown that, in the case of 5BS_RS, the local heterochromatin islands determined by the satellite DNA (119.2) and transposable element arrays, as well as the dissimilarity caused by large insertions/deletions (chromosome rearrangements) between 5BSs aestivum/dicoccoides, are likely the key determinants of recombination suppression in the region. Two major and two minor segments with significant loss of similarity were recognized within the 5BL_RS region. It was shown that the loss of similarity, which can lead to suppression of recombination in the 5BL_RS region, is caused by chromosomal rearrangements, driven by the activity of mobile genetic elements (both DNA transposons and long terminal repeat retrotransposons) and their divergence during evolution. It was noted that the regions marked as 5BS_RS and 5BL_RS are associated with chromosomal rearrangements identified earlier by С-banding analysis of intraspecific polymorphism of tetraploid emmer wheat. The revealed divergence in 5BS_RS and 5BL_RS may be a consequence of interspecific hybridization, plant genetic adaptation, or both.

Global Patterns of subgenome evolution in organelle-targeted genes of six allotetraploid angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary “mismatches” between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

The LONELY GUY gene family: from mosses to wheat, the key to the formation of active cytokinins in plants

LONELY GUY (LOG) was first identified in a screen of rice mutants with defects in meristem maintenance. In plants, LOG codes for cytokinin riboside 5'-monophosphate phosphoribohydrolase, which converts inactive cytokinin nucleotides directly to the active free bases. Many enzymes with the PGGxGTxxE motif have been misannotated as lysine decarboxylases; conversely not all enzymes containing this motif are cytokinin-specific LOGs. As LOG mutants clearly impact yield in rice, we investigated the LOG gene family in bread wheat. By interrogating the wheat (Triticum aestivum) genome database, we show that wheat has multiple LOGs. The close alignment of TaLOG1, TaLOG2 and TaLOG6 with the X-ray structures of two functional Arabidopsis thaliana LOGs allows us to infer that the wheat LOGs 1-11 are functional LOGs. Using RNA-seq data sets, we assessed TaLOG expression across 70 tissue types, their responses to various stressors, the pattern of cis-regulatory elements (CREs) and intron/exon patterns. TaLOG gene family members are expressed variously across tissue types. When the TaLOG CREs are compared with those of the cytokinin dehydrogenases (CKX) and glucosyltransferases (CGT), there is close alignment of CREs between TaLOGs and TaCKXs reflecting the key role of CKX in maintaining cytokinin homeostasis. However, we suggest that the main homeostatic mechanism controlling cytokinin levels in response to biotic and abiotic challenge resides in the CGTs, rather than LOG or CKX. However, LOG transgenics and identified mutants in rice variously impact yield, providing interesting avenues for investigation in wheat.

Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome

Common wheat (Triticum aestivum, BBAADD) is a major staple food crop worldwide. The diploid progenitors of the A and D subgenomes have been unequivocally identified; that of B, however, remains ambiguous and controversial but is suspected to be related to species of Aegilops, section Sitopsis. Here, we report the assembly of chromosome-level genome sequences of all five Sitopsis species, namely Aegilops bicornis, Ae. longissima, Ae. searsii, Ae. sharonensis, and Ae. speltoides, as well as the partial assembly of the Amblyopyrum muticum (synonym Aegilops mutica) genome for phylogenetic analysis. Our results reveal that the donor of the common wheat B subgenome is a distinct, and most probably extinct, diploid species that diverged from an ancestral progenitor of the B lineage to which the still extant Ae. speltoides and Am. muticum belong. In addition, we identified interspecific genetic introgressions throughout the evolution of the Triticum/Aegilops species complex. The five Sitopsis species have various assembled genome sizes (4.11-5.89 Gb) with high proportions of repetitive sequences (85.99%-89.81%); nonetheless, they retain high collinearity with other genomes or subgenomes of species in the Triticum/Aegilops complex. Differences in genome size were primarily due to independent post-speciation amplification of transposons. We also identified a set of Sitopsis genes pertinent to important agronomic traits that can be harnessed for wheat breeding. These newly assembled genome resources provide a new roadmap for evolutionary and genetic studies of the Triticum/Aegilops complex, as well as for wheat improvement.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

The wheat dioxygenase BX6 is involved in the formation of benzoxazinoids in planta and contributes to plant defense against insect herbivores

Benzoxazinoids are plant specialized metabolites with defense properties, highly abundant in wheat (Triticum), one of the world's most important crops. The goal of our study was to characterize dioxygenase BX6 genes in tetraploid and hexaploid wheat genotypes and to elucidate their effects on defense against herbivores. Phylogenetic analysis revealed four BX6 genes in the hexaploid wheat T. aestivum, but only one ortholog was found in the tetraploid (T. turgidum) wild emmer wheat and the cultivated durum wheat. Transcriptome sequencing of durum wheat plants, damaged by either aphids or caterpillars, revealed that several BX genes, including TtBX6, were upregulated upon caterpillar feeding, relative to the undamaged control plants. A virus-induced gene silencing approach was used to reduce the expression of BX6 in T. aestivum plants, which exhibited both reduced transcript levels and reduced accumulation of different benzoxazinoids. To elucidate the effect of BX6 on plant defense, bioassays with different herbivores feeding on BX6-silenced leaves were conducted. The results showed that plants with silenced BX6 were more susceptible to aphids and the two-spotted spider mite than the control. Overall, our study indicates that wheat BX6 is involved in benzoxazinoid formation in planta and contributes to plant resistance against insect herbivores.

Interfered chromosome pairing at high temperature promotes meiotic instability in autotetraploid Arabidopsis

Changes in environmental temperature affect multiple meiotic processes in flowering plants. Polyploid plants derived from whole-genome duplication (WGD) have enhanced genetic plasticity and tolerance to environmental stress but face challenges in organizing and segregating doubled chromosome sets. In this study, we investigated the impact of increased environmental temperature on male meiosis in autotetraploid Arabidopsis (Arabidopsis thaliana). Under low to mildly increased temperatures (5°C-28°C), irregular chromosome segregation universally occurred in synthetic autotetraploid Columbia-0 (Col-0). Similar meiotic lesions occurred in autotetraploid rice (Oryza sativa L.) and allotetraploid canola (Brassica napus cv Westar), but not in evolutionarily derived hexaploid wheat (Triticum aestivum). At extremely high temperatures, chromosome separation and tetrad formation became severely disordered due to univalent formation caused by the suppression of crossing-over. We found a strong correlation between tetravalent formation and successful chromosome pairing, both of which were negatively correlated with temperature elevation, suggesting that increased temperature interferes with crossing-over predominantly by impacting homolog pairing. We also showed that loading irregularities of axis proteins ASY1 and ASY4 co-localize on the chromosomes of the syn1 mutant and the heat-stressed diploid and autotetraploid Col-0, revealing that heat stress affects the lateral region of synaptonemal complex (SC) by impacting the stability of the chromosome axis. Moreover, we showed that chromosome axis and SC in autotetraploid Col-0 are more sensitive to increased temperature than those in diploid Arabidopsis. Taken together, our data provide evidence suggesting that WGD negatively affects the stability and thermal tolerance of meiotic recombination in newly synthetic autotetraploid Arabidopsis.

Genome-Wide Diversity of MADS-Box Genes in Bread Wheat is Associated with its Rapid Global Adaptability

MADS-box gene family members play multifarious roles in regulating the growth and development of crop plants and hold enormous promise for bolstering grain yield potential under changing global environments. Bread wheat (Triticum aestivum L.) is a key stable food crop around the globe. Until now, the available information concerning MADS-box genes in the wheat genome has been insufficient. Here, a comprehensive genome-wide analysis identified 300 high confidence MADS-box genes from the publicly available reference genome of wheat. Comparative phylogenetic analyses with Arabidopsis and rice MADS-box genes classified the wheat genes into 16 distinct subfamilies. Gene duplications were mainly identified in subfamilies containing unbalanced homeologs, pointing towards a potential mechanism for gene family expansion. Moreover, a more rapid evolution was inferred for M-type genes, as compared with MIKC-type genes, indicating their significance in understanding the evolutionary history of the wheat genome. We speculate that subfamily-specific distal telomeric duplications in unbalanced homeologs facilitate the rapid adaptation of wheat to changing environments. Furthermore, our in-silico expression data strongly proposed MADS-box genes as active guardians of plants against pathogen insurgency and harsh environmental conditions. In conclusion, we provide an entire complement of MADS-box genes identified in the wheat genome that could accelerate functional genomics efforts and possibly facilitate bridging gaps between genotype-to-phenotype relationships through fine-tuning of agronomically important traits.

Genome-wide characterization and functional analysis of class III peroxidase gene family in soybean reveal regulatory roles of GsPOD40 in drought tolerance

Class III peroxidases (PODs) are plant-specific glycoproteins, that play essential roles in various plant physiological processes and defence responses. To date, scarce information is available about the POD gene family in soybean. Hence, the present study is the first comprehensive report about the genome-wide characterization of GmPOD gene family in soybean (Glycine max L.). Here, we identified a total of 124 GmPOD genes in soybean, that are unevenly distributed across the genome. Phylogenetic analysis classified them into six distinct sub-groups (A-F), with one soybean specific subgroup. Exon-intron and motif analysis suggested the existence of structural and functional diversity among the sub-groups. Duplication analysis identified 58 paralogous gene pairs; segmental duplication and positive/Darwinian selection were observed as the major factors involved in the evolution of GmPODs. Furthermore, RNA-seq analysis revealed that 23 out of a total 124 GmPODs showed differential expression between drought-tolerant and drought-sensitive genotypes under stress conditions; however, two of them (GmPOD40 and GmPOD42) revealed the maximum deregulation in all contrasting genotypes. Overexpression (OE) lines of GsPOD40 showed considerably higher drought tolerance compared to wild type (WT) plants under stress treatment. Moreover, the OE lines showed enhanced photosynthesis and enzymatic antioxidant activities under drought stress, resulting in alleviation of ROS induced oxidative damage. Hence, the GsPOD40 enhanced drought tolerance in soybean by regulating the key physiological and biochemical pathways involved in the defence response. Lastly, the results of our study will greatly assist in further functional characterization of GsPODs in plant growth and stress tolerance in soybean.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

The origin and evolution of a two-component system of paralogous genes encoding the centromeric histone CENH3 in cereals

BACKGROUND

The cereal family Poaceae is one of the largest and most diverse angiosperm families. The central component of centromere specification and function is the centromere-specific histone H3 (CENH3). Some cereal species (maize, rice) have one copy of the gene encoding this protein, while some (wheat, barley, rye) have two. We applied a homology-based approach to sequenced cereal genomes, in order to finally trace the mutual evolution of the structure of the CENH3 genes and the nearby regions in various tribes.

RESULTS

We have established that the syntenic group or the CENH3 locus with the CENH3 gene and the boundaries defined by the CDPK2 and bZIP genes first appeared around 50 Mya in a common ancestor of the subfamilies Bambusoideae, Oryzoideae and Pooideae. This locus came to Pooideae with one copy of CENH3 in the most ancient tribes Nardeae and Meliceae. The βCENH3 gene as a part of the locus appeared in the tribes Stipeae and Brachypodieae around 35-40 Mya. The duplication was accompanied by changes in the exon-intron structure. Purifying selection acts mostly on αCENH3s, while βCENH3s form more heterogeneous structures, in which clade-specific amino acid motifs are present. In barley species, the βCENH3 gene assumed an inverted orientation relative to αCENH3 and the CDPK2 gene was substituted with LHCB-l. As the evolution and domestication of plant species went on, the locus was growing in size due to an increasing distance between αCENH3 and βCENH3 because of a massive insertion of the main LTR-containing retrotransposon superfamilies, gypsy and copia, without any evolutionary preference on either of them. A comparison of the molecular structure of the locus in the A, B and D subgenomes of the hexaploid wheat T. aestivum showed that invasion by mobile elements and concomitant rearrangements took place in an independent way even in evolutionarily close species.

CONCLUSIONS

The CENH3 duplication in cereals was accompanied by changes in the exon-intron structure of the βCENH3 paralog. The observed general tendency towards the expansion of the CENH3 locus reveals an amazing diversity of ways in which different species implement the scenario described in this paper.

Transcriptome based identification and validation of heat stress transcription factors in wheat progenitor species Aegilops speltoides

Wheat, one of the major cereal crops worldwide, get adversely affected by rising global temperature. We have identified the diploid B genome progenitor of wheat, Aegilops speltoides (SS), as a potential donor for heat stress tolerance. Therefore, the present work was planned to study the total transcriptome profile of heat stress-tolerant Ae. speltoides accession pau3809 (AS3809) and compare with that of tetraploid and hexaploid wheat cultivars PDW274 and PBW725, respectively. The comparative transcriptome was utilized to identify and validate heat stress transcription factors (HSFs), the key genes involved in imparting heat stress tolerance. Transcriptome analysis led to the identification of a total of 74 K, 68 K, and 76 K genes in AS3809, PDW274, and PBW725, respectively. There was a high uniformity of GO profiles under the biological, molecular, and cellular functions across the three wheat transcriptomes, suggesting the conservation of gene function. Twelve HSFs having the highest FPKM value were identified in the AS3809 transcriptome data, while six of these HSFs namely HSFA3, HSFA5, HSFA9, HSFB2a, HSFB2b, and HSFC1b, were validated with qRT PCR. These six HSFs were identified as an important component of thermotolerance in AS3809 as evident from their comparative higher expression under heat stress.

Nuclear-cytoplasmic balance: whole genome duplications induce elevated organellar genome copy number

The plant genome is partitioned across three distinct subcellular compartments: the nucleus, mitochondria, and plastids. Successful coordination of gene expression among these organellar genomes and the nuclear genome is critical for plant function and fitness. Whole genome duplication (WGD) events in the nucleus have played a major role in the diversification of land plants and are expected to perturb the relative copy number (stoichiometry) of nuclear, mitochondrial, and plastid genomes. Thus, elucidating the mechanisms whereby plant cells respond to the cytonuclear stoichiometric imbalance that follows WGDs represents an important yet underexplored question in understanding the evolutionary consequences of genome doubling. We used droplet digital PCR to investigate the relationship between nuclear and organellar genome copy numbers in allopolyploids and their diploid progenitors in both wheat and Arabidopsis. Polyploids exhibit elevated organellar genome copy numbers per cell, largely preserving the cytonuclear stoichiometry observed in diploids despite the change in nuclear genome copy number. To investigate the timescale over which cytonuclear stoichiometry may respond to WGD, we also estimated the organellar genome copy number in Arabidopsis synthetic autopolyploids and in a haploid-induced diploid line. We observed corresponding changes in organellar genome copy number in these laboratory-generated lines, indicating that at least some of the cellular response to cytonuclear stoichiometric imbalance is immediate following WGD. We conclude that increases in organellar genome copy numbers represent a common response to polyploidization, suggesting that maintenance of cytonuclear stoichiometry is an important component in establishing polyploid lineages.

Characterizing chloroplast genomes and inferring maternal divergence of the Triticum-Aegilops complex

The Triticum (wheat)–Aegilops (goatgrass) complex has been extensively studied, but the evolutionary history of polyploid wheats has not been fully elucidated. The chloroplast (cp) with maternal inheritance and homoplasy can simplify the sequence-based evolutionary inferences, but informative inferences would require a complete and accurate cp genome sequence. In this study, 16 cp genomes representing five Aegilops and 11 Triticum species and subspecies were sequenced, assembled and annotated, yielding five novel circular cp genome sequences. Analyzing the assembled cp genomes revealed no marked differences in genome structure and gene arrangement across the assayed species. A polymorphism analysis of 72 published cp genome sequences representing 10 Aegilops and 15 Triticum species and subspecies detected 1183 SNPs and 1881 SSRs. More than 80% SNPs detected resided on the downstream and upstream gene regions and only 2.78% or less SNPs were predicted to be deleterious. The largest nucleotide diversity was observed in the short single-copy genomic region. Relatively weak selection pressure on cp coding genes was detected. Different phylogenetic analyses confirmed that the maternal divergence of the Triticum–Aegilops complex had three deep lineages each representing a diploid species with nuclear A, B, or D genome. Dating the maternal divergence yielded age estimates of divergence that matched well with those reported previously. The divergence between emmer and bread wheats occurred at 8200–11,200 years ago. These findings are useful for further genomic studies, provide insight into cp genome evolvability and allow for better understanding of the maternal divergence of the Triticum–Aegilops complex.

Chromosome number is key to longevity of polyploid lineages

Polyploidy is ubiquitous and often recursive in plant lineages, most frequently resulting in extinction but occasionally associated with great evolutionary success. However, instead of chromosome numbers exponentially increasing due to recurrent polyploidy, most angiosperm species have fewer than 14 chromosome pairs. Following genome duplication, diploidisation can render one copy of essential genes nonfunctional without fitness cost. In isolated subpopulations, alternate (homoeologous) gene copies can be lost, creating incompatibilities that reduce fitness of hybrids between subpopulations, constraining exchange of favourable genetic changes and reducing species fitness. When multiple sets of incompatible genes are genetically linked, their deleterious effects are not independent. The effective number of independently acting sets of incompatible loci in hybrids is limited by chromosome number and recombination. Therefore, species with many chromosomes are subject to a higher fitness penalty during diploidisation. Karyotypic changes, especially fusions, that reduce gene flow are normally fitness disadvantages, but during the diploidisation process, can increase fitness by reducing mixing of differentially diploidised alleles. Fitness penalties caused by diploidisation favour accelerated karyotypic change, with each change increasing barriers to gene flow, contributing to speciation. Lower chromosome numbers and increased chromosome fusions confer advantages to surviving the diploidisation process following polyploid formation, by independent mechanisms.

Genetic Diversity, Linkage Disequilibrium and Population Structure of Bulgarian Bread Wheat Assessed by Genome-Wide Distributed SNP Markers: From Old Germplasm to Semi-Dwarf Cultivars

Genetic diversity and population structure are key resources for breeding purposes and genetic studies of important agronomic traits in crops. In this study, we described SNP-based genetic diversity, linkage disequilibrium and population structure in a panel of 179 bread wheat advanced cultivars and old accessions from Bulgaria, using an optimized wheat 25K Infinium iSelect array. Out of 19,019 polymorphic SNPs, 17,968 had a known chromosome position on the A (41%), B (42%) and D (11%) genome, and 6% were not assigned to any chromosome. Homoeologous group 4, in particular chromosome 4D, was the least polymorphic. In the total population, the Nei’s gene diversity was within the range 0.1–0.5, and the polymorphism information content ranged from 0.1 to 0.4. Significant differences between the old and modern collections were revealed with respect to the linkage disequilibrium (LD): the average values for LD (r²), the percentage of the locus pairs in LD and the LD decay were 0.64, 16% and 3.3 for the old germplasm, and 0.43, 30% and 4.1 for the modern releases, respectively. Structure and k-means clustering algorithm divided the panel into three groups. The old accessions formed a distinct subpopulation. The cluster analysis further distinguished the modern releases according to the geographic region and genealogy. Gene exchange was evidenced mainly between the subpopulations of contemporary cultivars. The achieved understanding of the genetic diversity and structure of the Bulgarian wheat population and distinctiveness of the old germplasm could be of interest for breeders developing cultivars with improved characteristics. The obtained knowledge about SNP informativeness and the LD estimation are worthwhile for selecting markers and for considering the composition of a population in association mapping studies of traits of interest.

Ploidy influences wheat mesophyll cell geometry, packing and leaf function

Leaf function is influenced by leaf structure, which is itself related not only to the spatial arrangement of constituent mesophyll cells, but also their size and shape. In this study, we used confocal microscopy to image leaves of Triticum genotypes varying in ploidy level to extract 3D information on individual mesophyll cell size and geometry. Combined with X-ray Computed Tomography and gas exchange analysis, the effect of changes in wheat mesophyll cell geometry upon leaf structure and function were investigated. Mesophyll cell size and shape were found to have changed during the course of wheat evolution. An unexpected linear relationship between mesophyll cell surface area and volume was discovered, suggesting anisotropic scaling of mesophyll cell geometry with increasing ploidy. Altered mesophyll cell size and shape were demonstrated to be associated with changes in mesophyll tissue architecture. Under experimental growth conditions, CO2 assimilation did not vary with ploidy, but stomatal conductance was lower in hexaploid plants, conferring a greater instantaneous water-use efficiency. We propose that as wheat mesophyll cells have become larger with increased ploidy, this has been accompanied by changes in cell geometry and packing which limit water loss while maintaining carbon assimilation.

Extensive structural variation in the Bowman-Birk inhibitor family in common wheat (Triticum aestivum L.)

BACKGROUND

Bowman-Birk inhibitors (BBI) are a family of serine-type protease inhibitors that modulate endogenous plant proteolytic activities during different phases of development. They also inhibit exogenous proteases as a component of plant defense mechanisms, and their overexpression can confer resistance to phytophagous herbivores and multiple fungal and bacterial pathogens. Dicot BBIs are multifunctional, with a "double-headed" structure containing two separate inhibitory loops that can bind and inhibit trypsin and chymotrypsin proteases simultaneously. By contrast, monocot BBIs have a non-functional chymotrypsin inhibitory loop, although they have undergone internal duplication events giving rise to proteins with multiple BBI domains.

RESULTS

We used a Hidden Markov Model (HMM) profile-based search to identify 57 BBI genes in the common wheat (Triticum aestivum L.) genome. The BBI genes are unevenly distributed, with large gene clusters in the telomeric regions of homoeologous group 1 and 3 chromosomes that likely arose through a series of tandem gene duplication events. The genomes of wheat progenitors also contain contiguous clusters of BBI genes, suggesting this family underwent expansion before the domestication of common wheat. However, the BBI gene family varied in size among different cultivars, showing this family remains dynamic. Because of these expansions, the BBI gene family is larger in wheat than other monocots such as maize, rice and Brachypodium. We found BBI proteins in common wheat with intragenic homologous duplications of cysteine-rich functional domains, including one protein with four functional BBI domains. This diversification may expand the spectrum of target substrates. Expression profiling suggests that some wheat BBI proteins may be involved in regulating endogenous proteases during grain development, while others were induced in response to biotic and abiotic stresses, suggesting a role in plant defense.

CONCLUSIONS

Genome-wide characterization reveals that the BBI gene family in wheat is subject to a high rate of homologous tandem duplication and deletion events, giving rise to a diverse set of encoded proteins. This information will facilitate the functional characterization of individual wheat BBI genes to determine their role in wheat development and stress responses, and their potential application in breeding.

Identification of individual components of a commercial wheat germ acid phosphatase preparation

Wheat germ acid phosphatase (WGAP) is a commercial preparation of partially purified protein commonly used in laboratory settings for non-specific enzymatic dephosphorylation. It is known that these preparations contain multiple phosphatase isozymes and are still relatively crude. This study therefore aimed to identify the protein components of a commercial preparation of wheat germ acid phosphatase using mass spectroscopy and comparative genomics. After one post-purchase purification step, the most prevalent fifteen proteins in the mixture included heat shock proteins, beta-amylases, glucoseribitol dehydrogenases, enolases, and an aminopeptidase. While not among the most abundant components, eight unique dephosphorylation enzymes were also present including three purple acid phosphatases. Furthermore, it is shown that some of these correspond to previously isolated isozymes; one of which has been also previously shown by transcriptome data to be overexpressed in wheat seeds. In summary, this study identified the major components of WGAP including phosphatases and hypothesizes the most active components towards a better understanding of this commonly used laboratory tool.

Differences between diploid donors are the main contributing factor for subgenome asymmetry measured in either gene ratio or relative diversity in allopolyploids

Subgenome asymmetry (SA) has routinely been attributed to different responses between the subgenomes of a polyploid to various stimuli during evolution. Here, we compared subgenome differences in gene ratio and relative diversity between artificial and natural genotypes of several allopolyploid species. Surprisingly, consistent differences were not detected between these two types of polyploid genotypes, although they differ in times exposed to evolutionary selection. The estimated ratio of shared genes between a subgenome and its diploid donor was invariably higher for the artificial allopolyploid genotypes than those for the natural genotypes, which is expected as it is now well-known that many genes in a species are not shared among all individuals. As the exact diploid parent for a given subgenome is unknown, the estimated ratios of shared genes for the natural genotypes would also include difference among individual genotypes of the diploid donor species. Further, we detected the presence of SA in genotypes before the completion of the polyploidization events as well as in those which were not formed via polyploidization. These results indicate that SA may, to a large degree, reflect differences between its diploid donors or that changes occurred during polyploid evolution are defined by their donor genomes.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Artificial selection in breeding extensively enriched a functional allelic variation in TaPHS1 for pre-harvest sprouting resistance in wheat

Pre-harvest sprouting (PHS) causes significant losses in wheat yield and quality worldwide. Previously, we cloned a PHS resistance gene, TaPHS1, and identified two causal mutations for reduced seed dormancy (SD) and increased PHS susceptibility. Here we identified a novel allelic variation of C to T transition in 3'-UTR of TaPHS1, which associated with reduced SD and PHS resistance. The T allele occurred in wild wheat progenitors and was likely the earliest functional mutation in TaPHS1 for PHS susceptibility. Allele frequency analysis revealed low frequency of the T allele in wild diploid and tetraploid wheat progenitors, but very high frequency in modern wheat cultivars and breeding lines, indicating that artificial selection quickly enriched the T allele during modern breeding. The T allele was significantly associated with short SD in both T. aestivum and T. durum, the two most cultivated species of wheat. This variation together with previously reported functional sequence variations co-regulated TaPHS1 expression levels and PHS resistance in different germplasms. Haplotype analysis of the four functional variations identified the best PHS resistance haplotype of TaPHS1. The resistance haplotype can be used in marker-assisted selection to transfer TaPHS1 to new wheat cultivars.

Coevolution in Hybrid Genomes: Nuclear-Encoded Rubisco Small Subunits and Their Plastid-Targeting Translocons Accompanying Sequential Allopolyploidy Events in Triticum

The Triticum/Aegilops complex includes hybrid species resulting from homoploid hybrid speciation and allopolyploid speciation. Sequential allotetra- and allohexaploidy events presumably result in two challenges for the hybrids, which involve 1) cytonuclear stoichiometric disruptions caused by combining two diverged nuclear genomes with the maternal inheritance of the cytoplasmic organellar donor; and 2) incompatibility of chimeric protein complexes with diverged subunits from nuclear and cytoplasmic genomes. Here, we describe coevolution of nuclear rbcS genes encoding the small subunits of Rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase) and nuclear genes encoding plastid translocons, which mediate recognition and translocation of nuclear-encoded proteins into plastids, in allopolyploid wheat species. We demonstrate that intergenomic paternal-to-maternal gene conversion specifically occurred in the genic region of the homoeologous rbcS3 gene from the D-genome progenitor of wheat (abbreviated as rbcS3D) such that it encodes a maternal-like or B-subgenome-like SSU3D transit peptide in allohexaploid wheat but not in allotetraploid wheat. Divergent and limited interaction between SSU3D and the D-subgenomic TOC90D translocon subunit is implicated to underpin SSU3D targeting into the chloroplast of hexaploid wheat. This implicates early selection favoring individuals harboring optimal maternal-like organellar SSU3D targeting in hexaploid wheat. These data represent a novel dimension of cytonuclear evolution mediated by organellar targeting and transportation of nuclear proteins.

Application of Genomics Tools in Wheat Breeding to Attain Durable Rust Resistance

Wheat is an important source of dietary protein and calories for the majority of the world's population. It is one of the largest grown cereal in the world occupying over 215 M ha. Wheat production globally is challenged by biotic stresses such as pests and diseases. Of the 50 diseases of wheat that are of economic importance, the three rust diseases are the most ubiquitous causing significant yield losses in the majority of wheat production environments. Under severe epidemics they can lead to food insecurity threats amid the continuous evolution of new races of the pathogens, shifts in population dynamics and their virulence patterns, thereby rendering several effective resistance genes deployed in wheat breeding programs vulnerable. This emphasizes the need to identify, characterize, and deploy effective rust-resistant genes from diverse sources into pre-breeding lines and future wheat varieties. The use of genetic resistance has been marked as eco-friendly and to curb the further evolution of rust pathogens. Deployment of multiple rust resistance genes including major and minor genes in wheat lines could enhance the durability of resistance thereby reducing pathogen evolution. Advances in next-generation sequencing (NGS) platforms and associated bioinformatics tools have revolutionized wheat genomics. The sequence alignment of the wheat genome is the most important landmark which will enable genomics to identify marker-trait associations, candidate genes and enhanced breeding values in genomic selection (GS) studies. High throughput genotyping platforms have demonstrated their role in the estimation of genetic diversity, construction of the high-density genetic maps, dissecting polygenic traits, and better understanding their interactions through GWAS (genome-wide association studies) and QTL mapping, and isolation of R genes. Application of breeder's friendly KASP assays in the wheat breeding program has expedited the identification and pyramiding of rust resistance alleles/genes in elite lines. The present review covers the evolutionary trends of the rust pathogen and contemporary wheat varieties, and how these research strategies galvanized to control the wheat killer genus Puccinia. It will also highlight the outcome and research impact of cost-effective NGS technologies and cloning of rust resistance genes amid the public availability of common and tetraploid wheat reference genomes.

Reduced chromatin accessibility underlies gene expression differences in homologous chromosome arms of diploid Aegilops tauschii and hexaploid wheat

Background

Polyploidy is centrally important in the evolution and domestication of plants because it leads to major genomic changes, such as altered patterns of gene expression, which are thought to underlie the emergence of new traits. Despite the common occurrence of these globally altered patterns of gene expression in polyploids, the mechanisms involved are not well understood.

Results

Using a precisely defined framework of highly conserved syntenic genes on hexaploid wheat chromosome 3DL and its progenitor 3 L chromosome arm of diploid Aegilops tauschii, we show that 70% of these gene pairs exhibited proportionately reduced gene expression, in which expression in the hexaploid context of the 3DL genes was ∼40% of the levels observed in diploid Ae tauschii. Several genes showed elevated expression during the later stages of grain development in wheat compared with Ae tauschii. Gene sequence and methylation differences probably accounted for only a few cases of differences in gene expression. In contrast, chromosome-wide patterns of reduced chromatin accessibility of genes in the hexaploid chromosome arm compared with its diploid progenitor were correlated with both reduced gene expression and the imposition of new patterns of gene expression.

Conclusions

Our pilot-scale analyses show that chromatin compaction may orchestrate reduced gene expression levels in the hexaploid chromosome arm of wheat compared to its diploid progenitor chromosome arm.

The Battle to Sequence the Bread Wheat Genome: A Tale of the Three Kingdoms

In the year 2018, the world witnessed the finale of the race to sequence the genome of the world's most widely grown crop, the common wheat. Wheat has been known to bear a notoriously large and complicated genome of a polyploidy nature. A decade competition to sequence the wheat genome initiated with a single consortium of multiple countries, taking a conventional strategy similar to that for sequencing Arabidopsis and rice, became ferocious over time as both sequencing technologies and genome assembling methodologies advanced. At different stages, multiple versions of genome sequences of the same variety (e.g., Chinese Spring) were produced by several groups with their special strategies. Finally, 16 years after the rice genome was finished and 9 years after that of maize, the wheat research community now possesses its own reference genome. Armed with these genomics tools, wheat will reestablish itself as a model for polyploid plants in studying the mechanisms of polyploidy evolution, domestication, genetic and epigenetic regulation of homoeolog expression, as well as defining its genetic diversity and breeding on the genome level. The enhanced resolution of the wheat genome should also help accelerate development of wheat cultivars that are more tolerant to biotic and/or abiotic stresses with better quality and higher yield.

Preferential Subgenome Elimination and Chromosomal Structural Changes Occurring in Newly Formed Tetraploid Wheat-Aegilops ventricosa Amphiploid (AABBDvDvNvNv)

Artificial allopolyploids derived from the genera Triticum and Aegilops have been used as genetic resources for wheat improvement and are a classic example of evolution via allopolyploidization. In this study, we investigated chromosomes and subgenome transmission behavior in the newly formed allopolyploid of wheat group via multicolor Fluorescence in situ hybridization (mc-FISH), using pSc119.2, pTa535, and (GAA)₇ as probe combinations, to enabled us to precisely identify individual chromosomes in 381 S₃ and S₄ generations plants derived from reciprocal crosses between Ae. ventricosa (D^vD^vN^vN^v) and T. turgidum (AABB). A higher rate of aneuploidy, constituting 66.04–86.41% individuals, was observed in these two early generations. Of the four constituent subgenomes, D^v showed the highest frequency of elimination, followed by N^v and B, while A was the most stable. In addition, structural chromosomal changes occurred ubiquitously in the selfed progenies of allopolyploids. Among the constituent subgenomes, B showed the highest number of aberrations. In terms of chromosomal dynamics, there was no significant association between the chromosomal behavior model and the cytoplasm, with the exception of chromosomal loss in the D^v subgenome. The chromosome loss frequency in the D^v subgenome was significantly higher in the T. turgidum × Ae. ventricosa cross than in the Ae. ventricosa × T. turgidum cross. This result indicates that, although the D subgenome showed great instability, allopolyploids containing D subgenome could probably be maintained after a certain hybridization in which the D subgenome donor was used as the maternal parent at its onset stage. Our findings provide valuable information pertaining to the behavior patterns of subgenomes during allopolyploidization. Moreover, the allopolyploids developed here could be used as potential resources for the genetic improvement of wheat.