Decrease in purifying selection pressures on wheat homoeologous genes: tetraploidization versus hexaploidization.

A series of polyploidizations in higher-order polyploids is the main event affecting gene content in a genome. Each polyploidization event can lead to massive functional divergence because of the subsequent decrease in selection pressure on duplicated genes; however, the causal relationship between multiple rounds of polyploidization and the functional divergence of duplicated genes is poorly understood. We focused on the Triticum-Aegilops complex lineage and compared selection pressure before and after tetraploidization and hexaploidization events. Although both events led to decreased selection pressure on homoeologous gene pairs (compared with diploids and tetraploids), the initial tetraploidization had a greater impact on selection pressure on homoeologous gene pairs than did subsequent hexaploidization. Consistent with this, selection pressure on expression patterns for the initial event relaxed more than those for the subsequent event. Surprisingly, the decreased selection pressure on these homoeologous genes was independent of the existence of in-paralogs within the same subgenome. Wheat homoeologous pairs had different evolutionary consequences compared with orthologs related to other mechanisms (ancient allopolyploidization, ancient autopolyploidization, and small-scale duplication). Furthermore, tetraploidization and hexaploidization also seemed to have different evolutionary consequences. This suggests that homoeologous genes retain unique functions, including functions that are unlikely to be preserved in genes generated by the other duplication mechanisms. We found that their unique functions differed between tetraploidization and hexaploidization (e.g., reproductive and chromosome segregation processes). These findings imply that the substantial number of gene pairs resulting from multiple allopolyploidization events, especially initial tetraploidization, may have been a unique source of functional divergence.

Complex but Clear Allopolyploid Pattern of Subtribe Tussilagininae (Asteraceae: Senecioneae) Revealed by Robust Phylogenomic Evidence, with Development of a Novel Homeolog-Sorting Pipeline.

Polyploidy is a significant mechanism in eukaryotic evolution and is particularly prevalent in the plant kingdom. However, our knowledge about this phenomenon and its effects on evolution remains limited. A major obstacle to the study of polyploidy is the great difficulty in untangling the origins of allopolyploids. Due to the drastic genome changes and the erosion of allopolyploidy signals caused by the combined effects of hybridization and complex post-polyploid diploidization processes, resolving the origins of allopolyploids has long been a challenging task. Here we revisit this issue with the interesting case of subtribe Tussilagininae (Asteraceae: Senecioneae) and by developing HomeoSorter, a new pipeline for network inferences by phasing homeologs to parental subgenomes. The pipeline is based on the basic idea of a previous study but with major changes to address the scaling problem and implement some new functions. With simulated data, we demonstrate that HomeoSorter works efficiently on genome-scale data and has high accuracy in identifying polyploid patterns and assigning homeologs. Using HomeoSorter, the maximum pseudo-likelihood model of Phylonet, and genome-scale data, we further address the complex origin of Tussilagininae, a speciose group (ca. 45 genera and 710 species) characterized by having high base chromosome numbers (mainly x = 30, 40). In particular, the inferred patterns are strongly supported by the chromosomal evidence. Tussilagininae is revealed to comprise two large groups with successive allopolyploid origins: Tussilagininae s.s. (mainly x = 30) and the Gynoxyoid group (x = 40). Two allopolyploidy events first give rise to Tussilagininae s.s., with the first event occurring between the ancestor of subtribe Senecioninae (x = 10) and a lineage (highly probably with x = 10) related to the Brachyglottis alliance, and the resulting hybrid lineage crossing with the ancestor of Chersodoma (x = 10) and leading to Tussilagininae s.s. Then, after early diversification, the Central American group (mainly x = 30) of Tussilagininae s.s., is involved in a third allopolyploidy event with, again, the Chersodoma lineage and produces the Gynoxyoid group. Our study highlights the value of HomeoSorter and the homeolog-sorting approach in polyploid phylogenetics. With rich species diversity and clear evolutionary patterns, Tussilagininae s.s. and the Gynoxyoid group are also excellent models for future investigations of polyploidy.

Deciphering peanut complex genomes paves a way to understand its origin and domestication.

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex.

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Nonadditive gene expression in polyploids.

Allopolyploidy involves hybridization and duplication of divergent parental genomes and provides new avenues for gene expression. The expression levels of duplicated genes in polyploids can show deviation from parental additivity (the arithmetic average of the parental expression levels). Nonadditive expression has been widely observed in diverse polyploids and comprises at least three possible scenarios: (a) The total gene expression level in a polyploid is similar to that of one of its parents (expression-level dominance); (b) total gene expression is lower or higher than in both parents (transgressive expression); and (c) the relative contribution of the parental copies (homeologs) to the total gene expression is unequal (homeolog expression bias). Several factors may result in expression nonadditivity in polyploids, including maternal-paternal influence, gene dosage balance, cis- and/or trans-regulatory networks, and epigenetic regulation. As our understanding of nonadditive gene expression in polyploids remains limited, a new generation of investigators should explore additional phenomena (i.e., alternative splicing) and use other high-throughput "omics" technologies to measure the impact of nonadditive expression on phenotype, proteome, and metabolome.

Homeolog expression quantification methods for allopolyploids.

Genome duplication with hybridization, or allopolyploidization, occurs in animals, fungi and plants, and is especially common in crop plants. There is an increasing interest in the study of allopolyploids because of advances in polyploid genome assembly; however, the high level of sequence similarity in duplicated gene copies (homeologs) poses many challenges. Here we compared standard RNA-seq expression quantification approaches used currently for diploid species against subgenome-classification approaches which maps reads to each subgenome separately. We examined mapping error using our previous and new RNA-seq data in which a subgenome is experimentally added (synthetic allotetraploid Arabidopsis kamchatica) or reduced (allohexaploid wheat Triticum aestivum versus extracted allotetraploid) as ground truth. The error rates in the two species were very similar. The standard approaches showed higher error rates (>10% using pseudo-alignment with Kallisto) while subgenome-classification approaches showed much lower error rates (<1% using EAGLE-RC, <2% using HomeoRoq). Although downstream analysis may partly mitigate mapping errors, the difference in methods was substantial in hexaploid wheat, where Kallisto appeared to have systematic differences relative to other methods. Only approximately half of the differentially expressed homeologs detected using Kallisto overlapped with those by any other method in wheat. In general, disagreement in low-expression genes was responsible for most of the discordance between methods, which is consistent with known biases in Kallisto. We also observed that there exist uncertainties in genome sequences and annotation which can affect each method differently. Overall, subgenome-classification approaches tend to perform better than standard approaches with EAGLE-RC having the highest precision.

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

BACKGROUND: Polyploid 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. RESULTS: We 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. CONCLUSION: Two 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.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species.

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Homeologs of Brassica SOC1, a central regulator of flowering time, are differentially regulated due to partitioning of evolutionarily conserved transcription factor binding sites in promoters.

Evolution of Brassica genome post-polyploidization reveals asymmetrical genome fractionation and copy number variation. Herein, we describe the impact of promoter divergence among SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) homeologs on expression and function in Brassica spp. SOC1, a regulated floral pathway integrator, is conserved as 3 redundant homeologs in diploid Brassicas. Even with high sequence identity within coding regions (92.8-100%), the spatio-temporal expression patterns of 9 SOC1 homologs in B. juncea and B. nigra indicates regulatory divergence. While LF and MF2 SOC1 homeologs are upregulated during floral transition, MF1 is barely expressed. Also, MF2 homeolog levels do not decline post-flowering, unlike LF. To investigate the underlying source of divergence, we analyzed the sequence and phylogeny of all reported (22) and isolated (21) upstream regions of Brassica SOC1. Full length upstream regions (4712-19189 bp) reveal 5 ubiquitously conserved ancestral Blocks, harboring binding sites of 18 TFs (TFBSs) characterized in Arabidopsis thaliana. The orthologs of these TFBSs are differentially conserved among Brassica SOC1 homeologs, imparting expression divergence. No crucial TFBSs are exclusively lost from LF_SOC1 promoter, while MF1_SOC1 has lost NF-Y binding site crucial for SOC1 activation by CONSTANS. MF2_SOC1 homeologs have lost important TFBSs (SEP3, AP1 and SMZ), responsible for SOC1 repression post-flowering. BjuAALF_SOC1 promoter (proximal 2 kb) shows ubiquitous reporter expression in B. juncea cv. Varuna transgenics, while BjuAAMF1_SOC1 promoter shows absence of reporter expression, validating the impact of TFBS divergence. Conservation of the original primary protein sequence is discovered in B. rapa homeologs (46) of 18 TFs. Co-regulation pattern of these TFs appeared similar for B. rapa LF and MF2 SOC1 homeologs; MF1 shows significant variation. Strong regulatory association is recorded for AP1, AP2, SEP3, FLC and CONSTANS/NF-Y, highlighting their importance in homeolog-specific SOC1 regulation. Correlation of B. juncea AP1, AP2 and FLC expression with SOC1 homeologs also complies with the TFBS differences. We thus conclude that redundant SOC1 loci contribute differentially to cumulative expression of SOC1 due to divergent selection of ancestral TFBSs.

A likelihood ratio test for changes in homeolog expression bias.

BACKGROUND: Gene duplications are a major source of raw material for evolution and a likely contributor to the diversity of life on earth. Duplicate genes (i.e., homeologs, in the case of a whole genome duplication) may retain their ancestral function, sub- or neofunctionalize, or be lost entirely. A primary way that duplicate genes evolve new functions is by altering their expression patterns. Comparing the expression patterns of duplicate genes gives clues as to whether any of these evolutionary processes have occurred. RESULTS: We develop a likelihood ratio test for the analysis of the expression ratios of duplicate genes across two conditions (e.g., tissues). We demonstrate an application of this test by comparing homeolog expression patterns of 1448 homeologous gene pairs using RNA-seq data generated from leaves and petals of an allotetraploid monkeyflower (Mimulus luteus). We assess the sensitivity of this test to different levels of homeolog expression bias and compare the method to several alternatives. CONCLUSIONS: The likelihood ratio test derived here is a direct, transparent, and easily implemented method for detecting changes in homeolog expression bias that outperforms alternative approaches. While our method was derived with homeolog analysis in mind, this method can be used to analyze changes in the ratio of expression levels between any two genes in any two conditions.

Conservatism and variability of gene expression profiles among homeologous transcription factors in Xenopus laevis.

Xenopus laevis has an allotetraploid genome of 3.1Gb, in contrast to the diploid genome of a closely related species, Xenopus tropicalis. Here, we identified 412 genes (189 homeolog pairs, one homeologous gene cluster pair, and 28 singletons) encoding transcription factors (TFs) in the X. laevis genome by comparing them with their orthologs from X. tropicalis. Those genes include the homeobox gene family (Mix/Bix, Lhx, Nkx, Paired, POU, and Vent), Sox, Fox, Pax, Dmrt, Hes, GATA, T-box, and some clock genes. Most homeolog pairs for TFs are retained in two X. laevis subgenomes, named L and S, at higher than average rates (87.1% vs 60.2%). Among the 28 singletons, 82.1% were deleted from chromosomes of the S subgenome, a rate similar to the genome-wide average (82.1% vs 74.6%). Interestingly, nkx2-1, nkx2-8, and pax9, which reside consecutively in a postulated functional gene cluster, were deleted from the S chromosome, suggesting cluster-level gene regulation. Transcriptome correlation analysis demonstrated that TF homeolog pairs tend to have more conservative developmental expression profiles than most other types of genes. In some cases, however, either of the homeologs may show strongly different spatio-temporal expression patterns, suggesting neofunctionalization, subfunctionalization, or nonfunctionalization after allotetraploidization. Analyses of otx1 suggests that homeologs with much lower expression levels have undergone greater amino acid sequence diversification. Our comprehensive study implies that TF homeologs are highly conservative after allotetraploidization, possibly because the DNA sequences that they bind were also duplicated, but in some cases, they differed in expression levels or became singletons due to dosage-sensitive regulation of their target genes.

Sequence and functional characterization of MIRNA164 promoters from Brassica shows copy number dependent regulatory diversification among homeologs.

The impact of polyploidy on functional diversification of cis-regulatory elements is poorly understood. This is primarily on account of lack of well-defined structure of cis-elements and a universal regulatory code. To the best of our knowledge, this is the first report on characterization of sequence and functional diversification of paralogous and homeologous promoter elements associated with MIR164 from Brassica. The availability of whole genome sequence allowed us to identify and isolate a total of 42 homologous copies of MIR164 from diploid species-Brassica rapa (A-genome), Brassica nigra (B-genome), Brassica oleracea (C-genome), and allopolyploids-Brassica juncea (AB-genome), Brassica carinata (BC-genome) and Brassica napus (AC-genome). Additionally, we retrieved homologous sequences based on comparative genomics from Arabidopsis lyrata, Capsella rubella, and Thellungiella halophila, spanning ca. 45 million years of evolutionary history of Brassicaceae. Sequence comparison across Brassicaceae revealed lineage-, karyotype, species-, and sub-genome specific changes providing a snapshot of evolutionary dynamics of miRNA promoters in polyploids. Tree topology of cis-elements associated with MIR164 was found to re-capitulate the species and family evolutionary history. Phylogenetic shadowing identified transcription factor binding sites (TFBS) conserved across Brassicaceae, of which, some are already known as regulators of MIR164 expression. Some of the TFBS were found to be distributed in a sub-genome specific (e.g., SOX specific to promoter of MIR164c from MF2 sub-genome), lineage-specific (YABBY binding motif, specific to C. rubella in MIR164b), or species-specific (e.g., VOZ in A. thaliana MIR164a) manner which might contribute towards genetic and adaptive variation. Reporter activity driven by promoters associated with MIR164 paralogs and homeologs was majorly in agreement with known role of miR164 in leaf shaping, regulation of lateral root development and senescence, and one previously un-described novel role in trichome. The impact of polyploidy was most profound when reporter activity across three MIR164c homeologs were compared that revealed negligible overlap, whereas reporter activity among two homeologs of MIR164a displays significant overlap. A copy number dependent cis-regulatory divergence thus exists in MIR164 genes in Brassica juncea. The full extent of regulatory diversification towards adaptive strategies will only be known when future endeavors analyze the promoter function under duress of stress and hormonal regimes.

Conserved but Attenuated Parental Gene Expression in Allopolyploids: Constitutive Zinc Hyperaccumulation in the Allotetraploid Arabidopsis kamchatica.

Allopolyploidization combines parental genomes and often confers broader species distribution. However, little is known about parentally transmitted gene expression underlying quantitative traits following allopolyploidization because of the complexity of polyploid genomes. The allopolyploid species Arabidopsis kamchatica is a natural hybrid of the zinc hyperaccumulator Arabidopsis halleri and of the nonaccumulator Arabidopsis lyrata We found that A. kamchatica retained the ability to hyperaccumulate zinc from A. halleri and grows in soils with both low and high metal content. Hyperaccumulation of zinc by A. kamchatica was reduced to about half of A. halleri, but is 10-fold greater than A. lyrata Homeologs derived from A. halleri had significantly higher levels of expression of genes such as HEAVY METAL ATPASE4 (HMA4), METAL TRANSPORTER PROTEIN1 and other metal ion transporters than those derived from A. lyrata, which suggests cis-regulatory differences. A. kamchatica has on average about half the expression of these genes compared with A. halleri due to fixed heterozygosity inherent in allopolyploids. Zinc treatment significantly changed the ratios of expression of 1% of homeologous pairs, including genes putatively involved in metal homeostasis. Resequencing data showed a significant reduction in genetic diversity over a large genomic region (290 kb) surrounding the HMA4 locus derived from the A. halleri parent compared with the syntenic A. lyrata-derived region, which suggests different evolutionary histories. We also estimated that three A. halleri-derived HMA4 copies are present in A. kamchatica Our findings support a transcriptomic model in which environment-related transcriptional patterns of both parents are conserved but attenuated in the allopolyploids.

Comprehensive analyses of hox gene expression in Xenopus laevis embryos and adult tissues.

From whole genome sequencing of an allotetraploid frog, Xenopus laevis, two homeologous sets (L and S) of four Hox clusters A through D (HoxA.L/S, HoxB.L/S, HoxC.L/S, and HoxD.L/S) and 13 paralogous groups (PGs) with 76 genes were identified, allowing us to carry out the first comprehensive analyses of hox gene expression in vertebrates. Expression of all hox genes during development and in adult tissues was analyzed by RNA-sequencing. The expression levels of most hox genes were similar between homeologs, but in some pairs, large differences were observed and several of these were confirmed by RT-PCR and whole mount in situ hybridization experiments. These results indicate that subfunctionalization of hox genes may have occurred since allotetraploidization. Furthermore, comprehensive analysis of hox gene expression during early development did not agree with the hypothesis of temporal collinearity especially in genes belonging to PG2 to PG10.

Recurrent sequence exchange between homeologous grass chromosomes.

All grass species evolved from an ancestor that underwent a whole-genome duplication (WGD) approximately 70 million years ago. Interestingly, the short arms of rice chromosomes 11 and 12 (and independently their homologs in sorghum) were found to be much more similar to each other than other homeologous regions within the duplicated genome. Based on detailed analysis of rice chromosomes 11 and 12 and their homologs in seven grass species, we propose a mechanism that explains the apparently 'younger' age of the duplication in this region of the genome, assuming a small number of reciprocal translocations at the chromosome termini. In each case the translocations were followed by unbalanced transmission and subsequent lineage sorting of the involved chromosomes to offspring. Molecular dating of these translocation events also allowed us to date major chromosome 'fusions' in the evolutionary lineages that led to Brachypodium and Triticeae. Furthermore, we provide evidence that rice is exceptional regarding the evolution of chromosomes 11 and 12, inasmuch as in other species the process of sequence exchange between homeologous chromosomes ceased much earlier than in rice. We presume that random events rather than selective forces are responsible for the observed high similarity between the short arm ends of rice chromosomes 11 and 12.

Genome-wide identification and homeolog-specific expression analysis of the SnRK2 genes in Brassica napus guard cells.

Sucrose non-fermenting-1-related protein kinase 2 (SnRK2) proteins constitute a small plant-specific serine/threonine kinase family involved in abscisic acid (ABA) signaling and plant responses to biotic and abiotic stresses. Although SnRK2s have been well-studied in Arabidopsis thaliana, little is known about SnRK2s in Brassica napus. Here we identified 30 putative sequences encoding 10 SnRK2 proteins in the B. napus genome and the expression profiles of a subset of 14 SnRK2 genes in guard cells of B. napus. In agreement with its polyploid origin, B. napus maintains both homeologs from its diploid parents. The results of quantitative real-time PCR (qRT-PCR) and reanalysis of RNA-Seq data showed that certain BnSnRK2 genes were commonly expressed in leaf tissues in different varieties of B. napus. In particular, qRT-PCR results showed that 12 of the 14 BnSnRK2s responded to drought stress in leaves and in ABA-treated guard cells. Among them, BnSnRK2.4 and BnSnRK2.6 were of interest because of their robust responsiveness to ABA treatment and drought stress. Notably, BnSnRK2 genes exhibited up-regulation of different homeologs, particularly in response to abiotic stress. The homeolog expression bias in BnSnRK2 genes suggests that parental origin of genes might be responsible for efficient regulation of stress responses in polyploids. This work has laid a foundation for future functional characterization of the different BnSnKR2 homeologs in B. napus and its parents, especially their functions in guard cell signaling and stress responses.

Subgenome-specific assembly of vitamin E biosynthesis genes and expression patterns during seed development provide insight into the evolution of oat genome.

Vitamin E is essential for humans and thus must be a component of a healthy diet. Among the cereal grains, hexaploid oats (Avena sativa L.) have high vitamin E content. To date, no gene sequences in the vitamin E biosynthesis pathway have been reported for oats. Using deep sequencing and orthology-guided assembly, coding sequences of genes for each step in vitamin E synthesis in oats were reconstructed, including resolution of the sequences of homeologs. Three homeologs, presumably representing each of the three oat subgenomes, were identified for the main steps of the pathway. Partial sequences, likely representing pseudogenes, were recovered in some instances as well. Pairwise comparisons among homeologs revealed that two of the three putative subgenome-specific homeologs are almost identical for each gene. Synonymous substitution rates indicate the time of divergence of the two more similar subgenomes from the distinct one at 7.9-8.7 MYA, and a divergence between the similar subgenomes from a common ancestor 1.1 MYA. A new proposed evolutionary model for hexaploid oat formation is discussed. Homeolog-specific gene expression was quantified during oat seed development and compared with vitamin E accumulation. Homeolog expression largely appears to be similar for most of genes; however, for some genes, homoeolog-specific transcriptional bias was observed. The expression of HPPD, as well as certain homoeologs of VTE2 and VTE4, is highly correlated with seed vitamin E accumulation. Our findings expand our understanding of oat genome evolution and will assist efforts to modify vitamin E content and composition in oats.

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.

Transcriptome shock in an interspecific F1 triploid hybrid of Oryza revealed by RNA sequencing.

Interspecific hybridization is a driving force in evolution and speciation of higher plants. Interspecific hybridization often induces immediate and saltational changes in gene expression, a phenomenon collectively termed "transcriptome shock". Although transcriptome shock has been reported in various plant and animal taxa, the extent and pattern of shock-induced expression changes are often highly idiosyncratic, and hence entails additional investigations. Here, we produced a set of interspecific F1 triploid hybrid plants between Oryza sativa, ssp. japonica (2n = 2x = 24, genome AA) and the tetraploid form of O. punctata (2n = 4x = 48, genome, BBCC), and conducted RNA-seq transcriptome profiling of the hybrids and their exact parental plants. We analyzed both homeolog expression bias and overall gene expression level difference in the hybrids relative to the in silico "hybrids" (parental mixtures). We found that approximately 16% (2,541) of the 16,112 expressed genes in leaf tissue of the F1 hybrids showed nonadditive expression, which were specifically enriched in photosynthesis-related pathways. Interestingly, changes in the maternal homeolog expression, including non-stochastic silencing, were the major causes for altered homeolog expression partitioning in the F1 hybrids. Our findings have provided further insights into the transcriptome response to interspecific hybridization and heterosis.