Unveiling kiwifruit TCP genes: evolution, functions, and expression insights

The TEOSINTE-BRANCHED1/CYCLOIDEA/PROLEFERATING-CELL-FACTORS (TCP) gene family is a plant-specific transcriptional factor family involved in leaf morphogenesis and senescence, lateral branching, hormone crosstalk, and stress responses. To date, a systematic study on the identification and characterization of the TCP gene family in kiwifruit has not been reported. Additionally, the function of kiwifruit TCPs in regulating kiwifruit responses to the ethylene treatment and bacterial canker disease pathogen (Pseudomonas syringae pv. actinidiae, Psa) has not been investigated. Here, we identified 40 and 26 TCP genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes, respectively. The synteny analysis of AcTCPs illustrated that whole-genome duplication accounted for the expansion of the TCP family in Ac. Phylogenetic, conserved domain, and selection pressure analysis indicated that TCP family genes in Ac and Ae had undergone different evolutionary patterns after whole-genome duplication (WGD) events, causing differences in TCP gene number and distribution. Our results also suggested that protein structure and cis-element architecture in promoter regions of TCP genes have driven the function divergence of duplicated gene pairs. Three and four AcTCP genes significantly affected kiwifruit responses to the ethylene treatment and Psa invasion, respectively. Our results provided insight into general characters, evolutionary patterns, and functional diversity of kiwifruit TCPs.

Evolutionary analysis of gene ages across TADs associates chromatin topology with whole-genome duplications

Topologically associated domains (TADs) are interaction subnetworks of chromosomal regions in 3D genomes. TAD boundaries frequently coincide with genome breaks while boundary deletion is under negative selection, suggesting that TADs may facilitate genome rearrangements and evolution. We show that genes co-localize by evolutionary age in humans and mice, resulting in TADs having different proportions of younger and older genes. We observe a major transition in the age co-localization patterns between the genes born during vertebrate whole-genome duplications (WGDs) or before and those born afterward. We also find that genes recently duplicated in primates and rodents are more frequently essential when they are located in old-enriched TADs and interact with genes that last duplicated during the WGD. Therefore, the evolutionary relevance of recent genes may increase when located in TADs with established regulatory networks. Our data suggest that TADs could play a role in organizing ancestral functions and evolutionary novelty.

Sustained ERK signaling promotes G2 cell cycle exit and primes cells for whole-genome duplication

Whole-genome duplication (WGD) is a frequent event in cancer evolution that fuels chromosomal instability. WGD can result from mitotic errors or endoreduplication, yet the molecular mechanisms that drive WGD remain unclear. Here, we use live single-cell analysis to characterize cell-cycle dynamics upon aberrant Ras-ERK signaling. We find that sustained ERK signaling in human cells leads to reactivation of the APC/C in G2, resulting in tetraploid G0-like cells that are primed for WGD. This process is independent of DNA damage or p53 but dependent on p21. Transcriptomics analysis and live-cell imaging showed that constitutive ERK activity promotes p21 expression, which is necessary and sufficient to inhibit CDK activity and which prematurely activates the anaphase-promoting complex (APC/C). Finally, either loss of p53 or reduced ERK signaling allowed for endoreduplication, completing a WGD event. Thus, sustained ERK signaling-induced G2 cell cycle exit represents an alternative path to WGD.

Extensive Loss and Gain of Conserved Noncoding Elements During Early Teleost Evolution

Conserved noncoding elements in vertebrates are enriched around transcription factor loci associated with development. However, loss and rapid divergence of conserved noncoding elements has been reported in teleost fish, albeit taking only few genomes into consideration. Taking advantage of the recent increase in high-quality teleost genomes, we focus on studying the evolution of teleost conserved noncoding elements, carrying out targeted genomic alignments and comparisons within the teleost phylogeny to detect conserved noncoding elements and reconstruct the ancestral teleost conserved noncoding elements repertoire. This teleost-centric approach confirms previous observations of extensive vertebrate conserved noncoding elements loss early in teleost evolution, but also reveals massive conserved noncoding elements gain in the teleost stem-group over 300 million years ago. Using synteny-based association to link conserved noncoding elements to their putatively regulated target genes, we show the most teleost gained conserved noncoding elements are found in the vicinity of orthologous loci involved in transcriptional regulation and embryonic development that are also associated with conserved noncoding elements in other vertebrates. Moreover, teleost and vertebrate conserved noncoding elements share a highly similar motif and transcription factor binding site vocabulary. We suggest that early teleost conserved noncoding element gains reflect a restructuring of the ancestral conserved noncoding element repertoire through both extreme divergence and de novo emergence. Finally, we support newly identified pan-teleost conserved noncoding elements have potential for accurate resolution of teleost phylogenetic placements in par with coding sequences, unlike ancestral only elements shared with spotted gar. This work provides new insight into conserved noncoding element evolution with great value for follow-up work on phylogenomics, comparative genomics, and the study of gene regulation evolution in teleosts.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Neopolyploidy-induced changes in giant duckweed (Spirodela polyrhiza) alter herbivore preference and performance and plant population performance

PREMISE

Polyploidy is a widespread mutational process in angiosperms that may alter population performance of not only plants but also their interacting species. Yet, knowledge of whether polyploidy affects plant-herbivore dynamics is scarce. Here, we tested whether aphid herbivores exhibit preference for diploid or neopolyploid plants, whether polyploidy impacts plant and herbivore performance, and whether these interactions depend on the plant genetic background.

METHODS

Using independently synthesized neotetraploid strains paired with their diploid progenitors of greater duckweed (Spirodela polyrhiza), we evaluated the effect of neopolyploidy on duckweed's interaction with the water-lily aphid (Rhopalosiphum nymphaeae). Using paired-choice experiments, we evaluated feeding preference of the herbivore. We then evaluated the consequences of polyploidy on aphid and plant performance by measuring population growth over multiple generations.

RESULTS

Aphids preferred neopolyploids when plants were provided at equal abundances but not at equal surface areas, suggesting the role of plant population surface area in driving this preference. Additionally, neopolyploidy increased aphid population performance, but this result was dependent on the plant's genetic lineage. Lastly, the impact of herbivory on neopolyploid vs. diploid duckweed varied greatly with genetic lineage, where neopolyploids appeared to be variably tolerant compared to diploids, sometimes mirroring the effect on herbivore performance.

CONCLUSIONS

By experimentally testing the impacts of polyploidy on trophic species interactions, we showed that polyploidization can impact the preference and performance of herbivores on their plant hosts. These results have significant implications for the establishment and persistence of plants and herbivores in the face of plant polyploidy.

A chromosome-scale assembly for 'd'Anjou' pear

Cultivated pear consists of several Pyrus species with Pyrus communis (European pear) representing a large fraction of worldwide production. As a relatively recently domesticated crop and perennial tree, pear can benefit from genome-assisted breeding. Additionally, comparative genomics within Rosaceae promises greater understanding of evolution within this economically important family. Here, we generate a fully phased chromosome-scale genome assembly of P. communis 'd'Anjou.' Using PacBio HiFi and Dovetail Omni-C reads, the genome is resolved into the expected 17 chromosomes, with each haplotype totaling nearly 540 Megabases and a contig N50 of nearly 14 Mb. Both haplotypes are highly syntenic to each other and to the Malus domestica 'Honeycrisp' apple genome. Nearly 45,000 genes were annotated in each haplotype, over 90% of which have direct RNA-seq expression evidence. We detect signatures of the known whole-genome duplication shared between apple and pear, and we estimate 57% of d'Anjou genes are retained in duplicate derived from this event. This genome highlights the value of generating phased diploid assemblies for recovering the full allelic complement in highly heterozygous crop species.

Chromosome-level genome of the transformable northern wattle, Acacia crassicarpa

The genus Acacia is a large group of woody legumes containing an enormous amount of morphological diversity in leaf shape. This diversity is at least in part the result of an innovation in leaf development where many Acacia species are capable of developing leaves of both bifacial and unifacial morphologies. While not unique in the plant kingdom, unifaciality is most commonly associated with monocots, and its developmental genetic mechanisms have yet to be explored beyond this group. In this study, we identify an accession of Acacia crassicarpa with high regeneration rates and isolate a clone for genome sequencing. We generate a chromosome-level assembly of this readily transformable clone, and using comparative analyses, confirm a whole-genome duplication unique to Caesalpinoid legumes. This resource will be important for future work examining genome evolution in legumes and the unique developmental genetic mechanisms underlying unifacial morphogenesis in Acacia.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

Integrative Analysis of Oleosin Genes Provides Insights into Lineage-Specific Family Evolution in Brassicales

Oleosins (OLEs) are a class of small but abundant structural proteins that play essential roles in the formation and stabilization of lipid droplets (LDs) in seeds of oil crops. Despite the proposal of five oleosin clades (i.e., U, SL, SH, T, and M) in angiosperms, their evolution in eudicots has not been well-established. In this study, we employed Brassicales, an economically important order of flowering plants possessing the lineage-specific T clade, as an example to address this issue. Three to 10 members were identified from 10 species representing eight plant families, which include Caricaceae, Moringaceae, Akaniaceae, Capparaceae, and Cleomaceae. Evolutionary and reciprocal best hit-based homologous analyses assigned 98 oleosin genes into six clades (i.e., U, SL, SH, M, N, and T) and nine orthogroups (i.e., U1, U2, SL, SH1, SH2, SH3, M, N, and T). The newly identified N clade represents an ancient group that has already appeared in the basal angiosperm Amborella trichopoda, which are constitutively expressed in the tree fruit crop Carica papaya, including pulp and seeds of the fruit. Moreover, similar to Clade N, the previously defined M clade is actually not Lauraceae-specific but an ancient and widely distributed group that diverged before the radiation of angiosperm. Compared with A. trichopoda, lineage-specific expansion of the family in Brassicales was largely contributed by recent whole-genome duplications (WGDs) as well as the ancient γ event shared by all core eudicots. In contrast to the flower-preferential expression of Clade T, transcript profiling revealed an apparent seed/embryo/endosperm-predominant expression pattern of most oleosin genes in Arabidopsis thaliana and C. papaya. Moreover, the structure and expression divergence of paralogous pairs was frequently observed, and a good example is the lineage-specific gain of an intron. These findings provide insights into lineage-specific family evolution in Brassicales, which facilitates further functional studies in nonmodel plants such as C. papaya.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

PloiDB: the plant ploidy database

See also the Commentary on this article by Spoelhof et al., 240: 909–911.

Plant neopolyploidy and genetic background differentiate the microbiome of duckweed across a variety of natural freshwater sources

Whole-genome duplication has long been appreciated for its role in driving phenotypic novelty in plants, often altering the way organisms interface with the abiotic environment. Only recently, however, have we begun to investigate how polyploidy influences interactions of plants with other species, despite the biotic niche being predicted as one of the main determinants of polyploid establishment. Nevertheless, we lack information about how polyploidy affects the diversity and composition of the microbial taxa that colonize plants, and whether this is genotype-dependent and repeatable across natural environments. This information is a first step towards understanding whether the microbiome contributes to polyploid establishment. We, thus, tested the immediate effect of polyploidy on the diversity and composition of the bacterial microbiome of the aquatic plant Spirodela polyrhiza using four pairs of diploids and synthetic autotetraploids. Under controlled conditions, axenic plants were inoculated with pond waters collected from 10 field sites across a broad environmental gradient. Autotetraploids hosted 4%-11% greater bacterial taxonomic and phylogenetic diversity than their diploid progenitors. Polyploidy, along with its interactions with the inoculum source and genetic lineage, collectively explained 7% of the total variation in microbiome composition. Furthermore, polyploidy broadened the core microbiome, with autotetraploids having 15 unique bacterial taxa in addition to the 55 they shared with diploids. Our results show that whole-genome duplication directly leads to novelty in the plant microbiome and importantly that the effect is dependent on the genetic ancestry of the polyploid and generalizable over many environmental contexts.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

Chromosome-level reference genome assembly provides insights into the evolution of Pennisetum alopecuroides

Pennisetum alopecuroides is an important forage grass resource, which plays a vital role in ecological environment improvement. Therefore, the acquisition of P. alopecuroides genome resources is conducive to the study of the adaptability of Pennisetum species in ecological remediation and forage breeding development. Here we assembled a P. alopecuroides cv. 'Liqiu' genome at the chromosome level with a size of approximately 845.71 Mb, contig N50 of 84.83Mb, and genome integrity of 99.13% as assessed by CEGMA. A total of 833.41-Mb sequences were mounted on nine chromosomes by Hi-C technology. In total, 60.66% of the repetitive sequences and 34,312 genes were predicted. The genomic evolution analysis showed that P. alopecuroides cv. 'Liqiu' was isolated from Setaria 7.53-13.80 million years ago and from Cenchrus 5.33-8.99 million years ago, respectively. The whole-genome event analysis showed that P. alopecuroides cv. 'Liqiu' underwent two whole-genome duplication (WGD) events in the evolution process, and the duplication events occurred at a similar time to that of Oryza sativa and Setaria viridis. The completion of the genome sequencing of P. alopecuroides cv. 'Liqiu' provides data support for mining high-quality genetic resources of P. alopecuroides and provides a theoretical basis for the origin and evolutionary characteristics of Pennisetum.

Whole Genome Duplication and Gene Evolution in the Hyperdiverse Venomous Gastropods

The diversity of venomous organisms and the toxins they produce have been increasingly investigated, but taxonomic bias remains important. Neogastropods, a group of marine predators representing almost 22% of the known gastropod diversity, evolved a wide range of feeding strategies, including the production of toxins to subdue their preys. However, whether the diversity of these compounds is at the origin of the hyperdiversification of the group and how genome evolution may correlate with both the compounds and species diversities remain understudied. Among the available gastropods genomes, only eight, with uneven quality assemblies, belong to neogastropods. Here, we generated chromosome-level assemblies of two species belonging to the Tonnoidea and Muricoidea superfamilies (Monoplex corrugatus and Stramonita haemastoma). The two obtained high-quality genomes had 3 and 2.2 Gb, respectively, and 92-89% of the total assembly conformed 35 pseudochromosomes in each species. Through the analysis of syntenic blocks, Hox gene cluster duplication, and synonymous substitutions distribution pattern, we inferred the occurrence of a whole genome duplication event in both genomes. As these species are known to release venom, toxins were annotated in both genomes, but few of them were found in homologous chromosomes. A comparison of the expression of ohnolog genes (using transcriptomes from osphradium and salivary glands in S. haemastoma), where both copies were differentially expressed, showed that most of them had similar expression profiles. The high quality of these genomes makes them valuable reference in their respective taxa, facilitating the identification of genome-level processes at the origin of their evolutionary success.

Emergence of the Synucleins

This study establishes the origin and evolutionary history of the synuclein genes. A combination of phylogenetic analyses of the synucleins from twenty-two model species, characterization of local synteny similarities among humans, sharks and lampreys, and statistical comparisons among lamprey and human chromosomes, provides conclusive evidence for the current diversity of synuclein genes arising from the whole-genome duplications (WGDs) that occurred in vertebrates. An ancestral synuclein gene was duplicated in a first WGD, predating the diversification of all living vertebrates. The two resulting genes are still present in agnathan vertebrates. The second WGD, specific to the gnathostome lineage, led to the emergence of the three classical synuclein genes, SNCA, SNCB and SNCG, which are present in all jawed vertebrate lineages. Additional WGDs have added new genes in both agnathans and gnathostomes, while some gene losses have occurred in particular species. The emergence of synucleins through WGDs prevented these genes from experiencing dosage effects, thus avoiding the potential detrimental effects associated with individual duplications of genes that encode proteins prone to aggregation. Additional insights into the structural and functional features of synucleins are gained through the analysis of the highly divergent synuclein proteins present in chondrichthyans and agnathans.

Polyploidy impacts population growth and competition with diploids: multigenerational experiments reveal key life-history trade-offs

Ecological theory predicts that early generation polyploids ('neopolyploids') should quickly go extinct owing to the disadvantages of rarity and competition with their diploid progenitors. However, polyploids persist in natural habitats globally. This paradox has been addressed theoretically by recognizing that reproductive assurance of neopolyploids and niche differentiation can promote establishment. Despite this, the direct effects of polyploidy at the population level remain largely untested despite establishment being an intrinsically population-level process. We conducted population-level experiments where life-history investment in current and future growth was tracked in four lineage pairs of diploids and synthetic autotetraploids of the aquatic plant Spirodela polyrhiza. Population growth was evaluated with and without competition between diploids and neopolyploids across a range of nutrient treatments. Although neopolyploid populations produce more biomass, they reach lower population sizes and have reduced carrying capacities when growing alone or in competition across all nutrient treatments. Thus, contrary to individual-level studies, our population-level data suggest that neopolyploids are competitively inferior to diploids. Conversely, neopolyploid populations have greater investment in dormant propagule production than diploids. Our results show that neopolyploid populations should not persist based on current growth dynamics, but high potential future growth may allow polyploids to establish in subsequent seasons.

Evolutionary loss of the ß1-adrenergic receptor in salmonids

Whole-genome duplications (WGDs) have been at the heart of the diversification of ß-adrenergic receptors (ß-ARs) in vertebrates. Non-teleost jawed vertebrates typically possess three ß-AR genes: adrb1 (ß1-AR), adrb2 (ß2-AR), and adrb3 (ß3-AR), originating from the ancient 2R (two rounds) WGDs. Teleost fishes, owing to the teleost-specific WGD, have five ancestral adrb paralogs (adrb1, adrb2a, adrb2b, adrb3a and adrb3b). Salmonids are particularly intriguing from an evolutionary perspective as they experienced an additional WGD after separating from other teleosts. Moreover, adrenergic regulation in salmonids, especially rainbow trout, has been intensively studied for decades. However, the repertoire of adrb genes in salmonids has not been yet characterized. An exhaustive genome survey of diverse salmonids, spanning five genera, complemented by phylogenetic sequence analysis, revealed each species has seven adrb paralogs: two adrb2a, two adrb2b, two adrb3a and one adrb3b. Surprisingly, salmonids emerge as the first known jawed vertebrate lineage to lack adrb1. adrb1 is nevertheless highly expressed in the hearts of non-salmonid teleosts, indicating that the wealth of data on adrenergic regulation in salmonids should be generalised to other teleost fishes with caution. It is hypothesised that the loss of adrb1 could have been viable because of the evolutionary radiation of adrb2 and adrb3 genes attributable to the salmonid WGD.

The chromosome-level genome of Eucommia ulmoides provides insights into sex differentiation and α-linolenic acid biosynthesis

Eucommia ulmoides Oliver is a typical dioecious plant endemic to China that has great medicinal and economic value. Here, we report a high-quality chromosome-level female genome of E. ulmoides obtained by PacBio and Hi-C technologies. The size of the female genome assembly was 1.01 Gb with 17 pseudochromosomes and 31,665 protein coding genes. In addition, Hi-C technology was used to reassemble the male genome released in 2018. The reassembled male genome was 1.24 Gb with the superscaffold N50 (48.30 Mb), which was increased 25.69 times, and the number of predicted genes increased by 11,266. Genome evolution analysis indicated that E. ulmoides has undergone two whole-genome duplication events before the divergence of female and male, including core eudicot γ whole-genome triplication event (γ-WGT) and a recent whole genome duplication (WGD) at approximately 27.3 million years ago (Mya). Based on transcriptome analysis, EuAP3 and EuAG may be the key genes involved in regulating the sex differentiation of E. ulmoides. Pathway analysis showed that the high expression of ω-3 fatty acid desaturase coding gene EU0103017 was an important reason for the high α-linolenic acid content in E. ulmoides. The genome of female and male E. ulmoides presented here is a valuable resource for the molecular biological study of sex differentiation of E. ulmoides and also will provide assistance for the breeding of superior varieties.

Neutral processes underlying the macro eco-evolutionary dynamics of mixed-ploidy systems

Polyploidy, i.e. the occurrence of multiple sets of chromosomes, is regarded as an important phenomenon in plant ecology and evolution, with all flowering plants likely having a polyploid ancestry. Owing to genome shock, minority cytotype exclusion and reduced fertility, polyploids emerging in diploid populations are expected to face significant challenges to successful establishment. Their establishment and persistence are often explained by possible fitness or niche differences that would relieve the competitive pressure with diploid progenitors. Experimental evidence for such advantages is, however, not unambiguous, and considerable niche overlap exists among most polyploid species and their diploid counterparts. Here, we develop a neutral spatially explicit eco-evolutionary model to understand whether neutral processes can explain the eco-evolutionary patterns of polyploids. We present a general mechanism for polyploid establishment by showing that sexually reproducing organisms assemble in space in an iterative manner, reducing frequency-dependent mating disadvantages and overcoming potential reduced fertility issues. Moreover, we construct a mechanistic theoretical framework that allows us to understand the long-term evolution of mixed-ploidy populations and show that our model is remarkably consistent with recent phylogenomic estimates of species extinctions in the Brassicaceae family.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Autopolyploids of Arabidopsis thaliana are more phenotypically plastic than their diploid progenitors

BACKGROUND AND AIMS

Polyploids are often hypothesized to have increased phenotypic plasticity compared with their diploid progenitors, but recent work suggests that the relationship between whole-genome duplication (WGD) and plasticity is not so straightforward. Impacts of WGD on plasticity are moderated by other evolutionary processes in nature, which has impeded generalizations regarding the effects of WGD alone. We assessed shifts in phenotypic plasticity and mean trait values accompanying WGD, as well as the adaptive consequences of these shifts.

METHODS

To isolate WGD effects, we compared two diploid lineages of Arabidopsis thaliana wiht corresponding autotetraploids grown across different salt and nutrient conditions in a growth chamber.

KEY RESULTS

For the few cases in which diploids and polyploids differed in plasticity, polyploids were more plastic, consistent with hypotheses that WGD increases plasticity. Under stress, increased plasticity was often adaptive (associated with higher total seed mass), but in other cases plasticity was unrelated to fitness. Mean trait values and plasticity were equally likely to be affected by WGD, but the adaptive consequences of these shifts were often context dependent or lineage specific. For example, polyploids had extended life spans, a shift that was adaptive in one polyploid lineage under amenable conditions but was maladaptive in the other lineage under stress.

CONCLUSIONS

Our work shows that increased phenotypic plasticity can result from WGD alone, independent of other evolutionary processes. We find that the effects of WGD can differ depending on the genotype of the progenitor and the environmental context. Though our experiment was limited to two genotypes of a single species, these findings support the idea that WGD can indeed increase plasticity.

Massive effects on chromatin after ploidy rearrangement in doubled haploids

This article comments on:

Piskorz EW, Xu L, Ma Y, Jiang H. 2023. Doubled-haploid induction generates extensive differential DNA methylation in Arabidopsis. Journal of Experimental Botany 74, 835–847

Doubled-haploid induction generates extensive differential DNA methylation in Arabidopsis

Doubled haploids (DHs) fix traits from hybrids in one generation. DH induction includes two changes in ploidy levels typically associated with variation in DNA methylation. However, DNA methylation patterns in DH plants and their biological significance are largely unknown. We generated three DH lines in Arabidopsis thaliana by crossing a haploid inducer with the accession Col-0, thus removing tissue culture and hybridization as a variable. DH induction produced thousands of differentially DNA methylated regions (DMRs), most of which were stochastic. Both haploidization and colchicine-induced genome duplication produced DMRs; the former mainly yielded DMRs at non-CG contexts, whereas the latter affected differential gene body methylation. Spontaneous genome doubling of haploid plants also induced DMRs in greater numbers than self-propagation. Our results provide the first evidence that haploid induction and genome doubling result in differential DNA methylation, offering a novel approach to induce epialleles.

Cyclin E-induced replicative stress drives p53-dependent whole-genome duplication

Whole-genome duplication (WGD) is a frequent event in cancer evolution and an important driver of aneuploidy. The role of the p53 tumor suppressor in WGD has been enigmatic: p53 can block the proliferation of tetraploid cells, acting as a barrier to WGD, but can also promote mitotic bypass, a key step in WGD via endoreduplication. In wild-type (WT) p53 tumors, WGD is frequently associated with activation of the E2F pathway, especially amplification of CCNE1, encoding cyclin E1. Here, we show that elevated cyclin E1 expression causes replicative stress, which activates ATR- and Chk1-dependent G2 phase arrest. p53, via its downstream target p21, together with Wee1, then inhibits mitotic cyclin-dependent kinase activity sufficiently to activate APC/C^Cdh1 and promote mitotic bypass. Cyclin E expression suppresses p53-dependent senescence after mitotic bypass, allowing cells to complete endoreduplication. Our results indicate that p53 can contribute to cancer evolution through the promotion of WGD.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Where whole-genome duplication is most beneficial: Adaptation of mangroves to a wide salinity range between land and sea

Whole-genome duplication (WGD) is believed to increase the chance of adaptation to a new environment. This conjecture may apply particularly well to new environments that are not only different but also more variable than ancestral habitats. One such prominent environment is the interface between land and sea, which has been invaded by woody plants, collectively referred as mangroves, multiple times. Here, we use two distantly related mangrove species (Avicennia marina and Rhizophora apiculata) to explore the effects of WGD on the adaptive process. We found that a high proportion of duplicated genes retained after WGD have acquired derived differential expression in response to salt gradient treatment. The WGD duplicates differentially expressed in at least one copy usually (>90%) diverge from their paralogues' expression profiles. Furthermore, both species evolved in parallel to have one paralogue expressed at a high level in both fresh water and hypersaline conditions but at a lower level at medium salinity. The pattern contrasts with the conventional view of monotone increase/decrease as salinity increases. Differentially expressed copies have thus probably acquired a new role in salinity tolerance. Our results indicate that the WGD duplicates may have evolved to function collaboratively in coping with different salinity levels, rather than specializing in the intermediate salinity optimal for mangrove plants. In conclusion, WGD and the retained duplicates appear to be an effective solution for adaptation to new and unstable environments.

Phylogenetic, Structural and Functional Evolution of the LHC Gene Family in Plant Species

Light-harvesting chlorophyll a/b-binding (LHC) superfamily proteins play a vital role in photosynthesis. Although the physiological and biochemical functions of LHC genes have been well-characterized, the structural evolution and functional differentiation of the products need to be further studied. In this paper, we report the genome-wide identification and phylogenetic analysis of LHC genes in photosynthetic organisms. A total of 1222 non-redundant members of the LHC family were identified from 42 species. According to the phylogenetic clustering of their homologues with Arabidopsis thaliana, they can be divided into four subfamilies. In the subsequent evolution of land plants, a whole-genome replication (WGD) event was the driving force for the evolution and expansion of the LHC superfamily, with its copy numbers rapidly increasing in angiosperms. The selection pressure of photosystem II sub-unit S (PsbS) and ferrochelatase (FCII) families were higher than other subfamilies. In addition, the transcriptional expression profiles of LHC gene family members in different tissues and their expression patterns under exogenous abiotic stress conditions significantly differed, and the LHC genes are highly expressed in mature leaves, which is consistent with the conclusion that LHC is mainly involved in the capture and transmission of light energy in photosynthesis. According to the expression pattern and copy number of LHC genes in land plants, we propose different evolutionary trajectories in this gene family. This study provides a basis for understanding the molecular evolutionary characteristics and evolution patterns of plant LHCs.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

A Revised Perspective on the Evolution of Troponin I and Troponin T Gene Families in Vertebrates

The troponin (Tn) complex, responsible for the Ca2+ activation of striated muscle, is composed of three interacting protein subunits: TnC, TnI, and TnT, encoded by TNNC, TNNI, and TNNT genes. TNNI and TNNT are sister gene families, and in mammals the three TNNI paralogs (TNNI1, TNNI2, TNNI3), which encode proteins with tissue-specific expression, are each in close genomic proximity with one of the three TNNT paralogs (TNNT2, TNNT3, TNNT1, respectively). It has been widely presumed that all vertebrates broadly possess genes of these same three classes, although earlier work has overlooked jawless fishes (cyclostomes) and cartilaginous fishes (chimeras, rays, and sharks), which are distantly related to other jawed vertebrates. With a new phylogenetic and synteny analysis of a diverse array of vertebrates including these taxonomic groups, we define five distinct TNNI classes (TNNI1-5), with TNNI4 and TNNI5 being only present in non-amniote vertebrates and typically found in tandem, and four classes of TNNT (TNNT1-4). These genes are located in four genomic loci that were generated by the 2R whole-genome duplications. TNNI3, encoding "cardiac TnI" in tetrapods, was independently lost in cartilaginous and ray-finned fishes. Instead, ray-finned fishes predominantly express TNNI1 in the heart. TNNI5 is highly expressed in shark hearts and contains a N-terminal extension similar to that of TNNI3 found in tetrapod hearts. Given that TNNI3 and TNNI5 are distantly related, this supports the hypothesis that the N-terminal extension may be an ancestral feature of vertebrate TNNI and not an innovation unique to TNNI3, as has been commonly believed.

How to survive in the world's third poplar: Insights from the genome of the highest altitude woody plant, Hippophae tibetana (Elaeagnaceae)

Hippophae tibetana (Tibetan sea-buckthorn) is one of the highest distributed woody plants in the world (3,000-5,200 meters a.s.l.). It is characterized by adaptation to extreme environment and important economic values. Here, we combined PacBio Hifi platform and Hi-C technology to assemble a 1,452.75 Mb genome encoding 33,367 genes with a Contig N50 of 74.31 Mb, and inferred its sexual chromosome. Two Hippophae-specific whole-genome duplication events (18.7-21.2 million years ago, Ma; 28.6-32.4 Ma) and long terminal repeats retroelements (LTR-RTs) amplifications were detected. Comparing with related species at lower altitude, Ziziphus jujuba (<1, 700 meters a.s.l.), H. tibetana had some significantly rapid evolving genes involved in adaptation to high altitude habitats. However, comparing with Hippophae rhamnoides (<3, 700 meters a.s.l.), no rapid evolving genes were found except microtubule and microtubule-based process genes, H. tibetana has a larger genome, with extra 2, 503 genes (7.5%) and extra 680.46 Mb transposable elements (TEs) (46.84%). These results suggest that the changes in the copy number and regulatory pattern of genes play a more important role for H. tibetana adapting to more extreme and variable environments at higher altitude by more TEs and more genes increasing genome variability and expression plasticity. This suggestion was supported by two findings: nitrogen-fixing genes of H. tibetana having more copies, and intact TEs being significantly closer genes than fragmentary TEs. This study provided new insights into the evolution of alpine plants.

Chromosome-scale Amaranthus tricolor genome provides insights into the evolution of the genus Amaranthus and the mechanism of betalain biosynthesis

Amaranthus tricolor is a vegetable and ornamental amaranth, with high lysine, dietary fibre and squalene content. The red cultivar of A. tricolor possesses a high concentration of betalains, which has been used as natural food colorants. Here, we constructed the genome of A. tricolor, the first reference genome for the subgenus Albersia, combining PacBio HiFi, Nanopore ultra-long and Hi-C data. The contig N50 size was 906 kb, and 99.58% of contig sequence was anchored to the 17 chromosomes, totalling 520 Mb. We annotated 27,813 protein-coding genes with an average 1.3 kb coding sequence and 5.3 exons. We inferred that A. tricolor underwent a whole-genome duplication (WGD) and that the WGD shared by amaranths occurred in the last common ancestor of subfamily Amaranthoideae. Moreover, we comprehensively identified candidate genes in betalain biosynthesis pathway. Among them, DODAα1 and CYP76ADα1, located in one topologically associated domain (TAD) of an active (A) compartment on chromosome 16, were more highly expressed in red leaves than in green leaves, and DODAα1 might be the rate-limiting enzyme gene in betalains biosynthesis. This study presents new genome resources and enriches our understanding of amaranth evolution, betalains production, facilitating molecular breeding improvements and the understanding of C4 plants evolution.

The chromosome-level genome of female ginseng (Angelica sinensis) provides insights into molecular mechanisms and evolution of coumarin biosynthesis

Coumarins are natural products with important medicinal values, and include simple coumarins, furanocoumarins and pyranocoumarins. Female ginseng (Angelica sinensis) is a renowned herb with abundant coumarins, originated in China and known for the treatment of female ailments for thousands of years. The molecular basis of simple coumarin biosynthesis in A. sinensis and the evolutionary history of the genes involved in furanocoumarin biosynthesis are largely unknown. Here, we generated the first chromosome-scale genome of A. sinensis. It has a genome size of 2.37 Gb, which was generated by combining PacBio and Hi-C sequencing technologies. The genome was predicted to contain 43 202 protein-coding genes dispersed mainly on 11 pseudochromosomes. We not only provided evidence for whole-genome duplication (WGD) specifically occurring in the Apioideae subfamily, but also demonstrated the vital role of tandem duplication for phenylpropanoid biosynthesis in A. sinensis. Combined analyses of transcriptomic and metabolomic data revealed key genes and candidate transcription factors regulating simple coumarin biosynthesis. Furthermore, phylogenomic synteny network analyses suggested prenyltransferase genes involved in furanocoumarin biosynthesis evolved independently in the Moraceae, Fabaceae, Rutaceae and Apiaceae after ζ and ε WGD. Our work sheds light on coumarin biosynthesis, and provides a benchmark for accelerating genetic research and molecular breeding in A. sinensis.

The Chromosome-level genome of Aesculus wilsonii provides new insights into terpenoid biosynthesis and Aesculus evolution

Aesculus L. (buckeye and horse chestnut) are woody plant species with important horticultural and medicinal values. Aesculus seeds are widely used as biomedicine and cosmetic ingredients due to their saponins. We report a chromosomal-scale genome of Aesculus wilsonii. Sequences amounting to a total of 579.01 Mb were assembled into 20 chromosomes. More than half of the genome (54.46%) were annotated as repetitive sequences, and 46,914 protein-coding genes were predicted. In addition to the widespread gamma event with core eudicots, a unique whole-genome duplication (WGD) event (17.69 Mya) occurred in Aesculus after buckeye differentiated from longan. Due to WGD events and tandem duplications, the related synthetic genes of triterpene saponins unique to Aesculus increased significantly. Combined with transcriptome characterization, the study preliminarily resolved the biosynthetic pathway of triterpenoid saponins like aescin in A. wilsonii genome. Analyses of the resequencing of 104 buckeye accessions revealed clear relationship between the geographic distribution and genetic differentiation of buckeye trees in China. We found that the buckeye species found in southern Shaanxi is A. wilsonii rather than A. chinensis. Population dynamics analysis further suggests that the population size and evolution of existing buckeye species have been influenced by climate fluctuations during the Pleistocene and recent domestication events. The genome of A. wilsonii and population genomics of Aesculus provide a resource for future research on Hippocastanaceae. These findings will contribute to the utilization and diversity protection of Aesculus.

De novo genome assembly of the medicinal plant Gentiana macrophylla provides insights into the genomic evolution and biosynthesis of iridoids

Gentiana macrophylla is a perennial herb in the Gentianaceae family, whose dried roots are used in traditional Chinese medicine. Here, we assembled a chromosome-level genome of G. macrophylla using a combination of Nanopore, Illumina, and Hi-C scaffolding approaches. The final genome size was ~1.79 Gb (contig N50 = 720.804 kb), and 98.89% of the genome sequences were anchored on 13 pseudochromosomes (scaffold N50 = 122.73 Mb). The genome contained 55,337 protein-coding genes, and 73.47% of the assemblies were repetitive sequences. Genome evolution analysis indicated that G. macrophylla underwent two rounds of whole-genome duplication after the core eudicot γ genome triplication event. We further identified candidate genes related to the biosynthesis of iridoids, and the corresponding gene families mostly expanded in G. macrophylla. In addition, we found that root-specific genes are enriched in pathways involved in defense responses, which may greatly improve the biological adaptability of G. macrophylla. Phylogenomic analyses showed a sister relationship of asterids and rosids, and all Gentianales species formed a monophyletic group. Our study contributes to the understanding of genome evolution and active component biosynthesis in G. macrophylla and provides important genomic resource for the genetic improvement and breeding of G. macrophylla.

Polyploidy increases storage but decreases structural stability in Arabidopsis thaliana

Whole-genome duplication, leading to polyploidy and endopolyploidy, is widespread throughout the tree of life.^1-3 Both polyploidy and endopolyploidy can increase cell size via nucleotypic effects, but the phenotypic consequences of increased cell size at the tissue and whole-organism levels are less well understood.^1-4 We quantified the consequences of autopolyploidy and endopolyploidy in nine diploid accessions of Arabidopsis thaliana, representing a gradient in endopolyploidy, to their corresponding experimentally synthesized neo-tetraploid and neo-octoploid cytotypes. The increase in cell size following genome duplication increased plant storage capacity, which increased tolerance of resource limitation, but also incurred biomechanical costs because of a reduction in the amount of cell wall per unit tissue volume. Our findings also show that the functional consequences of autopolyploidy can vary with accession identity, and the presence of this variation suggests that there is potential for adaptation following whole-genome duplication.

Idahoa and Subularia: Hidden polyploid origins of two enigmatic genera of crucifers

PREMISE

The monotypic Idahoa (I. scapigera) and the bispecific Subularia (S. aquatica and S. monticola) belong to Brassicaceae with unclear phylogenetic relationships and no tribal assignment. To fill this knowledge gap, we investigated these species and their closest relatives by combining cytogenomic and phylogenomic methods.

METHODS

We used whole plastome sequences in maximum likelihood and Bayesian inference analyses. We tested the phylogenetic informativeness of shared genomic repeats. We combined nuclear gene tree reconciliation and comparative chromosome painting (CCP) to examine the occurrence of past whole-genome duplications (WGDs).

RESULTS

The plastid data set corroborated the sister relationship between Idahoa and Subularia within the crucifer Lineage V but failed to resolve consistent topologies using both inference methods. The shared repetitive sequences provided conflicting pwhylogenetic signals. CCP analysis unexpectedly revealed that Idahoa (2n = 16) has a diploidized mesotetraploid genome, whereas two Subularia species (2n = 28 and 30) have diploidized mesoctoploid genomes. Several ancient allopolyploidy events have also been detected in closely related taxa (Chamira circaeoides, Cremolobeae, Eudemeae, and Notothlaspideae).

CONCLUSIONS

Our results suggest that the contentious phylogenetic placement of Idahoa and Subularia is best explained by two WGDs involving one or more shared parental genomes. The newly identified mesopolyploid genomes highlight the challenges of studying plant clades with complex polyploidy histories and provide a better framework for understanding genome evolution in the crucifer family.

Meiosis in allopolyploid Arabidopsis suecica

Polyploidy is a major force shaping eukaryote evolution but poses challenges for meiotic chromosome segregation. As a result, first-generation polyploids often suffer from more meiotic errors and lower fertility than established wild polyploid populations. How established polyploids adapt their meiotic behaviour to ensure genome stability and accurate chromosome segregation remains an active research question. We present here a cytological description of meiosis in the model allopolyploid species Arabidopsis suecica (2n = 4x = 26). In large part meiosis in A. suecica is diploid-like, with normal synaptic progression and no evidence of synaptic partner exchanges. Some abnormalities were seen at low frequency, including univalents at metaphase I, anaphase bridges and aneuploidy at metaphase II; however, we saw no evidence of crossover formation occurring between non-homologous chromosomes. The crossover number in A. suecica is similar to the combined number reported from its diploid parents Arabidopsis thaliana (2n = 2x = 10) and Arabidopsis arenosa (2n = 2x = 16), with an average of approximately 1.75 crossovers per chromosome pair. This contrasts with naturally evolved autotetraploid A. arenosa, where accurate chromosome segregation is achieved by restricting crossovers to approximately 1 per chromosome pair. Although an autotetraploid donor is hypothesized to have contributed the A. arenosa subgenome to A. suecica, A. suecica harbours diploid A. arenosa variants of key meiotic genes. These multiple lines of evidence suggest that meiosis in the recently evolved allopolyploid A. suecica is essentially diploid like, with meiotic adaptation following a very different trajectory to that described for autotetraploid A. arenosa.

Comparative genomics of orobanchaceous species with different parasitic lifestyles reveals the origin and stepwise evolution of plant parasitism

Orobanchaceae is the largest family of parasitic plants, containing autotrophic and parasitic plants with all degrees of parasitism. This makes it by far the best family for studying the origin and evolution of plant parasitism. Here we provide three high-quality genomes of orobanchaceous plants, the autotrophic Lindenbergia luchunensis and the holoparasitic plants Phelipanche aegyptiaca and Orobanche cumana. Phylogenomic analysis of these three genomes together with those previously published and the transcriptomes of other orobanchaceous species created a robust phylogenetic framework for Orobanchaceae. We found that an ancient whole-genome duplication (WGD; about 73.48 million years ago), which occurred earlier than the origin of Orobanchaceae, might have contributed to the emergence of parasitism. However, no WGD events occurred in any lineage of orobanchaceous parasites except for Striga after divergence from their autotrophic common ancestor, suggesting that, in contrast with previous speculations, WGD is not associated with the emergence of holoparasitism. We detected evident convergent gene loss in all parasites within Orobanchaceae and between Orobanchaceae and dodder Cuscuta australis. The gene families in the orobanchaceous parasites showed a clear pattern of recent gains and expansions. The expanded gene families are enriched in functions related to the development of the haustorium, suggesting that recent gene family expansions may have facilitated the adaptation of orobanchaceous parasites to different hosts. This study illustrates a stepwise pattern in the evolution of parasitism in the orobanchaceous parasites and will facilitate future studies on parasitism and the control of parasitic plants in agriculture.

Polyploidization: Consequences of genome doubling on the evolutionary potential of populations

Whole-genome duplication is common in plants and is considered to have a broad range of effects on individuals' phenotypes and genomes and to be an important driver of plant adaptation and speciation. Despite their increased capacity to cope with challenging environments, polyploid lineages are generally as prone to extinction, and sometimes more prone, than their diploid progenitors. Although several explanations have been proposed to explain the short- and long-term disadvantages of polyploidy on the survival probability of populations, the consequences of whole-genome doubling on the heritable variance remain poorly studied. Whole-genome doubling can have major effects not only on the genetics, but also on the ecology and life history of the populations. Modifications of other properties of populations can reverse the effects of polyploidization per se on heritable variance. In this synthesis, I summarize the empirical and theoretical knowledge about the multifarious consequences of genome doubling on the heritable variance of quantitative traits and on the evolutionary potential of polyploid populations compared to their diploid progenitors. I propose several ways to decipher the consequences of whole-genome doubling on survival probability and to study the further consequences of shifting the ecological niche and life-history traits of a population. I also highlight some practical considerations for comparing the heritable variance of a trait among different cytotypes. Such investigations appear to be timely and necessary to understand more about the paradoxical aspects of polyploidization and to understand the evolutionary potential of polyploid lineages in a global warming context.

High-quality reference genome sequences of two Cannaceae species provide insights into the evolution of Cannaceae

Canna edulis Ker-Gawl and Canna indica L. are species belonging to the Cannaceae family and both have a very high economic value. Here, we aimed to assemble genomes of C. edulis and C. indica at the chromosome level to generate a reference genome for the Cannaceae family. We also comparatively analyzed the genomes of C. edulis and C. indica and examined the molecular mechanisms responsible for the remarkable differences in plant characteristics in C. edulis varieties. Our results indicated that genome-wide duplication events had recently occurred in C. edulis and C. indica. The comparative analysis of the genomes of C. edulis and C. indica revealed that C. edulis exhibited a remarkable level of replication of genes in the starch and sucrose metabolic pathways, especially during sucrose hydrolysis. This finding is consistent with the fact that the starch content of the C. edulis tuber is higher than that of C. indica. Simplified genome re-sequencing revealed the population structure of 241 C. edulis genes, and a genome-wide association study of leaf traits revealed the location of key genes related to leaf color and morphology. These findings extend our understanding of Cannaceade at the molecular level, and provide an effective theoretical basis for further study and utilization of Cannaceae plants.

The genome of Corydalis reveals the evolution of benzylisoquinoline alkaloid biosynthesis in Ranunculales

Species belonging to the order Ranunculales have attracted much attention because of their phylogenetic position as a sister group to all other eudicot lineages and their ability to produce unique yet diverse benzylisoquinoline alkaloids (BIAs). The Papaveraceae family in Ranunculales is often used as a model system for studying BIA biosynthesis. Here, we report the chromosome-level genome assembly of Corydalis tomentella, a species of Fumarioideae, one of the two subfamilies of Papaveraceae. Based on comparisons of sequenced Ranunculalean species, we present clear evidence of a shared whole-genome duplication (WGD) event that has occurred before the divergence of Ranunculales but after its divergence from other eudicot lineages. The C. tomentella genome enabled us to integrate isotopic labeling and comparative genomics to reconstruct the BIA biosynthetic pathway for both sanguinarine biosynthesis shared by papaveraceous species and the cavidine biosynthesis that is specific to Corydalis. Also, our comparative analysis revealed that gene duplications, especially tandem gene duplications, underlie the diversification of BIA biosynthetic pathways in Ranunculales. In particular, tandemly duplicated berberine bridge enzyme-like genes appear to be involved in cavidine biosynthesis. In conclusion, our study of the C. tomentella genome provides important insights into the occurrence of WGDs during the early evolution of eudicots, as well as into the evolution of BIA biosynthesis in Ranunculales.

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

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