Reciprocal silencing, transcriptional bias and functional divergence of homeologs in polyploid cotton (gossypium)

Biochemical and morpho-physiological insights revealed low moisture stress adaptation mechanisms in cotton (Gossypium hirsutum L.)

Cotton (Gossypium hirsutum L.) is a multipurpose crop. Abiotic stresses, especially extreme heat and drought, limit crop growth and thus reduce cotton yield by about 50%. In this study, 30 cotton genotypes were tested against low moisture stress in a pot experiment in triplicates along with control under wire house conditions. At the 3-4 leaf stage, different morpho-physiological and biochemical parameters were measured in order to select the low moisture stress-tolerant genotypes. For the selection of the best performing genotypes, Multi-Trait Genotype-Ideotype Distance Index (MGIDI) was used for the ranking of genotypes on the basis of multiple indices. For biochemical traits, 09 (TPC, TF, TSP, MDA, SOD, POD, CAT, APX, and Proline) out of 24 showed significant genotypic effects and were used for MGIDI. Eight genotypes (N-812 N-1296 N-696 N-377 N-121-896 N-T86, and N-3496) were observed to be best performing than others at 25% selection pressure (SI = 25%). For morpho-physiological traits, 14 out of 15 showed significant genotypic effects and used for MGIDI. Ten genotypes (N-1237 N-812 N-1296 N-696 N-9078 N-377 N-512 N-121 N-375, and N-896) were observed to be best performing at 35% selection pressure (SI = 35%). Six genotypes, i.e. N-812-1296 N-696 N-377 N-121, and N-896 were found common in both MGIDI analysis. In conclusion, three genotypes, i.e. N-696, N-896, and N-T86 proved to be most resilient to low moisture stress. Develop protocols, identified genotypes and markers that can be used for development of climate-smart cotton genotypes.

Integrative analyses of long and short-read RNA sequencing reveal the spliced isoform regulatory network of seedling growth dynamics in upland cotton

The polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh_A02G0645_N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome.

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

BnaABF3 and BnaMYB44 regulate the transcription of zeaxanthin epoxidase genes in carotenoid and abscisic acid biosynthesis

Zeaxanthin epoxidase (ZEP) is a key enzyme that catalyzes the conversion of zeaxanthin to violaxanthin in the carotenoid and abscisic acid (ABA) biosynthesis pathways. The rapeseed (Brassica napus) genome has 4 ZEP (BnaZEP) copies that are suspected to have undergone subfunctionalization, yet the 4 genes' underlying regulatory mechanisms remain unknown. Here, we genetically confirmed the functional divergence of the gene pairs BnaA09.ZEP/BnaC09.ZEP and BnaA07.ZEP/BnaC07.ZEP, which encode enzymes with tissue-specific roles in carotenoid and ABA biosynthesis in flowers and leaves, respectively. Molecular and transgenic experiments demonstrated that each BnaZEP pair is transcriptionally regulated via ABA-responsive element-binding factor 3 s (BnaABF3s) and BnaMYB44s as common and specific regulators, respectively. BnaABF3s directly bound to the promoters of all 4 BnaZEPs and activated their transcription, with overexpression of individual BnaABF3s inducing BnaZEP expression and ABA accumulation under drought stress. Conversely, loss of BnaABF3s function resulted in lower expression of several genes functioning in carotenoid and ABA metabolism and compromised drought tolerance. BnaMYB44s specifically targeted and repressed the expression of BnaA09.ZEP/BnaC09.ZEP but not BnaA07.ZEP/BnaC07.ZEP. Overexpression of BnaA07.MYB44 resulted in increased carotenoid content and an altered carotenoid profile in petals. Additionally, RNA-seq analysis indicated that BnaMYB44s functions as a repressor in phenylpropanoid and flavonoid biosynthesis. These findings provide clear evidence for the subfunctionalization of duplicated genes and contribute to our understanding of the complex regulatory network involved in carotenoid and ABA biosynthesis in B. napus.

An insight into the gene expression evolution in Gossypium species based on the leaf transcriptomes

BACKGROUND

Gene expression pattern is associated with biological phenotype and is widely used in exploring gene functions. Its evolution is also crucial in understanding species speciation and divergence. The genus Gossypium is a bona fide model for studying plant evolution and polyploidization. However, the evolution of gene expression during cotton species divergence has yet to be extensively discussed.

RESULTS

Based on the seedling leaf transcriptomes, this work analyzed the transcriptomic content and expression patterns across eight cotton species, including six diploids and two natural tetraploids. Our findings indicate that, while the biological function of these cotton transcriptomes remains largely conserved, there has been significant variation in transcriptomic content during species divergence. Furthermore, we conducted a comprehensive analysis of expression distances across cotton species. This analysis lends further support to the use of G. arboreum as a substitute for the A-genome donor of natural cotton polyploids. Moreover, our research highlights the evolution of stress-responsive pathways, including hormone signaling, fatty acid degradation, and flavonoid biosynthesis. These processes appear to have evolved under lower selection pressures, presumably reflecting their critical role in the adaptations of the studied cotton species to diverse environments.

CONCLUSIONS

In summary, this study provided insights into the gene expression variation within the genus Gossypium and identified essential genes/pathways whose expression evolution was closely associated with the evolution of cotton species. Furthermore, the method of characterizing genes and pathways under unexpected high or slow selection pressure can also serve as a new strategy for gene function exploration.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae)

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Subgenome evolutionary dynamics in allotetraploid ferns: insights from the gene expression patterns in the allotetraploid species Phegopteris decursivepinnata (Thelypteridacea, Polypodiales)

Allopolyploidization often leads to disruptive conflicts among more than two sets of subgenomes, leading to genomic modifications and changes in gene expression. Although the evolutionary trajectories of subgenomes in allopolyploids have been studied intensely in angiosperms, the dynamics of subgenome evolution remain poorly understood in ferns, despite the prevalence of allopolyploidization. In this study, we have focused on an allotetraploid fern-Phegopteris decursivepinnata-and its diploid parental species, P. koreana (K) and P. taiwaniana (T). Using RNA-seq analyses, we have compared the gene expression profiles for 9,540 genes among parental species, synthetic F1 hybrids, and natural allotetraploids. The changes in gene expression patterns were traced from the F1 hybrids to the natural allopolyploids. This study has revealed that the expression patterns observed in most genes in the F1 hybrids are largely conserved in the allopolyploids; however, there were substantial differences in certain genes between these groups. In the allopolyploids compared with the F1 hybrids, the number of genes showing a transgressive pattern in total expression levels was increased. There was a slight reduction in T-dominance and a slight increase in K-dominance, in terms of expression level dominance. Interestingly, there is no obvious bias toward the T- or K-subgenomes in the number and expression levels overall, showing the absence of subgenome dominance. These findings demonstrated the impacts of the substantial transcriptome change after hybridization and the moderate modification during allopolyploid establishment on gene expression in ferns and provided important insights into subgenome evolution in polyploid ferns.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Domestication over Speciation in Allopolyploid Cotton Species: A Stronger Transcriptomic Pull

Cotton has been domesticated independently four times for its fiber, but the genomic targets of selection during each domestication event are mostly unknown. Comparative analysis of the transcriptome during cotton fiber development in wild and cultivated materials holds promise for revealing how independent domestications led to the superficially similar modern cotton fiber phenotype in upland (G. hirsutum) and Pima (G. barbadense) cotton cultivars. Here we examined the fiber transcriptomes of both wild and domesticated G. hirsutum and G. barbadense to compare the effects of speciation versus domestication, performing differential gene expression analysis and coexpression network analysis at four developmental timepoints (5, 10, 15, or 20 days after flowering) spanning primary and secondary wall synthesis. These analyses revealed extensive differential expression between species, timepoints, domestication states, and particularly the intersection of domestication and species. Differential expression was higher when comparing domesticated accessions of the two species than between the wild, indicating that domestication had a greater impact on the transcriptome than speciation. Network analysis showed significant interspecific differences in coexpression network topology, module membership, and connectivity. Despite these differences, some modules or module functions were subject to parallel domestication in both species. Taken together, these results indicate that independent domestication led G. hirsutum and G. barbadense down unique pathways but that it also leveraged similar modules of coexpression to arrive at similar domesticated phenotypes.

A pan-grass transcriptome reveals patterns of cellular divergence in crops

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

Genome-wide characterization of ubiquitin-conjugating enzyme gene family explores its genetic effects on the oil content and yield of Brassica napus

Ubiquitin-conjugating enzyme (UBC) is a critical part of the ubiquitin-proteasome pathway and plays crucial roles in growth, development and abiotic stress response in plants. Although UBC genes have been detected in several plant species, characterization of this gene family at the whole-genome level has not been conducted in Brassica napus. In the present study, 200 putative BnUBCs were identified in B. napus, which were clustered into 18 subgroups based on phylogenetic analysis. BnUBCs within each subgroup showed relatively conserved gene architectures and motifs. Moreover, the gene expression patterns in various tissues as well as the identification of cis-acting regulatory elements in BnUBC promoters suggested further investigation of their potential functions in plant growth and development. Furthermore, three BnUBCs were predicted as candidate genes for regulating agronomic traits related to oil content and yield through association mapping. In conclusion, this study provided a wealth of information on the UBC family in B. napus and revealed their effects on oil content and yield, which will aid future functional research and genetic breeding of B. napus.

Expression pattern of resynthesized allotetraploid Capsella is determined by hybridization, not whole-genome duplication

Polyploidization, the process leading to the increase in chromosome sets, is a major evolutionary transition in plants. Whole-genome duplication (WGD) within the same species gives rise to autopolyploids, whereas allopolyploids result from a compound process with two distinct components: WGD and interspecific hybridization. To dissect the instant effects of WGD and hybridization on gene expression and phenotype, we created a series of synthetic hybrid and polyploid Capsella plants, including diploid hybrids, autotetraploids of both parental species, and two kinds of resynthesized allotetraploids with different orders of WGD and hybridization. Hybridization played a major role in shaping the relative expression pattern of the neo-allopolyploids, whereas WGD had almost no immediate effect on relative gene expression pattern but, nonetheless, still affected phenotypes. No transposable element-mediated genomic shock scenario was observed in either neo-hybrids or neo-polyploids. Finally, WGD and hybridization interacted and the distorting effects of WGD were less strong in hybrids. Whole-genome duplication may even improve hybrid fertility. In summary, while the initial relative gene expression pattern in neo-allotetraploids was almost entirely determined by hybridization, WGD only had trivial effects on relative expression patterns, both processes interacted and had a strong impact on physical attributes and meiotic behaviors.

A transcriptomic-guided strategy used in identification of a wheat rust pathogen target and modification of the target enhanced host resistance to rust pathogens

Transcriptional reprogramming is an essential feature of plant immunity and is governed by transcription factors (TFs) and co-regulatory proteins associated with discrete transcriptional complexes. On the other hand, effector proteins from pathogens have been shown to hijack these vast repertoires of plant TFs. Our current knowledge of host genes' role (including TFs) involved in pathogen colonization is based on research employing model plants such as Arabidopsis and rice with minimal efforts in wheat rust interactions. In this study, we begun the research by identifying wheat genes that benefit rust pathogens during infection and editing those genes to provide wheat with passive resistance to rust. We identified the wheat MYC4 transcription factor (TF) located on chromosome 1B (TaMYC4-1B) as a rust pathogen target. The gene was upregulated only in susceptible lines in the presence of the pathogens. Down-regulation of TaMYC4-1B using barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) in the susceptible cultivar Chinese Spring enhanced its resistance to the stem rust pathogen. Knockout of the TaMYC4-1BL in Cadenza rendered new resistance to races of stem, leaf, and stripe rust pathogens. We developed new germplasm in wheat via modifications of the wheat TaMYC4-1BL transcription factor.

Homoeolog gene expression analysis reveals novel expression biases in upland hybrid cotton under intraspecific hybridization

Hybridization is useful to enhance the yield potential of agronomic crops in the world. Cotton has genome doubling due to the allotetraploid process and hybridization in coordination with duplicated genome can produce more yield and adaptability. Therefore, the expression of homoeologous gene pairs between hybrids and inbred parents is vital to characterize the genetic source of heterosis in cotton. Investigation results of homoeolog gene pairs between two contrasting hybrids and their respective inbred parents identified 36853 homoeolog genes in hybrids. It was observed both high and low hybrids had similar trends in homoeolog gene expression patterns in each tissue under study. An average of 96% of homoeolog genes had no biased expression and their expressions were derived from the equal contribution of both parents. Besides, very few homoeolog genes (an average of 1%) showed no biased or novel expression in both hybrids. The functional analysis described secondary metabolic pathways had a majority of novel biased homoeolog genes in hybrids. These results contribute preliminary knowledge about how hybridization affects expression patterns of homoeolog gene pairs in upland cotton hybrids. Our study also highlights the functional genomics of metabolic genes to explore the genetic mechanism of heterosis in cotton.

Genome-Wide Identification and Analysis of Ariadne Gene Family Reveal Its Genetic Effects on Agronomic Traits of Brassica napus

E3 ligases promote protein ubiquitination and degradation, which regulate every aspect of eukaryotic life. The Ariadne (ARI) proteins of RBR (ring between ring fingers) protein subfamily has been discovered as a group of potential E3 ubiquitin ligases. Only a few available research studies show their role in plant adaptations processes against the external environment. Presently, the functions of ARI proteins are largely unknown in plants. Therefore, in this study, we performed genome-wide analysis to identify the ARI gene family and explore their potential importance in B. napus. A total of 39 ARI genes were identified in the B. napus genome and were classified into three subfamilies (A, B and C) based on phylogenetic analysis. The protein–protein interaction networks and enrichment analysis indicated that BnARI genes could be involved in endoreduplication, DNA repair, proteasome assembly, ubiquitination, protein kinase activity and stress adaptation. The transcriptome data analysis in various tissues provided us an indication of some BnARI genes’ functional importance in tissue development. We also identified potential BnARI genes that were significantly responsive towards the abiotic stresses. Furthermore, eight BnARI genes were identified as candidate genes for multiple agronomic traits through association mapping analysis in B. napus; among them, BnaA02g12100D, which is the ortholog of AtARI8, was significantly associated with ten agronomic traits. This study provided useful information on BnARI genes, which could aid targeted functional research and genetic improvement for breeding in B. napus.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Wheat-Thinopyrum Substitution Lines Imprint Compensation Both From Recipients and Donors

Even frequently used in wheat breeding, we still have an insufficient understanding of the biology of the products via distant hybridization. In this study, a transcriptomic analysis was performed for six Triticum aestivum-Thinopyrum elongatum substitution lines in comparison with the host plants. All the six disomic substitution lines showed much stronger "transcriptomic-shock" occurred on alien genomes with 57.43-69.22% genes changed expression level but less on the recipient genome (2.19-8.97%). Genome-wide suppression of alien genes along chromosomes was observed with a high proportion of downregulated genes (39.69-48.21%). Oppositely, the wheat recipient showed genome-wide compensation with more upregulated genes, occurring on all chromosomes but not limited to the homeologous groups. Moreover, strong co-upregulation of the orthologs between wheat and Thinopyrum sub-genomes was enriched in photosynthesis with predicted chloroplastic localization, which indicates that the compensation happened not only on wheat host genomes but also on alien genomes.

Genome-Wide Characterization of Serine/Arginine-Rich Gene Family and Its Genetic Effects on Agronomic Traits of Brassica napus

Serine/arginine-rich (SR) proteins are indispensable factors for RNA splicing, and they play important roles in development and abiotic stress responses. However, little information on SR genes in Brassica napus is available. In this study, 59 SR genes were identified and classified into seven subfamilies: SR, SCL, RS2Z, RSZ, RS, SR45, and SC. In each subfamily, the genes showed relatively conserved structures and motifs, but displayed distinct expression patterns in different tissues and under abiotic stress, which might be caused by the varied cis-acting regulatory elements among them. Transcriptome datasets from Pacbio/Illumina platforms showed that alternative splicing of SR genes was widespread in B. napus and the majority of paralogous gene pairs displayed different splicing patterns. Protein-protein interaction analysis indicated that SR proteins were involved in the regulation of the whole lifecycle of mRNA, from synthesis to decay. Moreover, the association mapping analysis suggested that 12 SR genes were candidate genes for regulating specific agronomic traits, which indicated that SR genes could affect the development and hence influence the important agronomic traits of B. napus. In summary, this study provided elaborate information on SR genes in B. napus, which will aid further functional studies and genetic improvement of agronomic traits in B. napus.

Genetic Conservation of CBS Domain Containing Protein Family in Oryza Species and Their Association with Abiotic Stress Responses

Crop Wild Relatives (CWRs) form a comprehensive gene pool that can answer the queries related to plant domestication, speciation, and ecological adaptation. The genus ‘Oryza’ comprises about 27 species, of which two are cultivated, while the remaining are wild. Here, we have attempted to understand the conservation and diversification of the genes encoding Cystathionine β-synthase (CBS) domain-containing proteins (CDCPs) in domesticated and CWRs of rice. Few members of CDCPs were previously identified to be stress-responsive and associated with multiple stress tolerance in rice. Through genome-wide analysis of eleven rice genomes, we identified a total of 36 genes encoding CDCPs in O. longistaminata, 38 in O. glaberrima, 39 each in O. rufipogon, O. glumaepatula, O. brachyantha, O. punctata, and O. sativa subsp. japonica, 40 each in O. barthii and O. meridionalis, 41 in O. nivara, and 42 in O. sativa subsp. indica. Gene duplication analysis as well as non-synonymous and synonymous substitutions in the duplicated gene pairs indicated that this family is shaped majorly by the negative or purifying selection pressure through the long-term evolution process. We identified the presence of two additional hetero-domains, namely TerCH and CoatomerE (specifically in O. sativa subsp. indica), which were not reported previously in plant CDCPs. The in silico expression analysis revealed some of the members to be responsive to various abiotic stresses. Furthermore, the qRT-PCR based analysis identified some members to be highly inducive specifically in salt-tolerant genotype in response to salinity. The cis-regulatory element analysis predicted the presence of numerous stress as well as a few phytohormone-responsive elements in their promoter region. The data presented in this study would be helpful in the characterization of these CDCPs from rice, particularly in relation to abiotic stress tolerance.

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.

Genetic Divergence of Lineage-Specific Tandemly Duplicated Gene Clusters in Four Diploid Potato Genotypes

Potato (Solanum tuberosum L.) is the most important non-grain food crop. Tandem duplication significantly contributes to genome evolution. The objectives of this study were to (i) identify tandemly duplicated genes and compare their genomic distributions across potato genotypes, (ii) investigate the bias in functional specificities, (iii) explore the relationships among coding sequence, promoter and expression divergences associated with tandemly duplicated genes, (iv) examine the role of tandem duplication in generating and expanding lineage-specific gene families, (v) investigate the evolutionary forces affecting tandemly duplicated genes, and (vi) assess the similarities and differences with respect to above mentioned aspects between cultivated genotypes and their wild-relative. In this study, we used well-annotated and chromosome-scale de novo genome assemblies of multiple potato genotypes. Our results showed that tandemly duplicated genes are abundant and dispersed through the genome. We found that several functional specificities, such as disease resistance, stress-tolerance, and biosynthetic pathways of tandemly duplicated genes were differentially enriched across multiple potato genomes. Our results indicated the existence of a significant correlation among expression, promoter, and protein divergences in tandemly duplicated genes. We found about one fourth of tandemly duplicated gene clusters as lineage-specific among multiple potato genomes, and these tended to localize toward centromeres and revealed distinct selection signatures and expression patterns. Furthermore, our results showed that a majority of duplicated genes were retained through sub-functionalization followed by genetic redundancy, while only a small fraction of duplicated genes was retained though neo-functionalization. The lineage-specific expansion of gene families by tandem duplication coupled with functional bias might have significantly contributed to potato's genotypic diversity, and, thus, to adaption to environmental stimuli.

Transcriptional loss of domestication-driven cytoskeletal GhPRF1 gene causes defective floral and fiber development in cotton (Gossypium)

Constitutive- and fiber-specific RNAi of GhPRF1 gene illustrated strong correlation between domestication-driven profilin genes and floral/fiber architecture in cotton. During morpho-transformation of short-fuzz of wild cotton into the elongating spinnable fibers under the millennia of human selection, actin-polymerizing cytoskeletal profilin genes had undergone significant sequence alterations and spatiotemporal shift in their transcription levels. To comprehend the expression dynamics of profilin genes with their phenotypic implications, transgenic expression modulation of cotton profilin 1 (GhPRF1) gene was performed in the constitutive- and fiber-specific manner in Coker 310FR cotton cultivar. The constitutive GhPRF1-RNAi lines (35S:GhPRF1-RNAi) exhibited distorted 'monadelphous' staminal-tube, reduced pollen-viability and poorly developed fibers, whereas floral and fiber development of fiber-specific GhPRF1-RNAi lines showed no abnormalities. Moreover, the fiber-specific GhPRF1 overexpression lines (FBP7:GhPRF1-Ox) showed increased emergence of fiber-initials on the ovule surface, on the contrary to no fiber-initials in fiber-specific RNAi lines (FBP7:GhPRF1-RNAi). Interestingly, the average seed weight and fiber weight of FBP7:GhPRF1-Ox lines increased > 60% and > 38%, respectively, compared with FBP7:GhPRF1-RNAi lines and untransformed control seeds. On a molecular basis, the aberrant floral and fiber development of 35S:GhPRF1-RNAi lines was largely associated with sugar metabolism and hormone-signaling mechanisms. These observations illustrated the strong correlation between domestication-driven GhPRF genes, and floral/fiber development in cotton. Also, the enhanced agronomic traits in GhPRF1-Ox lines of cotton empowered us to recognize their imperative roles, and their future deployment for the sustainable cotton crop improvement.

Homoploid F1 hybrids and segmental allotetraploids of japonica and indica rice subspecies show similar and enhanced tolerance to nitrogen deficiency than parental lines

It remains unclear whether the merger of two divergent genomes by hybridization at the homoploid level or coupled with whole-genome duplication (WGD; allopolyploidy) can result in plants having better tolerance to stress conditions. In this study, we compared phenotypic performance and gene expression in the two diploid subspecies of rice (Oryza sativa subsp. japonica and indica), their reciprocal F1 hybrids, and in segmental allotetraploids under normal and nitrogen (N)-deficient conditions. We found that F1 hybrids and tetraploids showed higher and similar levels of tolerance to N deficiency than either parent. In parallel, total expression levels of 18 relevant functional genes were less perturbed by N deficiency in the F1 hybrids and tetraploids than in the parents. This was consistent with stable intrinsic partitioning of allelic/homoeologous expression defined by parental legacy in the homoploid F1 hybrids/tetraploids between the two conditions. The results suggest that genetic additivity at both the homoploid and allopolyploidy level might lead to similar beneficial phenotypic responses to nitrogen stress compared with the parents. The lack of synergistic responses to N limitation concomitant with WGD, relative to that exhibited by F1 hybrids, adds new empirical evidence in support of the emerging hypothesis that hybridization by itself can play a significant role in plant adaptive evolution in times of stress.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

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.

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.

Global insights into duplicated gene expression and alternative splicing in polyploid Brassica napus under heat, cold, and drought stress

Polyploidy has been a prevalent process during plant evolution and it has made a major impact on the structure and evolution of plant genomes. Many important crop plants are polyploid. There is considerable interest in expression patterns of duplicated genes in polyploids. Alternative splicing (AS) is a fundamental aspect of gene expression that produces multiple final transcript types from a single type of mRNAs. The effects of abiotic stress conditions on AS in polyploids has received little attention. We conducted a global transcriptome analysis of Brassica napus, an allotetraploid derived from B. rapa (A_T ) and B. oleracea (C_T ), by RNA-Seq of plants subjected to cold, heat, and drought stress treatments. Analyses of 27,360 pairs of duplicated genes revealed overall A_T subgenome biases in gene expression and C_T subgenome biases in the extent of alternative splicing under all three stress treatments. More genes increased in expression than decreased in response to the stresses. Negative correlations were found between expression levels and AS frequency for each type of AS. Cold stress produced the greatest changes in gene expression and AS. Cold-induced AS changes were more likely to be shared with those generated by drought than by heat stress. We used homeologs of FLC and CCA1 as case studies to show the dynamics of how duplicates in a polyploid respond to cold stress. Our results suggest that divergence in gene expression and AS patterns between duplicated genes may increase the flexibility of polyploids when responding to abiotic stressors.

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.

Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton

BACKGROUND

Utilization of heterosis has greatly improved the productivity of many crops worldwide. Understanding the potential molecular mechanism about how hybridization produces superior yield in upland cotton is critical for efficient breeding programs.

RESULTS

In this study, high, medium, and low hybrids varying in the level of yield heterosis were screened based on field experimentation of different years and locations. Phenotypically, high hybrid produced a mean of 14% more seed cotton yield than its better parent. Whole-genome RNA sequencing of these hybrids and their four inbred parents was performed using different tissues of the squaring stage. Comparative transcriptomic differences in each hybrid parent triad revealed a higher percentage of differentially expressed genes (DEGs) in each tissue. Expression level dominance analysis identified majority of hybrids DEGs were biased towards parent like expressions. An array of DEGs involved in ATP and protein binding, membrane, cell wall, mitochondrion, and protein phosphorylation had more functional annotations in hybrids. Sugar metabolic and plant hormone signal transduction pathways were most enriched in each hybrid. Further, these two pathways had most mapped DEGs on known seed cotton yield QTLs. Integration of transcriptome, QTLs, and gene co-expression network analysis discovered genes Gh_A03G1024, Gh_D08G1440, Gh_A08G2210, Gh_A12G2183, Gh_D07G1312, Gh_D08G1467, Gh_A03G0889, Gh_A08G2199, and Gh_D05G0202 displayed a complex regulatory network of many interconnected genes. qRT-PCR of these DEGs was performed to ensure the accuracy of RNA-Seq data.

CONCLUSIONS

Through genome-wide comparative transcriptome analysis, the current study identified nine key genes and pathways associated with biological process of yield heterosis in upland cotton. Our results and data resources provide novel insights and will be useful for dissecting the molecular mechanism of yield heterosis in cotton.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Expression patterns and functional divergence of homologous genes accompanied by polyploidization in cotton (Gossypium hirsutum L.)

Naturally allotetraploid cotton has been widely used as an ideal model to investigate gene expression remodeling as a consequence of polyploidization. However, the global gene pattern variation during early fiber development was unknown. In this study, through RNA-seq technology, we comprehensively investigated the expression patterns of homologous genes between allotetraploid cotton (G. hirsutum) and its diploid progenitors (G. arboreum and G. raimondii) at the fiber early development stage. In tetraploid cotton, genes showed expression level dominance (ELD) bias toward the A genome. This phenomenon was explained by the up-/downregulation of the homologs from the nondominant progenitor (D genome). Gene ontology (GO) enrichment results indicated that the ELD-A genes might be a prominent cause responsible for fiber property change through regulating the fatty acid biosynthesis/metabolism and microtubule procession, and the ELD-D genes might be involved in transcription regulation and stress inducement. In addition, the number and proportion of completely A- and D-subfunctionalized gene were similar at different fiber development stages. However, for neofunctionalization, the number and proportion of reactivated D-derived genes were greater than those of A at 3 and 5 DPA. Eventually, we found that some homologous genes belonging to several specific pathways might create novel asymmetric transcripts between two subgenomes during polyploidization and domestication process, further making the fiber property meet the human demands. Our study identified determinate pathways and their involved genes between allotetraploid cotton and their progenitors at early fiber development stages, providing new insights into the mechanism of cotton fiber evolution.

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

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

Early consequences of allopolyploidy alter floral evolution in Nicotiana (Solanaceae)

BACKGROUND

Polyploidy has played a major role in angiosperm evolution. Previous studies have examined polyploid phenotypes in comparison to their extant progenitors, but not in context of predicted progenitor phenotypes at allopolyploid origin. In addition, differences in the trends of polyploid versus diploid evolution have not been investigated. We use ancestral character-state reconstructions to estimate progenitor phenotype at allopolyploid origin to determine patterns of polyploid evolution leading to morphology of the extant species. We also compare trends in diploid versus allopolyploid evolution to determine if polyploidy modifies floral evolutionary patterns.

RESULTS

Predicting the ancestral phenotype of a nascent allopolyploid from reconstructions of diploid phenotypes at the time of polyploid formation generates different phenotype predictions than when extant diploid phenotypes are used, the outcome of which can alter conclusions about polyploid evolution; however, most analyses yield the same results. Using ancestral reconstructions of diploid floral phenotypes indicate that young polyploids evolve shorter, wider corolla tubes, but older polyploids and diploids do not show any detectable evolutionary trends. Lability of the traits examined (floral shape, corolla tube length, and corolla tube width) differs across young and older polyploids and diploids. Corolla length is more evolutionarily labile in older polyploids and diploids. Polyploids do not display unique suites of floral characters based on both morphological and color traits, but some suites of characters may be evolving together and seem to have arisen multiple times within Nicotiana, perhaps due to the influence of pollinators.

CONCLUSIONS

Young polyploids display different trends in floral evolution (shorter, wider corolla tubes, which may result in more generalist pollination) than older polyploids and diploids, suggesting that patterns of divergence are impacted by the early consequences of allopolyploidy, perhaps arising from genomic shock and/or subsequent genome stabilization associated with diploidization. Convergent evolution in floral morphology and color in Nicotiana can be consistent with pollinator preferences, suggesting that pollinators may have shaped floral evolution in Nicotiana.

Genome-wide analysis of the HAK potassium transporter gene family reveals asymmetrical evolution in tobacco (Nicotiana tabacum)

Being an essential mineral nutrient, potassium (K+) plays numerous important roles in plant growth and development and determines the yield and quality of crop products. The cellular level of K+ is controlled to a large extent by the K+ transporter, which belongs to the KT/HAK/KUP (HAK) family. However, little is known about these genes in tobacco. In this study, we surveyed the tobacco genome and identified 41 putative NtHAK genes (NtHAKS1-NtHAKS21 and NtHAKT1-NtHAKT20). Investigation of the cis-elements in upstream regions of these NtHAK genes suggests that members of this family respond to environmental cues and phytohormones. Expression data mining reveals that NtHAK genes showed clear sub-genome dominance. In all, these results will provide molecular insights into K+ transporter research in tobacco.

Allele-specific expression variation at different ploidy levels in Squalius alburnoides

Allopolyploid plants are long known to be subject to a homoeolog expression bias of varying degree. The same phenomenon was only much later suspected to occur also in animals based on studies of single selected genes in an allopolyploid vertebrate, the Iberian fish Squalius alburnoides. Consequently, this species became a good model for understanding the evolution of gene expression regulation in polyploid vertebrates. Here, we analyzed for the first time genome-wide allele-specific expression data from diploid and triploid hybrids of S. alburnoides and compared homoeolog expression profiles of adult livers and of juveniles. Co-expression of alleles from both parental genomic types was observed for the majority of genes, but with marked homoeolog expression bias, suggesting homoeolog specific reshaping of expression level patterns in hybrids. Complete silencing of one allele was also observed irrespective of ploidy level, but not transcriptome wide as previously speculated. Instead, it was found only in a restricted number of genes, particularly ones with functions related to mitochondria and ribosomes. This leads us to hypothesize that allelic silencing may be a way to overcome intergenomic gene expression interaction conflicts, and that homoeolog expression bias may be an important mechanism in the achievement of sustainable genomic interactions, mandatory to the success of allopolyploid systems, as in S. alburnoides.

Cis-trans controls and regulatory novelty accompanying allopolyploidization

Allopolyploidy is a prevalent process in plants, having important physiological, ecological and evolutionary consequences. Transcriptomic responses to genomic merger and doubling have been demonstrated in many allopolyploid systems, encompassing a diversity of phenomena including homoeolog expression bias, genome dominance, expression-level dominance and revamping of co-expression networks. Notwithstanding the foregoing, there remains a need to develop a conceptual framework that will stimulate a deeper understanding of these diverse phenomena and their mechanistic interrelationships. Here we introduce considerations relevant to this framework with a focus on cis-trans interactions among duplicated genes and alleles in hybrids and allopolyploids. By extending classic allele-specific expression analysis to the allopolyploid level, we distinguish the distinct effects of progenitor regulatory interactions from the novel intergenomic interactions that arise from genome merger and allopolyploidization. This perspective informs experiments designed to reveal the molecular genetic basis of gene regulatory control, and will facilitate the disentangling of genetic from epigenetic and higher-order effects that impact gene expression. Finally, we suggest that the extended cis-trans model may help conceptually unify several presently disparate hallmarks of allopolyploid evolution, including genome-wide expression dominance and biased fractionation, and lead to a new level of understanding of phenotypic novelty accompanying polyploidy.

Homeologous Epistasis in Wheat: The Search for an Immortal Hybrid

Hybridization between related species results in the formation of an allopolyploid with multiple subgenomes. These subgenomes will each contain complete, yet evolutionarily divergent, sets of genes. Like a diploid hybrid, allopolyploids will have two versions, or homeoalleles, for every gene. Partial functional redundancy between homeologous genes should result in a deviation from additivity. These epistatic interactions between homeoalleles are analogous to dominance effects, but are fixed across subgenomes through self pollination. An allopolyploid can be viewed as an immortalized hybrid, with the opportunity to select and fix favorable homeoallelic interactions within inbred varieties. We present a subfunctionalization epistasis model to estimate the degree of functional redundancy between homeoallelic loci and a statistical framework to determine their importance within a population. We provide an example using the homeologous dwarfing genes of allohexaploid wheat, Rht-1, and search for genome-wide patterns indicative of homeoallelic subfunctionalization in a breeding population. Using the IWGSC RefSeq v1.0 sequence, 23,796 homeoallelic gene sets were identified and anchored to the nearest DNA marker to form 10,172 homeologous marker sets. Interaction predictors constructed from products of marker scores were used to fit the homeologous main and interaction effects, as well as estimate whole genome genetic values. Some traits displayed a pattern indicative of homeoallelic subfunctionalization, while other traits showed a less clear pattern or were not affected. Using genomic prediction accuracy to evaluate importance of marker interactions, we show that homeologous interactions explain a portion of the nonadditive genetic signal, but are less important than other epistatic interactions.

Functional diversifications of GhERF1 duplicate genes after the formation of allotetraploid cotton

Whole genome duplication, a prevalent force of evolution in plants, results in massive genome restructuring in different organisms. Roles of the resultant duplicated genes are poorly understood, both functionally and evolutionarily. In the present study, differentially expressed ethylene responsive factors (GhERF1s), anchored on Chr-A07 and Chr-D07, were isolated from a high-yielding cotton hybrid (XZM2) and its parents. The GhERF1 was located in the B3 subgroup of the ethylene responsive factors subfamily involved in conferring tolerance to abiotic stress. Nucleotide sequence analysis of 524 diverse accessions, together with quantitative real-time polymerase chain reaction analysis, elucidated that de-functionalization of GhERF1-7A occurred due to one base insertion following formation of the allotetraploid cotton. Our quantitative trait loci and association mapping analyses highlighted a role for GhERF1-7A in conferring high boll number per plant in modern cotton cultivars. Overexpression of GhERF1-7A in transgenic Arabidopsis resulted in a substantial increase in the number of siliques and total seed yield. Neo-functionalization of GhERF1-7A was also observed in modern cultivars rather than in races and/or landraces, further supporting its role in the development of high-yielding cotton cultivars. Both de- and neo-functionalization occurred in one of the duplicate genes, thus providing new genomic insight into the evolution of allotetraploid cotton species.

Evolution and Expression Divergence of the CYP78A Subfamily Genes in Soybean

Gene expression divergence is an important evolutionary driving force for the retention of duplicate genes. In this study, we identified three CYP78A subfamily genes in soybean, GmCYP78A70, GmCYP78A57 and GmCYP78A72, which experienced different duplication events. GmCYP78A70 was mainly expressed in leaf tissue and the vegetative phase, whereas GmCYP78A57 was mainly expressed in floral tissue and seed, i.e., the reproductive phase. Expression of GmCYP78A72 could be detected in all the tissues and phases mentioned above. The expression levels of GmCYP78A70 and GmCYP78A57 in different soybean cultivars showed positive correlations with leaf size and 100-seed weight, respectively. The population genetics analysis indicated that the three genes had experienced different selective pressures during domestication and improved breeding of soybean. Deciphering the function of this subfamily of genes may well prove useful to breeders for improving soybean’s agronomic traits.

Transcriptional regulation of anthocyanin biosynthesis in a high-anthocyanin resynthesized Brassica napus cultivar

BACKGROUND

Anthocyanins are plant secondary metabolites with key roles in attracting insect pollinators and protecting against biotic and abiotic stresses. They have potential health-promoting effects as part of the human diet. Anthocyanin biosynthesis has been elucidated in many species, enabling the development of anthocyanin-enriched fruits, vegetables, and grains; however, few studies have investigated Brassica napus anthocyanin biosynthesis.

RESULTS

We developed a high-anthocyanin resynthesized B. napus line, Rs035, by crossing anthocyanin-rich B. rapa (A genome) and B. oleracea (C genome) lines, followed by chromosome doubling. We identified and characterized 73 and 58 anthocyanin biosynthesis genes in silico in the A and C genomes, respectively; these genes showed syntenic relationships with 41 genes in Arabidopsis thaliana and B. napus. Among the syntenic genes, twelve biosynthetic and six regulatory genes showed transgressively higher expression in Rs035, and eight structural genes and one regulatory gene showed additive expression. We identified three early-, four late-biosynthesis pathways, three transcriptional regulator genes, and one transporter as putative candidates enhancing anthocyanin accumulation in Rs035. Principal component analysis and Pearson's correlation coefficients corroborated the contribution of these genes to anthocyanin accumulation.

CONCLUSIONS

Our study lays the foundation for producing high-anthocyanin B. napus cultivars. The resynthesized lines and the differentially expressed genes we have identified could be used to transfer the anthocyanin traits to other commercial rapeseed lines using molecular and conventional breeding.

Genomic Analysis and Prediction within a US Public Collaborative Winter Wheat Regional Testing Nursery

The development of inexpensive, whole-genome profiling enables a transition to allele-based breeding using genomic prediction models. These models consider alleles shared between lines to predict phenotypes and select new lines based on estimated breeding values. This approach can leverage highly unbalanced datasets that are common to breeding programs. The Southern Regional Performance Nursery (SRPN) is a public nursery established by the USDA-ARS in 1931 to characterize performance and quality of near-release wheat ( L.) varieties from breeding programs in the US Central Plains. New entries are submitted annually and can be re-entered only once. The trial is grown at >30 locations each year and lines are evaluated for grain yield, disease resistance, and agronomic traits. Overall genetic gain is measured across years by including common check cultivars for comparison. We have generated whole-genome profiles via genotyping-by-sequencing (GBS) for 939 SPRN entries dating back to 1992 to explore the potential use of the nursery as a genomic selection (GS) training population (TP). The GS prediction models across years (average = 0.33) outperformed year-to-year phenotypic correlation for yield ( = 0.27) for a majority of the years evaluated, suggesting that genomic selection has the potential to outperform low heritability selection on yield in these highly variable environments. We also examined the predictability of programs using both program-specific and whole-set TPs. Generally, the predictability of a program was similar with both approaches. These results suggest that wheat breeding programs can collaboratively leverage the immense datasets that are generated from regional testing networks.

Rapid and pervasive development- and tissue-specific homeolog expression partitioning in newly formed inter-subspecific rice segmental allotetraploids

Background

In diverse plant taxa, whole-genome duplication (WGD) events are major sources of phenotypic novelty. Studies of gene expression in synthetic polyploids have shown immediate expression and functional partitioning of duplicated genes among different tissues. Many studies of the tissue-specific homeolog expression partitioning have focused on allopolyploids that have very different parental genomes, while few studies have focused on autopolyploids or allopolyploids that have similar parental genomes.

Results

In this study, we used a set of reciprocal F1 hybrids and synthetic tetraploids constructed from subspecies (japonica and indica) of Asian rice (Oryza sativa L.) as a model to gain insights into the expression partitioning of homeologs among tissues in a developmental context. We assayed the tissue-specific silencing (TSS) of the parental homeologs of 30 key genes in the hybrids and tetraploids relative to the in vitro “hybrids” (parental mixes) using Sequenom MassARRAY. We found that the parental mix and synthetic tetraploids had higher frequencies of homeolog TSS than the F1, revealing an instantaneous role of WGD on homeolog expression partitioning.

Conclusions

Our observations contradicted those of previous studies in which newly formed allopolyploids had a low TSS frequency, similar to that of F1 hybrids, suggesting that the impact of WGD on homeolog expression requires a longer time to manifest. In addition, we found that the TSS frequency in the tetraploids varied at different growth stages and that roots had a much higher frequency of TSS than leaves, which indicated that developmental and metabolic traits may influence the expression states of duplicated genes in newly formed plant polyploids.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5150-7) contains supplementary material, which is available to authorized users.

Gene expression in response to glyphosate treatment in fleabane (Conyza bonariensis) - glyphosate death response and candidate resistance genes

BACKGROUND

This study takes a whole-transcriptome approach to assess gene expression changes in response to glyphosate treatment in glyphosate-resistant fleabane. We assessed gene expression changes in both susceptible and resistant lines so that the glyphosate death response could be quantified, and constitutively expressed candidate resistance genes identified. There are three copies of the glyphosate target site (5-enolpyruvylshikimate-3-phosphate; EPSPS) gene in Conyza and because Conyza bonariensis is allohexaploid, there is a baseline nine copies of the gene in any individual.

RESULTS

Many genes were differentially expressed in response to glyphosate treatment. Known resistance mutations are present in EPSPS2 but they are present in a glyphosate-susceptible line as well as resistant lines and therefore not sufficient to confer resistance. EPSPS1 is expressed four times more than EPSPS2, further reducing the overall contribution of these mutations.

CONCLUSION

We demonstrate that glyphosate resistance in C. bonariensis is not the result of EPSPS mutations or overexpression, but due to a non-target-site mechanism. A large number of genes are affected by glyphosate treatment. We present a list of candidate non-target-site-resistance (NTSR) genes in fleabane for future studies into these mechanisms. © 2017 Society of Chemical Industry.

Genome-wide analysis of the MADS-box gene family in polyploid cotton (Gossypium hirsutum) and in its diploid parental species (Gossypium arboreum and Gossypium raimondii)

The MADS-box gene family encodes transcription factors that share a highly conserved domain known to bind to DNA. Members of this family control various processes of development in plants, from root formation to fruit ripening. In this work, a survey of diploid (Gossypium raimondii and Gossypium arboreum) and tetraploid (Gossypium hirsutum) cotton genomes found a total of 147, 133 and 207 MADS-box genes, respectively, distributed in the MIKC, Mα, Mβ, Mγ, and Mδ subclades. A comparative phylogenetic analysis among cotton species, Arabidopsis, poplar and grapevine MADS-box homologous genes allowed us to evaluate the evolution of each MADS-box lineage in cotton plants and identify sequences within well-established subfamilies. Chromosomal localization and phylogenetic analysis revealed that G. raimondii and G. arboreum showed a conserved evolution of the MIKC subclade and a distinct pattern of duplication events in the Mα, Mγ and Mδ subclades. Additionally, G. hirsutum showed a combination of its parental subgenomes followed by a distinct evolutionary history including gene gain and loss in each subclade. qPCR analysis revealed the expression patterns of putative homologs in the AP1, AP3, AGL6, SEP4, AGL15, AG, AGL17, TM8, SVP, SOC and TT16 subfamilies of G. hirsutum. The identification of putative cotton orthologs is discussed in the light of evolution and gene expression data from other plants. This analysis of the MADS-box genes in Gossypium species opens an avenue to understanding the origin and evolution of each gene subfamily within diploid and polyploid species and paves the way for functional studies in cotton species.

Allopolyploidization in Cucumis contributes to delayed leaf maturation with repression of redundant homoeologous genes

The important role of polyploidy in plant evolution is widely recognized. However, many questions remain to be explored to address how polyploidy affects the phenotype of the plant. To shed light on the phenotypic and molecular impacts of allopolyploidy, we investigated the leaf development of a synthesized allotetraploid (Cucumis × hytivus), with an emphasis on chlorophyll development. Delayed leaf maturation was identified in C. × hytivus, based on delayed leaf expansion, initial chlorophyll deficiency in the leaves and disordered sink-source transition. Anatomical observations also revealed disturbed chloroplast development in C. ×hytivus. The determination of chlorophyll biosynthesis intermediates suggested that the chlorophyll biosynthesis pathway of C. × hytivus is blocked at the site at which uroporphyrinogen III is catalysed to coproporphyrinogen III. Three chlorophyll biosynthesis-related genes, HEMA1, HEME2 and POR, were significantly repressed in C. × hytivus. Sequence alignment showed both synonymous and non-synonymous substitutions in the HEMA1, HEME2 and POR genes of the parents. Cloning of the chlorophyll biosynthetic genes suggested the retention of homoeologs. In addition, a chimeric clone of the HEMA1 gene that consisted of homologous genes from the parents was identified in C. × hytivus. Overall, our results showed that allopolyploidization in Cucumis has resulted in disturbed chloroplast development and reduced chlorophyll biosynthesis caused by the repressed expression of duplicated homologous genes, which further led to delayed leaf maturation in the allotetraploid, C. × hytivus. The preferential retention/loss of certain types of genes and non-reciprocal homoeologous recombination were also supported in the present study, which provides new insights into the impact of allopolyploidy.

Gene Co-expression Network Analysis Suggests the Existence of Transcriptional Modules Containing a High Proportion of Transcriptionally Differentiated Homoeologs in Hexaploid Wheat

Genome duplications aid in the formation of novel molecular networks through regulatory differentiation of the duplicated genes and facilitate adaptation to environmental change. Hexaploid wheat, Triticum aestivum, contains three homoeologous chromosome sets, the A-, B-, and D-subgenomes, which evolved through interspecific hybridization and subsequent whole-genome duplication. The divergent expression patterns of the homoeologs in hexaploid wheat suggest that they have undergone transcriptional and/or functional differentiation during wheat evolution. However, the distribution of transcriptionally differentiated homoeologs in gene regulatory networks and their related biological functions in hexaploid wheat are still largely unexplored. Therefore, we retrieved 727 publicly available wheat RNA-sequencing (RNA-seq) datasets from various tissues, developmental stages, and conditions, and identified 10,415 expressed homoeologous triplets. Examining the co-expression modules in the wheat transcriptome, we found that 66% of the expressed homoeologous triplets possess all three homoeologs grouped in the same co-expression modules. Among these, 15 triplets contain co-expressed homoeologs with differential expression levels between homoeoalleles across ≥ 95% of the 727 RNA-seq datasets, suggesting a consistent trend of homoeolog expression bias. In addition, we identified 2,831 differentiated homoeologs that showed gene expression patterns that deviated from those of the other two homoeologs. We found that seven co-expression modules contained a high proportion of such differentiated homoeologs, which accounted for ≥ 20% of the genes in each module. We also found that five of the co-expression modules are abundantly composed of genes involved in biological processes such as chloroplast biogenesis, RNA metabolism, putative defense response, putative posttranscriptional modification, and lipid metabolism, thereby suggesting that, the differentiated homoeologs might highly contribute to these biological functions in the gene network of hexaploid wheat.

BnDGAT1s Function Similarly in Oil Deposition and Are Expressed with Uniform Patterns in Tissues of Brassica napus

As an allotetraploid oilcrop, Brassica napus contains four duplicated Acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) genes, which catalyze one of the rate-limiting steps in triacylglycerol (TAG) biosynthesis in plants. While all four BnDGAT1s have been expressed functionally in yeast, their expression patterns in different germplasms and tissues and also consequent contribution to seed oil accumulation in planta remain to be elucidated. In this study, the coding regions of the four BnDGAT1s were expressed in an Arabidopsis dgat1 mutant. All four BnDGAT1s showed similar effects on oil content and fatty acid composition, a result which is different from that observed in previous studies of their expression in yeast. Expression patterns of BnDGAT1s were analyzed in developing seeds of 34 B. napus inbred lines and in different tissues of 14 lines. Different expression patterns were observed for the four BnDGAT1s, which suggests that they express independently or randomly in different germplasm sources. Higher expression of BnDGAT1s was correlated with higher seed oil content lines. Tissue-specific analyses showed that the BnDGAT1s were expressed in a uniform pattern in different tissues. Our results suggest that it is important to maintain expression of the four BnDGAT1s for maximum return on oil content.

Analysis of Three Sugarcane Homo/Homeologous Regions Suggests Independent Polyploidization Events of Saccharum officinarum and Saccharum spontaneum

Whole genome duplication has played an important role in plant evolution and diversification. Sugarcane is an important crop with a complex hybrid polyploid genome, for which the process of adaptation to polyploidy is still poorly understood. In order to improve our knowledge about sugarcane genome evolution and the homo/homeologous gene expression balance, we sequenced and analyzed 27 BACs (Bacterial Artificial Chromosome) of sugarcane R570 cultivar, containing the putative single-copy genes LFY (seven haplotypes), PHYC (four haplotypes), and TOR (seven haplotypes). Comparative genomic approaches showed that these sugarcane loci presented a high degree of conservation of gene content and collinearity (synteny) with sorghum and rice orthologous regions, but were invaded by transposable elements (TE). All the homo/homeologous haplotypes of LFY, PHYC, and TOR are likely to be functional, because they are all under purifying selection (dN/dS ≪ 1). However, they were found to participate in a nonequivalently manner to the overall expression of the corresponding gene. SNPs, indels, and amino acid substitutions allowed inferring the S. officinarum or S. spontaneum origin of the TOR haplotypes, which further led to the estimation that these two sugarcane ancestral species diverged between 2.5 and 3.5 Ma. In addition, analysis of shared TE insertions in TOR haplotypes suggested that two autopolyploidization may have occurred in the lineage that gave rise to S. officinarum, after its divergence from S. spontaneum.

Segmental allotetraploidy generates extensive homoeologous expression rewiring and phenotypic diversity at the population level in rice

Allopolyploidization, that is, concomitant merging and doubling of two or more divergent genomes in a common nucleus/cytoplasm, is known to instantly alter genomewide transcriptome dynamics, a phenomenon referred to as "transcriptomic shock." However, the immediate effects of transcriptomic alteration in generating phenotypic diversity at the population level remain underinvestigated. Here, we employed the MassARRAY-based Sequenom platform to assess and compare orthologous, allelic and homoeologous gene expression status in two tissues (leaf and root) of a set of randomly chosen individuals from populations of parental rice subspecies (indica and japonica), in vitro "hybrids" (parental mixes), reciprocal F1 hybrids and reciprocal tetraploids at the 5th-selfed generation (S5). We show that hybridization and whole genome duplication (WGD) have opposing effects on allelic and homoeologous expression in the F1 hybrids and tetraploids, respectively. Whereas hybridization exerts strong attenuating effects on allelic expression differences in diploid hybrids, WGD augments the intrinsic parental differences and generates extensive and variable homoeolog content which triggers diversification in expression patterning among the tetraploid plants. Coupled with the vast phenotypic diversity observed among the tetraploid individuals, our results provide experimental evidence in support of the notion that allopolyploidy catalyses rapid phenotypic diversification in higher plants. Our data further suggest that largely stochastic homoeolog content reshuffling rather than alteration in total expression level may be an important feature of evolution in young segmental allopolyploids, which underlies rapid expression diversity at the population level.

A comparative genomic and transcriptomic analysis at the level of isolated root hair cells reveals new conserved root hair regulatory elements

KEY MESSAGE

A comparative transcriptomic and genomic analysis between Arabidopsis thaliana and Glycine max root hair genes reveals the evolution of the expression of plant genes after speciation and whole genome duplication. Our understanding of the conservation and divergence of the expression patterns of genes between plant species is limited by the quality of the genomic and transcriptomic resources available. Specifically, the transcriptomes generated from plant organs are the reflection of the contribution of the different cell types composing the samples weighted by their relative abundances in the sample. These contributions can vary between plant species leading to the generation of datasets which are difficult to compare. To gain a deeper understanding of the evolution of gene transcription in and between plant species, we performed a comparative transcriptomic and genomic analysis at the level of one single plant cell type, the root hair cell, and between two model plants: Arabidopsis (Arabidopsis thaliana) and soybean (Glycine max). These two species, which diverged 90 million years ago, were selected as models based on the large amount of genomic and root hair transcriptomic information currently available. Our analysis revealed in detail the transcriptional divergence and conservation between soybean paralogs (i.e., the soybean genome is the product of two successive whole genome duplications) and between Arabidopsis and soybean orthologs in this single plant cell type. Taking advantage of this evolutionary study, we combined bioinformatics, molecular, cellular and microscopic tools to characterize plant promoter sequences and the discovery of two root hair regulatory elements (RHE1 and RHE2) consistently and specifically active in plant root hair cells.