Coordinated and fine-scale control of homoeologous gene expression in allotetraploid cotton

Function deficiency of GhOMT1 causes anthocyanidins over-accumulation and diversifies fibre colours in cotton (Gossypium hirsutum)

Naturally coloured cotton (NCC) fibres need little or no dyeing process in textile industry to low-carbon emission and are environment-friendly. Proanthocyanidins (PAs) and their derivatives were considered as the main components causing fibre coloration and made NCCs very popular and healthy, but the monotonous fibre colours greatly limit the wide application of NCCs. Here a G. hirsutum empurpled mutant (HS2) caused by T-DNA insertion is found to enhance the anthocyanidins biosynthesis and accumulate anthocyanidins in the whole plant. HPLC and LC/MS-ESI analysis confirmed the anthocyanidins methylation and peonidin, petunidin and malvidin formation are blocked. The deficiency of GhOMT1 in HS2 was associated with the activation of the anthocyanidin biosynthesis and the altered components of anthocyanidins. The transcripts of key genes in anthocyanidin biosynthesis pathway are significantly up-regulated in HS2, while transcripts of the genes for transport and decoration were at similar levels as in WT. To investigate the potential mechanism of GhOMT1 deficiency in cotton fibre coloration, HS2 mutant was crossed with NCCs. Surprisingly, offsprings of HS2 and NCCs enhanced PAs biosynthesis and increased PAs levels in their fibres from the accumulated anthocyanidins through up-regulated GhANR and GhLAR. As expected, multiple novel lines with improved fibre colours including orange red and navy blue were produced in their generations. Based on this work, a new strategy for breeding diversified NCCs was brought out by promoting PA biosynthesis. This work will help shed light on mechanisms of PA biosynthesis and bring out potential molecular breeding strategy to increase PA levels in NCCs.

Characterization of floral morphoanatomy and identification of marker genes preferentially expressed during specific stages of cotton flower development

Characterization of anther and ovule developmental programs and expression analyses of stage-specific floral marker genes in Gossypium hirsutum allowed to build a comprehensive portrait of cotton flower development before fiber initiation. Gossypium hirsutum is the most important cotton species that is cultivated worldwide. Although cotton reproductive development is important for fiber production, since fiber is formed on the epidermis of mature ovules, cotton floral development remains poorly understood. Therefore, this work aims to characterize the cotton floral morphoanatomy by performing a detailed description of anther and ovule developmental programs and identifying stage-specific floral marker genes in G. hirsutum. Using light microscopy and scanning electron microscopy, we analyzed anther and ovule development during 11 stages of flower development. To better characterize the ovule development in cotton, we performed histochemical analyses to evaluate the accumulation of phenolic compounds, pectin, and sugar in ovule tissues. After identification of major hallmarks of floral development, three key stages were established in G. hirsutum floral development: in stage 1 (early-EF), sepal, petal, and stamen primordia were observed; in stage 2 (intermediate-IF), primordial ovules and anthers are present, and the differentiating archesporial cells were observed, marking the beginning of microsporogenesis; and in stage 6 (late-LF), flower buds presented initial anther tapetum degeneration and microspore were released from the tetrad, and nucellus and both inner and outer integuments are developing. We used transcriptome data of cotton EF, IF and LF stages to identify floral marker genes and evaluated their expression by real-time quantitative PCR (qPCR). Twelve marker genes were preferentially expressed in a stage-specific manner, including the putative homologs for AtLEAFY, AtAPETALA 3, AtAGAMOUS-LIKE 19 and AtMALE STERILITY 1, which are crucial for several aspects of reproductive development, such as flower organogenesis and anther and petal development. We also evaluated the expression profile of B-class MADS-box genes in G. hirsutum floral transcriptome (EF, IF, and LF). In addition, we performed a comparative analysis of developmental programs between Arabidopsis thaliana and G. hirsutum that considered major morphoanatomical and molecular processes of flower, anther, and ovule development. Our findings provide the first detailed analysis of cotton flower development.

Virus-Induced Gene Silencing Identifies an Important Role of the TaRSR1 Transcription Factor in Starch Synthesis in Bread Wheat

The function of a wheat starch regulator 1 (TaRSR1) in regulating the synthesis of grain storage starch was determined using the barley stripe mosaic virus-virus induced gene-silencing (BSMV-VIGS) method in field experiments. Chlorotic stripes appeared on the wheat spikes infected with barley stripe mosaic virus-virus induced gene-silencing- wheat starch regulator 1 (BSMV-VIGS-TaRSR1) at 15 days after anthesis, at which time the transcription levels of the TaRSR1 gene significantly decreased. Quantitative real-time PCR was also used to measure the transcription levels of 26 starch synthesis-related enzyme genes in the grains of BSMV-VIGS-TaRSR1-silenced wheat plants at 20, 27, and 31 days after anthesis. The results showed that the transcription levels of some starch synthesis-related enzyme genes were markedly induced at different sampling time points: TaSSI, TaSSIV, TaBEIII, TaISA1, TaISA3, TaPHOL, and TaDPE1 genes were induced at each of the three sampling time points and TaAGPS1-b, TaAGPL1, TaAGPL2, TaSSIIb, TaSSIIc, TaSSIIIb, TaBEI, TaBEIIa, TaBEIIb, TaISA2, TaPHOH, and TaDPE2 genes were induced at one sampling time point. Moreover, both the grain starch contents, one thousand kernel weights, grain length and width of BSMV-VIGS-TaRSR1-infected wheat plants significantly increased. These results suggest that TaRSR1 acts as a negative regulator and plays an important role in starch synthesis in wheat grains by temporally regulating the expression of specific starch synthesis-related enzyme genes.

Polyploidy: Pitfalls and paths to a paradigm

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

Transcriptome analysis in Coffea eugenioides, an Arabica coffee ancestor, reveals differentially expressed genes in leaves and fruits

Studies in diploid parental species of polyploid plants are important to understand their contributions to the formation of plant and species evolution. Coffea eugenioides is a diploid species that is considered to be an ancestor of allopolyploid Coffea arabica together with Coffea canephora. Despite its importance in the evolutionary history of the main economic species of coffee, no study has focused on C. eugenioides molecular genetics. RNA-seq creates the possibility to generate reference transcriptomes and identify coding genes and potential candidates related to important agronomic traits. Therefore, the main objectives were to obtain a global overview of transcriptionally active genes in this species using next-generation sequencing and to analyze specific genes that were highly expressed in leaves and fruits with potential exploratory characteristics for breeding and understanding the evolutionary biology of coffee. A de novo assembly generated 36,935 contigs that were annotated using eight databases. We observed a total of ~5000 differentially expressed genes between leaves and fruits. Several genes exclusively expressed in fruits did not exhibit similarities with sequences in any database. We selected ten differentially expressed unigenes in leaves and fruits to evaluate transcriptional profiles using qPCR. Our study provides the first gene catalog for C. eugenioides and enhances the knowledge concerning the mechanisms involved in the C. arabica homeologous. Furthermore, this work will open new avenues for studies into specific genes and pathways in this species, especially related to fruit, and our data have potential value in assisted breeding applications.

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.

A new synthetic allotetraploid (A1A1G2G2) between Gossypium herbaceum and G. australe: bridging for simultaneously transferring favorable genes from these two diploid species into upland cotton

Gossypium herbaceum, a cultivated diploid cotton species (2n = 2x = 26, A₁A₁), has favorable traits such as excellent drought tolerance and resistance to sucking insects and leaf curl virus. G. australe, a wild diploid cotton species (2n = 2x = 26, G₂G₂), possesses numerous economically valuable characteristics such as delayed pigment gland morphogenesis (which is conducive to the production of seeds with very low levels of gossypol as a potential food source for humans and animals) and resistance to insects, wilt diseases and abiotic stress. Creating synthetic allotetraploid cotton from these two species would lay the foundation for simultaneously transferring favorable genes into cultivated tetraploid cotton. Here, we crossed G. herbaceum (as the maternal parent) with G. australe to produce an F₁ interspecific hybrid and doubled its chromosome complement with colchicine, successfully generating a synthetic tetraploid. The obtained tetraploid was confirmed by morphology, cytology and molecular markers and then self-pollinated. The S₁ seedlings derived from this tetraploid gradually became flavescent after emergence of the fifth true leaf, but they were rescued by grafting and produced S₂ seeds. The rescued S₁ plants were partially fertile due to the existence of univalents at Metaphase I of meiosis, leading to the formation of unbalanced, nonviable gametes lacking complete sets of chromosomes. The S₂ plants grew well and no flavescence was observed, implying that interspecific incompatibility, to some extent, had been alleviated in the S₂ generation. The synthetic allotetraploid will be quite useful for polyploidy evolutionary studies and as a bridge for transferring favorable genes from these two diploid species into Upland cotton through hybridization.

Persistence of subgenomes in paleopolyploid cotton after 60 my of evolution

The importance of whole-genome multiplication (WGM) in plant evolution has long been recognized. In flowering plants, WGM is both ubiquitous and in many lineages cyclical, each round followed by substantial gene loss (fractionation). This process may be biased with respect to duplicated chromosomes, often with overexpression of genes in less fractionated relative to more fractionated regions. This bias is hypothesized to arise through downregulation of gene expression through silencing of local transposable elements (TEs). We assess differences in gene expression between duplicated regions of the paleopolyploid cotton genome and demonstrate that the rate of fractionation is negatively correlated with gene expression. We examine recent hypotheses regarding the source of fractionation bias and show that TE-mediated, positional downregulation is absent in the modern cotton genome, seemingly excluding this phenomenon as the primary driver of biased gene loss. Nevertheless, the paleo subgenomes of diploid cotton are still distinguishable with respect to TE content, targeting of 24-nt-small interfering RNAs and GC content, despite approximately 60 My of evolution. We propose that repeat content per se and differential recombination rates may drive biased fractionation following WGM. These data highlight the likely importance of ancient genomic fractionation biases in shaping modern crop genomes.

Polyploidy and novelty: Gottlieb's legacy

Nearly four decades ago, Roose & Gottlieb (Roose & Gottlieb 1976 Evolution 30, 818-830. (doi:10.2307/2407821)) showed that the recently derived allotetraploids Tragopogon mirus and T. miscellus combined the allozyme profiles of their diploid parents (T. dubius and T. porrifolius, and T. dubius and T. pratensis, respectively). This classic paper addressed the link between genotype and biochemical phenotype and documented enzyme additivity in allopolyploids. Perhaps more important than their model of additivity, however, was their demonstration of novelty at the biochemical level. Enzyme multiplicity-the production of novel enzyme forms in the allopolyploids-can provide an extensive array of polymorphism for a polyploid individual and may explain, for example, the expanded ranges of polyploids relative to their diploid progenitors. In this paper, we extend the concept of evolutionary novelty in allopolyploids to a range of genetic and ecological features. We observe that the dynamic nature of polyploid genomes-with alterations in gene content, gene number, gene arrangement, gene expression and transposon activity-may generate sufficient novelty that every individual in a polyploid population or species may be unique. Whereas certain combinations of these features will undoubtedly be maladaptive, some unique combinations of newly generated variation may provide tremendous evolutionary potential and adaptive capabilities.

Polyploidy and small RNA regulation of cotton fiber development

Cotton is not only the most important source of renewal textile fibers, but also an excellent model for studying cell fate determination and polyploidy effects on gene expression and evolution of domestication traits. The combination of A and D-progenitor genomes into allotetraploid cotton induces intergenomic interactions and epigenetic effects, leading to the unequal expression of homoeologous genes. Small RNAs regulate the expression of transcription and signaling factors related to cellular growth, development and adaptation. An example is miRNA-mediated preferential degradation of homoeologous mRNAs encoding MYB-domain transcription factors that are required for the initiation of leaf trichomes in Arabidopsis and of seed fibers in cotton. This example of coevolution between small RNAs and their homoeologous targets could shape morphological traits such as fibers during the selection and domestication of polyploid crops.

Contribution of subgenomes to the transcriptome and their intertwined regulation in the allopolyploid Coffea arabica grown at contrasted temperatures

Polyploidy has occurred throughout the evolutionary history of plants and led to diversification and plant ecological adaptation. Functional plasticity of duplicate genes is believed to play a major role in the environmental adaptation of polyploids. In this context, we characterized genome-wide homoeologous gene expression in Coffea arabica, a recent allopolyploid combining two subgenomes that derive from two closely related diploid species, and investigated its variation in response to changing environment. The transcriptome of leaves of C. arabica cultivated at different growing temperatures suitable for one or the other parental species was examined using RNA-sequencing. The relative contribution of homoeologs to gene expression was estimated for 9959 and 10,628 genes in warm and cold conditions, respectively. Whatever the growing conditions, 65% of the genes showed equivalent levels of homoeologous gene expression. In 92% of the genes, relative homoeologous gene expression varied < 10% between growing temperatures. The subgenome contributions to the transcriptome appeared to be only marginally altered by the different conditions (involving intertwined regulations of homeologs) suggesting that C. arabica's ability to tolerate a broader range of growing temperatures than its diploid parents does not result from differential use of homoeologs.

Molecular characterisation of two homoeologous elicitor-responsive lipin genes in cotton

The identification and molecular characterisation of two lipin-like gene copies (GhLIPN) in cotton, Gossypium hirsutum, an allotetraploid derived from two progenitor diploid Gossypium species, is described. Sequence analyses of the GhLIPN copies, designated GhLIPN-1 and -2, revealed that they contain 11 exons, separated by ten introns. They each have a 2,643 bp open reading frame that encodes 880 aa proteins, and share a 97.7 and 95.5 % sequence similarity at the translated nucleotide and amino acid level, respectively. The GhLIPN genes have a distinct domain architecture consisting of an archetypical N-terminal lipin domain, followed by a haloacid dehalogenase (HAD) domain towards the C-terminus. A Southern blot did not distinguish between the two gene copies, which suggests that they may be homoeologs rather than paralogs. GhLIPN-2 is more similar to a homoeologous sequence from G. raimondii, representing the ancestral D-genome, compared to GhLIPN-1 that matches G. herbaceum and that represents the A-genome. Our data indicate that GhLIPN-1 and GhLIPN-2 are homoeologs that derive from the A- and the D-diploid genomes, respectively. The promoter sequences of GhLIPN-1 and -2 differ by 56 %, as a result of multiple indels. In silico analysis of the promoter regions revealed that both genes contain numerous putative defence-related and elicitor-responsive cis-elements that support a role for GhLIPN in defence responses. Relative quantification real-time PCR confirmed the up-regulation in response to a cell-wall-derived V. dahliae elicitor, which supported the association of GhLIPN with defence signalling. The results add a new dimension to the proposed roles of lipins in plants by suggesting that lipins may have a role in defence signalling.

Diversity of sequences and expression patterns among alleles of a sugarcane loading stem gene

Modern sugarcane cultivars are highly polyploid and aneuploid hybrids, which are propagated as clones. Their complex genome structure comprises 100-130 chromosomes and 10-13 hom(e)ologous copies of most loci. There is preliminary evidence of very high heterozygosity, with implications for genetic improvement approaches ranging from marker-assisted selection to transgenics. Here, we report that sugarcane cultivar Q200 has at least nine alleles at the Loading Stem Gene (ScLSG) locus. Exon-intron structure is identical and the predicted protein products show at least 92 % identity, across sugarcane alleles and the Sorghum homologue Sb07g027880. There is substantial variation in the 5' UTR and promoter regions including numerous allele-specific nucleotide polymorphisms, insertions and deletions. We developed an allele-specific qRT-PCR method to undertake the first compelling test of allele-specific expression in polyploid sugarcane. Seven alleles distinguished by this method all showed peak expression in the sucrose-loading zone of the stem, but there was apparent variability in expression patterns across other tissues. The ScLSG2 and ScLSG5 alleles appear promising for specificity of expression in stems, relative to leaf, meristem, emerging shoot and root tissues. Within the stem, there was activity in parenchyma, vascular and rind tissues. This expression pattern is of interest in basic research and biotechnology aimed at enhanced sucrose content, engineering value-added products, and manipulation of stem biomass composition.

Co-expression of soybean Dicer-like genes in response to stress and development

Regulation of gene transcription and post-transcriptional processes is critical for proper development, genome integrity, and stress responses in plants. Many genes involved in the key processes of transcriptional and post-transcriptional regulation have been well studied in model diploid organisms. However, gene and genome duplication may alter the function of the genes involved in these processes. To address this question, we assayed the stress-induced transcription patterns of duplicated gene pairs involved in RNAi and DNA methylation processes in the paleopolyploid soybean. Real-time quantitative PCR and Sequenom MassARRAY expression assays were used to profile the relative expression ratios of eight gene pairs across eight different biotic and abiotic stress conditions. The transcriptional responses to stress for genes involved in DNA methylation, RNAi processing, and miRNA processing were compared. The strongest evidence for pairwise co-expression in response to stresses was exhibited by non-paralogous Dicer-like (DCL) genes GmDCL2a-GmDCL3a and GmDCL1b-GmDCL2b, most profoundly in root tissues. Among homoeologous or paralogous DCL genes, the Dicer-like 2 (DCL2) gene pair exhibited the strongest response to stress and most conserved co-expression pattern. This was surprising because the DCL2 duplication event is more ancient than the other DCL duplications. Possible mechanisms that may be driving the DCL2 co-expression are discussed.

Genomic asymmetry in allopolyploid plants: wheat as a model

The evolvement of duplicated gene loci in allopolyploid plants has become the subject of intensive studies. Most duplicated genes remain active in neoallopolyploids contributing either to a favourable effect of an extra gene dosage or to the build-up of positive inter-genomic interactions when genes or regulation factors on homoeologous chromosomes are divergent. However, in a small number of loci (about 10%), genes of only one genome are active, while the homoeoalleles on the other genome(s) are either eliminated or partially or completely suppressed by genetic or epigenetic means. For several traits, the retention of controlling genes is not random, favouring one genome over the other(s). Such genomic asymmetry is manifested in allopolyploid wheat by the control of various morphological and agronomical traits, in the production of rRNA and storage proteins, and in interaction with pathogens. It is suggested that the process of cytological diploidization leading to exclusive intra-genomic meiotic pairing and, consequently, to complete avoidance of inter-genomic recombination, has two contrasting effects. Firstly, it provides a means for the fixation of positive heterotic inter-genomic interactions and also maintains genomic asymmetry resulting from loss or silencing of genes. The possible mechanisms and evolutionary advantages of genomic asymmetry are discussed.

Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton

BACKGROUND

Modern allotetraploid cotton contains an "A" and "D" genome from an ancestral polyploidy event that occurred approximately 1-2 million years ago. Diploid A- and D-genome species can be compared to the A- and D-genomes found within these allotetraploids to make evolutionary inferences about polyploidy. In this paper we present a comprehensive EST assembly derived from diploid and model allotetraploid cottons and demonstrate several evolutionary inferences regarding genic evolution that can be drawn from these data.

RESULTS

We generated a set of cotton expressed sequence tags (ESTs), comprising approximately 4.4 million Sanger and next-generation (454) transcripts supplemented by approximately 152 million Illumina reads from diploid and allotetraploid cottons. From the EST alignments we inferred 259,192 genome-specific single nucleotide polymorphisms (SNPs). Molecular evolutionary analyses of protein-coding regions demonstrate that the rate of nucleotide substitution has increased among both allotetraploid genomes relative to the diploids, and that the ratio of nonsynonymous to synonymous substitutions has increased in one of the two polyploid lineages we sampled. We also use these SNPs to show that a surprisingly high percentage of duplicate genes (~7 %) show a signature of non-independent evolution in the allotetraploid nucleus, having experienced one or more episodes of nonreciprocal homoeologous recombination (NRHR).

CONCLUSIONS

In this study we characterize the functional and mutational properties of the cotton transcriptome, produce a large genome-specific SNP database, and detect illegitimate genetic exchanges between duplicate genomes sharing a common allotetraploid nucleus. Our findings have important implications for our understanding of the consequences of polyploidy and duplicate gene evolution. We demonstrate that cotton genes have experienced an increased rate of molecular evolution following duplication by polyploidy, and that polyploidy has enabled considerable levels of nonreciprocal exchange between homoeologous genes.

Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton

BACKGROUND

Modern allotetraploid cotton contains an "A" and "D" genome from an ancestral polyploidy event that occurred approximately 1-2 million years ago. Diploid A- and D-genome species can be compared to the A- and D-genomes found within these allotetraploids to make evolutionary inferences about polyploidy. In this paper we present a comprehensive EST assembly derived from diploid and model allotetraploid cottons and demonstrate several evolutionary inferences regarding genic evolution that can be drawn from these data.

RESULTS

We generated a set of cotton expressed sequence tags (ESTs), comprising approximately 4.4 million Sanger and next-generation (454) transcripts supplemented by approximately 152 million Illumina reads from diploid and allotetraploid cottons. From the EST alignments we inferred 259,192 genome-specific single nucleotide polymorphisms (SNPs). Molecular evolutionary analyses of protein-coding regions demonstrate that the rate of nucleotide substitution has increased among both allotetraploid genomes relative to the diploids, and that the ratio of nonsynonymous to synonymous substitutions has increased in one of the two polyploid lineages we sampled. We also use these SNPs to show that a surprisingly high percentage of duplicate genes (~7 %) show a signature of non-independent evolution in the allotetraploid nucleus, having experienced one or more episodes of nonreciprocal homoeologous recombination (NRHR).

CONCLUSIONS

In this study we characterize the functional and mutational properties of the cotton transcriptome, produce a large genome-specific SNP database, and detect illegitimate genetic exchanges between duplicate genomes sharing a common allotetraploid nucleus. Our findings have important implications for our understanding of the consequences of polyploidy and duplicate gene evolution. We demonstrate that cotton genes have experienced an increased rate of molecular evolution following duplication by polyploidy, and that polyploidy has enabled considerable levels of nonreciprocal exchange between homoeologous genes.

Global transgenerational gene expression dynamics in two newly synthesized allohexaploid wheat (Triticum aestivum) lines

Background

Alteration in gene expression resulting from allopolyploidization is a prominent feature in plants, but its spectrum and extent are not fully known. Common wheat (Triticum aestivum) was formed via allohexaploidization about 10,000 years ago, and became the most important crop plant. To gain further insights into the genome-wide transcriptional dynamics associated with the onset of common wheat formation, we conducted microarray-based genome-wide gene expression analysis on two newly synthesized allohexaploid wheat lines with chromosomal stability and a genome constitution analogous to that of the present-day common wheat.

Results

Multi-color GISH (genomic in situ hybridization) was used to identify individual plants from two nascent allohexaploid wheat lines between Triticum turgidum (2n = 4x = 28; genome BBAA) and Aegilops tauschii (2n = 2x = 14; genome DD), which had a stable chromosomal constitution analogous to that of common wheat (2n = 6x = 42; genome BBAADD). Genome-wide analysis of gene expression was performed for these allohexaploid lines along with their parental plants from T. turgidum and Ae. tauschii, using the Affymetrix Gene Chip Wheat Genome-Array. Comparison with the parental plants coupled with inclusion of empirical mid-parent values (MPVs) revealed that whereas the great majority of genes showed the expected parental additivity, two major patterns of alteration in gene expression in the allohexaploid lines were identified: parental dominance expression and non-additive expression. Genes involved in each of the two altered expression patterns could be classified into three distinct groups, stochastic, heritable and persistent, based on their transgenerational heritability and inter-line conservation. Strikingly, whereas both altered patterns of gene expression showed a propensity of inheritance, identity of the involved genes was highly stochastic, consistent with the involvement of diverse Gene Ontology (GO) terms. Nonetheless, those genes showing non-additive expression exhibited a significant enrichment for vesicle-function.

Conclusions

Our results show that two patterns of global alteration in gene expression are conditioned by allohexaploidization in wheat, that is, parental dominance expression and non-additive expression. Both altered patterns of gene expression but not the identity of the genes involved are likely to play functional roles in stabilization and establishment of the newly formed allohexaploid plants, and hence, relevant to speciation and evolution of T. aestivum.

Homeologous gene expression in response to growing temperature in a recent Allopolyploid (Coffea arabica L.)

Allopolyploidy is considered as a major factor contributing to speciation, diversification, and plant ecological adaptation. In particular, the expression of duplicate genes (homeologs) can be altered leading to functional plasticity and to phenotypic novelty. This study investigated the influence of growing temperatures on homeologous gene expression in Coffea arabica L., a recent allopolyploid involving 2 closely related diploid parental species. The relative expression of homeologs of 13 genes all located in the same genomic region was analyzed using an SNP ratio quantification method based on dideoxy-terminated sequences of cDNA amplicons. The relative expression of homeologous genes varied depending on the gene, the organ, and the growing condition. Nevertheless, expression of both homeologs was always detected (i.e., no silencing). Although the growing conditions were suitable for one or other of the parental species, neither subgenome appeared preferentially expressed. Furthermore, relative homeologous expression showed moderate variations across organs and conditions and appeared uncorrelated between adjacent genes. These results indicate the absence of signs of subfunctionalization suggesting C. arabica has not undergone noticeable diploidization. Furthermore, these results suggest that the expression of homeologous genes in C. arabica is regulated by a shared trans-regulation mechanism acting similarly on the 2 subgenomes and that the observed biases in the relative homeolog expression may result from cis fine-scale factors.

Genomically biased accumulation of seed storage proteins in allopolyploid cotton

Allopolyploidy is an important process during plant evolution that results in the reunion of two divergent genomes into a common nucleus. Many of the immediate as well as longer-term genomic and epigenetic responses to polyploidy have become appreciated. To investigate the modifications of gene expression at the proteome level caused by allopolyploid formation, we conducted a comparative analysis of cotton seed proteomes from the allopolyploid Gossypium hirsutum (AD genome) and its model A-genome and D-genome diploid progenitors. An unexpectedly high level of divergence among the three proteomes was found, with about one-third of all protein forms being genome specific. Comparative analysis showed that there is a higher degree of proteomic similarity between the allopolyploid and its D-genome donor than its A-genome donor, reflecting a biased accumulation of seed proteins in the allopolyploid. Protein identification and genetic characterization of high-abundance proteins revealed that two classes of seed storage proteins, vicilins and legumins, compose the major component of cotton seed proteomes. Analyses further indicate differential regulation or modification of homoeologous gene products, as well as novel patterns in the polyploid proteome that may result from the interaction between homoeologous gene products. Our findings demonstrate that genomic merger and doubling have consequences that extend beyond the transcriptome into the realm of the proteome and that unequal expression of proteins from diploid parental genomes may occur in allopolyploids.

Large changes in anatomy and physiology between diploid Rangpur lime (Citrus limonia) and its autotetraploid are not associated with large changes in leaf gene expression

Very little is known about the molecular origin of the large phenotypic differentiation between genotypes arising from somatic chromosome set doubling and their diploid parents. In this study, the anatomy and physiology of diploid (2x) and autotetraploid (4x) Rangpur lime (Citrus limonia Osbeck) seedlings has been characterized. Growth of 2x was more vigorous than 4x although leaves, stems, and roots of 4x plants were thicker and contained larger cells than 2x that may have a large impact on cell-to-cell water exchanges. Leaf water content was higher in 4x than in 2x. Leaf transcriptome expression using a citrus microarray containing 21 081 genes revealed that the number of genes differentially expressed in both genotypes was less than 1% and the maximum rate of gene expression change within a 2-fold range. Six up-regulated genes in 4x were targeted to validate microarray results by real-time reverse transcription-PCR. Five of these genes were apparently involved in the response to water deficit, suggesting that, in control conditions, the genome expression of citrus autotetraploids may act in a similar way to diploids under water-deficit stress condition. The sixth up-regulated gene which codes for a histone may also play an important role in regulating the transcription of growth processes. These results show that the large phenotypic differentiation in 4x Rangpur lime compared with 2x is not associated with large changes in genome expression. This suggests that, in 4x Rangpur lime, subtle changes in gene expression may be at the origin of the phenotypic differentiation of 4x citrus when compared with 2x.

Characterization of duplicate gene evolution in the recent natural allopolyploid Tragopogon miscellus by next-generation sequencing and Sequenom iPLEX MassARRAY genotyping

Tragopogon miscellus (Asteraceae) is an evolutionary model for the study of natural allopolyploidy, but until now has been under-resourced as a genetic model. Using 454 and Illumina expressed sequence tag sequencing of the parental diploid species of T. miscellus, we identified 7782 single nucleotide polymorphisms that differ between the two progenitor genomes present in this allotetraploid. Validation of a sample of 98 of these SNPs in genomic DNA using Sequenom MassARRAY iPlex genotyping confirmed 92 SNP markers at the genomic level that were diagnostic for the two parental genomes. In a transcriptome profile of 2989 SNPs in a single T. miscellus leaf, using Illumina sequencing, 69% of SNPs showed approximately equal expression of both homeologs (duplicate homologous genes derived from different parents), 22% showed apparent differential expression and 8.5% showed apparent silencing of one homeolog in T. miscellus. The majority of cases of homeolog silencing involved the T. dubius SNP homeolog (164/254; 65%) rather than the T. pratensis homeolog (90/254). Sequenom analysis of genomic DNA showed that in a sample of 27 of the homeologs showing apparent silencing, 23 (85%) were because of genomic homeolog loss. These methods could be applied to any organism, allowing efficient and cost-effective generation of genetic markers.

Transcriptomic changes following recent natural hybridization and allopolyploidy in the salt marsh species Spartina x townsendii and Spartina anglica (Poaceae)

Allopolyploidy results from two events: the merger of divergent genomes and genome duplication. Both events have important functional consequences for the evolution and adaptation of newly formed allopolyploid species. In spite of the significant progress made in recent years, few studies have decoupled the effects of hybridization from genome duplication in the observed patterns of expression changes accompanying allopolyploidy in natural conditions. We used Agilent rice oligomicroarrays to explore gene expression changes following allopolyploidy in Spartina that includes a classic example of recent allopolyploid speciation: S. anglica formed during the 19th century following genome duplication of the hybrid S. x townsendii. Our data indicate important, but different, effects of hybridization and genome duplication in the expression patterns of the hybrid and allopolyploid. Deviation from parental additivity was most important following hybridization and was accompanied by maternal expression dominance, although transgressively expressed genes were also encountered. Maternal dominance was attenuated following genome duplication in S. anglica, but this species exhibits an increased number of transgressively overexpressed genes. These results reflect the decoupled effects of the 'genomic shock' following hybridization and genome redundancy on the genetic, epigenetic and regulatory mechanisms characterizing transcriptomic evolution in allopolyploids.

Tissue-specific silencing of homoeologs in natural populations of the recent allopolyploid Tragopogon mirus

Recent years have seen rapid advances in our knowledge of the transcriptomic consequences of allopolyploidy, primarily through the study of polyploid crops and model systems. However, few studies have distinguished between homoeologs and between tissues, and still fewer have examined young natural allopolyploid populations of independent origin, whose parental species are still present in the same location. Here, we examined the expression of 13 homoeolog pairs in seven tissues of 10 plants of allotetraploid Tragopogon mirus from two natural populations formed by independent polyploidizations between Tragopogon dubius and Tragopogon porrifolius c. 40 generations ago. We compare these with patterns of expression in the diploid parental species from the same locality. Of the 910 assays in T. mirus, 576 (63%) showed expression of both homoeologs, 63 (7%) showed no expression of either homoeolog, 186 (20%) showed nonexpression of one homoeolog across all tissues of a plant, and 72 (8%) showed non-expression of a homoeolog in a particular tissue within a plant. We found two cases of reciprocal tissue-specific expression between homoeologs, potentially indicative of subfunctionalization. Our study shows that tissue-specific silencing, and even apparent subfunctionalization, can arise rapidly in the early generations of natural allopolyploidy.