Evolution of the BBAA component of bread wheat during its history at the allohexaploid level

Yield increment and transcriptome response caused by blue light treatment in Hericium coralloides

BACKGROUND

Hericium coralloides is a traditional edible and medicinal mushroom. Light is a key factor in forming fruiting bodies of fungi; however, the effects of different light on the yield and morphogenesis of H. coralloides are still unknown. Therefore, the morphology, yield, and transcriptome of H. coralloides under blue, red, and white light conditions were investigated.

RESULT

Fruiting bodies under blue light exhibited superior morphological traits, such as milky white color, larger size, elongated stalks, and higher spine count, leading to higher yields. Different light treatments led to dramatic transcriptome changes ranging from 10,827 differentially expressed genes (DEGs) induced by blue light in Blue-4d to 11,375 DEGs induced by red light in Red-4d and accounted for 64.56% to 67.81% of all expressed genes. This massive amount of light-responsive genes has never been reported in fungi. Gene Ontology analysis showed that light affected nearly all aspects of life in H. coralloides; suggesting that the influence of light on fungi may have been underestimated. Blue light-induced yield increment may be achieved by specifically upregulating the growth-related processes such as DNA replication, chromosomes, and cell division.

CONCLUSIONS

This study offers preliminary insights into the potential role of blue light in modulating gene expression and yield stimulation in H. coralloides, potentially improving cultivation practice.

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Ploidy variation induces butterfly effect on chromatin topology in wheat

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

The SPL transcription factor TaSPL6 negatively regulates drought stress response in wheat

Environmental stresses, such as heat and drought, severely affect plant growth and development, and reduce wheat yield and quality globally. Squamosa promoter binding protein-like (SPL) proteins are plant-specific transcription factors that play a critical role in regulating plant responses to diverse stresses. In this study, we cloned and characterized TaSPL6, a wheat orthologous gene of rice OsSPL6. Three TaSPL6 homoeologs are located on the long arms of chromosomes 4A, 5B, and 5D, respectively, and share more than 98% sequence identity with each other. The TaSPL6 genes were preferentially expressed in roots, and their expression levels were downregulated in wheat seedlings subjected to heat, dehydration, and abscisic acid treatments. Subcellular localization experiments showed that TaSPL6 was localized in the nucleus. Overexpression of TaSPL6-A in wheat resulted in enhanced sensitivity to drought stress. The transgenic lines exhibited higher leaf water loss, malondialdehyde and reactive oxygen species (ROS) content, and lower antioxidant enzyme activities after drought treatment than wild-type plants. Gene silencing of TaSPL6 enhanced the drought tolerance of wheat, as reflected by better growth status. Additionally, RNA-seq and qRT-PCR analyses revealed that TaSPL6-A functions by decreasing the expression of a number of genes involved in stress responses. These findings suggest that TaSPL6 acts as a negative regulator of drought stress responses in plants, which may have major implications for understanding and enhancing crop tolerance to environmental stresses.

Developmental Instability and Gene Dysregulation in an Extracted Tetraploid from Hexaploid Wheat

The BBAA subgenomes of hexaploid common wheat can be 'extracted' to constitute a viable and self-reproducing novel tetraploid wheat, termed extracted tetraploid wheat (ETW). Prior studies have shown ETW manifesting phenotypic abnormalities and alteration in gene expression and epigenetic modifications. No population level investigation has been conducted, leaving the issue unclear regarding whether developmental stability, an essential property evolved in all natural organisms, might have been undermined in ETW. Here, we measured variations in five morphological traits and somatic chromosomal stability in populations of ETW and of its hexaploid donor, a resynthesized hexaploid and a natural tetraploid wheat. We observed phenotypic defects in ETW. Meanwhile, we documented much greater within-population variations in ETW than in the other wheat genotypes, most probably due to disrupted developmental stability in ETW. Also, somatic structural chromosome variations were detected only in ETW. Comparative transcriptome analyses indicated that the disrupted developmental stability of ETW is likely linked to massive dysregulation of genome-wide gene expression rather than to genetic mutations. Population network analysis of gene expression implicated intrinsic connectivity among the variable traits, while gene set enrichment analysis provided possible links between dysregulated gene expression and interlaced trait variation.

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.

Transcriptional Interactions of Single B-Subgenome Chromosome with C-Subgenome in B. oleracea-nigra Additional Lines

Serial monosomic alien addition lines (MAALs) provide an ideal system to elucidate the transcriptomic interactions between the alien chromosomes and recipient genome under aneuploidy. Herein, five available Brassica oleracea-nigra MAALs (CCB1, CCB4, CCB5, CCB6, CCB8), their derived B. oleracea plants (non-MAALs), and two parents were analyzed for their gene expressions by using high-throughput technology. Compared to parental B. oleracea, all MAALs showed various numbers of DEGs, but CCB8 gave much higher DEGs; the number of downregulated DEGs was slightly higher than the number of upregulated ones, except for in relation to CCB8. All derived B. oleracea plants also gave certain numbers of DEGs, despite these being much lower than in the respective MAALs. Compared to B. nigra, in all five MAALs more DEGs were downregulated than upregulated. Trans-effects were likely more prevailing than cis-effects, and these DEGs were predominantly associated with material transport by dysregulating the cellular component. Meanwhile, the orthologous genes on alien chromosomes could only play a feeble compensatory role for those gene pairs in C-subgenome, and different levels of the expressed genes had a greater tendency towards downregulation. These results revealed transcriptional aneuploidy response patterns between two genomes and suggested that cis- and trans-mechanisms synergistically regulated alien gene transcriptions after distant hybridization.

Comparative analysis of Fusarium crown rot resistance in synthetic hexaploid wheats and their parental genotypes

BACKGROUND

Fusarium crown rot (FCR) is a chronic disease of cereals worldwide. Compared with tetraploid wheat, hexaploid wheat is more resistant to FCR infection. The underlying reasons for the differences are still not clear. In this study, we compared FCR responses of 10 synthetic hexaploid wheats (SHWs) and their tetraploid and diploid parents. We then performed transcriptome analysis to uncover the molecular mechanism of FCR on these SHWs and their parents.

RESULTS

We observed higher levels of FCR resistance in the SHWs compared with their tetraploid parents. The transcriptome analysis suggested that multiple defense pathways responsive to FCR infection were upregulated in the SHWs. Notably, phenylalanine ammonia lyase (PAL) genes, involved in lignin and salicylic acid (SA) biosynthesis, exhibited a higher level of expression to FCR infection in the SHWs. Physiological and biochemical analysis validated that PAL activity and SA and lignin contents of the stem bases were higher in SHWs than in their tetraploid parents.

CONCLUSION

Overall, these findings imply that improved FCR resistance in SHWs compared with their tetraploid parents is probably related to higher levels of response on PAL-mediated lignin and SA biosynthesis pathways.

Development of homozygous tetraploid potato and whole genome doubling-induced the enrichment of H3K27ac and potentially enhanced resistance to cold-induced sweetening in tubers

Polyploid plants typically display advantages on some agronomically important traits over their diploid counterparts. Extensive studies have shown genetic, transcriptomic, and epigenetic dynamics upon polyploidization in multiple plant species. However, few studies have unveiled those alternations imposed only by ploidy level, without any interference from heterozygosity. Cultivated potato is highly heterozygous. Thus, in this study, we developed two homozygous autotetraploid lines and one homozygous diploid line in parallel from a homozygous diploid potato. We confirmed their ploidy levels using chloroplast counting and karyotyping. Oligo-FISH and genome re-sequencing validated that these potato lines are nearly homozygous. We investigated variations in phenotypes, transcription, and histone modifications between two ploidies. Both autotetraploid lines produced larger but fewer tubers than the diploid line. Interestingly, each autotetraploid line displayed ploidy-related differential expression for various genes. We also discovered a genome-wide enrichment of H3K27ac in genic regions upon whole-genome doubling (WGD). However, such enrichment was not associated with the differential gene expression between two ploidies. The tetraploid lines may exhibit better resistance to cold-induced sweetening (CIS) than the diploid line in tubers, potentially regulated through the expression of CIS-related key genes, which seems to be associated with the levels of H3K4me3 in cold-stored tubers. These findings will help to understand the impacts of autotetraploidization on dynamics of phenotypes, transcription, and histone modifications, as well as on CIS-related genes in response to cold storage.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Heat Stress Tolerance 2 confers basal heat stress tolerance in allohexaploid wheat (Triticum aestivum L.)

Heat stress substantially reduces the yield potential of wheat (Triticum aestivum L.), one of the most widely cultivated staple crops, and greatly threatens global food security in the context of global warming. However, few studies have explored the heat stress tolerance (HST)-related genetic resources in wheat. Here, we identified and fine-mapped a wheat HST locus, TaHST2, which is indispensable for HST in both the vegetative and reproductive stages of the wheat life cycle. The studied pair of near isogenic lines (NILs) exhibited diverse morphologies under heat stress, based on which we mapped TaHST2 to a 485 kb interval on chromosome arm 4DS. Under heat stress, TaHST2 confers a superior conversion rate from soluble sugars to starch in wheat grains, resulting in faster grain filling and a higher yield potential. A further exploration of genetic resources indicated that TaHST2 underwent strong artificial selection during wheat domestication, suggesting it is an essential locus for basal HST in wheat. Our findings provide deeper insights into the genetic basis of wheat HST and might be useful for global efforts to breed heat-stress-tolerant cultivars.

Pentaploidization Enriches the Genetic Diversity of Wheat by Enhancing the Recombination of AB Genomes

Allohexaploidization and continuous introgression play a key role in the origin and evolution of bread wheat. The genetic bottleneck of bread wheat resulting from limited germplasms involved in the origin and modern breeding may be compensated by gene flow from tetraploid wheat through introgressive hybridization. The inter-ploidy hybridization between hexaploid and tetraploid wheat generates pentaploid hybrids first, which absorbed genetic variations both from hexaploid and tetraploid wheat and have great potential for re-evolution and improvement in bread wheat. Therefore, understanding the effects of the pentaploid hybrid is of apparent significance in our understanding of the historic introgression and in informing breeding. In the current study, two sets of F2 populations of synthetic pentaploid wheat (SPW1 and SPW2) and synthetic hexaploid wheat (SHW1 and SHW2) were created to analyze differences in recombination frequency (RF) of AB genomes and distorted segregation of polymorphic SNP markers through SNP genotyping. Results suggested that (1) the recombination of AB genomes in the SPW populations was about 3- to 4-fold higher than that in the SHW populations, resulting from the significantly (P < 0.01) increased RF between adjacent and linked SNP loci, especially the variations that occurred in a pericentromeric region which would further enrich genetic diversity; (2) the crosses of hexaploid × tetraploid wheat could be an efficient way to produce pentaploid derivatives than the crosses of tetraploid × hexaploid wheat according to the higher germination rate found in the former crosses; (3) the high proportion of distorted segregation loci that skewed in favor of the female parent genotype/allele in the SPW populations might associate with the fitness and survival of the offspring. Based on the presented data, we propose that pentaploid hybrids should increasingly be used in wheat breeding. In addition, the contribution of gene flow from tetraploid wheat to bread wheat mediated by pentaploid introgressive hybridization also was discussed in the re-evolution of bread wheat.

Conservation and Divergence of SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family between Wheat and Rice

The SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family affects plant architecture, panicle structure, and grain development, representing key genes for crop improvements. The objective of the present study is to utilize the well characterized SPLs' functions in rice to facilitate the functional genomics of TaSPL genes. To achieve these goals, we combined several approaches, including genome-wide analysis of TaSPLs, comparative genomic analysis, expression profiling, and functional study of TaSPL3 in rice. We established the orthologous relationships of 56 TaSPL genes with the corresponding OsSPLs, laying a foundation for the comparison of known SPL functions between wheat and rice. Some TaSPLs exhibited different spatial-temporal expression patterns when compared to their rice orthologs, thus implicating functional divergence. TaSPL2/6/8/10 were identified to respond to different abiotic stresses through the combination of RNA-seq and qPCR expression analysis. Additionally, ectopic expression of TaSPL3 in rice promotes heading dates, affects leaf and stem development, and leads to smaller panicles and decreased yields per panicle. In conclusion, our work provides useful information toward cataloging of the functions of TaSPLs, emphasized the conservation and divergence between TaSPLs and OsSPLs, and identified the important SPL genes for wheat improvement.

Identification and characterization of QTL for spike morphological traits, plant height and heading date derived from the D genome of natural and resynthetic allohexaploid wheat

QHd.cau-7D.1 for heading date was delimited into the physical interval of approximately 17.38 Mb harboring three CONSTANS-like zinc finger genes. Spike morphological traits, plant height and heading date play important roles in yield improvement of wheat. To reveal the genetic factors that controlling spike morphological traits, plant height and heading date on the D genome, we conducted analysis of quantitative traits locus (QTL) using 198 F_7:8 recombinant inbred lines (RILs) derived from a cross between the common wheat TAA10 and resynthesized allohexaploid wheat XX329 with similar AABB genomes. A total of 23 environmentally stable QTL on the D sub-genome for spike length (SL), fertile spikelet number per spike (FSN), sterile spikelet number per spike (SSN), total spikelet number per spike (TSN), spike compactness (SC), plant height (PHT) and heading date (HD) were detected, among which eight appeared to be novel QTL. Furthermore, QHd.cau-7D.1 and QPht.cau-7D.2 shared identical confidence interval and were delimited into the physical interval of approximately 17.38 Mb with 145 annotated genes, including three CONSTANS-like zinc finger genes (TraesCS7D02G209000, TraesCS7D02G213000 and TraesCS7D02G220300). This study will help elucidate the molecular mechanism of the seven traits (SL, FSN, SSN, TSN, SC, PHT and HD) and provide a potentially valuable resource for genetic improvement.

De novo whole-genome assembly of Chrysanthemum makinoi, a key wild chrysanthemum

Chrysanthemum is among the top 10 cut, potted, and perennial garden flowers in the world. Despite this, to date, only the genomes of two wild diploid chrysanthemums have been sequenced and assembled. Here, we present the most complete and contiguous chrysanthemum de novo assembly published so far, as well as a corresponding ab initio annotation. The cultivated hexaploid varieties are thought to originate from a hybrid of wild chrysanthemums, among which the diploid Chrysanthemum makinoi has been mentioned. Using a combination of Oxford Nanopore long reads, Pacific Biosciences long reads, Illumina short reads, Dovetail sequences, and a genetic map, we assembled 3.1 Gb of its sequence into nine pseudochromosomes, with an N50 of 330 Mb and a BUSCO complete score of 92.1%. Our ab initio annotation pipeline predicted 95,074 genes and marked 80.0% of the genome as repetitive. This genome assembly of C. makinoi provides an important step forward in understanding the chrysanthemum genome, evolution, and history.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry

BACKGROUND

Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry.

RESULTS

Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat.

CONCLUSIONS

We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Conservation and trans-regulation of histone modification in the A and B subgenomes of polyploid wheat during domestication and ploidy transition

BACKGROUND

Polyploidy has played a prominent role in the evolution of plants and many other eukaryotic lineages. However, how polyploid genomes adapt to the abrupt presence of two or more sets of chromosomes via genome regulation remains poorly understood. Here, we analyzed genome-wide histone modification and gene expression profiles in relation to domestication and ploidy transition in the A and B subgenomes of polyploid wheat.

RESULTS

We found that epigenetic modification patterns by two typical euchromatin histone markers, H3K4me3 and H3K27me3, for the great majority of homoeologous triad genes in A and B subgenomes were highly conserved between wild and domesticated tetraploid wheats and remained stable in the process of ploidy transitions from hexaploid to extracted tetraploid and then back to resynthesized hexaploid. However, a subset of genes was differentially modified during tetraploid and hexaploid wheat domestication and in response to ploidy transitions, and these genes were enriched for particular gene ontology (GO) terms. The extracted tetraploid wheat manifested higher overall histone modification levels than its hexaploid donor, and which were reversible and restored to normal levels in the resynthesized hexaploid. Further, while H3K4me3 marks were distally distributed along each chromosome and significantly correlated with subgenome expression as expected, H3K27me3 marks showed only a weak distal bias and did not show a significant correlation with gene expression.

CONCLUSIONS

Our results reveal overall high stability of histone modification patterns in the A and B subgenomes of polyploid wheat during domestication and in the process of ploidy transitions. However, modification levels of a subset of functionally relevant genes in the A and B genomes were trans-regulated by the D genome in hexaploid wheat.

Histone H3K27 dimethylation landscapes contribute to genome stability and genetic recombination during wheat polyploidization

Bread wheat (Triticum aestivum) is an allohexaploid that was formed via two allopolyploidization events. Growing evidence suggests histone modifications are involved in the response to 'genomic shock' and environmental adaptation during polyploid formation and evolution. However, the role of histone modifications, especially histone H3 lysine-27 dimethylation (H3K27me2), in genome evolution remains elusive. Here we analyzed H3K27me2 and H3K27me3 profiles in hexaploid wheat and its tetraploid and diploid relatives. Although H3K27me3 levels were relatively stable among wheat species with different ploidy levels, H3K27me2 intensities increased concurrent with increased ploidy levels, and H3K27me2 peaks were colocalized with massively amplified DTC transposons (CACTA family) in euchromatin, which may silence euchromatic transposons to maintain genome stability during polyploid wheat evolution. Consistently, the distribution of H3K27me2 is mutually exclusive with another repressive histone mark, H3K9me2, that mainly silences transposons in heterochromatic regions. Remarkably, the regions with low H3K27me2 levels (named H3K27me2 valleys) were associated with the formation of DNA double-strand breaks in genomes of wheat, maize (Zea mays) and Arabidopsis. Our results provide a comprehensive view of H3K27me2 and H3K27me3 distributions during wheat evolution, which support roles for H3K27me2 in silencing euchromatic transposons to maintain genome stability and in modifying genetic recombination landscapes. These genomic insights may empower breeding improvement of crops.

New evidence confirming the CD genomic constitutions of the tetraploid Avena species in the section Pachycarpa Baum

The tetraploid Avena species in the section Pachycarpa Baum, including A. insularis, A. maroccana, and A. murphyi, are thought to be involved in the evolution of hexaploid oats; however, their genome designations are still being debated. Repetitive DNA sequences play an important role in genome structuring and evolution, so understanding the chromosomal organization and distribution of these sequences in Avena species could provide valuable information concerning genome evolution in this genus. In this study, the chromosomal organizations and distributions of six repetitive DNA sequences (including three SSR motifs (TTC, AAC, CAG), one 5S rRNA gene fragment, and two oat A and C genome specific repeats) were investigated using non-denaturing fluorescence in situ hybridization (ND-FISH) in the three tetraploid species mentioned above and in two hexaploid oat species. Preferential distribution of the SSRs in centromeric regions was seen in the A and D genomes, whereas few signals were detected in the C genomes. Some intergenomic translocations were observed in the tetraploids; such translocations were also detected between the C and D genomes in the hexaploids. These results provide robust evidence for the presence of the D genome in all three tetraploids, strongly suggesting that the genomic constitution of these species is DC and not AC, as had been thought previously.

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.

Changes in Alternative Splicing in Response to Domestication and Polyploidization in Wheat

Alternative splicing (AS) occurs extensively in eukaryotes as an important mechanism for regulating transcriptome complexity and proteome diversity, but variation in the AS landscape in response to domestication and polyploidization in crops is unclear. Hexaploid wheat (AABBDD, Triticum aestivum) has undergone two separate allopolyploidization events, providing an ideal model for studying AS changes during domestication and polyploidization events. In this study, we performed high-throughput transcriptome sequencing of roots and leaves from wheat species with varied ploidies, including wild diploids (A^bA^b, Triticum boeoticum) and tetraploids (AABB, Triticum dicoccoides), domesticated diploids (A^mA^m, Triticum monococcum) and tetraploids (AABB, Triticum dicoccum), hexaploid wheat (AABBDD, T aestivum), as well as newly synthesized hexaploids together with their parents. Approximately 22.1% of genes exhibited AS, with the major AS type being intron retention. The number of AS events decreased after domestication in both diploids and tetraploids. Moreover, the frequency of AS occurrence tended to decrease after polyploidization, consistent with the functional sharing model that proposes AS and duplicated genes are complementary in regulating transcriptome plasticity in polyploid crops. In addition, the subgenomes exhibited biased AS responses to polyploidization, and ∼87.1% of homeologs showed AS partitioning in hexaploid wheat. Interestingly, substitution of the D-subgenome modified 42.8% of AS patterns of the A- and B-subgenomes, indicating subgenome interplay reprograms AS profiles at a genome-wide level, although the causal-consequence relationship requires further study. Conclusively, our study shows that AS variation occurs extensively after polyploidization and domestication in wheat species.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Dynamic and reversible DNA methylation changes induced by genome separation and merger of polyploid wheat

Background

Wheat is a powerful genetic model for studying polyploid evolution and crop domestication. Hexaploid bread wheat was formed by two rounds of interspecific hybridization and polyploidization, processes which are often accompanied by genetic and epigenetic changes, including DNA methylation. However, the extent and effect of such changes during wheat evolution, particularly from tetraploid-to-hexaploid wheat, are currently elusive.

Results

Here we report genome-wide DNA methylation landscapes in extracted tetraploid wheat (ETW, AABB), natural hexaploid wheat (NHW, AABBDD), resynthesized hexaploid wheat (RHW, AABBDD), natural tetraploid wheat (NTW, AABB), and diploid (DD). In the endosperm, levels of DNA methylation, especially in CHG (H=A, T, or C) context, were dramatically decreased in the ETW relative to natural hexaploid wheat; hypo-differentially methylated regions (DMRs) (850,832) were 24-fold more than hyper-DMRs (35,111). Interestingly, those demethylated regions in ETW were remethylated in the resynthesized hexaploid wheat after the addition of the D genome. In ETW, hypo-DMRs correlated with gene expression, and TEs were demethylated and activated, which could be silenced in the hexaploid wheat. In NHW, groups of TEs were dispersed in genic regions of three subgenomes, which may regulate the expression of TE-associated genes. Further, hypo-DMRs in ETW were associated with reduced H3K9me2 levels and increased expression of histone variant genes, suggesting concerted epigenetic changes after separation from the hexaploid.

Conclusion

Genome merger and separation provoke dynamic and reversible changes in chromatin and DNA methylation. These changes correlate with altered gene expression and TE activity, which may provide insights into polyploid genome and wheat evolution.

Adaptive strategy of allohexaploid wheat to long-term salinity stress

BACKGROUND

Most studies of crop salinity tolerance are conducted under short-term stress condition within one growth stage. Understanding of the mechanisms of crop response to long-term salinity stress (LSS) is valuable for achieving the improvement of crop salinity tolerance. In the current study, we exposed allohexaploid wheat seeds to LSS conditions from germination stage to young seedling stage for 30 days. To elucidate the adaptive strategy of allohexaploid wheat to LSS, we analyzed chloroplast ultrastructure, leaf anatomy, transcriptomic profiling and concentrations of plant hormones and organic compatible solutes, comparing stressed and control plants.

RESULTS

Transcriptomic profiling and biochemical analysis showed that energy partitioning between general metabolism maintenance and stress response may be crucial for survival of allohexaploid wheat under LSS. Under LSS, wheat appeared to shift energy from general maintenance to stress response through stimulating the abscisic acid (ABA) pathway and suppressing gibberellin and jasmonic acid pathways in the leaf. We further distinguished the expression status of the A, B, and D homeologs of any gene triad, and also surveyed the effects of LSS on homeolog expression bias for salinity-tolerant triads. We found that LSS had similar effects on expression of the three homeologs for most salinity-tolerant triads. However, in some of these triads, LSS induced different effects on the expression of the three homeologs.

CONCLUSIONS

The shift of the energy from general maintenance to stress response may be important for wheat LSS tolerance. LSS influences homeolog expression bias of salinity-tolerant triads.

Multiple founder events explain the genetic diversity and structure of the model allopolyploid grass Brachypodium hybridum in the Iberian Peninsula hotspot

BACKGROUND AND AIMS

It is accepted that contemporary allopolyploid species have originated recurrently, but very few cases have been documented using multiple natural formations of the same species. To extend our knowledge, we have investigated the multiple origins, genetic variation and structure of the allotetraploid grass Brachypodium hybridum with respect to its progenitor diploid species B. distachyon (D genome) and B. stacei (S genome). For this, our primary focus is the Iberian Peninsula, an evolutionary hotspot for the genus Brachypodium.

METHODS

We analysed 342 B. hybridum individuals from 36 populations using ten nuclear SSR loci and two plastid loci. The B. hybridum genetic profiles were compared with those previously reported for B. stacei and B. distachyon. In addition, phylogenetic analysis of the plastid data was performed for a reduced subset of individuals.

KEY RESULTS

The nuclear simple sequence repeat (SSR) genetic analysis detected medium to high genetic diversity, with a strong south-to-north genetic structure cline, and a high selfing rate in B. hybridum. Comparative genetic analysis showed a close relatedness of current B. hybridum D allelic profiles with those of B. distachyon, but a lack of similarity with those of B. stacei, suggesting another B. stacei source for the B. hybridum S alleles. Plastid analysis detected three different bidirectional allopolyploidization events: two involved distinct B. distachyon-like ancestors and one involved a B. stacei-like ancestor. The south-eastern Iberian Peninsula B. hybridum populations were more genetically diverse and could have originated from at least two hybridization events whereas north-eastern/north-western Iberian Peninsula B. hybridum populations were less diverse and may have derived from at least one hybridization event.

CONCLUSIONS

The genetic and evolutionary evidence supports the plausible in situ origin of the south-eastern and northern Iberian Peninsula B. hybridum allopolyploids from their respective local B. distachyon and unknown B. stacei ancestors. The untapped multiple origins and genetic variation detected in these B. hybridum populations opens the way to future evolutionary analysis of allopolyploid formation and genomic dominance and expression in the B. hybridum-B. distachyon-B. stacei grass model complex.

Karyotype mosaicism in early generation synthetic hexaploid wheats

It is known that both the number and the structure of somatic chromosomes can vary in early generation hexaploid wheats. The phenomenon is generally assumed to arise as a result of the meiotic instability characteristic of freshly created allopolyploids. Here, an analysis of the somatic karyotype of a set of 33 early generation synthetic hexaploid wheats has revealed that variation, taking the form of sub-chromosomal fragments and inter-chromosomal translocations, can also arise in somatic tissue. A possible explanation for the observations was that karyotypic instability in early generation hexaploid wheat probably occurs not just during sporogenesis, but also in somatic tissue. However, other factors such as the use of nitrous oxide during the experiments could also cause the chromosome variations, and additional experimentation would be required to determine the most likely.

Gene Expression Changes During the Allo-/Deallopolyploidization Process of Brassica napus

Gene expression changes due to allopolyploidization have been extensively studied in plants over the past few decades. Nearly all these studies focused on comparing the changes before and after genome merger. In this study, we used the uniquely restituted Brassica rapa (RBR, A_eA_e, 2n = 20) obtained from Brassica napus (A_nA_nC_nC_n, 2n = 38) to analyze the gene expression changes and its potential mechanism during the process of allo-/deallopolyploidization. RNA-seq-based transcriptome profiling identified a large number of differentially expressed genes (DEGs) between RBR and natural B. rapa (A_rA_r), suggesting potential effects of allopolyploidization/domestication of AA component of B. napus at the tetrapolyploid level. Meanwhile, it was revealed that up to 20% of gene expressions were immediately altered when compared with those in the A_n-subgenome. Interestingly, one fifth of these changes are in fact indicative of the recovery of antecedent gene expression alternations occurring since the origin of B. napus and showed association with homoeologous expression bias between A_n and C_n subgenomes. Enrichment of distinct gene ontology (GO) categories of the above sets of genes further indicated potential functional cooperation of the A_n and C_n subgenome of B. napus. Whole genome methylation analysis revealed a small number of DEGs were identified in the differentially methylated regions.

An extracted tetraploid wheat harbouring the BBAA component of common wheat shows anomalous shikimate and sucrose metabolism

BACKGROUND

The BBAA subgenomes of hexaploid common wheat are structurally intact, which makes it possible to extract the BBAA subgenomes to constitute a novel plant type, namely, extracted tetraploid wheat (ETW). ETW displays multiple abnormal phenotypes such as massively reduced biomass and abnormal spike development, compared to extant tetraploid wheat with a BBAA genome. The genetic, biochemical and physiological basis underlying the phenotypic abnormality of ETW remains unknown.

RESULTS

To explore the biochemical basis of these phenotypic abnormalities, we analysed the metabolomic and proteomic profiles and quantified 46 physiological traits of ETW in comparison with its common wheat donor (genome BBAADD), and a durum tetraploid wheat cultivar (genome BBAA). Among these three types of wheat, ETW showed a saliently different pattern of nutrient accumulation and seed quality, markedly lower concentrations of many metabolites involved in carbohydrate metabolism, and higher concentrations of many metabolites related to amino acids. Among the metabolites, changes in shikimate and sucrose were the most conspicuous. Higher levels of shikimate and lower levels of sucrose influence many metabolic processes including carbohydrate and amino acid metabolism, which may contribute to the phenotypic abnormalities. Gene expression assay showed downregulation of a shikimate degradation enzyme (5-enolpyruvylshikimate-3-phosphate synthase) coding gene and upregulation of several genes coding for the sucrose hydrolysis enzyme, which could explain the higher levels of shikimate and lower levels of sucrose, respectively.

CONCLUSIONS

Our results suggest that significant and irreversible biochemical changes have occurred in the BBAA subgenomes of common wheat during the course of its co-evolution with the DD subgenome at the hexaploid level.

The Capacity to Buffer and Sustain Imbalanced D-Subgenome Chromosomes by the BBAA Component of Hexaploid Wheat Is an Evolved Dominant Trait

Successful generation of pentaploid wheat (genome, BBAAD) via interspecific hybridization between tetraploid wheat (BBAA) and hexaploid wheat (BBAADD) holds great promise to mutually exchange desirable traits between the two cultivated wheat species, as well as providing a novel facet for evolutionary studies of polyploid wheat. Taking advantage of the viable and fertile nature of an extracted tetraploid wheat (ETW) with a BBAA genome that is virtually identical with the BBAA component of a hexaploid common wheat, and a synthetic hexaploid wheat, we constructed four pentaploid wheats with several distinct yet complementary features, of which harboring homozygous BBAA subgenomes is a common feature. By using a combined FISH/GISH method that enables diagnosing all individual wheat chromosomes, we precisely karyotyped a larger number of cohorts from the immediate progenies of each of the four pentaploid wheats. We found that the BBAA component of hexaploid common wheat possesses a significantly stronger capacity to buffer and sustain imbalanced D genome chromosomes and appears to harbor more structural chromosome variations than the BBAA genome of tetraploid wheat. We also document that this stronger capacity of the hexaploid BBAA subgenomes behaves as a genetically controlled dominant trait. Our findings bear implications to the known greater than expected level of genetic diversity in, and the remarkable adaptability of, hexaploid common wheat as a staple crop of global significance, as well as in using pentaploidy as intermediates for reciprocal introgression of useful traits between tetraploid and hexaploid wheat cultivars.

Genome-Wide Survey and Comparative Analysis of Long Terminal Repeat (LTR) Retrotransposon Families in Four Gossypium Species

Long terminal repeat (LTR) retrotransposon is the most abundant DNA component and is largely responsible for plant genome size variation. Although it has been studied in plant species, very limited data is available for cotton, the most important fiber and texture crop. In this study, we performed a comprehensive analysis of LTR retrotransposon families across four cotton species. In tetraploid Gossypium species, LTR retrotransposon families from the progenitor D genome had more copies in D-subgenome, and families from the progenitor A genome had more copies in A-subgenome. Some LTR retrotransposon families that insert after polyploid formation may still distribute the majority of its copies in one of the subgenomes. The data also shows that families of 10~200 copies are abundant and they have a great influence on the Gossypium genome size; on the contrary, a small number of high copy LTR retrotransposon families have less contribution to the genome size. Kimura distance distribution indicates that high copy number family is not a recent outbreak, and there is no obvious relationship between family copy number and the period of evolution. Further analysis reveals that each LTR retrotransposon family may have their own distribution characteristics in cotton.

Asymmetrical changes of gene expression, small RNAs and chromatin in two resynthesized wheat allotetraploids

Polyploidy occurs in some animals and all flowering plants, including important crops such as wheat. The consequences of polyploidy in crops remain elusive, partly because their progenitors are unknown. Using two resynthesized wheat allotetraploids S^l S^l AA and AADD with known diploid progenitors, we analyzed mRNA and small RNA transcriptomes in the endosperm, compared transcriptomes between endosperm and root in AADD, and examined chromatin changes in the allotetraploids. In the endosperm, there were more non-additively expressed genes in S^l S^l AA than in AADD. In AADD, non-additively expressed genes were developmentally regulated, and the majority (62-70%) were repressed. The repressed genes in AADD included a group of histone methyltransferase gene homologs, which correlated with reduced histone H3K9me2 levels and activation of various transposable elements in AADD. In S^l S^l AA, there was a tendency for expression dominance of S^l over A homoeologs, but the histone methyltransferase gene homologs were additively expressed, correlating with insignificant changes in histone H3K9me2 levels. Moreover, more 24-nucleotide small inferring RNAs (siRNAs) in the A subgenome were disrupted in AADD than in S^l S^l AA, which were associated with expression changes of siRNA-associated genes. Our results indicate that asymmetrical changes in siRNAs, chromatin modifications, transposons and gene expression coincide with unstable AADD genomes and stable S^l S^l AA genomes, which could help explain the evolutionary trajectories of wheat allotetraploids formed by different progenitors.

Transcriptome analysis reveals potential mechanisms for different grain size between natural and resynthesized allohexaploid wheats with near-identical AABB genomes

BACKGROUND

Common wheat is a typical allohexaploid species (AABBDD) derived from the interspecific crossing between allotetraploid wheat (AABB) and Aegilops tauschii (DD). Wide variation in grain size and shape observed among Aegilops tauschii can be retained in synthetic allohexaploid wheats, but the underlying mechanism remains enigmatic. Here, the natural and resynthesized allohexaploid wheats with near-identical AB genomes and different D genomes (TAA10 and XX329) were employed for analysis.

RESULTS

Significant differences in grain size and weight between TAA10 and XX329 were observed at the early stages of development, which could be mainly attributed to the higher growth rates of the pericarp and endosperm cells in XX329 compared to TAA10. Furthermore, comparative transcriptome analysis identified that 8891 of 69,711 unigenes (12.75%) were differentially expressed between grains at 6 days after pollination (DAP) of TAA10 and XX329, including 5314 up-regulated and 3577 down-regulated genes in XX329 compared to TAA10. The MapMan functional annotation and enrichment analysis revealed that the differentially expressed genes were significantly enriched in categories of cell wall, carbohydrate and hormone metabolism. Notably, consistent with the up-regulation of sucrose synthase genes in resynthesized relative to natural allohexaploid wheat, the resynthesized allohexaploid wheat accumulated much higher contents of glucose and fructose in 6-DAP grains than those of the natural allohexaploid wheat.

CONCLUSIONS

These data indicated that the genetic variation of the D genome induced drastic alterations of gene expression in grains of the natural and resynthesized allohexaploid wheats, which may contribute to the observed differences in grain size and weight.

Evolution of the Aux/IAA Gene Family in Hexaploid Wheat

The Aux/IAA (IAA) gene family, involved in the auxin signalling pathway, acts as an important regulator in plant growth and development. In this study, we explored the evolutionary trajectory of the IAA family in common wheat. The results showed ten pairs of paralogs among 34 TaIAA family members. Seven of the pairs might have undergone segmental duplication, and the other three pairs appear to have experienced tandem duplication. Except for TaIAA15-16, these duplication events occurred in the ancestral genomes before the divergence of Triticeae. After that point, two polyploidization events shaped the current TaIAA family consisting of three subgenomic copies. The structure or expression pattern of the TaIAA family begins to differentiate in the hexaploid genome, where TaIAAs in the D genome lost more genes (eight) and protein secondary structures (α1, α3 and β5) than did the other two genomes. Expression analysis showed that six members of the TaIAA family were not expressed, and members such as TaIAA8, 15, 16, 28 and 33 exhibited tissue-specific expression patterns. In addition, three of the ten pairs of paralogs (TaIAA5-12, TaIAA15-16 and TaIAA29-30) showed similar expression patterns, and another five paralog pairs displayed differential expression patterns. Phylogenetic analysis showed that paralog pairs with high rates of evolution (ω > ω ₀), particularly TaIAA15-16 and TaIAA29-30, experienced greater motif loss, with only zero to two interacting IAA proteins. In contrast, most paralogous genes with low ω, such as TaIAA5-12, had more complete motifs and higher degrees of interaction with other family members.

Global Analysis of Gene Expression in Response to Whole-Chromosome Aneuploidy in Hexaploid Wheat

Aneuploidy, a condition of unbalanced chromosome content, represents a large-effect mutation that bears significant relevance to human health and microbe adaptation. As such, extensive studies of aneuploidy have been conducted in unicellular model organisms and cancer cells. Aneuploidy also frequently is associated with plant polyploidization, but its impact on gene expression and its relevance to polyploid genome evolution/functional innovation remain largely unknown. Here, we used a panel of diverse types of whole-chromosome aneuploidy of hexaploid wheat (Triticum aestivum), all under the common genetic background of cv Chinese Spring, to systemically investigate the impact of aneuploidy on genome-, subgenome-, and chromosome-wide gene expression. Compared with prior findings in haploid or diploid aneuploid systems, we unravel additional and novel features of alteration in global gene expression resulting from the two major impacts of aneuploidy, cis- and trans-regulation, as well as dosage compensation. We show that the expression-altered genes map evenly along each chromosome, with no evidence for coregulating aggregated expression domains. However, chromosomes and subgenomes in hexaploid wheat are unequal in their responses to aneuploidy with respect to the number of genes being dysregulated. Strikingly, homeologous chromosomes do not differ from nonhomologous chromosomes in terms of aneuploidy-induced trans-acting effects, suggesting that the three constituent subgenomes of hexaploid wheat are largely uncoupled at the transcriptional level of gene regulation. Together, our findings shed new insights into the functional interplay between homeologous chromosomes and interactions between subgenomes in hexaploid wheat, which bear implications to further our understanding of allopolyploid genome evolution and efforts in breeding new allopolyploid crops.

Revisiting Pivotal-Differential Genome Evolution in Wheat

An interesting and possibly unique pattern of genome evolution following polyploidy can be observed among allopolyploids of the Triticum and Aegilops genera (wheat group). Most polyploids in this group are presumed to share a common unaltered (pivotal) subgenome (U, D, or A) together with one or two modified (differential) subgenomes, a status that has been referred to as 'pivotal-differential' genome evolution. In this review we discuss various mechanisms that could be responsible for this evolutionary pattern, as well as evidence for and against the putative evolutionary mechanisms involved. We suggest that, in light of recent advances in genome sequencing and related technologies in the wheat group, the time has come to reopen the investigation into pivotal-differential genome evolution.

Functional Conservation and Divergence among Homoeologs of TaSPL20 and TaSPL21, Two SBP-Box Genes Governing Yield-Related Traits in Hexaploid Wheat

Maintaining high and stable yields has become an increasing challenge in wheat breeding due to climate change. Although Squamosa-promoter binding protein (SBP)-box genes have important roles in plant development, very little is known about the actual biological functions of wheat SBP-box family members. Here, we dissect the functional conservation, divergence, and exploitation of homoeologs of two paralogous TaSPL wheat loci during domestication and breeding. TaSPL20 and TaSPL21 were highly expressed in the lemma and palea. Ectopic expressions of TaSPL20/21 in rice exhibited similar functions in terms of promoting panicle branching but had different functions during seed development. We characterized all six TaSPL20/21 genes located across the three homoeologous (A, B, and D) genomes. According to the functional analysis of naturally occurring variants in 20 environments, four favorable haplotypes were identified. Together, they reduced plant height by up to 27.5%, and TaSPL21-6D-HapII increased 1000-grain weight by 9.73%. Our study suggests that TaSPL20 and TaSPL21 homoeologs underwent diversification in function with each evolving its own distinctive characteristics. During domestication and breeding of wheat in China, favorable haplotypes of each set were selected and exploited to varying degrees due to their large effects on plant height and 1000-grain weight.

Chromosomal structural changes and microsatellite variations in newly synthesized hexaploid wheat mediated by unreduced gametes

Allohexaploid wheat was derived from interspecific hybridization, followed by spontaneous chromosome doubling. Newly synthesized hexaploid wheat by crossing Triticum turgidum and Aegilops tauschii provides a classical model to understand the mechanisms of allohexaploidization in wheat. However, immediate chromosome level variation and microsatellite level variation of newly synthesized hexaploid wheat have been rarely reported. Here, unreduced gametes were applied to develop synthesized hexaploid wheat, NA0928, population by crossing T. turgidum ssp. dicoccum MY3478 and Ae. tauschii SY41, and further S0-S3 generations of NA0928 were assayed by sequential cytological and microsatellite techniques. We demonstrated that plentiful chromosomal structural changes and microsatellite variations emerged in the early generations of newly synthesized hexaploid wheat population NA0928, including aneuploidy with whole-chromosome loss or gain, aneuploidy with telosome formation, chromosome-specific repeated sequence elimination (indicated by fluorescence in situ hybridization) and microsatellite sequence elimination (indicated by sequencing), and many kinds of variations have not been previously reported. Additionally, we reported a new germplasm, T. turgidum accession MY3478 with excellent unreduced gametes trait, and then succeeded to transfer powdery mildew resistance from Ae. tauschii SY41 to synthesized allohexaploid wheat population NA0928, which would be valuable resistance resources for wheat improvement.

The poor lonesome A subgenome of Brassica napus var. Darmor (AACC) may not survive without its mate

Constitutive genomes of allopolyploid species evolve throughout their life span. However, the consequences of long-term alterations on the interdependency between each original genome have not been established. Here, we attempted an approach corresponding to subgenome extraction from a previously sequenced natural allotetraploid, offering a unique opportunity to evaluate plant viability and structural evolution of one of its diploid components. We employed two different strategies to extract the diploid AA component of the Brassica napus variety 'Darmor' (AACC, 2n = 4x = 38) and we assessed the genomic structure of the latest AA plants obtained (after four to five rounds of selection), using a 60K single nucleotide polymorphism Illumina array. Only one strategy was successful and the diploid AA plants that were structurally characterized presented a lower proportion of the B. napus A subgenome extracted than expected. In addition, our analyses revealed that some genes lost in a polyploid context appeared to be compensated for plant survival, either by conservation of genomic regions from B. rapa, used in the initial cross, or by some introgressions from the B. napus C subgenome. We conclude that as little as c. 7500 yr of coevolution could lead to subgenome interdependency in the allotetraploid B. napus as a result of structural modifications.

Identification of QTL for Grain Size and Shape on the D Genome of Natural and Synthetic Allohexaploid Wheats with Near-Identical AABB Genomes

Grain size and shape associated with yield and milling quality are important traits in wheat domestication and breeding. To reveal the genetic factors on the D genome that control grain size and shape variation, we conducted analysis of quantitative trait loci (QTL) using the F₂ and F_2:3 populations derived from a common allohexaploid wheat line TAA10 and a synthetic allohexaploid wheat XX329, which have near-identical AABB genomes and different DD genomes. Based on genotyping using wheat 660K single nucleotide polymorphism (SNP) array, TAA10 and XX329 exhibited 96.55, 98.10, and 66.26% genetic similarities of A, B, and D genomes, respectively. Phenotypic evaluation revealed that XX329 had higher thousand grain weight (TGW), grain length, width, area and perimeter than TAA10 across all environments, and the grain yield per plot of XX329 increased by 17.43-30.36% compared with that of TAA10 in two environments. A total of nine environmentally stable QTL associated with grain size and shape were mapped on chromosomes 2D and 7D and verified using near isogenic lines (NILs), with the synthetic allohexaploid wheat XX329 contributing favorable alleles. Notably, a novel QTL QTgw.cau-2D controlling grain weight was first identified from the synthetic allohexaploid wheat, which may be a more desirable target for genetic improvement in wheat breeding. Collectively, these results provide further insights into the genetic factors that shaped the grain morphology during wheat evolution and domestication.

Extraction of the Constituent Subgenomes of the Natural Allopolyploid Rapeseed (Brassica napus L.)

As the dynamic nature of progenitor genomes accompanies the speciation by interspecific hybridization, the extraction of the constituent subgenome(s) from a natural allopolyploid species of long history and then restitution of the progenitor(s) provides the unique opportunity to study the genome evolution and interplay. Herein, the A subgenome from the allotetraploid oilseed rape (Brassica napus L., AACC) was extracted through inducing the preferential elimination of C-subgenome chromosomes in intertribal crosses and the progenitor B. rapa was restituted (RBR). Then by crossing and backcrossing RBR with B. napus donor, the C subgenome was in situ dissected by adding each of its nine chromosomes to the extracted A subgenome and establishing the whole set of monosonic alien addition lines (MAALs). RBR from spring-type B. napus genotype "Oro" expressed a phenotype resembling some type of B. rapa never observed before, but showed a winter-type flowering habit. This RBR had weaker growth vigor and suffered more seriously from biotic and abiotic stresses compared with Oro. The phenotypes specific for these MAALs showed the location of the related genes on the particular C-subgenome chromosomes. These MAALs exhibited obviously different frequencies in homeologous pairing and transmission of additional C-subgenome chromosomes, which were associated with the distinct degrees of their relatedness, and even with the possible genetic regulation for meiotic pairing evolved in B. napus Finally, large scaffolds undetermined for sequence assembly of B. napus were anchored to specific C-subgenome chromosomes using MAALs.

Molecular evolution of Wcor15 gene enhanced our understanding of the origin of A, B and D genomes in Triticum aestivum

The allohexaploid bread wheat originally derived from three closely related species with A, B and D genome. Although numerous studies were performed to elucidate its origin and phylogeny, no consensus conclusion has reached. In this study, we cloned and sequenced the genes Wcor15-2A, Wcor15-2B and Wcor15-2D in 23 diploid, 10 tetraploid and 106 hexaploid wheat varieties and analyzed their molecular evolution to reveal the origin of the A, B and D genome in Triticum aestivum. Comparative analyses of sequences in diploid, tetraploid and hexaploid wheats suggest that T. urartu, Ae. speltoides and Ae. tauschii subsp. strangulata are most likely the donors of the Wcor15-2A, Wcor15-2B and Wcor15-2D locus in common wheat, respectively. The Wcor15 genes from subgenomes A and D were very conservative without insertion and deletion of bases during evolution of diploid, tetraploid and hexaploid. Non-coding region of Wcor15-2B gene from B genome might mutate during the first polyploidization from Ae. speltoides to tetraploid wheat, however, no change has occurred for this gene during the second allopolyploidization from tetraploid to hexaploid. Comparison of the Wcor15 gene shed light on understanding of the origin of the A, B and D genome of common wheat.

Polyploidy: Pitfalls and paths to a paradigm

Investigators have long searched for a polyploidy paradigm-rules or principles that might be common following polyploidization (whole-genome duplication, WGD). Here we attempt to integrate what is known across the more thoroughly investigated polyploid systems on topics ranging from genetics to ecology. We found that while certain rules may govern gene retention and loss, systems vary in the prevalence of gene silencing vs. homeolog loss, chromosomal change, the presence of a dominant genome (in allopolyploids), and the relative importance of hybridization vs. genome doubling per se. In some lineages, aspects of polyploidization are repeated across multiple origins, but in other species multiple origins behave more stochastically in terms of genetic and phenotypic change. Our investigation also reveals that the path to synthesis is hindered by numerous gaps in our knowledge of even the best-known systems. Particularly concerning is the absence of linkage between genotype and phenotype. Moreover, most recent studies have focused on the genetic and genomic attributes of polyploidy, but rarely is there an ecological or physiological context. To promote a path to a polyploidy paradigm (or paradigms), we propose a major community goal over the next 10-20 yr to fill the gaps in our knowledge of well-studied polyploids. Before a meaningful synthesis is possible, more complete data sets are needed for comparison-systems that include comparable genetic, genomic, chromosomal, proteomic, as well as morphological, physiological, and ecological data. Also needed are more natural evolutionary model systems, as most of what we know about polyploidy continues to come from a few crop and genetic models, systems that often lack the ecological context inherent in natural systems and necessary for understanding the drivers of biodiversity.

Transcriptome asymmetry in synthetic and natural allotetraploid wheats, revealed by RNA-sequencing

Allopolyploidization has occurred frequently within the Triticum-Aegilops complex which provides a suitable system to investigate how allopolyploidization shapes the expression patterns of duplicated homeologs. We have conducted transcriptome-profiling of leaves and young inflorescences in wild and domesticated tetraploid wheats, Triticum turgidum ssp. dicoccoides (BBAA) and ssp. durum (BBAA), an extracted tetraploid (BBAA), and a synthetic tetraploid (S(l) S(l) AA) wheat together with its diploid parents, Aegilops longissima (S(l) S(l) ) and Triticum urartu (AA). The two diploid species showed tissue-specific differences in genome-wide ortholog expression, which plays an important role in transcriptome shock-mediated homeolog expression rewiring and hence transcriptome asymmetry in the synthetic tetraploid. Further changes of homeolog expression apparently occurred in natural tetraploid wheats, which led to novel transcriptome asymmetry between the two subgenomes. In particular, our results showed that extremely biased homeolog expression can occur rapidly upon the allotetraploidzation and this trend is further enhanced in the course of domestication and evolution of polyploid wheats. Our results suggest that allopolyploidization is accompanied by distinct phases of homeolog expression changes, with parental legacy playing major roles in the immediate rewiring of homeolog expression upon allopolyploidization, while evolution and domestication under allotetraploidy drive further homeolog-expression changes toward re-established subgenome expression asymmetry.

Altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat

Polyploidy is a major driving force in plant evolution and speciation. Phenotypic changes often arise with the formation, natural selection and domestication of polyploid plants. However, little is known about the consequence of hybridization and polyploidization on root hair development. Here, we report that root hair length of synthetic and natural allopolyploid wheats is significantly longer than those of their diploid progenitors, whereas no difference is observed between allohexaploid and allotetraploid wheats. The expression of wheat gene TaRSL4, an orthologue of AtRSL4 controlling the root hair development in Arabidopsis, was positively correlated with the root hair length in diploid and allotetraploid wheats. Moreover, transcript abundance of TaRSL4 homoeologue from A genome (TaRSL4-A) was much higher than those of other genomes in natural allopolyploid wheat. Notably, increased root hair length by overexpression of the TaRSL4-A in wheat led to enhanced shoot fresh biomass under nutrient-poor conditions. Our observations indicate that increased root hair length in allohexaploid wheat originated in the allotetraploid progenitors and altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat.

Genome-Wide Gene Expressions Respond Differently to A-subgenome Origins in Brassica napus Synthetic Hybrids and Natural Allotetraploid

The young allotetraploid Brassica napus (2n = 38, AACC) is one of models to study genomic responses to allopolyploidization. The extraction of AA component from natural B. napus and then restitution of progenitor B. rapa should provide a unique opportunity to reveal the genome interplay for gene expressions during the evolution. Herein, B. napus hybrids (2n = 19, AC) between the extracted and extant B. rapa (2n = 20, AA) and the same B. oleracea genotype (2n = 18, CC) were studied by RNA-seq and compared with natural B. napus donor, to reveal the gene expression changes from hybridization and domestication and the effects of A genome with different origins. Upon the initial merger of two diploid genomes, additive gene expression was prevalent in these two hybrids, for non-additively expressed genes only represented a small portion of total expressed genes. A high proportion of genes exhibited expression level dominance, with no preference to either of the parental genomes. Comparison of homoeolog expressions also showed no bias toward any genomes and the parental expression patterns were often maintained in the hybrids and natural allotetraploids. Although, the overall patterns of gene expression were highly conserved between two hybrids, the extracted B. rapa responded less and appeared more compatible for hybridization than the extant B. rapa. Our results suggested that expression level dominance and homoeolog expressions bias were balanced at the initial stage of genome merger, and such balance were largely maintained during the domestication of B. napus, despite the increased extent over time.

The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres

Gossypium hirsutum contributes the most production of cotton fibre, but G. barbadense is valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome of G. barbadense has not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome of G. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade, and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, the CesA gene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism of CesA genes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a “relay race”-like fashion. We anticipate that the G. barbadense genome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

Making the Bread: Insights from Newly Synthesized Allohexaploid Wheat

Bread wheat (or common wheat, Triticum aestivum) is an allohexaploid (AABBDD, 2n = 6x = 42) that arose by hybridization between a cultivated tetraploid wheat T. turgidum (AABB, 2n = 4x = 28) and the wild goatgrass Aegilops tauschii (DD, 2n = 2x = 14). Polyploidization provided niches for rigorous genome modification at cytogenetic, genetic, and epigenetic levels, rendering a broader spread than its progenitors. This review summarizes the latest advances in understanding gene regulation mechanisms in newly synthesized allohexaploid wheat and possible correlation with polyploid growth vigor and adaptation. Cytogenetic studies reveal persistent association of whole-chromosome aneuploidy with nascent allopolyploids, in contrast to the genetic stability in common wheat. Transcriptome analysis of the euploid wheat shows that small RNAs are driving forces for homoeo-allele expression regulation via genetic and epigenetic mechanisms. The ensuing non-additively expressed genes and those with expression level dominance to the respective progenitor may play distinct functions in growth vigor and adaptation in nascent allohexaploid wheat. Further genetic diploidization of allohexaploid wheat is not random. Regional asymmetrical gene distribution, rather than subgenome dominance, is observed in both synthetic and natural allohexaploid wheats. The combinatorial effects of diverged genomes, subsequent selection of specific gene categories, and subgenome-specific traits are essential for the successful establishment of common wheat.