Allopolyploid Species Research Articles

Wheat speciation and adaptation: perspectives from reticulate evolution.

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

Imprints of independent allopolyploid formations on patterns of gene expression in two sibling yarrow species (Achillea, Asteraceae)

BackgroundPolyploid species often originate recurrently. While this is well known, there is little information on the extent to which distinct allotetraploid species formed from the same parent species differ in gene expression. The tetraploid yarrow species Achillea alpina and A. wilsoniana arose independently from allopolyploidization between diploid A. acuminata and A. asiatica. The genetics and geography of these origins are clear from previous studies, providing a solid basis for comparing gene expression patterns of sibling allopolyploid species that arose independently.ResultsWe conducted comparative RNA-sequencing analyses on the two Achillea tetraploid species and their diploid progenitors to evaluate: 1) species-specific gene expression and coexpression across the four species; 2) patterns of inheritance of parental gene expression; 3) parental contributions to gene expression in the allotetraploid species, and homeolog expression bias. Diploid A. asiatica showed a higher contribution than diploid A. acuminata to the transcriptomes of both tetraploids and also greater homeolog bias in these transcriptomes, possibly reflecting a maternal effect. Comparing expressed genes in the two allotetraploids, we found expression of ca. 30% genes were species-specific in each, which were most enriched for GO terms pertaining to “defense response”. Despite species-specific and differentially expressed genes between the two allotetraploids, they display similar transcriptome changes in comparison to their diploid progenitors.ConclusionTwo independently originated Achillea allotetraploid species exhibited difference in gene expression, some of which must be related to differential adaptation during their post-speciation evolution. On the other hand, they showed similar expression profiles when compared to their progenitors. This similarity might be expected when pairs of merged diploid genomes in tetraploids are similar, as is the case in these two particular allotetraploids.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids.

Allopolyploids have globally higher fitness than their diploid progenitors; however, by comparison, most resynthesized allopolyploids have poor fertility and highly unstable genome. Elucidating the evolutionary processes promoting genome stabilization and fertility is thus essential to comprehend allopolyploid success. Using the Brassica model, we mimicked the speciation process of a nascent allopolyploid species by resynthesizing allotetraploid Brassica napus and systematically selecting for euploid individuals over eight generations in four independent allopolyploidization events with contrasted genetic backgrounds, cytoplasmic donors, and polyploid formation type. We evaluated the evolution of meiotic behavior and fertility and identified rearrangements in S1 to S9 lineages to explore the positive consequences of euploid selection on B.napus genome stability. Recurrent selection of euploid plants for eight generations drastically reduced the percentage of aneuploid progenies as early as the fourth generation, concomitantly with a decrease in number of newly fixed homoeologous rearrangements. The consequences of homoeologous rearrangements on meiotic behavior and seed number depended strongly on the genetic background and cytoplasm donor. The combined use of both self-fertilization and recurrent euploid selection allowed identification of genomic regions associated with fertility and meiotic behavior, providing complementary evidence to explain B.napus speciation success.

Stable, fertile lines produced by hybridization between allotetraploids Brassica juncea (AABB) and Brassica carinata (BBCC) have merged the A and C genomes.

Many flowering plant taxa contain allopolyploids that share one or more genomes in common. In the Brassica genus, crop species Brassica juncea and Brassica carinata share the B genome, with 2n=AABB and 2n=BBCC genome complements, respectively. Hybridization results in 2n=BBAC hybrids, but the fate of these hybrids over generations of self-pollination has never been reported. We produced and characterized B.juncea×B. carinata (2n=BBAC) interspecific hybrids over six generations of self-pollination under selection for high fertility using a combination of genotyping, fertility phenotyping, and cytogenetics techniques. Meiotic pairing behaviour improved from 68% bivalents in the F1 to 98% in the S5 /S6 generations, and initially low hybrid fertility also increased to parent species levels. The S5 /S6 hybrids contained an intact B genome (16 chromosomes) plus a new, stable A/C genome (18-20 chromosomes) resulting from recombination and restructuring of A and C-genome chromosomes. Our results provide the first experimental evidence that two genomes can come together to form a new, restructured genome in hybridization events between two allotetraploid species that share a common genome. This mechanism should be considered in interpreting phylogenies in taxa with multiple allopolyploid species.

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

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

Genome evolution in Oryza allopolyploids of various ages: Insights into the process of diploidization.

The prevalence and recurrence of whole-genome duplication in plants and its major role in evolution have been well recognized. Despite great efforts, many aspects of genome evolution, particularly the temporal progression of genomic responses to allopolyploidy and the underlying mechanisms, remain poorly understood. The rice genus Oryza consists of both recently formed and older allopolyploid species, representing an attractive system for studying the genome evolution after allopolyploidy. In this study, through screening BAC libraries and sequencing and annotating the targeted BAC clones, we generated orthologous genomic sequences surrounding the DEP1 locus, a major grain yield QTL in cultivated rice, from four Oryza polyploids of various ages and their likely diploid genome donors or close relatives. Based on sequenced DEP1 region and published data from three other genomic regions, we investigated the temporal evolutionary dynamics of four polyploid genomes at both genetic and expression levels. In the recently formed BBCC polyploid, Oryza minuta, genome dominance was not observed and its short-term responses to allopolyploidy are mainly manifested as a high proportion of homoeologous gene pairs showing unequal expression. This could partly be explained by parental legacy, rewiring of divergent regulatory networks and epigenetic modulation. Moreover, we detected an ongoing diploidization process in this genus, and suggest that the expression divergence driven by changes of selective constraint probably plays a big role in the long-term diploidization. These findings add novel insights into our understanding of genome evolution after allopolyploidy, and could facilitate crop improvements through hybridization and polyploidization.

Metal accumulation and its effect on leaf herbivory in an allopolyploid species Arabidopsis kamchatica inherited from a diploid hyperaccumulator A. halleri

AbstractExcessive amounts of metal ions in soil are toxic for most plant species, yet metal can also facilitate plant survival by elemental defense against herbivores and pathogens. Zinc and cadmium hyperaccumulation in Arabidopsis halleri is known to be effective for the defense against natural enemies. The allotetraploid species A. kamchatica, derived from A. halleri and a non‐hyperaccumulator A. lyrata, has a lower hyperaccumulation level of zinc than A. halleri, but its significance for elemental defense remains unknown. In this study, we evaluated the accumulation levels of zinc and cadmium in the allotetraploid compared with its diploid progenitors, and evaluated the contribution of metal treatments to anti‐herbivore resistance under field conditions. The accumulation level of zinc in A. kamchatica was intermediate between the progenitors, but that of cadmium was lower than in both diploid progenitors. The elemental defense of A. kamchatica and A. halleri was supported by a field experiment comparing the herbivory level between a control group and metal‐supplemented plants. Moreover, the effect of elemental defense was lower in A. kamchatica than in the hyperaccumulator progenitor A. halleri, which is consistent with the metal accumulation level. This result reveals that the allotetraploid plant inherited its hyperaccumulating ability from one progenitor as an advantageous trait but at an intermediate level.

DEFECTIVE ENDOSPERM-D1 (Dee-D1) is crucial for endosperm development in hexaploid wheat

Hexaploid wheat (Triticum aestivum L.) is a natural allopolyploid and provides a usable model system to better understand the genetic mechanisms that underlie allopolyploid speciation through the hybrid genome doubling. Here we aimed to identify the contribution of chromosome 1D in the development and evolution of hexaploid wheat. We identified and mapped a novel DEFECTIVE ENDOSPERM–D1 (Dee-D1) locus on 1DL that is involved in the genetic control of endosperm development. The absence of Dee-D1 leads to non-viable grains in distant crosses and alters grain shape, which negatively affects grain number and thousand-grain weight. Dee-D1 can be classified as speciation locus with a positive effect on the function of genes which are involved in endosperm development in hybrid genomes. The presence of Dee-D1 is necessary for the normal development of endosperm, and thus play an important role in the evolution and improvement of grain yield in hexaploid wheat.

Distinct nucleotide patterns among three subgenomes of bread wheat and their potential origins during domestication after allopolyploidization

BackgroundThe speciation and fast global domestication of bread wheat have made a great impact on three subgenomes of bread wheat. DNA base composition is an essential genome feature, which follows the individual-strand base equality rule and [AT]-increase pattern at the genome, chromosome, and polymorphic site levels among thousands of species. Systematic analyses on base compositions of bread wheat and its wild progenitors could facilitate further understanding of the evolutionary pattern of genome/subgenome-wide base composition of allopolyploid species and its potential causes.ResultsGenome/subgenome-wide base-composition patterns were investigated by using the data of polymorphic site in 93 accessions from worldwide populations of bread wheat, its diploid and tetraploid progenitors, and their corresponding reference genome sequences. Individual-strand base equality rule and [AT]-increase pattern remain in recently formed hexaploid species bread wheat at the genome, subgenome, chromosome, and polymorphic site levels. However, D subgenome showed the fastest [AT]-increase across polymorphic site from Aegilops tauschii to bread wheat than that on A and B subgenomes from wild emmer to bread wheat. The fastest [AT]-increase could be detected almost all chromosome windows on D subgenome, suggesting different mechanisms between D and other two subgenomes. Interestingly, the [AT]-increase is mainly contributed by intergenic regions at non-selective sweeps, especially the fastest [AT]-increase of D subgenome. Further transition frequency and sequence context analysis indicated that three subgenomes shared same mutation type, but D subgenome owns the highest mutation rate on high-frequency mutation type. The highest mutation rate on D subgenome was further confirmed by using a bread-wheat-private SNP set. The exploration of loci/genes related to the [AT] value of D subgenome suggests the fastest [AT]-increase of D subgenome could be involved in DNA repair systems distributed on three subgenomes of bread wheat.ConclusionsThe highest mutation rate is detected on D subgenome of bread wheat during domestication after allopolyploidization, leading to the fastest [AT]-increase pattern of D subgenome. The phenomenon may come from the joint action of multiple repair systems inherited from its wild progenitors.

The transition to flowering in winter rapeseed during vernalization.

Flowering time is a major determinant of adaptation, fitness and yield in the allopolyploid species rapeseed (Brassica napus). Despite being a close relative to Arabidopsis thaliana, little is known about the timing of floral transition and the genes that govern this process. Winter, semi-winter and spring type plants have important life history characteristics that differ in vernalization requirements for flowering and are important for growing rapeseed in different regions of the world. In this study, we investigated the timing of vernalization-driven floral transition in winter rapeseed and the effect of photoperiod and developmental age on flowering time and vernalization responsiveness. Microscopy and whole transcriptome analyses at the shoot apical meristems of plants grown under controlled conditions showed that floral transition is initiated within few weeks of vernalization. Certain Bna.SOC1 and Bna.SPL5 homeologs were among the induced genes, suggesting that they are regulating the timing of cold-induced floral transition. Moreover, the flowering response of plants with shorter pre-vernalization period correlated with a delayed expression of Bna.SOC1 and Bna.SPL5 genes. In essence, this study presents a detailed analysis of vernalization-driven floral transition and the aspects of juvenility and dormancy and their effect on flowering time in rapeseed.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms.

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

Genomic mosaicism due to homoeologous exchange generates extensive phenotypic diversity in nascent allopolyploids.

Allopolyploidy is an important process in plant speciation, yet newly formed allopolyploid species typically suffer from extreme genetic bottlenecks. One escape from this impasse might be homoeologous meiotic pairing, during which homoeologous exchanges (HEs) generate phenotypically variable progeny. However, the immediate genome-wide patterns and resulting phenotypic diversity generated by HEs remain largely unknown. Here, we analyzed the genome composition of 202 phenotyped euploid segmental allopolyploid individuals from the fourth selfed generation following chromosomal doubling of reciprocal F1 hybrids of crosses between rice subspecies, using whole-genome sequencing. We describe rampant occurrence of HEs that, by overcoming incompatibility or conferring superiority of hetero-cytonuclear interactions, generate extensive and individualized genomic mosaicism across the analyzed tetraploids. We show that the resulting homoeolog copy number alteration in tetraploids affects known-function genes and their complex genetic interactions, in the process creating extraordinary phenotypic diversity at the population level following a single initial hybridization. Our results illuminate the immediate genomic landscapes possible in a tetraploid genomic environment, and underscore HE as an important mechanism that fuels rapid phenotypic diversification accompanying the initial stages of allopolyploid evolution.

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

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

Natural polymorphisms in a pair of NSP2 homoeologs can cause loss of nodulation in peanut.

Microbial symbiosis in legumes is achieved through nitrogen-fixing root nodules, and these are important for sustainable agriculture. The molecular mechanisms underlying development of root nodules in polyploid legume crops are largely understudied. Through map-based cloning and QTL-seq approaches, we identified a pair of homoeologous GRAS transcription factor genes, Nodulation Signaling Pathway 2 (AhNSP2-B07 or Nb) and AhNSP2-A08 (Na), controlling nodulation in cultivated peanut (Arachis hypogaea L.), an allotetraploid legume crop, which exhibited non-Mendelian and Mendelian inheritance, respectively. The segregation of nodulation in the progeny of Nananbnb genotypes followed a 3:1 Mendelian ratio, in contrast to the 5:3~1:1 non-Mendelian ratio for nanaNbnb genotypes. Additionally, a much higher frequency of the nb allele (13%) than the na allele (4%) exists in the peanut germplasm collection, suggesting that Nb is less essential than Na in nodule organogenesis. Our findings reveal the genetic basis of naturally occurred non-nodulating peanut plants, which can be potentially used for nitrogen fixation improvement in peanut. Furthermore, the results have implications for and provide insights into the evolution of homoeologous genes in allopolyploid species.

Experimental and Field Data Support Range Expansion in an Allopolyploid Arabidopsis Owing to Parental Legacy of Heavy Metal Hyperaccumulation.

Empirical evidence is limited on whether allopolyploid species combine or merge parental adaptations to broaden habitats. The allopolyploid Arabidopsis kamchatica is a hybrid of the two diploid parents Arabidopsis halleri and Arabidopsis lyrata. A. halleri is a facultative heavy metal hyperaccumulator, and may be found in cadmium (Cd) and zinc (Zn) contaminated environments, as well as non-contaminated environments. A. lyrata is considered non-tolerant to these metals, but can be found in serpentine habitats. Therefore, the parents have adaptation to different environments. Here, we measured heavy metals in soils from native populations of A. kamchatica. We found that soil Zn concentration of nearly half of the sampled 40 sites was higher than the critical toxicity level. Many of the sites were near human construction, suggesting adaptation of A. kamchatica to artificially contaminated soils. Over half of the A. kamchatica populations had >1,000 μg g–1 Zn in leaf tissues. Using hydroponic treatments, most genotypes accumulated >3,000 μg g–1 Zn, with high variability among them, indicating substantial genetic variation in heavy metal accumulation. Genes involved in heavy metal hyperaccumulation showed an expression bias in the A. halleri-derived homeolog in widely distributed plant genotypes. We also found that two populations were found growing on serpentine soils. These data suggest that A. kamchatica can inhabit a range of both natural and artificial soil environments with high levels of ions that either of the parents specializes and that it can accumulate varying amount of heavy metals. Our field and experimental data provide a compelling example of combining genetic toolkits for soil adaptations to expand the habitat of an allopolyploid species.

Molecular Evolution and Expansion of the KUP Family in the Allopolyploid Cotton Species Gossypium hirsutum and Gossypium barbadense.

The comprehensive analysis of gene family evolution will elucidate the origin and evolution of gene families. The K+ uptake (KUP) gene family plays important roles in K+ uptake and transport, plant growth and development, and abiotic stress responses. However, the current understanding of the KUP family in cotton is limited. In this study, 51 and 53 KUPs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. These KUPs were divided into five KUP subfamilies, with subfamily 2 containing three groups. Different subfamilies had different member numbers, conserved motifs, gene structures, regulatory elements, and gene expansion and loss rates. A paleohexaploidization event caused the expansion of GhKUP and GbKUP in cotton, and duplication events in G. hirsutum and G. barbadense have happened in a common ancestor of Gossypium. Meanwhile, the KUP members of the two allopolyploid subgenomes of G. hirsutum and G. barbadense exhibited unequal gene proportions, gene structural diversity, uneven chromosomal distributions, asymmetric expansion rates, and biased gene loss rates. In addition, the KUP families of G. hirsutum and G. barbadense displayed evolutionary conservation and divergence. Taken together, these results illustrated the molecular evolution and expansion of the KUP family in allopolyploid cotton species.

The speciation and adaptation of the polyploids: a case study of the Chinese Isoetes L. diploid-polyploid complex

BackgroundThe Chinese Isoetes L. are distributed in a stairway pattern: diploids in the high altitude and polyploids in the low altitude. The allopolyploid I. sinensis and its diploid parents I. yunguiensis and I. taiwanensis is an ideal system with which to investigate the relationships between polyploid speciation and the ecological niches preferences.ResultsThere were two major clades in the nuclear phylogenetic tree, all of the populations of polyploid were simultaneously located in both clades. The chloroplast phylogenetic tree included two clades with different populations of the polyploid clustered with the diploids separately: I. yunguiensis with partial populations of the I. sinensis and I. taiwanensis with the rest populations of the I. sinensis. The crow node of the I. sinensis allopolyploid system was 4.43 Ma (95% HPD: 2.77–6.97 Ma). The divergence time between I. sinensis and I. taiwanensis was estimated to 0.65 Ma (95% HPD: 0.26–1.91 Ma). The narrower niche breadth in I.sinensis than those of its diploid progenitors and less niche overlap in the pairwise comparisons between the polyploid and its progenitors.ConclusionsOur results elucidate that I. yunguinensis and I. taiwanensis contribute to the speciation of I. sinensis, the diploid parents are the female parents of different populations. The change of altitude might have played an important role in allopolyploid speciation and the pattern of distribution of I. sinensis. Additionally, niche novelty of the allopolyploid population of I. sinensis has been detected, in accordance with the hypothesis that niche shift between the polyploids and its diploid progenitors is important for the establishment and persistence of the polyploids.

Long-read sequencing reveals widespread intragenic structural variants in a recent allopolyploid crop plant.

SummaryGenome structural variation (SV) contributes strongly to trait variation in eukaryotic species and may have an even higher functional significance than single‐nucleotide polymorphism (SNP). In recent years, there have been a number of studies associating large chromosomal scale SV ranging from hundreds of kilobases all the way up to a few megabases to key agronomic traits in plant genomes. However, there have been little or no efforts towards cataloguing small‐ (30–10 000 bp) to mid‐scale (10 000–30 000 bp) SV and their impact on evolution and adaptation‐related traits in plants. This might be attributed to complex and highly duplicated nature of plant genomes, which makes them difficult to assess using high‐throughput genome screening methods. Here, we describe how long‐read sequencing technologies can overcome this problem, revealing a surprisingly high level of widespread, small‐ to mid‐scale SV in a major allopolyploid crop species, Brassica napus. We found that up to 10% of all genes were affected by small‐ to mid‐scale SV events. Nearly half of these SV events ranged between 100 bp and 1000 bp, which makes them challenging to detect using short‐read Illumina sequencing. Examples demonstrating the contribution of such SV towards eco‐geographical adaptation and disease resistance in oilseed rape suggest that revisiting complex plant genomes using medium‐coverage long‐read sequencing might reveal unexpected levels of functional gene variation, with major implications for trait regulation and crop improvement.

Phylogenetic relationships and the maternal donor of Roegneria (Triticeae: Poaceae) based on three nuclear DNA sequences (ITS, Acc1, and Pgk1) and one chloroplast region (trnL‐F)

AbstractRoegneria is a polyploid perennial genus in the tribe Triticeae. Some species of Roegneria are morphologically similar to genus Elymus and have been classified in Elymus. To investigate the delimitation and phylogenetic relationships of Roegneria, nuclear (ITS, Acc1, and Pgk1) and chloroplast (trnL–trnF) DNA regions were sequenced for 38 allopolyploid species and 32 diploid species of Triticeae. Phylogenetic analyses of nuclear DNA revealed that all Roegneria species were included in the St and Y genome clades, and that the Y genome was closely related to the V and Xp genomes. The chloroplast DNA dataset showed that Roegneria species were grouped with Pseudoroegneria species. The Pseudoroegneria species from the Middle East (P. libanotica and P. tauri) and Central Asia (P. strigosa) were more closely related to Roegneria species. The results suggested that: (i) the species containing the St and Y genomes should be segregated from Elymus and treated as a distinct genus, Roegneria, based on the genomic constitution; (ii) P. libanotica, P. tauri, and/or P. strigosa potentially served as the maternal donor of the St genome in Roegneria; (iii) The Y genome of Roegneria originated from a diploid Y genome species, and the V and Xp genomes may have contributed to Y genome formation; (iv) among Roegneria species of previously uncertain genomic constitution, R. seriotina was tetraploid and possessed the StY genomes, E. calcicolus was hexaploid with the StYH genomic constitution and should be classified in Campeiostachys, R. glaucifolia possessed the StStY genomes, and R. tschimganica had the genomic constitution St1St2Y.

Parental legacy and regulatory novelty in Brachypodium diurnal transcriptomes accompanying their polyploidy

Polyploidy is a widespread phenomenon in eukaryotes that can lead to phenotypic novelty and has important implications for evolution and diversification. The modification of phenotypes in polyploids relative to their diploid progenitors may be associated with altered gene expression. However, it is largely unknown how interactions between duplicated genes affect their diurnal expression in allopolyploid species. In this study, we explored parental legacy and hybrid novelty in the transcriptomes of an allopolyploid species and its diploid progenitors. We compared the diurnal transcriptomes of representative Brachypodium cytotypes, including the allotetraploid Brachypodium hybridum and its diploid progenitors Brachypodium distachyon and Brachypodium stacei. We also artificially induced an autotetraploid B. distachyon. We identified patterns of homoeolog expression bias (HEB) across Brachypodium cytotypes and time-dependent gain and loss of HEB in B. hybridum. Furthermore, we established that many genes with diurnal expression experienced HEB, while their expression patterns and peak times were correlated between homoeologs in B. hybridum relative to B. distachyon and B. stacei, suggesting diurnal synchronization of homoeolog expression in B. hybridum. Our findings provide insight into the parental legacy and hybrid novelty associated with polyploidy in Brachypodium, and highlight the evolutionary consequences of diurnal transcriptional regulation that accompanied allopolyploidy.