Insights into the Ecology and Evolution of Polyploid Plants through Network Analysis

Friends without benefits: Extensive cytotype sympatry and polyploid persistence in an African geophyte

PREMISE

Polyploidy is a major factor in plant adaptation and speciation. Multiple mechanisms contribute to autopolyploid frequency within populations, but uncertainties remain regarding mechanisms that facilitate polyploid establishment and persistence. Here we aimed to document and predict cytotype distributions of Oxalis obliquifolia Steud. ex A. Rich. across Gauteng, South Africa, and test for evidence of possible mechanisms, including morphological, phenological, and reproductive traits, that may potentially facilitate polyploid persistence.

METHODS

Over 320 O. obliquifolia plants from 25 sites were cytotyped using flow cytometry, and DNA ploidy was confirmed using meiotic chromosome squashes. Cytotypes were mapped and correlations with abiotic variables assessed using ordinations. To assess morphological and phenological associations with cytotype, we grew multiple cytotypes in a common garden, measured phenotypic traits and compared them using linear models and discriminant analyses. Intercytotype reproductive isolation was assessed using crossing experiments, and AMOVAs based on ITS DNA sequences tested for cytogeographic structure.

RESULTS

Six cytotypes were identified, and most sites had multiple cytotypes. Abiotic variables were not predictive of cytotype distribution. A clear gigas effect was present. Differences in flower size and phenology suggested pollinator interactions could play a role in polyploid persistence. Intercytotype crosses produced seed at low frequency. DNA data suggested diploids and polyploids were largely reproductively isolated in situ, and polyploidization events were not frequent enough to explain high cytotype sympatry.

CONCLUSIONS

Diploids and polyploids are behaving as separate species, despite little observable niche differentiation and non-zero potential intercytotype seed set. Tests on biotic interactions and intercytotype F1 fitness may provide insights into diploid and polyploid coexistence.

Polyploidy and the evolution of phenotypic integration: Network analysis reveals relationships among anatomy, morphology, and physiology

Premise

Most traits are polygenic and most genes are pleiotropic, resulting in complex, integrated phenotypes. Polyploidy presents an excellent opportunity to explore the evolution of phenotypic integration as entire genomes are duplicated, allowing for new associations among traits and potentially leading to enhanced or reduced phenotypic integration. Despite the multivariate nature of phenotypic evolution, studies often rely on simplistic bivariate correlations that cannot accurately represent complex phenotypes or data reduction techniques that can obscure specific trait relationships.

Methods

We apply network modeling, a common gene co-expression analysis, to the study of phenotypic integration to identify multivariate patterns of phenotypic evolution, including anatomy and morphology (structural) and physiology (functional) traits in response to whole genome duplication in the genus Brassica.

Results

We identify four key structural traits that are overrepresented in the evolution of phenotypic integration. Seeding networks with key traits allowed us to identify structure-function relationships not apparent from bivariate analyses. In general, allopolyploids exhibited larger, more robust networks indicative of increased phenotypic integration compared to diploids.

Discussion

Phenotypic network analysis may provide important insights into the effects of selection on non-target traits, even when they lack direct correlations with the target traits. Network analysis may allow for more nuanced predictions of both natural and artificial selection.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

Genome-wide identification and comparative analysis of CLE family in rapeseed and its diploid progenitors

Crop genomics and breeding CLAVATA3/EMBRYO SURROUNDING REGION-RELATED (CLE) proteins belong to a small peptide family in plants. During plant development, CLE gene family members play a pivotal role in regulating cell-to-cell communication and stem cell maintenance. However, the evolutionary process and functional importance of CLEs are unclear in Brassicaceae. In this study, a total of 70 BnCLEs were identified in Brassica napus (2n = 4x = 38, A_nC_n): 32 from the A_n subgenome, 36 from the C_n subgenome, and 2 from the unanchored subgenome. Meanwhile, 29 BrCLE and 32 BoCLE genes were explored in Brassica rapa (2n = 2x = 20, A_r) and Brassica oleracea (2n = 2x = 18, C_o). Phylogenetic analysis revealed that 163 CLEs derived from three Brassica species and Arabidopsis thaliana can be divided into seven subfamilies. Homology and synteny analyses indicated whole-genome triplication (WGT) and segmental duplication may be the major contributors to the expansion of CLE family. In addition, RNA-seq and qPCR analysis indicated that 19 and 16 BnCLEs were more highly expressed in immature seeds and roots than in other tissues. Some CLE gene pairs exhibited different expression patterns in the same tissue, which indicated possible functional divergence. Furthermore, genetic variations and regional association mapping analysis indicated that 12 BnCLEs were potential genes for regulating important agronomic traits. This study provided valuable information to understand the molecular evolution and biological function of CLEs in B. napus and its diploid progenitors, which will be helpful for genetic improvement of high-yield breeding in B. napus.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Integrative analyses of targeted metabolome and transcriptome of Isatidis Radix autotetraploids highlighted key polyploidization-responsive regulators

Background

Isatidis Radix, the root of Isatis indigotica Fort. (Chinese woad) can produce a variety of efficacious compound with medicinal properties. The tetraploid I. indigotica plants exhibit superior phenotypic traits, such as greater yield, higher bioactive compounds accumulation and enhanced stress tolerance. In this study, a comparative transcriptomic and metabolomic study on Isatidis Radix autotetraploid and its progenitor was performed.

Results

Through the targeted metabolic profiling, 283 metabolites were identified in Isatidis Radix, and 70 polyploidization-altered metabolites were obtained. Moreover, the production of lignans was significantly increased post polyploidization, which implied that polyploidization-modulated changes in lignan biosynthesis. Regarding the transcriptomic shift, 2065 differentially expressed genes (DEGs) were identified as being polyploidy-responsive genes, and the polyploidization-altered DEGs were enriched in phenylpropanoid biosynthesis and plant hormone signal transduction. The further integrative analysis of polyploidy-responsive metabolome and transcriptome showed that 1584 DEGs were highly correlated with the 70 polyploidization-altered metabolites, and the transcriptional factors TFs-lignans network highlighted 10 polyploidy-altered TFs and 17 fluctuated phenylpropanoid pathway compounds.

Conclusions

These results collectively indicated that polyploidization contributed to the high content of active compounds in autotetraploid roots, and the gene–lignan pathway network analysis highlighted polyploidy–responsive key functional genes and regulators.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07980-w.

Intricate Distribution Patterns of Six Cytotypes of Allium oleraceum at a Continental Scale: Niche Expansion and Innovation Followed by Niche Contraction With Increasing Ploidy Level

The establishment and success of polyploids are thought to often be facilitated by ecological niche differentiation from diploids. Unfortunately, most studies compared diploids and polyploids, ignoring variation in ploidy level in polyploids. To fill this gap, we performed a large-scale study of 11,163 samples from 1,283 populations of the polyploid perennial geophyte Allium oleraceum with reported mixed-ploidy populations, revealed distribution ranges of cytotypes, assessed their niches and explored the pattern of niche change with increasing ploidy level. Altogether, six ploidy levels (3x-8x) were identified. The most common were pentaploids (53.6%) followed by hexaploids (22.7%) and tetraploids (21.6%). Higher cytotype diversity was found at lower latitudes than at higher latitudes (>52° N), where only tetraploids and pentaploids occurred. We detected 17.4% of mixed-ploidy populations, usually as a combination of two, rarely of three, cytotypes. The majority of mixed-ploidy populations were found in zones of sympatry of the participating cytotypes, suggesting they have arisen through migration (secondary contact zone). Using coarse-grained variables (climate, soil), we found evidence of both niche expansion and innovation in tetraploids related to triploids, whereas higher ploidy levels showed almost zero niche expansion, but a trend of increased niche unfilling of tetraploids. Niche unfilling in higher ploidy levels was caused by a contraction of niche envelopes toward lower continentality of the climate and resulted in a gradual decrease of niche breadth and a gradual shift in niche optima. Field-recorded data indicated wide habitat breadth of tetraploids and pentaploids, but also a pattern of increasing synanthropy in higher ploidy levels. Wide niche breadth of tetra- and pentaploids might be related to their multiple origins from different environmental conditions, higher "age", and retained sexuality, which likely preserve their adaptive potential. In contrast, other cytotypes with narrower niches are mostly asexual, probably originating from a limited range of contrasting environments. Persistence of local ploidy mixtures could be enabled by the perenniality of A. oleraceum and its prevalence of vegetative reproduction, facilitating the establishment and decreasing exclusion of minority cytotype due to its reproductive costs. Vegetative reproduction might also significantly accelerate colonization of new areas, including recolonization of previously glaciated areas.

An optimised chromatin immunoprecipitation (ChIP) method for starchy leaves of Nicotiana benthamiana to study histone modifications of an allotetraploid plant

All flowering plants have evolved through multiple rounds of polyploidy throughout the evolutionary process. Intergenomic interactions between subgenomes in polyploid plants are predicted to induce chromatin modifications such as histone modifications to regulate expression of gene homoeologs. Nicotiana benthamiana is an ancient allotetraploid plant with ecotypes collected from climatically diverse regions of Australia. Studying the chromatin landscape of this unique collection will likely shed light on the importance of chromatin modifications in gene regulation in polyploids as well its implications in adaptation of plants in environmentally diverse conditions. Generally, chromatin immunoprecipitation and high throughput DNA sequencing (ChIP-seq) is used to study chromatin modifications. However, due to the starchy nature of mature N. benthamiana leaves, previously published protocols were unsuitable. The higher amounts of starch in leaves that co-precipitated with nuclei hindered downstream processing of DNA. Here we present an optimised ChIP protocol for N. benthamiana leaves to facilitate comparison of chromatin modifications in two closely related ecotypes. Several steps of ChIP were optimised including tissue harvesting, nuclei isolation, nuclei storage, DNA shearing and DNA recovery. Commonly available antibodies targeting histone 3 lysine 4 trimethylation (H3K4me3) and histone 3 lysine 9 dimethylation (H3K9me2) histone modifications were used and success of ChIP was confirmed by PCR and next generation sequencing. Collectively, our optimised method is the first comprehensive ChIP method for mature starchy leaves of N. benthamiana to enable studies of chromatin landscape at the genome-wide scale.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

Polyploidy in liver development, homeostasis and disease

Polyploidy (or whole-genome duplication) is the condition of having more than two basic sets of chromosomes. Polyploidization is well tolerated in many species and can lead to specific biological functions. In mammals, programmed polyploidization takes place during development in certain tissues, such as the heart and placenta, and is considered a feature of differentiation. However, unscheduled polyploidization can cause genomic instability and has been observed in pathological conditions, such as cancer. Polyploidy of the liver parenchyma was first described more than 100 years ago. The liver is one of the few mammalian organs that display changes in polyploidy during homeostasis, regeneration and in response to damage. In the human liver, approximately 30% of hepatocytes are polyploid. The polyploidy of hepatocytes results from both nuclear polyploidy (an increase in the amount of DNA per nucleus) and cellular polyploidy (an increase in the number of nuclei per cell). In this Review, we discuss the regulation of polyploidy in liver development and pathophysiology. We also provide an overview of current knowledge about the mechanisms of hepatocyte polyploidization, its biological importance and the fate of polyploid hepatocytes during liver tumorigenesis.

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

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

[Hepatic polyploidy: Dr Jekyll or Mr Hyde]

Polyploidy (alias whole genome amplification) refers to organisms containing more than two basic sets of chromosomes. Polyploidy was first observed in plants more than a century ago, and it is known that such processes occur in many eukaryotes under a variety of circumstances. In mammals, the development of polyploid cells can contribute to tissue differentiation and therefore possibly a gain of function. Alternately, it can be associated with development of disease such as cancer. Polyploidy can occur because of cell fusion or abnormal cell division. Polyploidy is a common characteristic of the mammalian liver. Polyploidization occurs notably during liver development, but also in adults because of cellular stress. Recent progresses have unraveled the mechanisms and functional consequences of hepatocytes polyploidization during normal and pathological liver growth.

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.

Integrating Networks, Phylogenomics, and Population Genomics for the Study of Polyploidy

Duplication events are regarded as sources of evolutionary novelty, but our understanding of general trends for the long-term trajectory of additional genomic material is still lacking. Organisms with a history of whole genome duplication (WGD) offer a unique opportunity to study potential trends in the context of gene retention and/or loss, gene and network dosage, and changes in gene expression. In this review, we discuss the prevalence of polyploidy across the tree of life, followed by an overview of studies investigating genome evolution and gene expression. We then provide an overview of methods in network biology, phylogenomics, and population genomics that are critical for advancing our understanding of evolution post-WGD, highlighting the need for models that can accommodate polyploids. Finally, we close with a brief note on the importance of random processes in the evolution of polyploids with respect to neutral versus selective forces, ancestral polymorphisms, and the formation of autopolyploids versus allopolyploids.

Non-Additive Transcriptomic Responses to Inoculation with Rhizobia in a Young Allopolyploid Compared with Its Diploid Progenitors

Root nodule symbioses (nodulation) and whole genome duplication (WGD, polyploidy) are both important phenomena in the legume family (Leguminosae). Recently, it has been proposed that polyploidy may have played a critical role in the origin or refinement of nodulation. However, while nodulation and polyploidy have been studied independently, there have been no direct studies of mechanisms affecting the interactions between these phenomena in symbiotic, nodule-forming species. Here, we examined the transcriptome-level responses to inoculation in the young allopolyploid Glycine dolichocarpa (T2) and its diploid progenitor species to identify underlying processes leading to the enhanced nodulation responses previously identified in T2. We assessed the differential expression of genes and, using weighted gene co-expression network analysis (WGCNA), identified modules associated with nodulation and compared their expression between species. These transcriptomic analyses revealed patterns of non-additive expression in T2, with evidence of transcriptional responses to inoculation that were distinct from one or both progenitors. These differential responses elucidate mechanisms underlying the nodulation-related differences observed between T2 and the diploid progenitors. Our results indicate that T2 has reduced stress-related transcription, coupled with enhanced transcription of modules and genes implicated in hormonal signaling, both of which are important for nodulation.

The evolutionary significance of polyploidy

Polyploidy, or the duplication of entire genomes, has been observed in prokaryotic and eukaryotic organisms, and in somatic and germ cells. The consequences of polyploidization are complex and variable, and they differ greatly between systems (clonal or non-clonal) and species, but the process has often been considered to be an evolutionary 'dead end'. Here, we review the accumulating evidence that correlates polyploidization with environmental change or stress, and that has led to an increased recognition of its short-term adaptive potential. In addition, we discuss how, once polyploidy has been established, the unique retention profile of duplicated genes following whole-genome duplication might explain key longer-term evolutionary transitions and a general increase in biological complexity.

The role of polyploidy in shaping morphological diversity in natural populations of Phlox amabilis

PREMISE OF THE STUDY

Studies of natural populations of polyploids increasingly highlight complex patterns of variation in ploidy and geographic distribution of cytotypes. As our understanding of the complexity of polyploidy grows, our understanding of the morphological correlates of polyploidy should expand as well. Here we examine in what ways, and to what degree, polyploidy affects the overall phenotype of a species across its distribution when there are three ploidies and geographic complexity in cytotype distribution.

METHODS

We measured 31 morphological traits from stems, leaves, and flowers from up to 25 individuals from 11 sites across the distribution of Phlox amabilis. Chromosome counts and flow cytometry confirmed and expanded upon earlier research documenting diploid, tetraploid, and hexaploid populations, and also identified a site with two ploidies. Univariate and multivariate statistics were used to characterize the morphological effects of polyploidy.

KEY RESULTS

We detected significant associations between morphology and ploidy in 11 traits spread across vegetative and reproductive structures. Generally, diploid individuals differed from polyploid individuals to a greater extent, and in different ways, than tetraploid and hexaploid plants differed from each other. Multivariate morphometrics demonstrated that the two primary axes of overall variation are driven by morphological traits associated with polyploidy, and individuals of different ploidies can be discriminated with 95% success.

CONCLUSIONS

Polyploidy plays a major role in shaping overall morphological diversity in natural populations of P. amabilis.