Polyploidy: adaptation to the genomic environment

Expanding the toolkit for ploidy manipulation in Chlamydomonas reinhardtii

Whole-genome duplications, widely observed in plant lineages, have significant evolutionary and ecological impacts. Yet, our current understanding of the direct implications of ploidy shifts on short- and long-term plant evolution remains fragmentary, necessitating further investigations across multiple ploidy levels. Chlamydomonas reinhardtii is a valuable model organism with profound potential to study the impact of ploidy increase on the longer term in a laboratory environment. This is partly due to the ability to increase the ploidy level. We developed a strategy to engineer ploidy in C. reinhardtii using noninterfering, antibiotic, selectable markers. This approach allows us to induce higher ploidy levels in C. reinhardtii and is applicable to field isolates, which expands beyond specific auxotroph laboratory strains and broadens the genetic diversity of parental haploid strains that can be crossed. We implement flow cytometry for precise measurement of the genome size of strains of different ploidy. We demonstrate the creation of diploids, triploids, and tetraploids by engineering North American field isolates, broadening the application of synthetic biology principles in C. reinhardtii. However, our newly formed triploids and tetraploids show signs of rapid aneuploidization. Our study greatly facilitates the application of C. reinhardtii to study polyploidy, in both fundamental and applied settings.

Overexpression of GFP Fusions of Regulators of the SAC from Arabidopsis thaliana

Cell division is a fundamental biological process, essential for sustaining life on Earth. Accurate replication followed by uniform segregation of the genome is required to ensure cell division is sustainable and reduces the likelihood of aneuploidy. The cell cycle has various checkpoints to safeguard proper replication, for example, the spindle assembly checkpoint (SAC) which ensures that all chromosomes are correctly aligned and attached to the spindle, before the transition to anaphase. The precise function of the SAC and SAC components in plants is so far unclear. First, the high level of polyploidy in plants raises concerns about the efficacy of the SAC. Second, many plant SAC components are implicated in other cellular processes, such as MAD1, which has been implicated in the reproductive transition of Arabidopsis thaliana. Overexpression of GFP fusions of core SAC components provides a key route to establish the functions of the different SAC components in plants. Here we describe two methods for agrobacterium-mediated transformation of plants.

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological, and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb 'Qinjiao'. Herein, we present its first chromosome-level genome sequence assembly and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5 Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70 Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

Phased chromosome-scale genome assembly of an asexual, allopolyploid root-knot nematode reveals complex subgenomic structure

We present the chromosome-scale genome assembly of the allopolyploid root-knot nematode Meloidogyne javanica. We show that the M. javanica genome is predominantly allotetraploid, comprising two subgenomes, A and B, that most likely originated from hybridisation of two ancestral parental species. The assembly was annotated using full-length non-chimeric transcripts, comparison to reference databases, and ab initio prediction techniques, and the subgenomes were phased using ancestral k-mer spectral analysis. Subgenome B appears to show fission of chromosomal contigs, and while there is substantial synteny between subgenomes, we also identified regions lacking synteny that may have diverged in the ancestral genomes prior to or following hybridisation. This annotated and phased genome assembly forms a significant resource for understanding the origins and genetics of these globally important plant pathogens.

Incorporating eco-evolutionary information into species distribution models provides comprehensive predictions of species range shifts under climate change

As climate changes increasingly influence species distributions, ecosystem functions, and biodiversity, the urgency to understand how species' ranges shift under those changes is great. Species distribution models (SDMs) are vital approaches that can predict species distributions under changing climates. However, SDMs based on the species' current occurrences may underestimate the species' climatic tolerances. Integrating species' realized niches at different periods, also known as multi-temporal calibration, can provide an estimation closer to its fundamental niche. Based on this, we further proposed an integrated framework that combines eco-evolutionary data and SDMs (phylogenetically-informed SDMs) to provide comprehensive predictions of species range shifts under climate change. To evaluate our approach's performance, we applied it to a group of related species, the Chrysanthemum zawadskii species complex (Anthemidae, Asteracee). First, we investigated the niche differentiation between species and intraspecific lineages of the complex and estimated their rates of niche evolution. Next, using both standard SDMs and our phylogenetically-informed SDMs, we generated predictions of suitability areas for all species and lineages and compared the results. Finally, we reconstructed the historical range dynamics for the species of this complex. Our results showed that the species and intraspecific lineages of the complex had varying degrees of niche differentiation and different rates of niche evolution. Lineage-level SDMs can provide more realistic predictions for species with intraspecific differentiation than species-level models can. The phylogenetically-informed SDMs provided more complete environmental envelopes and predicted broader potential distributions for all species than the standard SDMs did. Range dynamics varied among the species that have different rates of niche evolution. Our framework integrating eco-evolutionary data and SDMs contributes to a better understanding of the species' responses to climate change and can help to make more targeted conservation efforts for the target species under climate change, particularly for rare species.

Haplotype-based phylogenetic analysis and population genomics uncover the origin and domestication of sweetpotato

The hexaploid sweetpotato (Ipomoea batatas) is one of the most important root crops worldwide. However, its genetic origin remains controversial, and its domestication history remains unknown. In this study, we used a range of genetic evidence and a newly developed haplotype-based phylogenetic analysis to identify two probable progenitors of sweetpotato. The diploid progenitor was likely closely related to Ipomoea aequatoriensis and contributed the B1 subgenome, IbT-DNA2, and the lineage 1 type of chloroplast genome to sweetpotato. The tetraploid progenitor of sweetpotato was most likely I. batatas 4x, which donated the B2 subgenome, IbT-DNA1, and the lineage 2 type of chloroplast genome. Sweetpotato most likely originated from reciprocal crosses between the diploid and tetraploid progenitors, followed by a subsequent whole-genome duplication. In addition, we detected biased gene exchanges between the subgenomes; the rate of B1 to B2 subgenome conversions was nearly three times higher than that of B2 to B1 subgenome conversions. Our analyses revealed that genes involved in storage root formation, maintenance of genome stability, biotic resistance, sugar transport, and potassium uptake were selected during the speciation and domestication of sweetpotato. This study sheds light on the evolution of sweetpotato and paves the way for improvement of this crop.

Separating phases of allopolyploid evolution with resynthesized and natural Capsella bursa-pastoris

Allopolyploidization is a frequent evolutionary transition in plants that combines whole-genome duplication (WGD) and interspecific hybridization. The genome of an allopolyploid species results from initial interactions between parental genomes and long-term evolution. Distinguishing the contributions of these two phases is essential to understanding the evolutionary trajectory of allopolyploid species. Here, we compared phenotypic and transcriptomic changes in natural and resynthesized Capsella allotetraploids with their diploid parental species. We focused on phenotypic traits associated with the selfing syndrome and on transcription-level phenomena such as expression-level dominance (ELD), transgressive expression (TRE), and homoeolog expression bias (HEB). We found that selfing syndrome, high pollen, and seed quality in natural allotetraploids likely resulted from long-term evolution. Similarly, TRE and most down-regulated ELD were only found in natural allopolyploids. Natural allotetraploids also had more ELD toward the self-fertilizing parental species than resynthesized allotetraploids, mirroring the establishment of the selfing syndrome. However, short-term changes mattered, and 40% of the cases of ELD in natural allotetraploids were already observed in resynthesized allotetraploids. Resynthesized allotetraploids showed striking variation of HEB among chromosomes and individuals. Homoeologous synapsis was its primary source and may still be a source of genetic variation in natural allotetraploids. In conclusion, both short- and long-term mechanisms contributed to transcriptomic and phenotypic changes in natural allotetraploids. However, the initial gene expression changes were largely reshaped during long-term evolution leading to further morphological changes.

Induction of tetraploids in Paper Mulberry (Broussonetia papyrifera (L.) L'Hér. ex Vent.) by colchicine

BACKGROUND

Broussonetia papyrifera (L.) L'Hér. ex Vent. has the characteristics of strong stress resistance, high crude protein content, and pruning tolerance. It is an ecological, economic, and medicinal plant. Polyploid plants usually perform better than their corresponding diploid plants in terms of nutrients, active substances, and stress resistance.

RESULTS

In this study, the leaves, calli, and seeds of diploid B. papyrifera were used for tetraploid induction by colchicine. The induction effect of colchicine on B. papyrifera was summarized through the early morphology, chromosome count and flow cytometry. It was concluded that the best induction effect (18.6%) was obtained when the leaves of B. papyrifera were treated in liquid MS (Murashige and Skoog) medium containing 450 mg·L-1 colchicine for 3 d. The comparative analysis of the growth characteristics of diploid and tetraploid B. papyrifera showed that tetraploid B. papyrifera has larger ground diameter, larger stomata, thicker palisade tissue and thicker sponge tissue than diploid B. papyrifera. In addition, the measurement of photosynthetic features also showed that tetraploids had higher chlorophyll content and higher photosynthetic rates.

CONCLUSION

This study showed that tetraploid B. papyrifera could be obtained by treating leaves, callus and seeds with liquid and solid colchicine, but the induction efficiency was different. Moreover, there were differences in stomata, leaf cell structure and photosynthetic features between tetraploid B. papyrifera and its corresponding diploid. The induced tetraploid B. papyrifera can provide a technical basis and breeding material for the creation of B. papyrifera germplasm resources in the future.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Artificial Diploid Escherichia coli by a CRISPR Chromosome-Doubling Technique

Synthetic biology has been represented by the creation of artificial life forms at the genomic scale. In this work, a CRISPR-based chromosome-doubling technique is designed to first construct an artificial diploid Escherichia coli cell. The stable single-cell diploid E. coli is isolated by both maximal dilution plating and flow cytometry, and confirmed with quantitative PCR, fluorescent in situ hybridization, and third-generation genome sequencing. The diploid E. coli has a greatly reduced growth rate and elongated cells at 4-5 µm. It is robust against radiation, and the survival rate after exposure to UV increased 40-fold relative to WT. As a novel life form, the artificial diploid E. coli is an ideal substrate for research fundamental questions in life science concerning polyploidy. And this technique may be applied to other bacteria.

Phylogeography and population structure of the global, wide host-range hybrid pathogen Phytophthora × cambivora

Invasive, exotic plant pathogens pose a major threat to native and agricultural ecosystems. Phytophthora × cambivora is an invasive, destructive pathogen of forest and fruit trees causing severe damage worldwide to chestnuts (Castanea), apricots, peaches, plums, almonds and cherries (Prunus), apples (Malus), oaks (Quercus), and beech (Fagus). It was one of the first damaging invasive Phytophthora species to be introduced to Europe and North America, although its origin is unknown. We determined its population genetic history in Europe, North and South America, Australia and East Asia (mainly Japan) using genotyping-by-sequencing. Populations in Europe and Australia appear clonal, those in North America are highly clonal yet show some degree of sexual reproduction, and those in East Asia are partially sexual. Two clonal lineages, each of opposite mating type, and a hybrid lineage derived from these two lineages, dominated the populations in Europe and were predominantly found on fagaceous forest hosts (Castanea, Quercus, Fagus). Isolates from fruit trees (Prunus and Malus) belonged to a separate lineage found in Australia, North America, Europe and East Asia, indicating the disease on fruit trees could be caused by a distinct lineage of P. × cambivora, which may potentially be a separate sister species and has likely been moved with live plants. The highest genetic diversity was found in Japan, suggesting that East Asia is the centre of origin of the pathogen. Further surveys in unsampled, temperate regions of East Asia are needed to more precisely identify the location and range of the centre of diversity.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Whole-genome duplication (WGD) occurs broadly and repeatedly across the history of eukaryotes and is recognized as a prominent evolutionary force, especially in plants. Immediately following WGD, most genes are present in two copies as paralogs. Due to this redundancy, one copy of a paralog pair commonly undergoes pseudogenization and is eventually lost. When speciation occurs shortly after WGD; however, differential loss of paralogs may lead to spurious phylogenetic inference resulting from the inclusion of pseudoorthologs–paralogous genes mistakenly identified as orthologs because they are present in single copies within each sampled species. The influence and impact of including pseudoorthologs versus true orthologs as a result of gene extinction (or incomplete laboratory sampling) are only recently gaining empirical attention in the phylogenomics community. Moreover, few studies have yet to investigate this phenomenon in an explicit coalescent framework. Here, using mathematical models, numerous simulated data sets, and two newly assembled empirical data sets, we assess the effect of pseudoorthologs on species tree estimation under varying degrees of incomplete lineage sorting (ILS) and differential gene loss scenarios following WGD. When gene loss occurs along the terminal branches of the species tree, alignment-based (BPP) and gene-tree-based (ASTRAL, MP-EST, and STAR) coalescent methods are adversely affected as the degree of ILS increases. This can be greatly improved by sampling a sufficiently large number of genes. Under the same circumstances, however, concatenation methods consistently estimate incorrect species trees as the number of genes increases. Additionally, pseudoorthologs can greatly mislead species tree inference when gene loss occurs along the internal branches of the species tree. Here, both coalescent and concatenation methods yield inconsistent results. These results underscore the importance of understanding the influence of pseudoorthologs in the phylogenomics era. [Coalescent method; concatenation method; incomplete lineage sorting; pseudoorthologs; single-copy gene; whole-genome duplication.]

Rapid and visual monitoring of alien sequences using crop wild relatives specific oligo-painting: The case of cucumber chromosome engineering

Wild species related to domesticated crops (crop wild relatives, or CWRs) represent a high level of genetic diversity that provides a practical gene pool for crop pre-breeding employed to address climate change and food demand challenges globally. Nevertheless, rapid identifying and visual tracking of alien chromosomes and sequences derived from CWRs have been a technical challenge for crop chromosome engineering. Here, a species-specific oligonucleotide (oligo) pool was developed by using the reference genome of Cucumis hystrix (HH, 2n = 2x = 24), a wild species carrying many favorable traits and interspecific compatibility with cultivated cucumber (C. sativus, CC, 2n = 2x = 14). These synthetic double-stranded oligo probes were applied to validate the assembly and characterize the chromosome architectures of C. hystrix, as well as to rapidly identify C. hystrix-chromosomes in diverse C. sativus-hystrix chromosome-engineered germplasms, including interspecific hybrid F1 (HC), synthetic allopolyploids (HHCC, CHC, and HCH) and alien additional lines (CC-H). Moreover, a ∼2Mb of C. hystrix-specific sequences, introduced into cultivated cucumber, were visualized by CWR-specific oligo-painting. These results demonstrate that the CWR-specific oligo-painting technique holds broad applicability for chromosome engineering of numerous crops, as it allows rapid identification of alien chromosomes, reliable detection of homoeologous recombination, and visual tracking of the introgression process. It is promising to achieve directed and high-precision crop pre-breeding combined with other breeding techniques, such as CRISPR/Cas9-mediated chromosome engineering.

Phylogenomic Analyses of Alismatales Shed Light into Adaptations to Aquatic Environments

Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptations to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptations, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the three main clades that was best explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light- or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and WGDs using phylogenomics.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

The Different Faces of Arabidopsis arenosa-A Plant Species for a Special Purpose

The following review article collects information on the plant species Arabidopsis arenosa. Thus far, A. arenosa has been known as a model species for autotetraploidy studies because, apart from diploid individuals, there are also tetraploid populations, which is a unique feature of this Arabidopsis species. In addition, A arenosa has often been reported in heavy metal-contaminated sites, where it occurs together with a closely related species A. halleri, a model plant hyperaccumulator of Cd and Zn. Recent studies have shown that several populations of A. arenosa also exhibit Cd and Zn hyperaccumulation. However, it is assumed that the mechanism of hyperaccumulation differs between these two Arabidopsis species. Nevertheless, this phenomenon is still not fully understood, and thorough research is needed. In this paper, we summarize the current state of knowledge regarding research on A. arenosa.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

Genomic insight into the divergence and adaptive potential of a forgotten landrace G. hirsutum L. purpurascens

Wild progenitors are an excellent source for strengthening the genetic basis and accumulation of desirable variation lost because of directional selection and adaptation in modern cultivars. Here, we re-evaluate a landrace of Gossypium hirsutum, formerly known as Gossypium purpurascens. Our study seeks to understand the genomic structure, variation, and breeding potential of this landrace, providing potential insights into the biogeographic history and genomic changes likely associated with domestication. A core set of accessions, including current varieties, obsolete accessions, G. purpurascens, and other geographical landraces, are subjected to genotyping along with multilocation phenotyping. Population fixation statistics suggests a marked differentiation between G. purpurascens and three other groups, emphasizing the divergent genomic behavior of G. purpurascens. Phylogenetic analysis establishes the primitive nature of G. purpurascens, identifying it as a vital source of functional variation, the inclusion of which in the upland cotton (cultivated G. hirsutum) gene pool may broaden the genetic basis of modern cultivars. Genome-wide association results indicate multiple loci associated with domestication regions corresponding to flowering and fiber quality. Moreover, the conserved nature of G. purpurascens can also provide insights into the evolutionary process of G. hirsutum.

Maternal karyogene and cytoplasmic genotype affect the induction efficiency of doubled haploid inducer in Brassica napus

BACKGROUND

Artificial synthesis of octoploid rapeseed double haploid (DH) induction lines Y3380 and Y3560 was made possible by interspecific hybridization and genome doubling techniques. Production of pure lines by DH induction provides a new way to achieve homozygosity earlier in B.napus. Previously, the mechanism of induction, and whether the induction has obvious maternal genotypic differences or not, are not known so far.

RESULTS

In this study, different karyogene and cytoplasmic genotype of B.napus were pollinated with the previously reported DH inducers e.g. Y3380 and Y3560. Our study presents a fine comparison of different cytoplasmic genotypes hybridization to unravel the mechanism of DH induction. Ploidy identification, fertility and SSR marker analysis of induced F₁ generation, revealed that ploidy and phenotype of the induced F₁ plants were consistent with that type of maternal, rather than paternal parent. The SNP chip analysis revealed that induction efficiency of DH inducers were affected by the karyogene when the maternal cytoplasmic genotypes were the same. However, DH induction efficiency was also affected by cytoplasmic genotype when the karyogenes were same, and the offspring of the ogura cytoplasm showed high frequency inducer gene hybridization or low-frequency infiltration.

CONCLUSION

The induction effect is influenced by the interaction between maternal karyogene and cytoplasmic genotype, and the results from the partial hybridization of progeny chromosomes indicate that the induction process may be attributed to the selective elimination of paternal chromosome. This study provides a basis for exploring the mechanism of DH inducer in B.napus, and provides new insights for utilization of inducers in molecular breeding.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Different strategies of Cd tolerance and accumulation in Arabidopsis halleri and Arabidopsis arenosa

Pseudometallophytes are commonly used to study the evolution of metal tolerance and accumulation traits in plants. Within the Arabidopsis genus, the adaptation of Arabidopsis halleri to metalliferous soils has been widely studied, which is not the case for the closely related species Arabidopsis arenosa. We performed an in-depth physiological comparison between the A. halleri and A. arenosa populations from the same polluted site, together with the geographically close non-metallicolous (NM) populations of both species. The ionomes, growth, photosynthetic parameters and pigment content were characterized in the plants that were growing on their native site and in a hydroponic culture under Cd treatments. In situ, the metallicolous (M) populations of both species hyperaccumulated Cd and Zn. The NM population of A. halleri hyperaccumulated Cd and Zn while the NM A. arenosa did not. In the hydroponic experiments, the NM populations of both species accumulated more Cd in their shoots than the M populations. Our research suggests that the two Arabidopsis species evolved different strategies of adaptation to extreme metallic environments that involve fine regulation of metal homeostasis, adjustment of the photosynthetic apparatus and accumulation of flavonols and anthocyanins.

Two cytotype niche shifts are of different magnitude in Solidago gigantea

PREMISE

Polyploidy may serve to contribute to range size if autopolyploid cytotypes are adapted to differing ecological conditions. This study aims to establish the geographic distribution of cytotypes within the giant goldenrod (Solidago gigantea), and to assess whether cytotypes exhibit differing ecological tolerances and morphology.

METHODS

A range-wide set of 629 Solidago gigantea individuals was obtained through field collecting, sampling from herbarium specimens, and incorporating existing chromosome counts. Cytotype of each unknown sample was estimated by observing allele numbers at twelve microsatellite loci, a strategy that was assessed by comparing estimated to known cytotype in 20 chromosome-counted samples. Abiotic ecological differentiation was assessed for two transitions: diploid-tetraploid and tetraploid-hexaploid. Morphological differentiation among cytotypes was assessed.

RESULTS

Microsatellite repeat variation accurately estimated cytotype in 85% of samples for which ploidy was known. Applying this approach to samples of unknown ploidy established that the three cytotypes are non-randomly distributed. Although niche modeling and MANOVA approaches identified significant differences in macro-climatic conditions for both cytotype transitions, the tetraploid to hexaploid transition was more substantial. Leaf length and width did not differ among cytotypes. Although leaf vestiture exhibited strong trends, no absolute differences were observed among cytotypes.

CONCLUSIONS

With the largest such study to date, we established niche transitions among giant goldenrod cytotypes of differing magnitudes. Collectively, this suggests that whole-genome duplication has contributed to Solidago gigantea's large range.

Comprehensive multiomics analysis reveals key roles of NACs in plant growth and development and its environmental adaption mechanism by regulating metabolite pathways

Abnormal environmental conditions induce polyploidization and exacerbate vulnerability to agricultural production. Polyploidization is a pivotal event for plant adaption to stress and the expansion of transcription factors. NACs play key roles in plant stress resistance and growth and development, but the adaptive mechanism of NACs during plant polyploidization remain to be explored. Here, we identified and analyzed NACs from 15 species and found that the expansion of NACs was contributed by polyploidization. The regulatory networks were systematically analyzed based on polyomics. NACs might influence plant phenotypes and were correlated with amino acids acting as nitrogen source, indicating that NACs play a vital role in plant development. More importantly, in quinoa and Arabidopsis thaliana, NACs enabled plants to resist stress by regulating flavonoid pathways, and the universality was further confirmed by the Arabidopsis population. Our study provides a cornerstone for future research into improvement of important agronomic traits by transcription factors in a changing global environment.

Parental legacy and regulatory novelty in Brachypodium diurnal transcriptomes accompanying their polyploidy

Polyploidy is a widespread phenomenon in eukaryotes that can lead to phenotypic novelty and has important implications for evolution and diversification. The modification of phenotypes in polyploids relative to their diploid progenitors may be associated with altered gene expression. However, it is largely unknown how interactions between duplicated genes affect their diurnal expression in allopolyploid species. In this study, we explored parental legacy and hybrid novelty in the transcriptomes of an allopolyploid species and its diploid progenitors. We compared the diurnal transcriptomes of representative Brachypodium cytotypes, including the allotetraploid Brachypodium hybridum and its diploid progenitors Brachypodium distachyon and Brachypodium stacei. We also artificially induced an autotetraploid B. distachyon. We identified patterns of homoeolog expression bias (HEB) across Brachypodium cytotypes and time-dependent gain and loss of HEB in B. hybridum. Furthermore, we established that many genes with diurnal expression experienced HEB, while their expression patterns and peak times were correlated between homoeologs in B. hybridum relative to B. distachyon and B. stacei, suggesting diurnal synchronization of homoeolog expression in B. hybridum. Our findings provide insight into the parental legacy and hybrid novelty associated with polyploidy in Brachypodium, and highlight the evolutionary consequences of diurnal transcriptional regulation that accompanied allopolyploidy.

Comparative Genomic Analyses and a Novel Linkage Map for Cisco (Coregonus artedi) Provide Insights into Chromosomal Evolution and Rediploidization Across Salmonids

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco (Coregonus artedi), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts.

Molecular footprints of selection effects and whole genome duplication (WGD) events in three blueberry species: detected by transcriptome dataset

BACKGROUND

Both selection effects and whole genome duplication played very important roles in plant speciation and evolution, and to decipher the corresponding molecular footprint has always been a central task of geneticists. Vaccinium is species rich genus that comprised of about 450 species, and blueberry is one of the most important species of Vaccinium genus, which is gaining popularity because of high healthful value. In this article, we aimed to decipher the molecular footprints of natural selection on the single copy genes and WGD events occur in the evolutionary history of blueberry species.

RESULTS

We identified 30,143, 29,922 and 28,891 putative protein coding sequences from 45,535, 42,914 and 43,630 unigenes assembled from the leaves' transcriptome assembly of 19 rabbiteye (T1), 13 southern highbush (T2) and 22 northern highbush (T3) blueberry cultivars. A total of 17, 21 and 27 single copy orthologs were found to undergone positive selection in T1 versus T2, T1 versus T3, and T2 versus T3, respectively, and these orthologs were enriched in metabolic pathways including "Terpenoid backbone biosynthesis", "Valine, leucine and isoleucine biosynthesis", "Butanoate metabolism", "C5-Branched dibasic acid metabolism" "Pantothenate and CoA biosynthesis". We also detected significant molecular footprints of a recent (about 9.04 MYA), medium (about 43.44 MYA) and an ancient (about 116.39 MYA) WGD events that occurred in the evolutionary history of three blueberry species.

CONCLUSION

Some important functional genes revealed positive selection effect in blueberry. At least three rounds of WGD events were detected in the evolutionary history of blueberry species. Our work provides insights about the genetic mechanism of adaptive evolution in blueberry and species radiation of Vaccinium in short geological scale time.

The hybrid protein interactome contributes to rice heterosis as epistatic effects

Heterosis is the phenomenon in which hybrid progeny exhibits superior traits in comparison with those of their parents. Genomic variations between the two parental genomes may generate epistasis interactions, which is one of the genetic hypotheses explaining heterosis. We postulate that protein-protein interactions specific to F₁ hybrids (F₁ -specific PPIs) may occur when two parental genomes combine, as the proteome of each parent may supply novel interacting partners. To test our assumption, an inter-subspecies hybrid interactome was simulated by in silico PPI prediction between rice japonica (cultivar Nipponbare) and indica (cultivar 9311). Four-thousand, six-hundred and twelve F₁ -specific PPIs accounting for 20.5% of total PPIs in the hybrid interactome were found. Genes participating in F₁ -specific PPIs tend to encode metabolic enzymes and are generally localized in genomic regions harboring metabolic gene clusters. To test the genetic effect of F₁ -specific PPIs in heterosis, genomic selection analysis was performed for trait prediction with additive, dominant and epistatic effects separately considered in the model. We found that the removal of single nucleotide polymorphisms associated with F₁ -specific PPIs reduced prediction accuracy when epistatic effects were considered in the model, but no significant changes were observed when additive or dominant effects were considered. In summary, genomic divergence widely dispersed between japonica and indica rice may generate F₁ -specific PPIs, part of which may accumulatively contribute to heterosis according to our computational analysis. These candidate F₁ -specific PPIs, especially for those involved in metabolic biosynthesis pathways, are worthy of experimental validation when large-scale protein interactome datasets are generated in hybrid rice in the future.

Is selfing a reproductive assurance promoting polyploid establishment? Reduced fitness, leaky self-incompatibility and lower inbreeding depression in neotetraploids

PREMISE

Newly formed polyploids face significant obstacles to persistence and population establishment because of fitness costs of intercytotype mating. Selfing provides the opportunity to escape mate limitation, enabling production of new individuals and increasing the likelihood of fixation of new polyploid lineages. Still, association between self-compatibility and polyploidy is not always clear. We compared self-incompatibility and inbreeding depression in neotetraploids and their diploid progenitor to explore the direct effects of whole genome duplications on self-incompatibility and the implications of ploidy-driven changes for polyploid establishment.

METHODS

Outcross and self-pollinations were performed in diploids and synthetic neotetraploids of Jasione maritima var. maritima, and reproductive success was measured through fruit and seed production and seed germination. Self- and outcross offspring were grown under controlled conditions, and plant performance was measured through several fitness parameters.

RESULTS

Neotetraploids showed an overall lower performance than diploids. Reproductive success was negatively affected by selfing in both cytotypes. However, greater variation in the expression of self-incompatibility was observed in neotetraploids; additionally, developmental and physiological parameters were not affected by selfing on neotetraploids, with an overall similar fitness of outcrossed and selfed individuals, resulting in lower inbreeding depression indexes.

CONCLUSIONS

Neotetraploids might have benefited from selfing at initial stages after their formation. Genome duplications resulted in leaky self-incompatibility, enabling the production of offspring under minority cytotype disadvantage with similar fitness as outcrossed offspring. Our results support theoretical assumptions that selfing might be important for neopolyploid establishment, although changes in self-incompatibility might not be abrupt.

Unraveling the Hexaploid Sweetpotato Inheritance Using Ultra-Dense Multilocus Mapping

The hexaploid sweetpotato (Ipomoea batatas (L.) Lam., 2n = 6x = 90) is an important staple food crop worldwide and plays a vital role in alleviating famine in developing countries. Due to its high ploidy level, genetic studies in sweetpotato lag behind major diploid crops significantly. We built an ultra-dense multilocus integrated genetic map and characterized the inheritance system in a sweetpotato full-sib family using our newly developed software, MAPpoly. The resulting genetic map revealed 96.5% collinearity between I. batatas and its diploid relative I. trifida We computed the genotypic probabilities across the whole genome for all individuals in the mapping population and inferred their complete hexaploid haplotypes. We provide evidence that most of the meiotic configurations (73.3%) were resolved in bivalents, although a small portion of multivalent signatures (15.7%), among other inconclusive configurations (11.0%), were also observed. Except for low levels of preferential pairing in linkage group 2, we observed a hexasomic inheritance mechanism in all linkage groups. We propose that the hexasomic-bivalent inheritance promotes stability to the allelic transmission in sweetpotato.

Unraveling the Hexaploid Sweetpotato Inheritance Using Ultra-Dense Multilocus Mapping

The hexaploid sweetpotato (Ipomoea batatas (L.) Lam., 2n = 6x = 90) is an important staple food crop worldwide and plays a vital role in alleviating famine in developing countries. Due to its high ploidy level, genetic studies in sweetpotato lag behind major diploid crops significantly. We built an ultra-dense multilocus integrated genetic map and characterized the inheritance system in a sweetpotato full-sib family using our newly developed software, MAPpoly. The resulting genetic map revealed 96.5% collinearity between I. batatas and its diploid relative I. trifida We computed the genotypic probabilities across the whole genome for all individuals in the mapping population and inferred their complete hexaploid haplotypes. We provide evidence that most of the meiotic configurations (73.3%) were resolved in bivalents, although a small portion of multivalent signatures (15.7%), among other inconclusive configurations (11.0%), were also observed. Except for low levels of preferential pairing in linkage group 2, we observed a hexasomic inheritance mechanism in all linkage groups. We propose that the hexasomic-bivalent inheritance promotes stability to the allelic transmission in sweetpotato.

Oligo-painting and GISH reveal meiotic chromosome biases and increased meiotic stability in synthetic allotetraploid Cucumis ×hytivus with dysploid parental karyotypes

BACKGROUND

Meiosis of newly formed allopolyploids frequently encounter perturbations induced by the merging of divergent and hybridizable genomes. However, to date, the meiotic properties of allopolyploids with dysploid parental karyotypes have not been studied in detail. The allotetraploid Cucumis ×hytivus (HHCC, 2n = 38) was obtained from interspecific hybridization between C. sativus (CC, 2n = 14) and C. hystrix (HH, 2n = 24) followed by chromosome doubling. The results of this study thus offer an excellent opportunity to explore the meiotic properties of allopolyploids with dysploid parental karyotypes.

RESULTS

In this report, we describe the meiotic properties of five chromosomes (C5, C7, H1, H9 and H10) and two genomes in interspecific hybrids and C. ×hytivus (the 4th and 14th inbred family) through oligo-painting and genomic in situ hybridization (GISH). We show that 1) only two translocations carrying C5-oligo signals were detected on the chromosomes C2 and C4 of one 14th individual by the karyotyping of eight 4th and 36 14th plants based on C5- and C7-oligo painting, and possible cytological evidence was observed in meiosis of the 4th generation; 2) individual chromosome have biases for homoeologous pairing and univalent formation in F₁ hybrids and allotetraploids; 3) extensive H-chromosome autosyndetic pairings (e.g., H-H, 25.5% PMCs) were observed in interspecific F₁ hybrid, whereas no C-chromosome autosyndetic pairings were observed (e.g. C-C); 4) the meiotic properties of two subgenomes have significant biases in allotetraploids: H-subgenome exhibits higher univalent and chromosome lagging frequencies than C-subgenome; and 5) increased meiotic stability in the S₁₄ generation compared with the S₄ generation, including synchronous meiosis behavior, reduced incidents of univalent and chromosome lagging.

CONCLUSIONS

These results suggest that the meiotic behavior of two subgenomes has dramatic biases in response to interspecific hybridization and allopolyploidization, and the meiotic behavior harmony of subgenomes is a key subject of meiosis evolution in C. ×hytivus. This study helps to elucidate the meiotic properties and evolution of nascent allopolyploids with the dysploid parental karyotypes.

DNA methylation repatterning accompanying hybridization, whole genome doubling and homoeolog exchange in nascent segmental rice allotetraploids

Allopolyploidization, which entails interspecific hybridization and whole genome duplication (WGD), is associated with emergent genetic and epigenetic instabilities that are thought to contribute to adaptation and evolution. One frequent genomic consequence of nascent allopolyploidization is homoeologous exchange (HE), which arises from compromised meiotic fidelity and generates genetically and phenotypically variable progenies. Here, we used a genetically tractable synthetic rice segmental allotetraploid system to interrogate genome-wide DNA methylation and gene expression responses and outcomes to the separate and combined effects of hybridization, WGD and HEs. Progenies of the tetraploid rice were genomically diverse due to genome-wide HEs that affected all chromosomes, yet they exhibited overall methylome stability. Nonetheless, regional variation of cytosine methylation states was widespread in the tetraploids. Transcriptome profiling revealed genome-wide alteration of gene expression, which at least in part associates with changes in DNA methylation. Intriguingly, changes of DNA methylation and gene expression could be decoupled from hybridity and sustained and amplified by HEs. Our results suggest that HEs, a prominent genetic consequence of nascent allopolyploidy, can exacerbate, diversify and perpetuate the effects of allopolyploidization on epigenetic and gene expression variation, and hence may contribute to allopolyploid evolution.

Next-generation transcriptome assembly and analysis: Impact of ploidy

Whole genome duplications (WGD) occur widely in plants, but the effects of these events impact all branches of life. WGD events have major evolutionary impacts, often leading to major structural changes within the chromosomes and massive changes in gene expression that facilitate rapid speciation and gene diversification. Even for species that currently have diploid genomes, the impact of ancestral duplication events is still present in the genomes, especially in the context of highly similar gene families that are retained from WGD. However, the impact of these ploidies on various bioinformatics workflows has not been studied well. In this review, we overview biological significance of polyploidy in different organisms. We describe the impact of having polyploid transcriptomes on bioinformatics analyses, especially focusing on transcriptome assembly and transcript quantification. We discuss the benefits of using simulated benchmarking data when we examine the performance of various methods. We also present an example strategy to generate simulated allopolyploid transcriptomes and RNAseq datasets and how these benchmark datasets can be used to assess the performance of transcript assembly and quantification methods. Our benchmarking study shows that all transcriptome assembly methods are affected by having polyploid genomes. Quantification accuracy is also impacted by polyploidy depending on the method. These simulated datasets can be adapted for testing, such as, read mapping, variant calling, and differential expression using biologically realistic conditions.

Long-Term Conservation of Ohnologs Through Partial Tetrasomy Following Whole-Genome Duplication in Salmonidae

Whole-genome duplications (WGDs) have occurred repeatedly and broadly throughout the evolutionary history of eukaryotes. However, the effects of WGD on genome function and evolution remain unclear. The salmonid WGD that occurred approximately 88 million years ago presents an excellent opportunity for studying the effects of WGD as ∼10-15% of each salmonid genome still exhibits tetrasomic inheritance. Herein, we utilized the rainbow trout (Oncorhynchus mykiss) genome assembly and brain transcriptome data to examine the fate of gene pairs (ohnologs) following the salmonid whole-genome duplication. We find higher sequence identity between ohnologs located within known tetrasomic regions than between ohnologs found in disomic regions, and that tetrasomically inherited ohnologs showed greater similarity in patterns of gene expression and per ohnolog were lower expressed, than disomically inherited ohnologs. Enrichment testing for Gene Ontology terms identified 49 over-represented terms in tetrasomically inherited ohnologs compared to disomic ohnologs. However, why these ohnologs are retained as tetrasomic is difficult to answer. It could be that we have identified salmonid specific "dangerous duplicates", that is, genes that cannot take on new roles following WGD. Alternatively, there may be adaptive advantages for retaining genes as functional duplicates in tetrasomic regions, as presumably, movement of these genes into disomic regions would affect both their sequence identity and their gene expression patterns.

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.

Transgenerationally Precipitated Meiotic Chromosome Instability Fuels Rapid Karyotypic Evolution and Phenotypic Diversity in an Artificially Constructed Allotetraploid Wheat (AADD)

Although a distinct karyotype with defined chromosome number and structure characterizes each biological species, it is intrinsically labile. Polyploidy or whole-genome duplication has played a pervasive and ongoing role in the evolution of all eukaryotes, and is the most dramatic force known to cause rapid karyotypic reconfiguration, especially at the initial stage. However, issues concerning transgenerational propagation of karyotypic heterogeneity and its translation to phenotypic diversity in nascent allopolyploidy, at the population level, have yet to be studied in detail. Here, we report a large-scale examination of transgenerationally propagated karyotypic heterogeneity and its phenotypic manifestation in an artificially constructed allotetraploid with a genome composition of AADD, that is, involving two of the three progenitor genomes of polyploid wheat. Specifically, we show that 1) massive organismal karyotypic heterogeneity is precipitated after 12 consecutive generations of selfing from a single euploid founder individual, 2) there exist dramatic differences in aptitudes between subgenomes and among chromosomes for whole-chromosome gain and/or loss and structural variations, 3) majority of the numerical and structural chromosomal variations are concurrent due to mutual contingency and possible functional constraint, 4) purposed and continuous selection and propagation for euploidy over generations did not result in enhanced karyotype stabilization, and 5) extent of karyotypic variation correlates with variability of phenotypic manifestation. Together, our results document that allopolyploidization catalyzes rampant and transgenerationally heritable organismal karyotypic heterogeneity that drives population-level phenotypic diversification, which lends fresh empirical support to the still contentious notion that whole-genome duplication enhances organismal evolvability.

Paleogenomics: reconstruction of plant evolutionary trajectories from modern and ancient DNA

How contemporary plant genomes originated and evolved is a fascinating question. One approach uses reference genomes from extant species to reconstruct the sequence and structure of their common ancestors over deep timescales. A second approach focuses on the direct identification of genomic changes at a shorter timescale by sequencing ancient DNA preserved in subfossil remains. Merged within the nascent field of paleogenomics, these complementary approaches provide insights into the evolutionary forces that shaped the organization and regulation of modern genomes and open novel perspectives in fostering genetic gain in breeding programs and establishing tools to predict future population changes in response to anthropogenic pressure and global warming.

Protocols for the Molecular Evolutionary Analysis of Membrane Protein Gene Duplicates

Gene duplication is an important process in the evolution of gene content in eukaryotic genomes. Understanding when gene duplicates contribute new molecular functions to genomes through molecular adaptation is one important goal in comparative genomics. In large gene families, however, characterizing adaptation and neofunctionalization across species is challenging, as models have traditionally quantified the timing of duplications without considering underlying gene trees. This protocol combines multiple approaches to detect adaptation in protein duplicates at a phylogenetic scale. We include a description of models for gene tree-species tree reconciliation that enable different types of inference, as well as a practical guide to their use. Although simulation-based approaches successfully detect shifts in the rate of duplication/retention, the conflation between the duplication and retention processes, the distinct trajectories of duplicates under non-, sub-, and neofunctionalization, as well as dosage effects offer hitherto unexplored analytical avenues. We introduce mathematical descriptions of these probabilities and offer a road map to computational implementation whose starting point is parsimony reconciliation. Sequence evolution information based on the ratio of nonsynonymous to synonymous nucleotide substitution rates (dN/dS) can be combined with duplicate survival probabilities to better predict the emergence of new molecular functions in retained duplicates. Together, these methods enable characterization of potentially adaptive candidate duplicates whose neofunctionalization may contribute to phenotypic divergence across species.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae)

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

Meiotic chromosome stability of a newly formed allohexaploid wheat is facilitated by selection under abiotic stress as a spandrel

Polyploidy is a prominent route to speciation in plants; however, this entails resolving the challenges of meiotic instability facing abrupt doubling of chromosome complement. This issue remains poorly understood. We subjected progenies of a synthetic hexaploid wheat, analogous to natural common wheat, but exhibiting extensive meiotic chromosome instability, to heat or salt stress. We selected stress-tolerant cohorts and generated their progenies under normal condition. We conducted fluorescent in situ hybridization/genomic in situ hybridization-based meiotic/mitotic analysis, RNA-Seq and virus-induced gene silencing (VIGS)-mediated assay of meiosis candidate genes. We show that heritability of stress tolerance concurred with increased euploidy frequency due to enhanced meiosis stability. We identified a set of candidate meiosis genes with altered expression in the stress-tolerant plants vs control, but the expression was similar to that of common wheat (cv Chinese Spring, CS). We demonstrate VIGS-mediated downregulation of individual candidate meiosis genes in CS is sufficient to confer an unstable meiosis phenotype mimicking the synthetic wheat. Our results suggest that heritable regulatory changes of preexisting meiosis genes may be hitchhiked as a spandrel of stress tolerance, which significantly improves meiosis stability in the synthetic wheat. Our findings implicate a plausible scenario that the meiosis machinery in hexaploid wheat may have already started to evolve at its onset stage.

Tandem Duplicate Genes in Maize Are Abundant and Date to Two Distinct Periods of Time

Tandem duplicate genes are proximally duplicated and as such occur in similar genomic neighborhoods. Using the maize B73 and PH207 de novo genome assemblies, we identified thousands of tandem gene duplicates that account for ∼10% of the annotated genes. These tandem duplicates have a bimodal distribution of ages, which coincide with ancient allopolyploidization and more recent domestication. Tandem duplicates are smaller on average and have a higher probability of containing LTR elements than other genes, suggesting origins in nonhomologous recombination. Within relatively recent tandem duplicate genes, ∼26% appear to be undergoing degeneration or divergence in function from the ancestral copy. Our results show that tandem duplicates are abundant in maize, arose in bursts throughout maize evolutionary history under multiple potential mechanisms, and may provide a substrate for novel phenotypic variation.

Immediate vs. evolutionary consequences of polyploidy on clonal reproduction in an autopolyploid plant

Background and Aims

Clonal reproduction in polyploids is expected to exceed that in diploids, due to either the immediate direct effects of whole-genome duplication (WGD) or selection during establishment. The timing of polyploidy effects on clonality are largely unknown despite its hypothesized influence on polyploid success. This study tests the direction and timing of divergence in clonal traits in diploid and polyploid Chamerion angustifolium.

Methods

Root bud production and biomass allocation patterns were compared between diploids and synthesized tetraploids (neotetraploids), and between neotetraploids and naturally occurring tetraploids grown in a common environment.

Key Results

Neotetraploids produced more root buds and fewer sexual structures than diploids and natural tetraploids; diploids and natural tetraploids had similar root bud numbers and sexual investment. The root bud:inflorescence biomass ratio was 71 % higher in neotetraploids than in natural tetraploids. Root bud location suggests that ramet density in neotetraploid genets could be higher than in diploid genets.

Conclusions

WGD immediately increases investment in asexual vs. sexual reproduction in C. angustifolium, potentially promoting within-cytotype mating and establishment for neopolyploids. However, evolutionary change after the polyploidization event negates the direct effects of WGD. Natural polyploids and diploids have similar root bud production and biomass allocation patterns, probably resulting from habitat- and ploidy-mediated selection on polyploids to become more like diploids. These results highlight the value of studying the effects of polyploidization in young vs. established polyploids.

Transient Activation of Apomixis in Sexual Neotriploids May Retain Genomically Altered States and Enhance Polyploid Establishment

Polyploid genomes evolve and follow a series of dynamic transfigurations along with adaptation and speciation. The initial formation of a new polyploid individual within a diploid population usually involves a triploid bridge, a two-step mechanism of cell fusions between ubiquitous (reduced) and rare (unreduced) gametes. The primary fusion event creates an intermediate triploid individual with unbalanced genome sets, a situation of genomic-shock characterized by gene expression dysregulation, high dosage sensitivity, disturbed cell divisions, and physiological and reproductive attributes drastically altered. This near-sterile neotriploid must produce (even) eupolyploids through secondary fusion events to restore genome steadiness, meiotic balance, and fertility required for the demographic establishment of a nascent lineage. Natural conditions locate several difficulties to polyploid establishment, including the production of highly unbalanced and rarely unreduced (euploid) gametes, frequency-dependent disadvantages (minority cytotype exclusion), severe fitness loss, and ecological competition with diploid parents. Persistence and adaptation of neopolyploids depend upon genetic and phenotypic novelty coupled to joint selective forces that preserve shock-induced genomic changes (subgenome homeolog partitioning) and drive meiotic (reproductive) stabilization and ecological diversification. Thus, polyploid establishment through the triploid bridge is a feasible but not ubiquitous process that requires a number of low-probability events and singular circumstances. Yet, frequencies of polyploids suggest that polyploid establishment is a pervasive process. To explain this disparity, and supported in experimental evidence, I propose that situations like hybridization and ploidy-state transitions associated to genomic shock and substantial developmental alterations can transiently activate apomixis as a mechanism to halt genomic instability and cancel factors restraining neopolyploid's sexual fertility, particularly in triploids. Apomixis -as a temporal alternative to sex- skip meiosis and syngamy, and thus can freeze genomic attributes, avoid unbalanced chromosomal segregation and increase the formation of unreduced euploid gametes, elude frequency-dependent reproductive disadvantages by parthenogenetic development of the embryo and permissive development of endosperm during seed formation, and increase the effective population size of the neopolyploid lineage favoring the formation rate of eupolyploids compared to aneuploids. The subsequent action of genome resilience mechanisms that alleviate transcriptomic shock and selection upon gene interactions might restore a stable meiosis and sexual fertility within few generations, as observed in synthetic polyploids. Alternatively, provided that resilience mechanisms fail, the neopolyploid might retain apomixis and hold genomically and transcriptionally altered states for many generations.