Jeremy E. Coate

Quantifying Whole Transcriptome Size, a Prerequisite for Understanding Transcriptome Evolution Across Species: An Example from a Plant Allopolyploid

Evolutionary biologists are increasingly comparing gene expression patterns across species. Due to the way in which expression assays are normalized, such studies provide no direct information about expression per gene copy (dosage responses) or per cell and can give a misleading picture of genes that are differentially expressed. We describe an assay for estimating relative expression per cell. When used in conjunction with transcript profiling data, it is possible to compare the sizes of whole transcriptomes, which in turn makes it possible to compare expression per cell for each gene in the transcript profiling data set. We applied this approach, using quantitative reverse transcriptase-polymerase chain reaction and high throughput RNA sequencing, to a recently formed allopolyploid and showed that its leaf transcriptome was approximately 1.4-fold larger than either progenitor transcriptome (70% of the sum of the progenitor transcriptomes). In contrast, the allopolyploid genome is 94.3% as large as the sum of its progenitor genomes and retains ≥93.5% of the sum of its progenitor gene complements. Thus, “transcriptome downsizing” is greater than genome downsizing. Using this transcriptome size estimate, we inferred dosage responses for several thousand genes and showed that the majority exhibit partial dosage compensation. Homoeologue silencing is nonrandomly distributed across dosage responses, with genes showing extreme responses in either direction significantly more likely to have a silent homoeologue. This experimental approach will add value to transcript profiling experiments involving interspecies and interploidy comparisons by converting expression per transcriptome to expression per genome, eliminating the need for assumptions about transcriptome size.

The legacy of diploid progenitors in allopolyploid gene expression patterns

Allopolyploidization (hybridization and whole-genome duplication) is a common phenomenon in plant evolution with immediate saltational effects on genome structure and gene expression. New technologies have allowed rapid progress over the past decade in our understanding of the consequences of allopolyploidy. A major question, raised by early pioneer of this field Leslie Gottlieb, concerned the extent to which gene expression differences among duplicate genes present in an allopolyploid are a legacy of expression differences that were already present in the progenitor diploid species. Addressing this question necessitates phylogenetically well-understood natural study systems, appropriate technology, availability of genomic resources and a suitable analytical framework, including a sufficiently detailed and generally accepted terminology. Here, we review these requirements and illustrate their application to a natural study system that Gottlieb worked on and recommended for this purpose: recent allopolyploids of Tragopogon (Asteraceae). We reanalyse recent data from this system within the conceptual framework of parental legacies on duplicate gene expression in allopolyploids. On a broader level, we highlight the intellectual connection between Gottliebs phrasing of this issue and the more contemporary framework of cis- versus trans-regulation of duplicate gene expression in allopolyploid plants.

Comparative Evolution of Photosynthetic Genes in Response to Polyploid and Nonpolyploid Duplication

The likelihood of duplicate gene retention following polyploidy varies by functional properties (e.g. gene ontologies or protein family domains), but little is known about the effects of whole-genome duplication on gene networks related by a common physiological process. Here, we examined the effects of both polyploid and nonpolyploid duplications on genes encoding the major functional groups of photosynthesis (photosystem I, photosystem II, the light-harvesting complex, and the Calvin cycle) in the cultivated soybean (Glycine max), which has experienced two rounds of whole-genome duplication. Photosystem gene families exhibit retention patterns consistent with dosage sensitivity (preferential retention of polyploid duplicates and elimination of nonpolyploid duplicates), whereas Calvin cycle and light-harvesting complex gene families do not. We observed similar patterns in barrel medic (Medicago truncatula), which shared the older genome duplication with soybean but has evolved independently for approximately 50 million years, and in Arabidopsis (Arabidopsis thaliana), which experienced two nested polyploidy events independent from the legume duplications. In both soybean and Arabidopsis, Calvin cycle gene duplicates exhibit a greater capacity for functional differentiation than do duplicates within the photosystems, which likely explains the greater retention of ancient, nonpolyploid duplicates and larger average gene family size for the Calvin cycle relative to the photosystems.

Lessons from Natural and Artificial Polyploids in Higher Plants

Polyploidy in higher plants is a major source of genetic novelty upon which selection may act to drive evolution, as evidenced by the widespread success of polyploid species in the wild. However, research into the effects of polyploidy can be confounded by the entanglement of several processes: genome duplication, hybridisation (allopolyploidy is frequent in plants) and subsequent evolution. The discovery of the chemical agent colchicine, which can be used to produce artificial polyploids on demand, has enabled scientists to unravel these threads and understand the complex genomic changes involved in each. We present here an overview of lessons learnt from studies of natural and artificial polyploids, and from comparisons between the 2, covering basic cellular and metabolic consequences through to alterations in epigenetic gene regulation, together with 2 in-depth case studies in Senecio and Glycine. See also the sister article focusing on animals by Arai and Fujimoto in this themed issue.

Extensive Translational Regulation of Gene Expression in an Allopolyploid (Glycine dolichocarpa)

This work shows that translational regulation of gene expression in a recently formed allopolyploid is widespread, reduces transcriptional differences between the polyploid and its diploid progenitors, and correlates with the retention of genes from an older polyploidy event. These findings suggest that translational regulation is significant in both early and long-term responses to polyploidy. All flowering plants have experienced repeated rounds of polyploidy (whole-genome duplication), which has in turn driven the evolution of novel phenotypes and ecological tolerances and been a major driver of speciation. The effects of polyploidy on gene expression have been studied extensively at the level of transcription and, to a much lesser extent, at the level of the steady state proteome, but not at the level of translation. We used polysome profiling by RNA-Seq to quantify translational regulation of gene expression in a recently formed (∼100,000 years ago) allotetraploid (Glycine dolichocarpa) closely related to the cultivated soybean (Glycine max). We show that there is a high level of concordance between the allopolyploid transcriptome and translatome overall but that at least one-quarter of the transcriptome is translationally regulated. We further show that translational regulation preferentially targets genes involved in transcription, translation, and photosynthesis, causes regional and possibly whole-chromosome shifts in expression bias between duplicated genes (homoeologs), and reduces transcriptional differences between the polyploid and its diploid progenitors, possibly attenuating misregulation resulting from genome merger and/or doubling. Finally, translational regulation correlates positively with long-term retention of homoeologs from a paleopolyploidy event, suggesting that it plays a significant role in polyploid evolution.

Transgressive physiological and transcriptomic responses to light stress in allopolyploid Glycine dolichocarpa (Leguminosae).

Allopolyploidy is often associated with increased photosynthetic capacity as well as enhanced stress tolerance. Excess light is a ubiquitous plant stress associated with photosynthetic light harvesting. We show that under chronic excess light, the capacity for non-photochemical quenching (NPQmax), a photoprotective mechanism, was higher in a recently formed natural allotetraploid (Glycine dolichocarpa, designated ‘T2’) than in its diploid progenitors (G. tomentella, ‘D3’; and G. syndetika, ‘D4’). This enhancement in NPQmax was due to an increase in energy-dependent quenching (qE) relative to D3, combined with an increase in zeaxanthin-dependent quenching (qZ) relative to D4. To explore the genetic basis for this phenotype, we profiled D3, D4 and T2 leaf transcriptomes and found that T2 overexpressed genes of the water–water cycle relative to both diploid progenitors, as well as genes involved in cyclic electron flow around photosystem I (CEF-PSI) and the xanthophyll cycle, relative to D4. Xanthophyll pigments have critical roles in NPQ, and the water–water cycle and CEF-PSI are non-photosynthetic electron transport pathways believed to facilitate NPQ formation. In the absence of CO2, T2 also exhibited greater quantum yield of photosystem II than either diploid, indicating a greater capacity for non-photosynthetic electron transport. We postulate that, relative to its diploid progenitors, T2 is able to achieve higher NPQmax due to an increase in xanthophyll pigments coupled with enhanced electron flow through the water–water cycle and CEF-PSI.

Anatomical, biochemical, and photosynthetic responses to recent allopolyploidy in Glycine dolichocarpa (Fabaceae)

PREMISE OF THE STUDY Previous studies have shown that polyploidy has pronounced effects on photosynthesis. Most of these studies have focused on synthetic or recently formed autopolyploids, and comparatively little is known about the integrated effects of natural allopolyploidy, which involves hybridity and genome doubling and often incorporates multiple genotypes through recurrent origins and lineage recombination. METHODS Glycine dolichocarpa (designated T2) is a natural allotetraploid with multiple origins. We quantified 21 anatomical, biochemical, and physiological phenotypes relating to photosynthesis in T2 and its diploid progenitors, G. tomentella (D3) and G. syndetika (D4). To assess how direction of cross affects these phenotypes, we included three T2 accessions having D3-like plastids (T2(D3)) and two accessions having D4-like plastids (T2(D4)). KEY RESULTS T2 accessions were transgressive (more extreme than any diploid accession) for 17 of 21 phenotypes, and species means differed significantly in T2 vs. both progenitors for four of 21 phenotypes (higher for guard cell length, electron transport capacity [J(max)] per palisade cell, and J(max) per mesophyll cell; lower for palisade cells per unit leaf area). Within T2, four of 21 parameters differed significantly between T2(D3) and T2(D4) (palisade cell volume; chloroplast number and volume per unit leaf area; and J(max) per unit leaf area). CONCLUSIONS T2 is characterized by transgressive photosynthesis-related phenotypes (including an ca. 2-fold increase in J(max) per cell), as well as by significant intraspecies variation correlating with plastid type. These data indicate prominent roles for both nucleotypic effects and cytoplasmic factors in photosynthetic responses to allopolyploidy.

The wild side of a major crop: soybean's perennial cousins from Down Under.

The accumulation of over 30 years of basic research on the biology, genetic variation, and evolution of the wild perennial relatives of soybean (Glycine max) provides a foundation to improve cultivated soybean. The cultivated soybean and its wild progenitor, G. soja, have a center of origin in eastern Asia and are the only two species in the annual subgenus Soja. Systematic and evolutionary studies of the ca. 30 perennial species of subgenus Glycine, native to Australia, have benefited from the availability of the G. max genomic sequence. The perennial species harbor many traits of interest to soybean breeders, among them resistance to major soybean pathogens such as cyst nematode and leaf rust. New species in the Australian subgenus continue to be described, due to the collection of new material and to insights gleaned through systematic studies of accessions in germplasm collections. Ongoing studies in perennial species focus on genomic regions that contain genes for key traits relevant to soybean breeding. These comparisons also include the homoeologous regions that are the result of polyploidy in the common ancestor of all Glycine species. Subgenus Glycine includes a complex of recently formed allopolyploids that are the focus of studies aimed at elucidating genomic, transcriptomic, physiological, taxonomic, morphological, developmental, and ecological processes related to polyploid evolution. Here we review what has been learned over the past 30 years and outline ongoing work on photosynthesis, nitrogen fixation, and floral biology, much of it drawing on new technologies and resources.

Variation in transcriptome size: are we getting the message?

The number of RNA molecules per cell (transcriptome size) is highly variable, differing among and within cell types depending on cell size, stage of the cell cycle, ploidy level, age, disease state, and growth condition. Such variation has been observed at the level of total RNA, ribosomal RNA, messenger RNA (mRNA), and the polyadenylated fraction of mRNA, and these distinct RNA species can also vary in abundance with respect to each other. This variation in transcriptome size has been largely ignored or overlooked, and in fact, standard data normalization procedures for transcript profiling experiments implicitly assume that mRNA transcriptome size is constant. Consequently, variation in transcriptome size has important technical implications for such experiments, as well as profound biological implications for the affected cells and underlying genomes. Here, we review what is known about transcriptome size variation, explore how ignoring this variation introduces systematic bias into standard expression profiling experiments, and present examples of how such biases have led to erroneous conclusions in expression studies of sex chromosome dosage compensation, cancer, Rett syndrome, embryonic development, aging, and polyploidy. We also discuss how quantifying transcriptome size will help to elucidate the selective forces underlying patterns of gene and genome evolution and review the evidence that cells exert tight control over transcriptome size in order to maintain cell size homeostasis and to optimize chemical reactions within the cell, such that loss of control over transcriptome size is associated with cancer and aging. Thus, transcriptome size is an important phenotype in its own right. Finally, we discuss strategies for quantifying transcriptome size and individual gene dosage responses in order to account for and better understand this important biological phenomenon.

Mining transcriptomic data to study the origins and evolution of a plant allopolyploid complex

Allopolyploidy combines two progenitor genomes in the same nucleus. It is a common speciation process, especially in plants. Deciphering the origins of polyploid species is a complex problem due to, among other things, extinct progenitors, multiple origins, gene flow between different polyploid populations, and loss of parental contributions through gene or chromosome loss. Among the perennial species of Glycine, the plant genus that includes the cultivated soybean (G. max), are eight allopolyploid species, three of which are studied here. Previous crossing studies and molecular systematic results from two nuclear gene sequences led to hypotheses of origin for these species from among extant diploid species. We use several phylogenetic and population genomics approaches to clarify the origins of the genomes of three of these allopolyploid species using single nucleotide polymorphism data and a guided transcriptome assembly. The results support the hypothesis that all three polyploid species are fixed hybrids combining the genomes of the two putative parents hypothesized on the basis of previous work. Based on mapping to the soybean reference genome, there appear to be no large regions for which one homoeologous contribution is missing. Phylogenetic analyses of 27 selected transcripts using a coalescent approach also are consistent with multiple origins for these allopolyploid species, and suggest that origins occurred within the last several hundred thousand years.