Intrinsic karyotype stability and gene copy number variations may have laid the foundation for tetraploid wheat formation

Genomic prediction of synthetic hexaploid wheat upon tetraploid durum and diploid Aegilops parental pools

Bread wheat (Triticum aestivum L.) is a globally important food crop, which was domesticated about 8-10,000 years ago. Bread wheat is an allopolyploid, and it evolved from two hybridization events of three species. To widen the genetic base in breeding, bread wheat has been re-synthesized by crossing durum wheat (Triticum turgidum ssp. durum) and goat grass (Aegilops tauschii Coss), leading to so-called synthetic hexaploid wheat (SHW). We applied the quantitative genetics tools of "hybrid prediction"-originally developed for the prediction of wheat hybrids generated from different heterotic groups - to a situation of allopolyploidization. Our use-case predicts the phenotypes of SHW for three quantitatively inherited global wheat diseases, namely tan spot (TS), septoria nodorum blotch (SNB), and spot blotch (SB). Our results revealed prediction abilities comparable to studies in 'traditional' elite or hybrid wheat. Prediction abilities were highest using a marker model and performing random cross-validation, predicting the performance of untested SHW (0.483 for SB to 0.730 for TS). When testing parents not necessarily used in SHW, combination prediction abilities were slightly lower (0.378 for SB to 0.718 for TS), yet still promising. Despite the limited phenotypic data, our results provide a general example for predictive models targeting an allopolyploidization event and a method that can guide the use of genetic resources available in gene banks.

Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus

BACKGROUND

Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown.

RESULTS

Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels.

CONCLUSIONS

Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids

KEY MESSAGE

We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

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.

Chromosome translocation affects multiple phenotypes, causes genome-wide dysregulation of gene expression, and remodels metabolome in hexaploid wheat

Chromosomal rearrangements (CRs) may occur in newly formed polyploids due to compromised meiotic fidelity. Moreover, CRs can be more readily tolerated in polyploids allowing their longer-term retention and hence potential spreading/fixation within a lineage. The direct functional consequences of CRs in plant polyploids remain unexplored. Here, we identified a heterozygous individual from a synthetic allohexaploid wheat in which the terminal parts of the long-arms of chromosomes 2D (approximately 193 Mb) and 4A (approximately 167 Mb) were reciprocally translocated. Five homogeneous translocation lines including both unbalanced and balanced types were developed by selfing fertilization of the founder mutant (RT [2DL; 4AL]-ter/1, reciprocal translocation). We investigated impacts of these translocations on phenotype, genome-wide gene expression and metabolome. We find that, compared with sibling wild-type, CRs in the form of both unbalanced and balanced translocations induced substantial changes of gene expression primarily via trans-regulation in the nascent allopolyploid wheat. The CRs also manifested clear phenotypic and metabolic consequences. In particular, the genetically balanced, stable reciprocal translocations lines showed immediate enhanced reproductive fitness relative to wild type. Our results underscore the profound impact of CRs on gene expression in nascent allopolyploids with wide-ranging phenotypic and metabolic consequences, suggesting CRs are an important source of genetic variation that can be exploited for crop breeding.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Homoeologous exchange enables rapid evolution of tolerance to salinity and hyper-osmotic stresses in a synthetic allotetraploid wheat

The link between polyploidy and enhanced adaptation to environmental stresses could be a result of polyploidy itself harbouring higher tolerance to adverse conditions, or polyploidy possessing higher evolvability than diploids under stress conditions. Natural polyploids are inherently unsuitable to disentangle these two possibilities. Using selfed progenies of a synthetic allotetraploid wheat AT3 (AADD) along with its diploid parents, Triticum urartu TMU38 (AA) and Aegilops tauschii TQ27 (DD), we addressed the foregoing issue under abiotic salinity and hyper-osmotic (drought-like) stress. Under short duration of both stresses, euploid plants of AT3 showed intermediate tolerance of diploid parents; under life-long duration of both stresses, tolerant individuals to either stress emerged from selfed progenies of AT3, but not from comparable-sized diploid parent populations. Tolerance to both stresses were conditioned by the same two homoeologous exchanges (HEs; 2DS/2AS and 3DL/3AL), and at least one HE needed to be at the homozygous state. Transcriptomic analyses revealed that hyper-up-regulation of within-HE stress responsive genes of the A sub-genome origin is likely responsible for the dual-stress tolerant phenotypes. Our results suggest that HE-mediated inter-sub-genome rearrangements can be an important mechanism leading to adaptive evolution in allopolyploids as well as a promising target for genetic manipulation in crop improvement.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Chromosomal instability and phenotypic variation in a specific lineage derived from a synthetic allotetraploid wheat

Newly formed plant allopolyploids usually have meiosis defect, resulting in chromosomal instability manifested as variation in chromosome number and/or structure. However, not all nascent allopolyploids are equally unstable. The wheat group (Aegilops/Triticum) contains 13 diploid species with distinct genome types. Many of these species can be artificially hybridized to produce viable but sterile inter-specific/intergeneric F1 hybrids, which can generate fertile synthetic allotetraploid wheats after whole genome doubling. Compared with synthetic allotetraploid wheats that contain genome combinations of AADD and S*S*DD (S* refers to related S genomes of a different species), those containing an S*S*AA genome are significantly more stable. However, robustness of the relative stability of S*S*AA genomes is unknown, nor are the phenotypic and fitness consequences during occurrences of secondary chromosomal instability. Here, we report a specific lineage originated from a single individual plant of a relatively stable synthetic allotetraploid wheat with genomes S ^l S ^l AA (S ^l and A subgenomes were from Ae. longissima and T. urartu, respectively) that showed a high degree of transgenerational chromosomal instability. Both numerical chromosome variation (NCV) and structural chromosome variation (SCV) occurred widely. While substantial differences in frequencies of both NCV and SCV were detected across the different chromosomes, only NCV frequencies were significantly different between the two subgenomes. We found that NCVs and SCVs occurred primarily due to perturbed meiosis, allowing formation of multivalents and univalents as well as homoeologous exchanges. Thus, the combination of NCVs and SCVs affected multiple phenotypic traits, particularly those related to reproductive fitness.

Salinity Tolerance in a Synthetic Allotetraploid Wheat (SlSlAA) Is Similar to Its Higher Tolerant Parent Aegilops longissima (SlSl) and Linked to Flavonoids Metabolism

Allotetraploidization between A and S (closely related to B) genome species led to the speciation of allotetraploid wheat (genome BBAA). However, the immediate metabolic outcomes and adaptive changes caused by the allotetraploidization event are poorly understood. Here, we investigated how allotetraploidization affected salinity tolerance using a synthetic allotetraploid wheat line (genome S^lS^lAA, labeled as 4x), its Aegilops longissima (genome S^lS^l, labeled as S^lS^l) and Triticum urartu (AA genome, labeled as AA) parents. We found that the degree of salinity tolerance of 4x was similar to its S^lS^l parent, and both were substantially more tolerant to salinity stress than AA. This suggests that the S^lS^l subgenome exerts a dominant effect for this trait in 4x. Compared with S^lS^l and 4x, the salinity-stressed AA plants did not accumulate a higher concentration of Na⁺ in leaves, but showed severe membrane peroxidation and accumulated a higher concentration of ROS (H₂O₂ and O₂ ^⋅⁣-) and a lesser concentration of flavonoids, indicating that ROS metabolism plays a key role in saline sensitivity. Exogenous flavonoid application to roots of AA plants significantly relieved salinity-caused injury. Our results suggest that the higher accumulation of flavonoids in S^lS^l may contribute to ROS scavenging and salinity tolerance, and these physiological properties were stably inherited by the nascent allotetraploid S^lS^lAA.

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.

Targeting Induced Local Lesions in the Wheat DEMETER and DRE2 Genes, Responsible for Transcriptional Derepression of Wheat Gluten Proteins in the Developing Endosperm

Wheat is a major source of energy and nutrition worldwide, but it is also a primary cause of frequent diet-induced health issues, specifically celiac disease, for which the only effective therapy so far is strict dietary abstinence from gluten-containing grains. Wheat gluten proteins are grouped into two major categories: high-molecular-weight glutenin subunits (HMWgs), vital for mixing and baking properties, and gliadins plus low-molecular-weight glutenin subunits (LMWgs) that contain the overwhelming majority of celiac-causing epitopes. We put forth a hypothesis that eliminating gliadins and LMWgs while retaining HMWgs might allow the development of reduced-immunogenicity wheat genotypes relevant to most gluten-sensitive individuals. This hypothesis stems from the knowledge that the molecular structures and regulatory mechanisms of the genes encoding the two groups of gluten proteins are quite different, and blocking one group's transcription, without affecting the other's, is possible. The genes for gliadins and LMWgs have to be de-methylated by 5-methylcytosine DNA glycosylase/lyase (DEMETER) and an iron-sulfur (Fe-S) cluster biogenesis enzyme (DRE2) early during endosperm development to permit their transcription. In this study, a TILLING (Targeting Induced Local Lesions IN Genomes) approach was undertaken to identify mutations in the homoeologous DEMETER (DME) and DRE2 genes in common and durum wheat. Lines with mutations in these genes were obtained that displayed reduced content of immunogenic gluten proteins while retaining essential baking properties. Although our data at first glance suggest new possibilities for treating celiac disease and are therefore of medical and agronomical interest, it also shows that inducing mutations in the DME and DRE2 genes analyzed here affected pollen viability and germination. Hence there is a need to develop other approaches in the future to overcome this undesired effect.

Genetic diversity, distribution and domestication history of the neglected GGAtAt genepool of wheat

We present a comprehensive survey of cytogenetic and genomic diversity of the GGA^tA^t genepool of wheat, thereby unlocking these plant genetic resources for wheat improvement. Wheat yields are stagnating around the world and new sources of genes for resistance or tolerances to abiotic traits are required. In this context, the tetraploid wheat wild relatives are among the key candidates for wheat improvement. Despite its potential huge value for wheat breeding, the tetraploid GGA^tA^t genepool is largely neglected. Understanding the population structure, native distribution range, intraspecific variation of the entire tetraploid GGA^tA^t genepool and its domestication history would further its use for wheat improvement. The paper provides the first comprehensive survey of genomic and cytogenetic diversity sampling the full breadth and depth of the tetraploid GGA^tA^t genepool. According to the results obtained, the extant GGA^tA^t genepool consists of three distinct lineages. We provide detailed insights into the cytogenetic composition of GGA^tA^t wheats, revealed group- and population-specific markers and show that chromosomal rearrangements play an important role in intraspecific diversity of T. araraticum. The origin and domestication history of the GGA^tA^t lineages is discussed in the context of state-of-the-art archaeobotanical finds. We shed new light on the complex evolutionary history of the GGA^tA^t wheat genepool and provide the basis for an increased use of the GGA^tA^t wheat genepool for wheat improvement. The findings have implications for our understanding of the origins of agriculture in southwest Asia.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Frequent numerical and structural chromosome changes in early generations of synthetic hexaploid wheat

Modern hexaploid wheat (Triticum aestivum L.; AABBDD) evolved from a hybrid of tetraploid wheat (closely related to Triticum turgidum L. ssp. durum (Desf.) Husn., AABB) and goatgrass (Aegilops tauschii Coss., DD). Variations in chromosome structure and ploidy played important roles in wheat evolution. How these variations occurred and their role in expanding the genetic diversity in modern wheat is mostly unknown. Synthetic hexaploid wheat (SHW) can be used to investigate chromosome variation that occurs during the early generations of existence. SHW lines derived by crossing durum wheat 'Langdon' with twelve Ae. tauschii accessions were analyzed using oligonucelotide probe multiplex fluorescence in situ hybridization (FISH) to metaphase chromosomes and SNP markers. Cluster analysis based on SNP markers categorized them into three groups. Among 702 plants from the S8 and S9 generations, 415 (59.12%) carried chromosome variations involving all 21 chromosomes but with different frequencies for each chromosome and sub-genome. Total chromosome variation frequencies varied between lines, but there was no significant difference among the three groups. The non-random chromosome variations in SHW lines detected in this research may be an indication that similar variations occurred in the early stages of wheat polyploidization and played important roles in wheat evolution.

High quality haplotype-resolved genome assemblies of Populus tomentosa Carr., a stabilized interspecific hybrid species widespread in Asia

Populus has a wide ecogeographical range spanning the Northern Hemisphere, and interspecific hybrids are common. Populus tomentosa Carr. is widely distributed and cultivated in the eastern region of Asia, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. Reference genomes are available for several Populus species, however, our goals were to produce a very high quality de novo chromosome-level genome assembly in P. tomentosa genome that could serve as a reference for evolutionary and ecological studies of hybrid speciation throughout the genus. Here, combining long-read sequencing and Hi-C scaffolding, we present a high-quality, haplotype-resolved genome assembly. The genome size was 740.2 Mb, with a contig N50 size of 5.47 Mb and a scaffold N50 size of 46.68 Mb, consisting of 38 chromosomes, as expected with the known diploid chromosome number (2n = 2x = 38). A total of 59,124 protein-coding genes were identified. Phylogenomic analyses revealed that P. tomentosa is comprised of two distinct subgenomes, which we deomonstrate is likely to have resulted from hybridization between Populus adenopoda as the female parent and Populus alba var. pyramidalis as the male parent, with an origin of approximately 3.93 Ma. Although highly colinear, significant structural variation was found between the two subgenomes. Our study provides a valuable resource for ecological genetics and forest biotechnology.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

High Chromosomal Stability and Immortalized Totipotency Characterize Long-Term Tissue Cultures of Chinese Ginseng (Panax Ginseng)

Chinese ginseng (Panax ginseng C. A. Meyer) is a highly cherished traditional Chinese medicine, with several confirmed medical effects and many more asserted health-boosting functions. Somatic chromosomal instability (CIN) is a hallmark of many types of human cancers and also related to other pathogenic conditions such as miscarriages and intellectual disabilities, hence, the study of this phenomenon is of wide scientific and translational medical significance. CIN also ubiquitously occurs in cultured plant cells, and is implicated as a major cause of the rapid decline/loss of totipotency with culture duration, which represents a major hindrance to the application of transgenic technologies in crop improvement. Here, we report two salient features of long-term cultured callus cells of ginseng, i.e., high chromosomal stability and virtually immortalized totipotency. Specifically, we document that our callus of ginseng, which has been subcultured for 12 consecutive years, remained highly stable at the chromosomal level and showed little decline in totipotency. We show that these remarkable features of cultured ginseng cells are likely relevant to the robust homeostasis of the transcriptional expression of specific genes (i.e., genes related to tissue totipotency and chromosomal stability) implicated in the manifestation of these two complex phenotypes. To our knowledge, these two properties of ginseng have not been observed in any animals (with respect to somatic chromosomal stability) and other plants. We posit that further exploration of the molecular mechanisms underlying these unique properties of ginseng, especially somatic chromosomal stability in protracted culture duration, may provide novel clues to the mechanistic understanding of the occurrence of CIN in human disease.

Evolutionary Contribution of Duplicated Genes to Genome Evolution in the Ginseng Species Complex

Genes duplicated by whole genome duplication (WGD) and small-scale duplication (SSD) have played important roles in adaptive evolution of all flowering plants. However, it still remains underinvestigated how the distinct models of duplication events and their contending evolutionary patterns have shaped the genome and epigenomes of extant plant species. In this study, we investigated the contribution of the WGD- and SSD-derived duplicate genes to the genome evolution of one diploid and three closely related allotetraploid Panax species based on genome, methylome, and proteome data sets. Our genome-wide comparative analyses revealed that although the ginseng species complex was recently diverged, they have evolved distinct overall patterns of nucleotide variation, cytosine methylation, and protein-level expression. In particular, genetic and epigenetic asymmetries observed in the recent WGD-derived genes are largely consistent across the ginseng species complex. In addition, our results revealed that gene duplicates generated by ancient WGD and SSD mechanisms exhibited distinct evolutionary patterns. We found the ancient WGD-derived genes (i.e., ancient collinear gene) are genetically more conserved and hypomethylated at the cytosine sites. In contrast, some of the SSD-derived genes (i.e., dispersal duplicated gene) showed hypermethylation and high variance in nucleotide variation pattern. Functional enrichment analyses of the duplicated genes indicated that adaptation-related traits (i.e., photosynthesis) created during the distant ancient WGDs are further strengthened by both the more recent WGD and SSD. Together, our findings suggest that different types of duplicated genes may have played distinct but relaying evolutionary roles in the polyploidization and speciation processes in the ginseng species complex.

Cytological and molecular characterizations of a novel 2A nullisomic line derived from a widely-grown wheat cultivar Zhoumai 18 conferring male sterility

A dwarf, multi-pistil and male sterile dms mutant was previously reported by us. However, the genetic changes in this dms are unclear. To examine the genetic changes, single nucleotide polymorphism (SNP) association, chromosome counting, and high-resolution chromosome fluorescence in situ hybridization (FISH) techniques were employed. By comparing tall plants (T) with dwarf plants (D) in the offspring of dms mutant plants, SNP association analysis indicated that most SNPs were on chromosome 2A. There were three types in offspring of dms plants, with 42, 41 and 40 chromosomes respectively. High-resolution chromosome painting analysis demonstrated that T plants had all 42 wheat chromosomes; the medium plants (M) had 41 chromosomes, lacking one chromosome 2A; while D plants had 40 wheat chromosomes, and lacked both 2A chromosomes. These data demonstrated that dms resulted from a loss of chromosome 2A. We identified 23 genes on chromosome 2A which might be involved in the development of stamens or pollen grains. These results lay a solid foundation for further analysis of the molecular mechanisms of wheat male sterility. Because D plants can be used as a female parent to cross with other wheat genotypes, dms is a unique germplasm for any functional study of chromosome 2A and wheat breeding specifically targeting genes on 2A.

Illegitimate Recombination Between Homeologous Genes in Wheat Genome

Polyploidies produce a large number of duplicated regions and genes in genomes, which have a long-term impact and stimulate genetic innovation. The high similarity between homeologous chromosomes, forming different subgenomes, or homologous regions after genome repatterning, may permit illegitimate DNA recombination. Here, based on gene colinearity, we aligned the (sub)genomes of common wheat (Triticum aestivum, AABBDD genotype) and its relatives, including Triticum urartu (AA), Aegilops tauschii (DD), and T. turgidum ssp. dicoccoides (AABB) to detect the homeologous (paralogous or orthologous) colinear genes within and between (sub)genomes. Besides, we inferred more ancient paralogous regions produced by a much ancient grass-common tetraploidization. By comparing the sequence similarity between paralogous and orthologous genes, we assumed abnormality in the topology of constructed gene trees, which could be explained by gene conversion as a result of illegitimate recombination. We found large numbers of inferred converted genes (>2,000 gene pairs) suggested long-lasting genome instability of the hexaploid plant, and preferential donor roles by DD genes. Though illegitimate recombination was much restricted, duplicated genes produced by an ancient whole-genome duplication, which occurred millions of years ago, also showed evidence of likely gene conversion. As to biological function, we found that ~40% catalytic genes in colinearity, including those involved in starch biosynthesis, were likely affected by gene conversion. The present study will contribute to understanding the functional and structural innovation of the common wheat genome.

The effects of hybridization and genome doubling in plant evolution via allopolyploidy

Polyploidy is a pervasive and recurring phenomenon across the tree of life, which occurred at variable time scales, ecological amplitudes and cell types, and is especially prominent in the evolutionary histories of plants. Importantly, many of the world's most important crops and noxious invasive weeds are recent polyploids. Polyploidy includes two major types, autopolyploidy, referring to doubling of a single species genome, and allopolyploidy referring to doubling of two or more merged genomes via biological hybridization of distinct but related species. The prevalence of both types of polyploidy implies that both genome doubling alone and doubling coupled with hybridization confer selective advantages over their diploid progenitors under specific circumstances. In cases of allopolyploidy, the two events, genome doubling and hybridization, have both advantages and disadvantages. Accumulated studies have established that, in allopolyploidy, some advantage(s) of doubling may compensate for the disadvantage(s) of hybridity and vice versa, although further study is required to validate generality of this trend. Some studies have also revealed a variety of non-Mendelian genetic and genomic consequences induced by doubling and hybridization separately or concertedly in nascent allopolyploidy; however, the significance of which to the immediate establishment and longer-term evolutionary success of allopolyploid species remain to be empirically demonstrated and ecologically investigated. This review aims to summarize recent advances in our understanding of the roles of hybridization and genome doubling, in separation and combination, in the evolution of allopolyploid genomes, as well as fruitful future research directions that are emerging from these studies.

Homoeologous exchanges occur through intragenic recombination generating novel transcripts and proteins in wheat and other polyploids

Recombination between homeologous chromosomes, also known as homeologous exchange (HE), plays a significant role in shaping genome structure and gene expression in interspecific hybrids and allopolyploids of several plant species. However, the molecular mechanisms that govern HEs are not well understood. Here, we studied HE events in the progeny of a nascent allotetraploid (genome AADD) derived from two diploid progenitors of hexaploid bread wheat using cytological and whole-genome sequence analyses. In total, 37 HEs were identified and HE junctions were mapped precisely. HEs exhibit typical patterns of homologous recombination hotspots, being biased toward low-copy, subtelomeric regions of chromosome arms and showing association with known recombination hotspot motifs. But, strikingly, while homologous recombination preferentially takes place upstream and downstream of coding regions, HEs are highly enriched within gene bodies, giving rise to novel recombinant transcripts, which in turn are predicted to generate new protein fusion variants. To test whether this is a widespread phenomenon, a dataset of high-resolution HE junctions was analyzed for allopolyploid Brassica, rice, Arabidopsis suecica, banana, and peanut. Intragenic recombination and formation of chimeric genes was detected in HEs of all species and was prominent in most of them. HE thus provides a mechanism for evolutionary novelty in transcript and protein sequences in nascent allopolyploids.

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.

Uncovering homeologous relationships between tetraploid Agropyron cristatum and bread wheat genomes using COS markers

Using COS markers, the study reveals homeologous relationships between tetraploid Agropyron cristatum and bread wheat to support alien introgression breeding of wheat. Crested wheatgrass (Agropyron cristatum L. Gaertn.) is a wild relative of wheat that possesses many genes that are potentially useful in wheat improvement. The species comprises a complex of diploid, tetraploid and hexaploid forms. In this study, wheat-A. cristatum chromosome, telosome and translocation lines were used to characterize syntenic relationships between tetraploid A. cristatum and bread wheat. Prior to mapping COS markers, the cytogenetic stock lines were characterized for fertility and by FISH and GISH for karyotype stability. Out of 328 COS markers selected for the study, 279 consistently amplified products in tetraploid A. cristatum, and, out of these, 139 were polymorphic between tetraploid crested wheatgrass and wheat. Sixty-nine markers were found to be suitable for the detection of tetraploid A. cristatum chromosomes 1P-6P in wheat, ranging from 6 to 17 markers per chromosome. BLASTn of the source ESTs resulted in significant hits for 67 markers on the wheat pseudomolecules. Generally, COS markers of the same homeologous group were detected on similar arms in both Agropyron and wheat. However, some intragenomic duplications and chromosome rearrangements were detected in tetraploid A. cristatum. These results provide new insights into the structure and evolution of the tetraploid A. cristatum genome and will facilitate the exploitation of the wild species for introgression breeding of bread wheat.

Enriching LMW-GS alleles and strengthening gluten properties of common wheat through wide hybridization with wild emmer

Two advanced lines (BAd7-209 and BAd7-213) with identical high-molecular-weight glutenin subunit composition were obtained via wide hybridization between low-gluten cultivar chuannong16 (CN16) and wild emmer D97 (D97). BAd7-209 was better than BAd7-213, and both of them were much better than CN16 in a dough quality test. We found that BAd7-209 had more abundant and higher expression levels of low-molecular-weight glutenin subunit (LMW-GS) proteins than those of BAd7-213. Twenty-nine novel LMW-GS genes at Glu-A3 locus were isolated from BAd7-209, BAd7-213 and their parents. We found that all 29 LMW-GS genes possessed the same primary structure shared by other known LMW-GSs. Twenty-seven genes encode LMW-m-type subunits, and two encode LMW-i-type subunits. BAd7-209 had a higher number of LMW-GS genes than BAd7-213, CN16, and D97. Two wild emmer genes MG574329 and MG574330 were present in the two advanced lines. Most of the LMW-m-type genes showed minor nucleotide variations between wide hybrids and their parents that could be induced through the wide hybridization process. Our results demonstrated that the wild emmer LMW-GS alleles could be feasibly transferred and integrated into common wheat background via wide hybridization and the potential value of the wild emmer LMW-GS alleles in breeding programs designed to improve wheat flour quality.

Identification of P genome chromosomes in Agropyron cristatum and wheat-A. cristatum derivative lines by FISH

Agropyron cristatum (L.) Gaertn. (P genome) is cultivated as pasture fodder and can provide many desirable genes for wheat improvement. With the development of genomics and fluorescence in situ hybridization (FISH) technology, probes for identifying plant chromosomes were also developed. However, there are few reports on A. cristatum chromosomes. Here, FISH with the repeated sequences pAcTRT1 and pAcpCR2 enabled the identification of all diploid A. cristatum chromosomes. An integrated idiogram of A. cristatum chromosomes was constructed based on the FISH patterns of five diploid A. cristatum individuals. Structural polymorphisms of homologous chromosomes were observed not only among different individuals but also within individuals. Moreover, seventeen wheat-A. cristatum introgression lines containing different P genome chromosomes were identified with pAcTRT1 and pAcpCR2 probes. The arrangement of chromosomes in diploid A. cristatum was determined by identifying correspondence between the P chromosomes in these genetically identified introgression lines and diploid A. cristatum chromosomes. The two probes were also effective for discriminating all chromosomes of tetraploid A. cristatum, and the differences between two tetraploid A. cristatum accessions were similar to the polymorphisms among individuals of diploid A. cristatum. Collectively, the results provide an effective means for chromosome identification and phylogenetic studies of P genome chromosomes.

Transgenerationally Precipitated Meiotic Chromosome Instability Fuels Rapid Karyotypic Evolution and Phenotypic Diversity in an Artificially Constructed Allotetraploid Wheat (AADD)

Although a distinct karyotype with defined chromosome number and structure characterizes each biological species, it is intrinsically labile. Polyploidy or whole-genome duplication has played a pervasive and ongoing role in the evolution of all eukaryotes, and is the most dramatic force known to cause rapid karyotypic reconfiguration, especially at the initial stage. However, issues concerning transgenerational propagation of karyotypic heterogeneity and its translation to phenotypic diversity in nascent allopolyploidy, at the population level, have yet to be studied in detail. Here, we report a large-scale examination of transgenerationally propagated karyotypic heterogeneity and its phenotypic manifestation in an artificially constructed allotetraploid with a genome composition of AADD, that is, involving two of the three progenitor genomes of polyploid wheat. Specifically, we show that 1) massive organismal karyotypic heterogeneity is precipitated after 12 consecutive generations of selfing from a single euploid founder individual, 2) there exist dramatic differences in aptitudes between subgenomes and among chromosomes for whole-chromosome gain and/or loss and structural variations, 3) majority of the numerical and structural chromosomal variations are concurrent due to mutual contingency and possible functional constraint, 4) purposed and continuous selection and propagation for euploidy over generations did not result in enhanced karyotype stabilization, and 5) extent of karyotypic variation correlates with variability of phenotypic manifestation. Together, our results document that allopolyploidization catalyzes rampant and transgenerationally heritable organismal karyotypic heterogeneity that drives population-level phenotypic diversification, which lends fresh empirical support to the still contentious notion that whole-genome duplication enhances organismal evolvability.

The Agropyron cristatum karyotype, chromosome structure and cross-genome homoeology as revealed by fluorescence in situ hybridization with tandem repeats and wheat single-gene probes

Fluorescence in situ hybridization with probes for 45 cDNAs and five tandem repeats revealed homoeologous relationships of Agropyron cristatum with wheat. The results will contribute to alien gene introgression in wheat improvement. Crested wheatgrass (Agropyron cristatum L. Gaertn.) is a wild relative of wheat and a promising source of novel genes for wheat improvement. To date, identification of A. cristatum chromosomes has not been possible, and its molecular karyotype has not been available. Furthermore, homoeologous relationship between the genomes of A. cristatum and wheat has not been determined. To develop chromosome-specific landmarks, A. cristatum genomic DNA was sequenced, and new tandem repeats were discovered. Their distribution on mitotic chromosomes was studied by fluorescence in situ hybridization (FISH), which revealed specific patterns for five repeats in addition to 5S and 45S ribosomal DNA and rye subtelomeric repeats pSc119.2 and pSc200. FISH with one tandem repeat together with 45S rDNA enabled identification of all A. cristatum chromosomes. To analyze the structure and cross-species homoeology of A. cristatum chromosomes with wheat, probes for 45 mapped wheat cDNAs covering all seven chromosome groups were localized by FISH. Thirty-four cDNAs hybridized to homoeologous chromosomes of A. cristatum, nine hybridized to homoeologous and non-homoeologous chromosomes, and two hybridized to unique positions on non-homoeologous chromosomes. FISH using single-gene probes revealed that the wheat-A. cristatum collinearity was distorted, and important structural rearrangements were observed for chromosomes 2P, 4P, 5P, 6P and 7P. Chromosomal inversions were found for pericentric region of 4P and whole chromosome arm 6PL. Furthermore, reciprocal translocations between 2PS and 4PL were detected. These results provide new insights into the genome evolution within Triticeae and will facilitate the use of crested wheatgrass in alien gene introgression into wheat.

Characterization of an Integrated Active Glu-1Ay Allele in Common Wheat from Wild Emmer and Its Potential Role in Flour Improvement

Glu-1Ay, one of six genes encoding a high molecular weight glutenin subunit (HMW-GS), is frequently silenced in hexaploid common wheat. Here, an active allele of Glu-1Ay was integrated from wild emmer wheat (Triticum turgidum ssp. dicoccoides) accession D97 into the common wheat (Triticum aestivum) cultivar Chuannong 16 via the repeated self-fertilization of the pentaploid interspecific hybrid, culminating in the selection of a line TaAy7-40 shown to express the wild emmer Glu-1Ay allele. The open reading frame of this allele was a 1830 bp long sequence, demonstrated by its heterologous expression in Escherichia coli to encode a 608-residue polypeptide. Its nucleotide sequence was 99.2% identical to that of the sequence within the wild emmer parent. The TaAy7-40 introgression line containing the active Glu-1Ay allele showed higher protein content, higher sodium dodecyl sulfate (SDS) sedimentation value, higher content of wet gluten in the flour, higher grain weight, and bigger grain size than Chuannong 16. The end-use quality parameters of the TaAy7-40 were superior to those of the medium gluten common wheat cultivars Mianmai 37 and Neimai 9. Thus, the active Glu-1Ay allele might be of potential value in breeding programs designed to improve wheat flour quality.

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.

Validation of an updated Associative Transcriptomics platform for the polyploid crop species Brassica napus by dissection of the genetic architecture of erucic acid and tocopherol isoform variation in seeds

An updated platform was developed to underpin association genetics studies in the polyploid crop species Brassica napus (oilseed rape). Based on 1.92 × 10¹² bases of leaf mRNAseq data, functional genotypes, comprising 355 536 single-nucleotide polymorphism markers and transcript abundance were scored across a genetic diversity panel of 383 accessions using a transcriptome reference comprising 116 098 ordered coding DNA sequence (CDS) gene models. The use of the platform for Associative Transcriptomics was first tested by analysing the genetic architecture of variation in seed erucic acid content, as high-erucic rapeseed oil is highly valued for a variety of applications in industry. Known loci were identified, along with a previously undetected minor-effect locus. The platform was then used to analyse variation for the relative proportions of tocopherol (vitamin E) forms in seeds, and the validity of the most significant markers was assessed using a take-one-out approach. Furthermore, the analysis implicated expression variation of the gene Bo2g050970.1, an orthologue of VTE4 (which encodes a γ-tocopherol methyl transferase converting γ-tocopherol into α-tocopherol) associated with the observed trait variation. The establishment of the first full-scale Associative Transcriptomics platform for B. napus enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

Mixed-Ploidy Species: Progress and Opportunities in Polyploid Research

Mixed-ploidy species harbor a unique form of genomic and phenotypic variation that influences ecological interactions, facilitates genetic divergence, and offers insights into the mechanisms of polyploid evolution. However, there have been few attempts to synthesize this literature. We review here research on the cytotype distribution, diversity, and dynamics of intensively studied mixed-ploidy species and consider the implications for understanding mechanisms of polyploidization such as cytotype formation, establishment, coexistence, and post-polyploid divergence. In general, mixed-ploidy species are unevenly represented among families: they exhibit high cytotype diversity, often within populations, and frequently comprise rare and odd-numbered ploidies. Odd-ploidies often occur in association with asexuality. We highlight research hypotheses and opportunities that take advantage of the unique properties of ploidy variation.

Surviving a Genome Collision: Genomic Signatures of Allopolyploidization in the Recent Crop Species Brassica napus

Polyploidization has played a major role in crop plant evolution, leading to advantageous traits that have been selected by humans. Here, we describe restructuring patterns in the genome of L., a recent allopolyploid species. Widespread segmental deletions, duplications, and homeologous chromosome exchanges were identified in diverse genome sequences from 32 natural and 20 synthetic accessions, indicating that homeologous exchanges are a major driver of postpolyploidization genome diversification. Breakpoints of genomic rearrangements are rich in microsatellite sequences that are known to interact with the meiotic recombination machinery. In both synthetic and natural , a subgenome bias was observed toward exchanges replacing larger chromosome segments from the C-subgenome by their smaller, homeologous A-subgenome segments, driving postpolyploidization genome size reduction. Selection in natural favored segmental deletions involving genes associated with immunity, reproduction, and adaptation. Deletions affecting mismatch repair system genes, which are assumed to control homeologous recombination, were also found to be under selection. Structural exchanges between homeologous subgenomes appear to be a major source of novel genetic diversity in de novo allopolyploids. Documenting the consequences of genomic collision by genomic resequencing gives insights into the adaptive processes accompanying allopolyploidization.

Karyotype Stability and Unbiased Fractionation in the Paleo-Allotetraploid Cucurbita Genomes

The Cucurbita genus contains several economically important species in the Cucurbitaceae family. Here, we report high-quality genome sequences of C. maxima and C. moschata and provide evidence supporting an allotetraploidization event in Cucurbita. We are able to partition the genome into two homoeologous subgenomes based on different genetic distances to melon, cucumber, and watermelon in the Benincaseae tribe. We estimate that the two diploid progenitors successively diverged from Benincaseae around 31 and 26 million years ago (Mya), respectively, and the allotetraploidization happened at some point between 26 Mya and 3 Mya, the estimated date when C. maxima and C. moschata diverged. The subgenomes have largely maintained the chromosome structures of their diploid progenitors. Such long-term karyotype stability after polyploidization has not been commonly observed in plant polyploids. The two subgenomes have retained similar numbers of genes, and neither subgenome is globally dominant in gene expression. Allele-specific expression analysis in the C. maxima × C. moschata interspecific F₁ hybrid and their two parents indicates the predominance of trans-regulatory effects underlying expression divergence of the parents, and detects transgressive gene expression changes in the hybrid correlated with heterosis in important agronomic traits. Our study provides insights into polyploid genome evolution and valuable resources for genetic improvement of cucurbit crops.

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.

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.

Wheat paleohistory created asymmetrical genomic evolution

Following the triplication reported in Brassiceae ∼10million years ago, and at the basis of rosids ∼100million years ago, bias in organization and regulation, known as subgenome dominance, has been reported between the three post-polyploidy compartments referenced to as less fractionated (LF), medium fractionated (MF1) and more fractionated (MF2), that have been proposed to derive from an hexaploidization event involving ancestors of 7-14-21 chromosomes. Modern bread wheat experienced similar paleohistory during the last half million year of evolution opening a new hypothesis where the wheat genome is at the earliest stages on the road of diploidization through subgenome dominance driving asymmetry in gene content, gene expression abundance, transposable element content as dynamics and epigenetic control between the A, B and D subgenomes.

The Impact of Open Pollination on the Structural Evolutionary Dynamics, Meiotic Behavior, and Fertility of Resynthesized Allotetraploid Brassica napus L

Allopolyploidy, which results from the merger and duplication of two divergent genomes, has played a major role in the evolution and diversification of flowering plants. The genomic changes that occur in resynthesized or natural neopolyploids have been extensively studied, but little is known about the effects of the reproductive mode in the initial generations that may precede its successful establishment. To truly reflect the early generations of a nascent polyploid, two resynthesized allotetraploid Brassica napus populations were obtained for the first time by open pollination. In these populations, we detected a much lower level of aneuploidy (third generation) compared with those previously published populations obtained by controlled successive selfing. We specifically studied 33 resynthesized B. napus individuals from our two open pollinated populations, and showed that meiosis was affected in both populations. Their genomes were deeply shuffled after allopolyploidization: up to 8.5 and 3.5% of the C and A subgenomes were deleted in only two generations. The identified deletions occurred mainly at the distal part of the chromosome, and to a significantly greater extent on the C rather than the A subgenome. Using Fluorescent In Situ Hybridization (BAC-FISH), we demonstrated that four of these deletions corresponded to fixed translocations (via homeologous exchanges). We were able to evaluate the size of the structural variations and their impact on the whole genome size, gene content, and allelic diversity. In addition, the evolution of fertility was assessed, to better understand the difficulty encountered by novel polyploid individuals before the putative formation of a novel stable species.

Isolation and application of P genome-specific DNA sequences of Agropyron Gaertn. in Triticeae

Different types of P genome sequences and markers were developed, which could be used to analyze the evolution of P genome in Triticeae and identify precisely wheat- A. cristatum introgression lines. P genome of Agropyron Gaertn. plays an important role in Triticeae and could provide many desirable genes conferring high yield, disease resistance, and stress tolerance for wheat genetic improvement. Therefore, it is significant to develop specific sequences and functional markers of P genome. In this study, 126 sequences were isolated from the degenerate oligonucleotide primed-polymerase chain reaction (DOP-PCR) products of microdissected chromosome 6PS. Forty-eight sequences were identified as P genome-specific sequences by dot-blot hybridization and DNA sequences analysis. Among these sequences, 22 displayed the characteristics of retrotransposons, nine and one displayed the characteristics of DNA transposons and tandem repetitive sequence, respectively. Fourteen of 48 sequences were determined to distribute on different regions of P genome chromosomes by fluorescence in situ hybridization, and the distributing regions were as following: all over P genome chromosomes, centromeres, pericentromeric regions, distal regions, and terminal regions. We compared the P genome sequences with other genome sequences of Triticeae and found that the similar sequences of the P genome sequences were widespread in Triticeae, but differentiation occurred to various extents. Additionally, thirty-four molecular markers were developed from the P genome sequences, which could be used for analyzing the evolutionary relationship among 16 genomes of 18 species in Triticeae and identifying P genome chromatin in wheat-A. cristatum introgression lines. These results will not only facilitate the study of structure and evolution of P genome chromosomes, but also provide a rapid detecting tool for effective utilization of desirable genes of P genome in wheat improvement.

Sequencing and comparative analyses of Aegilops tauschii chromosome arm 3DS reveal rapid evolution of Triticeae genomes

Bread wheat (Triticum aestivum, AABBDD) is an allohexaploid species derived from two rounds of interspecific hybridizations. A high-quality genome sequence assembly of diploid Aegilops tauschii, the donor of the wheat D genome, will provide a useful platform to study polyploid wheat evolution. A combined approach of BAC pooling and next-generation sequencing technology was employed to sequence the minimum tiling path (MTP) of 3176 BAC clones from the short arm of Ae. tauschii chromosome 3 (At3DS). The final assembly of 135 super-scaffolds with an N50 of 4.2 Mb was used to build a 247-Mb pseudomolecule with a total of 2222 predicted protein-coding genes. Compared with the orthologous regions of rice, Brachypodium, and sorghum, At3DS contains 38.67% more genes. In comparison to At3DS, the short arm sequence of wheat chromosome 3B (Ta3BS) is 95-Mb large in size, which is primarily due to the expansion of the non-centromeric region, suggesting that transposable element (TE) bursts in Ta3B likely occurred there. Also, the size increase is accompanied by a proportional increase in gene number in Ta3BS. We found that in the sequence of short arm of wheat chromosome 3D (Ta3DS), there was only less than 0.27% gene loss compared to At3DS. Our study reveals divergent evolution of grass genomes and provides new insights into sequence changes in the polyploid wheat genome.

Cytological identification of an Aegilops variabilis chromosome carrying stripe rust resistance in wheat

Aegilops variabilis (UUSvSv), an important sources for wheat improvement, originated from chromosome doubling of a natural hybrid between Ae. umbellulata (UU) with Ae. longissima (SlSl). The Ae. variabilis karyotype was poorly characterized by fluorescent in situ hybridization (FISH). The FISH probe combination of pSc119.2, pTa71 and pTa-713 identified each of the 14 pairs of Ae. variabilis chromosomes. Our FISH ideogram was further used to detect an Ae. variabilis chromosome carrying stripe rust resistance in the background of wheat lines developed from crosses of the stripe rust susceptible bread wheat cultivar Yiyuan 2 with a resistant Ae. variabilis accession. Among the 15 resistant BC1F7 lines, three were 2Sv + 4Sv addition lines (2n = 46) and 12 were 2Sv(2B) or 2Sv(2D) substitution lines that were confirmed with SSR markers. SSR marker gwm148 can be used to trace 2Sv in common wheat background. Chromosome 2Sv probably carries gametocidal(Gc) gene(s) since cytological instability and chromosome structural variations, including non-homologous translocations, were observed in some lines with this chromosome. Due to the effects of photoperiod genes, substitution lines 2Sv(2D) and 2Sv(2B) exhibited late heading with 2Sv(2D) lines being later than 2Sv(2B) lines. 2Sv(2D) substitution lines were also taller and exhibited higher spikelet numbers and longer spikes.

How diploidization turned a tetraploid into a pseudotriploid

PREMISE OF THE STUDY

Despite being highly fertile and occupying a large geographic region, the North American heartleaf bittercress (Cardamine cordifolia; Brassicaceae) has a puzzling triploid-like chromosome number (2n = 3x = 24). As most triploids are sterile, we embarked on a detailed analysis of the C. cordifolia genome to elucidate its origin and structure.

METHODS

Mitotic and meiotic chromosome complement of C. cordifolia was analyzed by comparative chromosome painting using chromosome-specific BAC contigs of Arabidopsis thaliana. Resulting chromosome patterns were documented by multicolor fluorescence microscopy and compared with known ancestral and extant Brassicaceae genomes.

KEY RESULTS

We discovered that C. cordifolia is not a triploid hybrid but a diploidized tetraploid with the prevalence of regular, diploid-like meiotic pairing. The ancestral tetraploid chromosome number (2n = 32) was reduced to a triploid-like number (2n = 24) through four terminal chromosome translocations.

CONCLUSIONS

The structure of the pseudotriploid C. cordifolia genome results from a stepwise diploidization process after whole-genome duplication. We showed that translocation-based descending dysploidy (from n = 16 to n = 12) was mediated by the formation of five new chromosomes. The genome of C. cordifolia represents the diploidization process in statu nascendi and provides valuable insights into mechanisms of postpolyploidy rediploidization in land plants. Our data further suggest that chromosome number alone does not need to be a reliable proxy of species' evolutionary past and that the same chromosome number may originate either by polyploidization (hybridization) or due to descending dysploidy.

Transcriptome shock invokes disruption of parental expression-conserved genes in tetraploid wheat

Allopolyploidy often triggers phenotypic novelty and gene expression remolding in the resulting polyploids. In this study, we employed multiple phenotypic and genetic approaches to investigate the nature and consequences of allotetraploidization between A- and S-subgenome of tetraploid wheat. Results showed that karyotype of the nascent allopolyploid plants (AT2) is stable but they showed clear novelty in multiple morphological traits which might have positively contributed to the initial establishment of the tetraploids. Further microarray-based transcriptome profiling and gene-specific cDNA-pyrosequencing have documented that transcriptome shock was exceptionally strong in AT2, but a substantial proportion of the induced expression changes was rapidly stabilized in early generations. Meanwhile, both additive and nonadditive expression genes showed extensive homeolog expression remodeling and which have led to the subgenome expression dominance in leaf and young inflorescence of AT2. Through comparing the homeolog-expressing patterns between synthetic and natural tetraploid wheats, it appears that the shock-induced expression changes at both the total expression level and subgenome homeolog partitioning are evolutionarily persistent. Together, our study shed new light on how gene expression changes have rapidly occurred at the initial stage following allotetraploidization, as well as their evolutionary relevance, which may have implications for wheat improvements.

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.

Allopolyploidy and the complex phylogenetic relationships within the Hordeum brachyantherum taxon

Hordeum brachyantherum Nevski includes two subspecies: the diploid (2×) subsp. californicum, and subsp. brachyantherum, which itself includes a tetraploid (4×) and a hexaploid (6×) cytotype. The phylogenetic relationships between these taxa and the origin of the polyploids remain controversial. To provide additional information to the many molecular phylogenetic analyses conducted within Hordeum, FISH-based karyotypes were produced for all subspecies/cytotypes within H. brachyantherum. Chromosomes of H. roshevitzii and H. marinum subsp. gussoneanum were also analysed since these species are potentially involved in the origin of the polyploids. For karyotyping, ten repetitive DNA sequences were screened to indentify repeats showing sufficient diversity in terms of copy number and localisation that they might serve as physical markers for distinguishing between each mitotic chromosome pair in all accessions. Genomic in situ hybridisation (GISH) was used to distinguish between subgenomes in polyploids. The karyotype maps allowed the assessment of the chromosomal diversity within species/cytotypes and the identification of possibly homoeologous chromosomes. The results show a wide divergence between the chromosomes of subsp. californicum and H. roshevitzii, and with their supposed derivatives in subsp. brachyantherum 4×. One of the three subgenomes of subsp. brachyantherum 6× is derived from subsp. gussoneanum with no genomic reorganisation (i.e., neither amplification nor loss of the repetitive DNA sequences analysed). It is generally accepted that subsp. brachyantherum 4× is the other progenitor of subsp. brachyantherum 6×, but the present results suggest this to be unlikely. The present findings thus show the cytogenetic diversity and genomic structure of H. brachyantherum, and reveal its complex evolutionary history, in which chromosomal diversification and allopolyploidy have played important roles.

Heritable alteration of DNA methylation induced by whole-chromosome aneuploidy in wheat

Aneuploidy causes changes in gene expression and phenotypes in all organisms studied. A previous study in the model plant Arabidopsis thaliana showed that aneuploidy-generated phenotypic changes can be inherited to euploid progenies and implicated an epigenetic underpinning of the heritable variations. Based on an analysis by amplified fragment length polymorphism and methylation-sensitive amplified fragment length polymorphism markers, we found that although genetic changes at the nucleotide sequence level were negligible, extensive changes in cytosine DNA methylation patterns occurred in all studied homeologous group 1 whole-chromosome aneuploid lines of common wheat (Triticum aestivum), with monosomic 1A showing the greatest amount of methylation changes. The changed methylation patterns were inherited by euploid progenies derived from the aneuploid parents. The aneuploidy-induced DNA methylation alterations and their heritability were verified at selected loci by bisulfite sequencing. Our data have provided empirical evidence supporting earlier suggestions that heritability of aneuploidy-generated, but aneuploidy-independent, phenotypic variations may have an epigenetic basis. That at least one type of aneuploidy - monosomic 1A - was able to cause significant epigenetic divergence of the aneuploid plants and their euploid progenies also lends support to recent suggestions that aneuploidy may have played an important and protracted role in polyploid genome evolution.

Genome wide identification of C1-2i zinc finger proteins and their response to abiotic stress in hexaploid wheat

The C1-2i wheat Q-type C2H2 zinc finger protein (ZFP) transcription factor subclass has been reported to play important roles in plant stress responses. This subclass of ZFPs has not been studied in hexaploid wheat (Triticum aestivum) and we aimed to identify all members of this subclass and evaluate their responses to different abiotic stresses causing oxidative stress. Exploiting the recently published wheat draft genome sequence, we identified 53 members (including homoeologs from A, B and D genomes) of the C1-2i wheat Q-type C2H2 ZFPs (TaZFPs) representing 21 genes. Evolution analysis revealed that 9 TaZFPs members are directly inherited from the parents Triticum urartu and Aegilops tauschii, while 15 diverged through neoploidization events. This TaZFP subclass is responsive to the oxidative stress generator H2O2 and to high light, drought stress and flooding. Most TaZFPs are responsive to H2O2 (37/53), high light (44/53), flooding (31/53) or drought (37/53); 32 TaZFPs were up-regulated by at least 3 stresses and 16 were responsive to all stresses tested. A large number of these TaZFPs were physically mapped on different wheat draft genome sequences with known markers useful for QTL mapping. Our results show that the C1-2i subclass of TaZFPs is associated with responses to different abiotic stresses and that most TaZFPs (30/53 or 57 %) are located on group 5 chromosomes known to be involved in environment adaptation. Detailed characterization of these novel wheat TaZFPs and their association to QTL or eQTL may help to design wheat cultivars with improved tolerance to abiotic stress.