Polyploidy and the petal transcriptome of Gossypium

Comparative Analysis of Floral Transcriptomes in Gossypium hirsutum (Malvaceae)

Organ-specific transcriptomes provide valuable insight into the genes involved in organ identity and developmental control. This study investigated transcriptomes of floral organs and subtending bracts in wild and domesticated Gossypium hirsutum, focusing on MADS-box genes critical for floral development. The expression profiles of A, B, C, D, and E class genes were analyzed, confirming their roles in floral organ differentiation. Hierarchical clustering revealed similar expression patterns between bracts and sepals, as well as between petals and stamens, while carpels clustered with developing cotton fibers, reflecting their shared characteristics. Beyond MADS-box genes, other transcription factors were analyzed to explore the genetic basis of floral development. While wild and domesticated cotton showed similar expression patterns for key genes, domesticated cotton exhibited significantly higher expression in carpels compared to wild cotton, which aligns with the increased number of ovules in the carpels of domesticated cotton. Functional enrichment analysis highlighted organ-specific roles: genes upregulated in bracts were enriched for photosynthesis-related GO terms, while diverse functions were enriched in floral organs, supporting their respective functions. Notably, A class genes were not significantly expressed in petals, deviating from the ABCDE model, which warrants further analysis. Lastly, the ABCDE class genes exhibited differential homoeolog expression bias toward each subgenome between two accessions, suggesting that the domestication process has influenced homoeolog utilization despite functional constraints in floral organogenesis.

In vitro Induction and Phenotypic Variations of Autotetraploid Garlic (Allium sativum L.) With Dwarfism

Garlic (Allium sativum L.) is a compelling horticultural crop with high culinary and therapeutic values. Commercial garlic varieties are male-sterile and propagated asexually from individual cloves or bulbils. Consequently, its main breeding strategy has been confined to the time-consuming and inefficient selection approach from the existing germplasm. Polyploidy, meanwhile, plays a prominent role in conferring plants various changes in morphological, physiological, and ecological properties. Artificial polyploidy induction has gained pivotal attention to generate new genotype for further crop improvement as a mutational breeding method. In our study, efficient and reliable in vitro induction protocols of autotetraploid garlic were established by applying different antimitotic agents based on high-frequency direct shoot organogenesis initiated from inflorescence explant. The explants were cultured on solid medium containing various concentrations of colchicine or oryzalin for different duration days. Afterward, the ploidy levels of regenerated plantlets with stable and distinguished characters were confirmed by flow cytometry and chromosome counting. The colchicine concentration at 0.2% (w/v) combined with culture duration for 20 days was most efficient (the autotetraploid induction rate was 21.8%) compared to the induction rate of 4.3% using oryzalin at 60 μmol L-1 for 20 days. No polymorphic bands were detected by simple sequence repeat analysis between tetraploid and diploid plantlets. The tetraploids exhibited a stable and remarkable dwarfness effect rarely reported in artificial polyploidization among wide range of phenotypic variations. There are both morphological and cytological changes including extremely reduced plant height, thickening and broadening of leaves, disappearance of pseudostem, density reduction, and augmented width of stomatal. Furthermore, the level of phytohormones, including, indole propionic acid, gibberellin, brassinolide, zeatin, dihydrozeatin, and methyl jasmonate, was significantly lower in tetraploids than those in diploid controls, except indole acetic acid and abscisic acid, which could partly explain the dwarfness in hormonal regulation aspect. Moreover, as the typical secondary metabolites of garlic, organosulfur compounds including allicin, diallyl disulfide, and diallyl trisulfide accumulated a higher content significantly in tetraploids. The obtained dwarf genotype of autotetraploid garlic could bring new perspectives for the artificial polyploids breeding and be implemented as a new germplasm to facilitate investigation into whole-genome doubling consequences.

In silico modelling and characterization of eight blast resistance proteins in resistant and susceptible rice cultivars

BACKGROUND

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) resistance genes are the largest class of plant resistance genes which play an important role in the plant defense response. These genes are better conserved than others and function as a recognition-based immune system in plants through their encoded proteins.

RESULTS

Here, we report the effect of Magnaporthe oryzae, the rice blast pathogen inoculation in resistant BR2655 and susceptible HR12 rice cultivars. Transcriptomic profiling was carried out to analyze differential gene expression in these two cultivars. A total of eight NBS-LRR uncharacterized resistance proteins (RP1, RP2, RP3, RP4, RP5, RP6, RP7, and RP8) were selected in these two cultivars for in silico modeling. Modeller 9.22 and SWISS-MODEL servers were used for the homology modeling of eight RPs. ProFunc server was utilized for the prediction of secondary structure and function. The CDvist Web server and Interpro scan server detected the motif and domains in eight RPs. Ramachandran plot of eight RPs confirmed that the modeled structures occupied favorable positions.

CONCLUSIONS

From the present study, computational analysis of these eight RPs may afford insights into their role, function, and valuable resource for studying the intricate details of the plant defense mechanism. Furthermore, the identification of resistance proteins is useful for the development of molecular markers linked to resistance genes.

Transcriptome analysis of diploid and triploid Populus tomentosa

Triploid Chinese white poplar (Populus tomentosa Carr., Salicaceae) has stronger advantages in growth and better stress resistance and wood quality than diploid P. tomentosa. Using transcriptome sequencing technology to identify candidate transcriptome-based markers for growth vigor in young tree tissue is of great significance for the breeding of P. tomentosa varieties in the future. In this study, the cuttings of diploid and triploid P. tomentosa were used as plant materials, transcriptome sequencing was carried out, and their tissue culture materials were used for RT-qPCR verification of the expression of genes. The results showed that 12,240 differentially expressed genes in diploid and triploid P. tomentosa transcripts were annotated and enriched into 135 metabolic pathways. The top six pathways that enriched the most significantly different genes were plant-pathogen interaction, phenylpropanoid biosynthesis, MAPK signalling pathway-plant, ascorbate and aldarate metabolism, diterpenoid biosynthesis, and the betalain biosynthesis pathway. Ten growth-related genes were selected from pathways of plant hormone signal transduction and carbon fixation in photosynthetic organisms for RT-qPCR verification. The expression levels of MDH and CYCD3 in tissue-cultured and greenhouse planted triploid P. tomentosa were higher than those in tissue-cultured diploid P. tomentosa, which was consist ent with the TMM values calculated by transcriptome.

Transcriptome Dynamics of the Inflorescence in Reciprocally Formed Allopolyploid Tragopogon miscellus (Asteraceae)

Polyploidy is an important evolutionary mechanism and is prevalent among land plants. Most polyploid species examined have multiple origins, which provide genetic diversity and may enhance the success of polyploids. In some polyploids, recurrent origins can result from reciprocal crosses between the same diploid progenitors. Although great progress has been made in understanding the genetic consequences of polyploidy, the genetic implications of reciprocal polyploidization remain poorly understood, especially in natural polyploids. Tragopogon (Asteraceae) has become an evolutionary model system for studies of recent and recurrent polyploidy. Allotetraploid T. miscellus has formed reciprocally in nature with resultant distinctive floral and inflorescence morphologies (i.e., short- vs. long-liguled forms). In this study, we performed comparative inflorescence transcriptome analyses of reciprocally formed T. miscellus and its diploid parents, T. dubius and T. pratensis. In both forms of T. miscellus, homeolog expression of ∼70% of the loci showed vertical transmission of the parental expression patterns (i.e., parental legacy), and ∼20% of the loci showed biased homeolog expression, which was unbalanced toward T. pratensis. However, 17.9% of orthologous pairs showed different homeolog expression patterns between the two forms of T. miscellus. No clear effect of cytonuclear interaction on biased expression of the maternal homeolog was found. In terms of the total expression level of the homeologs studied, 22.6% and 16.2% of the loci displayed non-additive expression in short- and long-liguled T. miscellus, respectively. Unbalanced expression level dominance toward T. pratensis was observed in both forms of T. miscellus. Significantly, genes annotated as being involved in pectin catabolic processes were highly expressed in long-liguled T. miscellus relative to the short-liguled form, and the majority of these differentially expressed genes were transgressively down-regulated in short-liguled T. miscellus. Given the known role of these genes in cell expansion, they may play a role in the differing floral and inflorescence morphologies of the two forms. In summary, the overall inflorescence transcriptome profiles are highly similar between reciprocal origins of T. miscellus. However, the dynamic homeolog-specific expression and non-additive expression patterns observed in T. miscellus emphasize the importance of reciprocal origins in promoting the genetic diversity of polyploids.

Towards the new normal: Transcriptomic convergence and genomic legacy of the two subgenomes of an allopolyploid weed (Capsella bursa-pastoris)

Allopolyploidy has played a major role in plant evolution but its impact on genome diversity and expression patterns remains to be understood. Some studies found important genomic and transcriptomic changes in allopolyploids, whereas others detected a strong parental legacy and more subtle changes. The allotetraploid C. bursa-pastoris originated around 100,000 years ago and one could expect the genetic polymorphism of the two subgenomes to follow similar trajectories and their transcriptomes to start functioning together. To test this hypothesis, we sequenced the genomes and the transcriptomes (three tissues) of allotetraploid C. bursa-pastoris and its parental species, the outcrossing C. grandiflora and the self-fertilizing C. orientalis. Comparison of the divergence in expression between subgenomes, on the one hand, and divergence in expression between the parental species, on the other hand, indicated a strong parental legacy with a majority of genes exhibiting a conserved pattern and cis-regulation. However, a large proportion of the genes that were differentially expressed between the two subgenomes, were also under trans-regulation reflecting the establishment of a new regulatory pattern. Parental dominance varied among tissues: expression in flowers was closer to that of C. orientalis and expression in root and leaf to that of C. grandiflora. Since deleterious mutations accumulated preferentially on the C. orientalis subgenome, the bias in expression towards C. orientalis observed in flowers indicates that expression changes could be adaptive and related to the selfing syndrome, while biases in the roots and leaves towards the C. grandiflora subgenome may be reflective of the differential genetic load.

Most plant species have a polyploid at some stage of their ancestry. Polyploidy, genome doubling through either multiple copies of a single species or through genomes of different species coming into the same nucleus, is therefore a crucial step in plant evolution. Understanding its impact on basic biological functions is thus a matter of interest. Shepherd’s purse (Capsella bursa-pastoris) is a major weed that appeared about 100,000 years ago through hybridization of two diploid species of the same genus. In the present project, we measured genetic diversity and analyzed gene expression patterns in flowers, roots, and leaves of C. bursa-pastoris individuals as well as in its two parental species, the outcrossing C. grandiflora and the self-fertilizing C. orientalis. Our data shows that, after 100,000 generations of evolution, the origin of the two subgenomes can still be seen: the genome inherited from C. grandiflora still differs from the one inherited from self-fertilizing C. orientalis. However, there are also signs that the two genomes have started to work together and are jointly regulated, and the way expression pattern varied across the three tissues indicates that the evolution of gene expression was adaptive.

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.

Divergent gene expression levels between diploid and autotetraploid Tolmiea relative to the total transcriptome, the cell, and biomass

PREMISE OF THE STUDY

Studies of gene expression and polyploidy are typically restricted to characterizing differences in transcript concentration. Using diploid and autotetraploid Tolmiea, we present an integrated approach for cross-ploidy comparisons that account for differences in transcriptome size and cell density and make multiple comparisons of transcript abundance.

METHODS

We use RNA spike-in standards in concert with cell size and density to identify and correct for differences in transcriptome size and compare levels of gene expression across multiple scales: per transcriptome, per cell, and per biomass.

KEY RESULTS

In total, ~17% of all loci were identified as differentially expressed (DEGs) between the diploid and autopolyploid species. The per-transcriptome normalization, the method researchers typically use, captured the fewest DEGs (58% of total DEGs) and failed to detect any DEGs not found by the alternative normalizations. When transcript abundance was normalized per biomass and per cell, ~66% and ~82% of the total DEGs were recovered, respectively. The discrepancy between per-transcriptome and per-cell recovery of DEGs occurs because per-transcriptome normalizations are concentration-based and therefore blind to differences in transcriptome size.

CONCLUSIONS

While each normalization enables valid comparisons at biologically relevant scales, a holistic comparison of multiple normalizations provides additional explanatory power not available from any single approach. Notably, autotetraploid loci tend to conserve diploid-like transcript abundance per biomass through increased gene expression per cell, and these loci are enriched for photosynthesis-related functions.

A Robust Methodology for Assessing Differential Homeolog Contributions to the Transcriptomes of Allopolyploids

Polyploidy has played a pivotal and recurring role in angiosperm evolution. Allotetraploids arise from hybridization between species and possess duplicated gene copies (homeologs) that serve redundant roles immediately after polyploidization. Although polyploidization is a major contributor to plant evolution, it remains poorly understood. We describe an analytical approach for assessing homeolog-specific expression that begins with de novo assembly of parental transcriptomes and effectively (i) reduces redundancy in de novo assemblies, (ii) identifies putative orthologs, (iii) isolates common regions between orthologs, and (iv) assesses homeolog-specific expression using a robust Bayesian Poisson-Gamma model to account for sequence bias when mapping polyploid reads back to parental references. Using this novel methodology, we examine differential homeolog contributions to the transcriptome in the recently formed allopolyploids Tragopogon mirus and T. miscellus (Compositae). Notably, we assess a larger Tragopogon gene set than previous studies of this system. Using carefully identified orthologous regions and filtering biased orthologs, we find in both allopolyploids largely balanced expression with no strong parental bias. These new methods can be used to examine homeolog expression in any tetrapolyploid system without requiring a reference genome.

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.

Genomic and transcriptomic alterations following intergeneric hybridization and polyploidization in the Chrysanthemum nankingense×Tanacetum vulgare hybrid and allopolyploid (Asteraceae)

Allopolyploid formation involves two major events: interspecific hybridization and polyploidization. A number of species in the Asteraceae family are polyploids because of frequent hybridization. The effects of hybridization on genomics and transcriptomics in Chrysanthemum nankingense×Tanacetum vulgare hybrids have been reported. In this study, we obtained allopolyploids by applying a colchicine treatment to a synthesized C. nankingense×T. vulgare hybrid. Sequence-related amplified polymorphism (SRAP), methylation-sensitive amplification polymorphism (MSAP), and high-throughput RNA sequencing (RNA-Seq) technologies were used to investigate the genomic, epigenetic, and transcriptomic alterations in both the hybrid and allopolyploids. The genomic alterations in the hybrid and allopolyploids mainly involved the loss of parental fragments and the gain of novel fragments. The DNA methylation level of the hybrid was reduced by hybridization but was restored somewhat after polyploidization. There were more significant differences in gene expression between the hybrid/allopolyploid and the paternal parent than between the hybrid/allopolyploid and the maternal parent. Most differentially expressed genes (DEGs) showed down-regulation in the hybrid/allopolyploid relative to the parents. Among the non-additive genes, transgressive patterns appeared to be dominant, especially repression patterns. Maternal expression dominance was observed specifically for down-regulated genes. Many methylase and methyltransferase genes showed differential expression between the hybrid and parents and between the allopolyploid and parents. Our data indicate that hybridization may be a major factor affecting genomic and transcriptomic changes in newly formed allopolyploids. The formation of allopolyploids may not simply be the sum of hybridization and polyploidization changes but also may be influenced by the interaction between these processes.

Differential transcript profiling alters regulatory gene expression during the development of Gossypium arboreum, G.stocksii and somatic hybrids

Polyploidy or genome doubling (i.e., the presence of two or more diploid parental genome sets within an organism) are very important in higher plants. Of particular interest are the mechanisms in the new microenvironment of the common nucleus, where doubled regulatory networks interact to generate a viable genetic system capable of regulating growth, development and responses to the environment. To determine the effects of whole genome merging and doubling on the global gene expression architecture of a new polyploid, derived from protoplast fusion of the A₁A₁ genome of Gossypium arboreum and the E₁E₁ genome of Gossypium stocksii, we monitored gene expression through cDNA-AFLP in the somatic hybrids (G. arboreum + G. stocksii). The genomic expression patterns of the somatic hybrids revealed that changes in expression levels mainly involved regulatory genes (31.8% of the gene expression profiles), and the AA and EE genomes contributed equally to genome-wide expression in the newly formed AAEE genome from additivity and dominance effects. These results provide a novel perspective on polyploid gene regulation and hint at the underlying genetic basis of allopolyploid adaption in the new microenvironmental nucleus.

Molecular characterisation and expression profiling of calcineurin B-like (CBL) genes in Chinese cabbage under abiotic stresses

Calcium signals act as a second messenger in plant responses to various abiotic stresses, which regulate a range of physiological processes. Calcium-binding proteins, like calcineurin B-like (CBL) proteins, belong to a unique group of calcium sensors that play a role in calcium signalling. However, their identities and functions are unknown in Chinese cabbage. In this study, 17 CBL genes were identified from the Brassica rapa L. (Chinese cabbage) database and Br135K microarray datasets. They were used to construct a phylogenetic tree with known CBL proteins of other species. Analysis of genomic distribution and evolution revealed different gene duplication in Chinese cabbage compared to Arabidopsis. The microarray expression analysis showed differential expression of BrCBL genes at various temperatures. Organ-specific expression was observed by RT-PCR, and qRT-PCR analyses revealed responsiveness of BrCBL genes to cold, drought and salt stresses. Our findings confirm that CBL genes are involved in calcium signalling and regulate responses to environmental stimuli, suggesting this family gene have crucial role to play in plant responses to abiotic stresses. The results facilitate selection of candidate genes for further functional characterisation. In addition, abiotic stress-responsive genes reported in this study might be exploited for marker-aided backcrossing of Chinese cabbage.

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 and metabolome of synthetic Solanum autotetraploids reveal key genomic stress events following polyploidization

Polyploids are generally classified as autopolyploids, derived from a single species, and allopolyploids, arising from interspecific hybridization. The former represent ideal materials with which to study the consequences of genome doubling and ascertain whether there are molecular and functional rules operating following polyploidization events. To investigate whether the effects of autopolyploidization are common to different species, or if species-specific or stochastic events are prevalent, we performed a comprehensive transcriptomic and metabolomic characterization of diploids and autotetraploids of Solanum commersonii and Solanum bulbocastanum. Autopolyploidization remodelled the transcriptome and the metabolome of both species. In S. commersonii, differentially expressed genes (DEGs) were highly enriched in pericentromeric regions. Most changes were stochastic, suggesting a strong genotypic response. However, a set of robustly regulated transcripts and metabolites was also detected, including purine bases and nucleosides, which are likely to underlie a common response to polyploidization. We hypothesize that autopolyploidization results in nucleotide pool imbalance, which in turn triggers a genomic shock responsible for the stochastic events observed. The more extensive genomic stress and the higher number of stochastic events observed in S. commersonii with respect to S. bulbocastanum could be the result of the higher nucleoside depletion observed in this species.

Transcriptomic changes following synthesis of a Populus full-sib diploid and allotriploid population with different heterozygosities driven by three types of 2n female gamete

Diploid gametes are usually applied to produce triploids of Populus [originating from first-division restitution (FDR), second-division restitution (SDR), and postmeiotic restitution (PMR) 2n eggs]. Three types of 2n gametes transmitted different parental heterozygosities in Populus. Failed spindle formation and no chromosomal separation to opposite poles during meiosis I mean that FDR 2n gametes carry nonsister chromatids that are potentially heterozygous. By contrast, SDR 2n gametes result from failed sister chromatid separation in meiosis II, and therefore, they carry sister chromatid that are potentially homozygous. Completely homozygous 2n gametes can arise from the PMR mechanism. The alteration of gene expression resulting from allopolyploidization is a prominent feature in plants. We compared gene expression in the full-sib progeny of three allotriploid Populus populations (triploid-F, triploid-S, and triploid-P) with that in its parent species, and their full-sib diploid F1 hybrid. Genome-wide expression level dominance was biased toward the maternal in the diploid F1 hybrid and three allotriploid populations, whereas our data indicated important, but different, effects of the transmission of different heterozygosity by 2n female gametes in the expression patterns of allopolyploids. Because of the higher level of heterozygosity, the triploids had higher rates of non-additive and transgressive expression patterns in the triploid-F than in triploid-S and triploid-P. Compared with diploid F1, about 30-fold more genes (251) were differently expressed in the triploid-F than in the triploid-S (9) and triploid-P (8), respectively. These findings indicate that hybridization and polyploidization have immediate and distinct effects on the large-scale patterns of gene expression, and different effects on the transmission of heterozygosity by three 2n female gametes.

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.

Analysis of [Gossypium capitis-viridis × (G.hirsutum × G.australe)2] Trispecific Hybrid and Selected Characteristics

Speciation is always a contentious and challenging issue following with the presence of gene flow. In Gossypium, there are many valuable resources and wild diploid cotton especially C and B genome species possess some excellent traits which cultivated cotton always lacks. In order to explore character transferring rule from wild cotton to upland tetraploid cotton, the [G. capitis-viridis × (G. hirsutum × G. australe)2] triple hybrid was synthesized by interspecies hybridization and chromosome doubling. Morphology comparisons were measured among this hybrid and its parents. It showed that trispecific hybrid F1 had some intermediate morphological characters like leaf style between its parents and some different characters from its parents, like crawl growth characteristics and two kind flower color. It is highly resistant to insects comparing with other cotton species by four year field investigation. By cytogenetic analysis, triple hybrid was further confirmed by meiosis behavior of pollen mother cells. Comparing with regular meiosis of its three parents, it was distinguished by the occurrence of polyads with various numbers of unbalanced microspores and finally generating various abnormal pollen grains. All this phenomenon results in the sterility of this hybrid. This hybrid was further identified by SSR marker from DNA molecular level. It showed that 98 selected polymorphism primers amplified effective bands in this hybrids and its parents. The genetic proportion of three parents in this hybrid is 47.8% from G. hirsutum, 14.3% from G. australe, 7.0% from G. capitis-viridis, and 30.9% recombination bands respectively. It was testified that wild genetic material has been transferred into cultivated cotton and this new germplasm can be incorporated into cotton breeding program.

Methods for mapping and categorization of DNA sequence reads from allopolyploid organisms

Genome read categorization determines the genome of origin for sequence reads from an allopolyploid organism. Different techniques have been used to perform read categorization, mostly based on homoeo-SNPs identified between extant diploid relatives of allopolyploids. We present a novel technique for read categorization implemented by the software PolyDog. We demonstrate its accuracy and improved categorization relative to other methods. We discuss the situations in which one method or another might be most appropriate.

Genome-wide analysis of the omega-3 fatty acid desaturase gene family in Gossypium

BACKGROUND

The majority of commercial cotton varieties planted worldwide are derived from Gossypium hirsutum, which is a naturally occurring allotetraploid produced by interspecific hybridization of A- and D-genome diploid progenitor species. While most cotton species are adapted to warm, semi-arid tropical and subtropical regions, and thus perform well in these geographical areas, cotton seedlings are sensitive to cold temperature, which can significantly reduce crop yields. One of the common biochemical responses of plants to cold temperatures is an increase in omega-3 fatty acids, which protects cellular function by maintaining membrane integrity. The purpose of our study was to identify and characterize the omega-3 fatty acid desaturase (FAD) gene family in G. hirsutum, with an emphasis on identifying omega-3 FADs involved in cold temperature adaptation.

RESULTS

Eleven omega-3 FAD genes were identified in G. hirsutum, and characterization of the gene family in extant A and D diploid species (G. herbaceum and G. raimondii, respectively) allowed for unambiguous genome assignment of all homoeologs in tetraploid G. hirsutum. The omega-3 FAD family of cotton includes five distinct genes, two of which encode endoplasmic reticulum-type enzymes (FAD3-1 and FAD3-2) and three that encode chloroplast-type enzymes (FAD7/8-1, FAD7/8-2, and FAD7/8-3). The FAD3-2 gene was duplicated in the A genome progenitor species after the evolutionary split from the D progenitor, but before the interspecific hybridization event that gave rise to modern tetraploid cotton. RNA-seq analysis revealed conserved, gene-specific expression patterns in various organs and cell types and semi-quantitative RT-PCR further revealed that FAD7/8-1 was specifically induced during cold temperature treatment of G. hirsutum seedlings.

CONCLUSIONS

The omega-3 FAD gene family in cotton was characterized at the genome-wide level in three species, showing relatively ancient establishment of the gene family prior to the split of A and D diploid progenitor species. The FAD genes are differentially expressed in various organs and cell types, including fiber, and expression of the FAD7/8-1 gene was induced by cold temperature. Collectively, these data define the genetic and functional genomic properties of this important gene family in cotton and provide a foundation for future efforts to improve cotton abiotic stress tolerance through molecular breeding approaches.

Genome-wide analysis of the omega-3 fatty acid desaturase gene family in Gossypium

BACKGROUND

The majority of commercial cotton varieties planted worldwide are derived from Gossypium hirsutum, which is a naturally occurring allotetraploid produced by interspecific hybridization of A- and D-genome diploid progenitor species. While most cotton species are adapted to warm, semi-arid tropical and subtropical regions, and thus perform well in these geographical areas, cotton seedlings are sensitive to cold temperature, which can significantly reduce crop yields. One of the common biochemical responses of plants to cold temperatures is an increase in omega-3 fatty acids, which protects cellular function by maintaining membrane integrity. The purpose of our study was to identify and characterize the omega-3 fatty acid desaturase (FAD) gene family in G. hirsutum, with an emphasis on identifying omega-3 FADs involved in cold temperature adaptation.

RESULTS

Eleven omega-3 FAD genes were identified in G. hirsutum, and characterization of the gene family in extant A and D diploid species (G. herbaceum and G. raimondii, respectively) allowed for unambiguous genome assignment of all homoeologs in tetraploid G. hirsutum. The omega-3 FAD family of cotton includes five distinct genes, two of which encode endoplasmic reticulum-type enzymes (FAD3-1 and FAD3-2) and three that encode chloroplast-type enzymes (FAD7/8-1, FAD7/8-2, and FAD7/8-3). The FAD3-2 gene was duplicated in the A genome progenitor species after the evolutionary split from the D progenitor, but before the interspecific hybridization event that gave rise to modern tetraploid cotton. RNA-seq analysis revealed conserved, gene-specific expression patterns in various organs and cell types and semi-quantitative RT-PCR further revealed that FAD7/8-1 was specifically induced during cold temperature treatment of G. hirsutum seedlings.

CONCLUSIONS

The omega-3 FAD gene family in cotton was characterized at the genome-wide level in three species, showing relatively ancient establishment of the gene family prior to the split of A and D diploid progenitor species. The FAD genes are differentially expressed in various organs and cell types, including fiber, and expression of the FAD7/8-1 gene was induced by cold temperature. Collectively, these data define the genetic and functional genomic properties of this important gene family in cotton and provide a foundation for future efforts to improve cotton abiotic stress tolerance through molecular breeding approaches.

Understanding the relationship between cotton fiber properties and non-cellulosic cell wall polysaccharides

A detailed knowledge of cell wall heterogeneity and complexity is crucial for understanding plant growth and development. One key challenge is to establish links between polysaccharide-rich cell walls and their phenotypic characteristics. It is of particular interest for some plant material, like cotton fibers, which are of both biological and industrial importance. To this end, we attempted to study cotton fiber characteristics together with glycan arrays using regression based approaches. Taking advantage of the comprehensive microarray polymer profiling technique (CoMPP), 32 cotton lines from different cotton species were studied. The glycan array was generated by sequential extraction of cell wall polysaccharides from mature cotton fibers and screening samples against eleven extensively characterized cell wall probes. Also, phenotypic characteristics of cotton fibers such as length, strength, elongation and micronaire were measured. The relationship between the two datasets was established in an integrative manner using linear regression methods. In the conducted analysis, we demonstrated the usefulness of regression based approaches in establishing a relationship between glycan measurements and phenotypic traits. In addition, the analysis also identified specific polysaccharides which may play a major role during fiber development for the final fiber characteristics. Three different regression methods identified a negative correlation between micronaire and the xyloglucan and homogalacturonan probes. Moreover, homogalacturonan and callose were shown to be significant predictors for fiber length. The role of these polysaccharides was already pointed out in previous cell wall elongation studies. Additional relationships were predicted for fiber strength and elongation which will need further experimental validation.

Doubling down on genomes: polyploidy and crop plants

Polyploidy, or whole genome multiplication, is ubiquitous among angiosperms. Many crop species are relatively recent allopolyploids, resulting from interspecific hybridization and polyploidy. Thus, an appreciation of the evolutionary consequences of (allo)polyploidy is central to our understanding of crop plant domestication, agricultural improvement, and the evolution of angiosperms in general. Indeed, many recent insights into plant biology have been gleaned from polyploid crops, including, but not limited to wheat, tobacco, sugarcane, apple, and cotton. A multitude of evolutionary processes affect polyploid genomes, including rapid and substantial genome reorganization, transgressive gene expression alterations, gene fractionation, gene conversion, genome downsizing, and sub- and neofunctionalization of duplicate genes. Often these genomic changes are accompanied by heterosis, robustness, and the improvement of crop yield, relative to closely related diploids. Historically, however, the genome-wide analysis of polyploid crops has lagged behind those of diploid crops and other model organisms. This lag is partly due to the difficulties in genome assembly, resulting from the genomic complexities induced by combining two or more evolutionarily diverged genomes into a single nucleus and by the significant size of polyploid genomes. In this review, we explore the role of polyploidy in angiosperm evolution, the domestication process and crop improvement. We focus on the potential of modern technologies, particularly next-generation sequencing, to inform us on the patterns and processes governing polyploid crop improvement and phenotypic change subsequent to domestication.

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.

Ancient gene duplicates in Gossypium (cotton) exhibit near-complete expression divergence

Whole genome duplication (WGD) is widespread in flowering plants and is a driving force in angiosperm diversification. The redundancy introduced by WGD allows the evolution of novel gene interactions and functions, although the patterns and processes of diversification are poorly understood. We identified ∼ 2,000 pairs of paralogous genes in Gossypium raimondii (cotton) resulting from an approximately 60 My old 5- to 6-fold ploidy increase. Gene expression analyses revealed that, in G. raimondii, 99.4% of the gene pairs exhibit differential expression in at least one of the three tissues (petal, leaf, and seed), with 93% to 94% exhibiting differential expression on a per-tissue basis. For 1,666 (85%) pairs, differential expression was observed in all tissues. These observations were mirrored in a time series of G. raimondii seed, and separately in leaf, petal, and seed of G. arboreum, indicating expression level diversification before species divergence. A generalized linear model revealed 92.4% of the paralog pairs exhibited expression divergence, with most exhibiting significant gene and tissue interactions indicating complementary expression patterns in different tissues. These data indicate massive, near-complete expression level neo- and/or subfunctionalization among ancient gene duplicates, suggesting these processes are essential in their maintenance over ∼ 60 Ma.

Comparative evolutionary and developmental dynamics of the cotton (Gossypium hirsutum) fiber transcriptome

The single-celled cotton (Gossypium hirsutum) fiber provides an excellent model to investigate how human selection affects phenotypic evolution. To gain insight into the evolutionary genomics of cotton domestication, we conducted comparative transcriptome profiling of developing cotton fibers using RNA-Seq. Analysis of single-celled fiber transcriptomes from four wild and five domesticated accessions from two developmental time points revealed that at least one-third and likely one-half of the genes in the genome are expressed at any one stage during cotton fiber development. Among these, ∼5,000 genes are differentially expressed during primary and secondary cell wall synthesis between wild and domesticated cottons, with a biased distribution among chromosomes. Transcriptome data implicate a number of biological processes affected by human selection, and suggest that the domestication process has prolonged the duration of fiber elongation in modern cultivated forms. Functional analysis suggested that wild cottons allocate greater resources to stress response pathways, while domestication led to reprogrammed resource allocation toward increased fiber growth, possibly through modulating stress-response networks. This first global transcriptomic analysis using multiple accessions of wild and domesticated cottons is an important step toward a more comprehensive systems perspective on cotton fiber evolution. The understanding that human selection over the past 5,000+ years has dramatically re-wired the cotton fiber transcriptome sets the stage for a deeper understanding of the genetic architecture underlying cotton fiber synthesis and phenotypic evolution.

Ever since Darwin biologists have recognized that comparative study of crop plants and their wild relatives offers a powerful framework for generating insights into the mechanisms that underlie evolutionary change. Here, we study the domestication process in cotton, Gossypium hirsutum, an allopolyploid species (containing two different genomes) which initially was domesticated approximately 5000 years ago, and which primarily is grown for its single-celled seed fibers. Strong directional selection over the millennia was accompanied by transformation of the short, coarse, and brown fibers of wild plants into the long, strong, and fine white fibers of the modern cotton crop plant. To explore the evolutionary genetics of cotton domestication, we conducted transcriptome profiling of developing cotton fibers from multiple accessions of wild and domesticated cottons. Comparative analysis revealed that the domestication process dramatically rewired the transcriptome, affecting more than 5,000 genes, and with a more evenly balanced usage of the duplicated copies arising from genome doubling. We identify many different biological processes that were involved in this transformation, including those leading to a prolongation of fiber elongation and a reallocation of resources toward increased fiber growth in modern forms. The data provide a rich resource for future functional analyses targeting crop improvement and evolutionary objectives.