Widespread Whole Genome Duplications Contribute to Genome Complexity and Species Diversity in Angiosperms

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

Genome engineering on size reduction and complexity simplification: A review

BACKGROUND

Genome simplification is an important topic in the field of life sciences that has attracted attention from its conception to the present day. It can help uncover the essential components of the genome and, in turn, shed light on the underlying operating principles of complex biological systems. This has made it a central focus of both basic and applied research in the life sciences. With the recent advancements in related technologies and our increasing knowledge of the genome, now is an opportune time to delve into this topic.

AIM OF REVIEW

Our review investigates the progress of genome simplification from two perspectives: genome size reduction and complexity simplification. In addition, we provide insights into the future development trends of genome simplification.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Reducing genome size requires eliminating non-essential elements as much as possible. This process has been facilitated by advances in genome manipulation and synthesis techniques. However, we still need a better and clearer understanding of living systems to reduce genome complexity. As there is a lack of quantitative and clearly defined standards for this task, we have opted to approach the topic from various perspectives and present our findings accordingly.

Diploid and tetraploid cytotypes of the flagship Cape species Dicerothamnus rhinocerotis (Asteraceae): variation in distribution, ecological niche, morphology and genetics

BACKGROUND AND AIMS

The Greater Cape Floristic Region is one of the world's biodiversity hotspots and is considered poor in polyploids. To test this assumption, ploidy variation was investigated in a widespread Cape shrub, Dicerothamnus rhinocerotis (renosterbos, Asteraceae). The aim was to elucidate the cytotype distribution and population composition across the species range, and to assess differences in morphology, environmental niches and genetics.

METHODS

Ploidy level and genome size were determined via flow cytometry and cytotype assignment was confirmed by chromosome counting. Restriction site-associated DNA sequencing (RADseq) analyses were used to infer genetic relationships. Cytotype climatic and environmental niches were compared using a range of environmental layers and a soil model, while morphological differences were examined using multivariate methods.

KEY RESULTS

The survey of 171 populations and 2370 individuals showed that the species comprises diploid and tetraploid cytotypes, no intermediates and only 16.8 % of mixed populations. Mean 2C values were 1.80-2.06 pg for diploids and 3.48-3.80 pg for tetraploids, with very similar monoploid genome sizes. Intra-cytotype variation showed a significant positive correlation with altitude and longitude in both cytotypes and with latitude in diploids. Although niches of both cytotypes were highly equivalent and similar, their optima and breadth were shifted due to differences mainly in isothermality and available water capacity. Morphometric analyses showed significant differences in the leaves and corolla traits, the number of florets per capitulum, and cypsela dimensions between the two cytotypes. Genetic analyses revealed four groups, three of them including both cytotypes.

CONCLUSIONS

Dicerothamnus rhinocerotis includes two distinct cytotypes that are genetically similar. While tetraploids arise several times independently within different genetic groups, morphological and ecological differences are evident between cytotypes. Our results open up new avenues for questions regarding the importance of ploidy in the megadiverse Cape flora, and exemplify the need for population-based studies focused on ploidy variation.

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.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

LocoGSE, a sequence-based genome size estimator for plants

Extensive research has focused on exploring the range of genome sizes in eukaryotes, with a particular emphasis on land plants, where significant variability has been observed. Accurate estimation of genome size is essential for various research purposes, but existing sequence-based methods have limitations, particularly for low-coverage datasets. In this study, we introduce LocoGSE, a novel genome size estimator designed specifically for low-coverage datasets generated by genome skimming approaches. LocoGSE relies on mapping the reads on single copy consensus proteins without the need for a reference genome assembly. We calibrated LocoGSE using 430 low-coverage Angiosperm genome skimming datasets and compared its performance against other estimators. Our results demonstrate that LocoGSE accurately predicts monoploid genome size even at very low depth of coverage (<1X) and on highly heterozygous samples. Additionally, LocoGSE provides stable estimates across individuals with varying ploidy levels. LocoGSE fills a gap in sequence-based plant genome size estimation by offering a user-friendly and reliable tool that does not rely on high coverage or reference assemblies. We anticipate that LocoGSE will facilitate plant genome size analysis and contribute to evolutionary and ecological studies in the field. Furthermore, at the cost of an initial calibration, LocoGSE can be used in other lineages.

PlantFUNCO: Integrative Functional Genomics Database Reveals Clues into Duplicates Divergence Evolution

Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Genome-wide analysis of the peanut CaM/CML gene family reveals that the AhCML69 gene is associated with resistance to Ralstonia solanacearum

BACKGROUND

Calmodulins (CaMs)/CaM-like proteins (CMLs) are crucial Ca2+-binding sensors that can decode and transduce Ca2+ signals during plant development and in response to various stimuli. The CaM/CML gene family has been characterized in many plant species, but this family has not yet been characterized and analyzed in peanut, especially for its functions in response to Ralstonia solanacearum. In this study, we performed a genome-wide analysis to analyze the CaM/CML genes and their functions in resistance to R. solanacearum.

RESULTS

Here, 67, 72, and 214 CaM/CML genes were identified from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. The genes were divided into nine subgroups (Groups I-IX) with relatively conserved exon‒intron structures and motif compositions. Gene duplication, which included whole-genome duplication, tandem repeats, scattered repeats, and unconnected repeats, produced approximately 81 pairs of homologous genes in the AhCaM/CML gene family. Allopolyploidization was the main reason for the greater number of AhCaM/CML members. The nonsynonymous (Ka) versus synonymous (Ks) substitution rates (less than 1.0) suggested that all homologous pairs underwent intensive purifying selection pressure during evolution. AhCML69 was constitutively expressed in different tissues of peanut plants and was involved in the response to R. solanacearum infection. The AhCML69 protein was localized in the cytoplasm and nucleus. Transient overexpression of AhCML69 in tobacco leaves increased resistance to R. solanacearum infection and induced the expression of defense-related genes, suggesting that AhCML69 is a positive regulator of disease resistance.

CONCLUSIONS

This study provides the first comprehensive analysis of the AhCaM/CML gene family and potential genetic resources for the molecular design and breeding of peanut bacterial wilt resistance.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Genome-wide analysis of the ABC gene family in almond and functional predictions during flower development, freezing stress, and salt stress

ABC (ATP-binding cassette) transporter proteins are one of the most extensive protein families known to date and are ubiquitously found in animals, plants, and microorganisms. ABCs have a variety of functions, such as plant tissue development regulation, hormone transport, and biotic and abiotic stress resistance. However, the gene characterization and function of the ABC gene family in almond (Prunus dulcis) have not been thoroughly studied. In this study, we identified 117 PdABC genes using the whole genome of 'Wanfeng' almond obtained by sequencing and explored their protein characterization. The PdABC family members were classified into eight subfamilies. The members of the same subfamily had conserved motifs but poorly conserved numbers of exons and introns and were unevenly distributed among the eight subfamilies and on the eight chromosomes. Expression patterns showed that PdABC family members were significantly differentially expressed during almond development, dormant freezing stress, and salt stress. We found that PdABC59 and PdABC77 had extremely high expression levels in pollen. PdABC63 and PdABC64 had high expression levels during almond petal development and multiple stages of flower development. PdABC98 was highly expressed in annual dormant branches after six temperature-freezing stress treatments. PdABC29, PdABC69, and PdABC98 were highly expressed under different concentrations of salt stress. This study preliminarily investigated the expression characteristics of ABC genes in different tissues of almond during flower development, freezing stress and salt stress, and the results will provide a reference for further in-depth research and breeding of almond in the future.

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.

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Phylogenomics insights into gene evolution, rapid species diversification, and morphological innovation of the apple tribe (Maleae, Rosaceae)

Maleae is one of the most widespread tribes of Rosaceae and includes several important fruit crops and ornamental plants. We used nuclear genes from 62 transcriptomes/genomes, including 26 newly generated transcriptomes, to reconstruct a well-supported phylogeny and study the evolution of fruit and leaf morphology and the possible effect of whole genome duplication (WGD). Our phylogeny recovered 11 well-supported clades and supported the monophyly of most genera (except Malus, Sorbus, and Pourthiaea) with at least two sampled species. A WGD was located to the most recent common ancestor (MRCA) of Maleae and dated to c. 54 million years ago (Ma) near the Early Eocene Climatic Optimum, supporting Gillenieae (x = 9) being a parental lineage of Maleae (x = 17) and including duplicate regulatory genes related to the origin of the fleshy pome fruit. Whole genome duplication-derived paralogs that are retained in specific lineages but lost in others are predicted to function in development, metabolism, and other processes. An upshift of diversification and innovations of fruit and leaf morphologies occurred at the MRCA of the Malinae subtribe, coinciding with the Eocene-Oligocene transition (c. 34 Ma), following a lag from the time of the WGD event. Our results provide new insights into the Maleae phylogeny, its rapid diversification, and morphological and molecular evolution.

Geographical and life-history traits associated with low and high species richness across angiosperm families

Introduction

The phenomenal expansion of angiosperms has prompted many investigations into the factors driving their diversification, but there remain significant gaps in our understanding of flowering plant species diversity.

Methods

Using the crown age of families from five studies, we used a maximum likelihood approach to classify families as having poor, predicted or high species richness (SR) using strict consensus criteria. Using these categories, we looked for associations between family SR and i) the presence of an inferred familial ancestral polyploidization event, ii) 23 life history and floral traits compiled from previously published datasets and papers, and iii) sexual system (dioecy) or genetically determined self-incompatibility (SI) mating system using an updated version of our own database and iv) geographic distribution using a new database describing the global distribution of plant species/families across realms and biomes and inferred range.

Results

We find that more than a third of angiosperm families (65%) had predicted SR, a large proportion (30.2%) were species poor, while few (4.8%) had high SR. Families with poor SR were less likely to have undergone an ancestral polyploidization event, exhibited deficits in diverse traits, and were more likely to have unknown breeding systems and to be found in only one or few biomes and realms, especially the Afrotropics or Australasia. On the other hand, families with high SR were more likely to have animal mediated pollination or dispersal, are enriched for epiphytes and taxa with an annual life history, and were more likely to harbour sporophytic SI systems. Mapping the global distribution of georeferenced taxa by their family DR, we find evidence of regions dominated by taxa from lineages with high vs low SR.

Discussion

These results are discussed within the context of the literature describing "depauperons" and the factors contributing to low and high biodiversity in angiosperm clades.

Integrating genomic and multiomic data for Angelica sinensis provides insights into the evolution and biosynthesis of pharmaceutically bioactive compounds

Angelica sinensis roots (Angelica roots) are rich in many bioactive compounds, including phthalides, coumarins, lignans, and terpenoids. However, the molecular bases for their biosynthesis are still poorly understood. Here, an improved chromosome-scale genome for A. sinensis var. Qinggui1 is reported, with a size of 2.16 Gb, contig N50 of 4.96 Mb and scaffold N50 of 198.27 Mb, covering 99.8% of the estimated genome. Additionally, by integrating genome sequencing, metabolomic profiling, and transcriptome analysis of normally growing and early-flowering Angelica roots that exhibit dramatically different metabolite profiles, the pathways and critical metabolic genes for the biosynthesis of these major bioactive components in Angelica roots have been deciphered. Multiomic analyses have also revealed the evolution and regulation of key metabolic genes for the biosynthesis of pharmaceutically bioactive components; in particular, TPSs for terpenoid volatiles, ACCs for malonyl CoA, PKSs for phthalide, and PTs for coumarin biosynthesis were expanded in the A. sinensis genome. These findings provide new insights into the biosynthesis of pharmaceutically important compounds in Angelica roots for exploration of synthetic biology and genetic improvement of herbal quality.

Exploration of the truncated cytosolic Hsp70 in plants - unveiling the diverse T1 lineage and the conserved T2 lineage

The 70-kDa heat shock proteins (Hsp70s) are chaperone proteins involved in protein folding processes. Truncated Hsp70 (Hsp70T) refers to the variant lacking a conserved C-terminal motif, which is crucial for co-chaperone interactions or protein retention. Despite their significance, the characteristics of Hsp70Ts in plants remain largely unexplored. In this study, we performed a comprehensive genome-wide analysis of 192 sequenced plant and green algae genomes to investigate the distribution and features of Hsp70Ts. Our findings unveil the widespread occurrence of Hsp70Ts across all four Hsp70 forms, including cytosolic, endoplasmic reticulum, mitochondrial, and chloroplast Hsp70s, with cytosolic Hsp70T being the most prevalent and abundant subtype. Cytosolic Hsp70T is characterized by two distinct lineages, referred to as T1 and T2. Among the investigated plant and green algae species, T1 genes were identified in approximately 60% of cases, showcasing a variable gene count ranging from one to several dozens. In contrast, T2 genes were prevalent across the majority of plant genomes, usually occurring in fewer than five gene copies per species. Sequence analysis highlights that the putative T1 proteins exhibit higher similarity to full-length cytosolic Hsp70s in comparison to T2 proteins. Intriguingly, the T2 lineage demonstrates a higher level of conservation within their protein sequences, whereas the T1 lineage presents a diverse range in the C-terminal and SBDα region, leading to categorization into four distinct subtypes. Furthermore, we have observed that T1-rich species characterized by the possession of 15 or more T1 genes exhibit an expansion of T1 genes into tandem gene clusters. The T1 gene clusters identified within the Laurales order display synteny with clusters found in a species of the Chloranthales order and another species within basal angiosperms, suggesting a conserved evolutionary relationship of T1 gene clusters among these plants. Additionally, T2 genes demonstrate distinct expression patterns in seeds and under heat stress, implying their potential roles in seed development and stress response.

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

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.

Production of Single-Cell Oil from a Local Isolate Bacillus subtilis Using Palm Fronds

This study, conducted at the Department of Biology, University of Anbar, Iraq, focuses on addressing the escalating issue of contamination and aims to acquire microbial oils to alleviate the global shortage in plant and animal oil production, utilizing environmental waste as a carbon source to reduce global pollution and select efficient local bacterial isolates of Bacillus subtilis for the production of single-cell oil (SCO) using local soil and environmental waste as a carbon source. Four isolates were selected as the best in producing single-cell oil, with the isolate with code C4 standing out as it recorded the highest production. It is worth noting that all these isolates belong to the bacteria type Bacillus subtilis. Palm fronds were found to be the most suitable environmental residue for SCO production compared to other waste materials (wheat straw and wheat bran). Submerged cultures were used to improve SCO production, with optimal conditions determined as pH 7, a temperature of 30°C, carbon source concentration of 3 g/100 ml, inoculum volume of 3 ml/100 ml, inoculum density of 20 × 107 cells, and an incubation period of 72 hours. The Soxhlet extraction method was used to obtain the oil, which was found to contain high percentages of unsaturated fatty acids, particularly linoleic acid (46.030%) and palmitoleic acid (16.579%). The oil was highly soluble in chloroform and ethanol but insoluble in water. The saponification test indicated the potential for soap production from the oil. This comprehensive research addresses the need for locally sourced and sustainable SCO production, offering insights into the selection of efficient bacterial isolates, the optimization of cultivation conditions, and the valuable properties of the resulting SCO. The significance of this study lies in the production of single-cell oil from soil-isolated Bacillus subtilis bacteria for use in food applications.

Deciphering peanut complex genomes paves a way to understand its origin and domestication

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

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.

The chromosome-scale reference genome and transcriptome analysis of Solanum torvum provides insights into resistance to root-knot nematodes

Solanum torvum (Swartz) (2n = 24) is a wild Solanaceae plant with high economic value that is used as a rootstock in grafting for Solanaceae plants to improve the resistance to a soil-borne disease caused by root-knot nematodes (RKNs). However, the lack of a high-quality reference genome of S. torvum hinders research on the genetic basis for disease resistance and application in horticulture. Herein, we present a chromosome-level assembly of genomic sequences for S. torvum combining PacBio long reads (HiFi reads), Illumina short reads and Hi-C scaffolding technology. The assembled genome size is ~1.25 Gb with a contig N50 and scaffold N50 of 38.65 Mb and 103.02 Mb, respectively as well as a BUSCO estimate of 98%. GO enrichment and KEGG pathway analysis of the unique S. torvum genes, including NLR and ABC transporters, revealed that they were involved in disease resistance processes. RNA-seq data also confirmed that 48 NLR genes were highly expressed in roots and fibrous roots and that three homologous NLR genes (Sto0288260.1, Sto0201960.1 and Sto0265490.1) in S. torvum were significantly upregulated after RKN infection. Two ABC transporters, ABCB9 and ABCB11 were identified as the hub genes in response to RKN infection. The chromosome-scale reference genome of the S. torvum will provide insights into RKN resistance.

Combined genomic and transcriptomic analysis reveals the contribution of tandem duplication genes to low-temperature adaptation in perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is an agronomically important cool-season grass species that is widely used as forage for ruminant animal production and cultivated in temperate regions for the establishment of lawns. However, the underlying genetic mechanism of the response of L. perenne to low temperature is still unclear. In the present study, we performed a comprehensive study and identified 3,770 tandem duplication genes (TDGs) in L. perenne, and evolutionary analysis revealed that L. perenne might have undergone a duplication event approximately 7.69 Mya. GO and KEGG pathway functional analyses revealed that these TDGs were mainly enriched in photosynthesis, hormone-mediated signaling pathways and responses to various stresses, suggesting that TDGs contribute to the environmental adaptability of L. perenne. In addition, the expression profile analysis revealed that the expression levels of TDGs were highly conserved and significantly lower than those of all genes in different tissues, while the frequency of differentially expressed genes (DEGs) from TDGs was much higher than that of DEGs from all genes in response to low-temperature stress. Finally, in-depth analysis of the important and expanded gene family indicated that the members of the ELIP subfamily could rapidly respond to low temperature and persistently maintain higher expression levels during all low temperature stress time points, suggesting that ELIPs most likely mediate low temperature responses and help to facilitate adaptation to low temperature in L. perenne. Our results provide evidence for the genetic underpinning of low-temperature adaptation and valuable resources for practical application and genetic improvement for stress resistance in L. perenne.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Genome-Wide Identification and Expression Analysis of Auxin Response Factor Gene Family in Linum usitatissimum

Auxin response factors (ARFs) are critical components of the auxin signaling pathway, and are involved in diverse plant biological processes. However, ARF genes have not been investigated in flax (Linum usitatissimum L.), an important oilseed and fiber crop. In this study, we comprehensively analyzed the ARF gene family and identified 33 LuARF genes unevenly distributed on the 13 chromosomes of Longya-10, an oil-use flax variety. Detailed analysis revealed wide variation among the ARF family members and predicted nuclear localization for all proteins. Nineteen LuARFs contained a complete ARF structure, including DBD, MR, and CTD, whereas the other fourteen lacked the CTD. Phylogenetic analysis grouped the LuARFs into four (I-V) clades. Combined with sequence analysis, the LuARFs from the same clade showed structural conservation, implying functional redundancy. Duplication analysis identified twenty-seven whole-genome-duplicated LuARF genes and four tandem-duplicated LuARF genes. These duplicated gene pairs' K_a/K_s ratios suggested a strong purifying selection pressure on the LuARF genes. Collinearity analysis revealed that about half of the LuARF genes had homologs in other species, indicating a relatively conserved nature of the ARFs. The promoter analysis identified numerous hormone- and stress-related elements, and the qRT-PCR experiment revealed that all LuARF genes were responsive to phytohormone (IAA, GA3, and NAA) and stress (PEG, NaCl, cold, and heat) treatments. Finally, expression profiling of LuARF genes in different tissues by qRT-PCR indicated their specific functions in stem or capsule growth. Thus, our findings suggest the potential functions of LuARFs in flax growth and response to an exogenous stimulus, providing a basis for further functional studies on these genes.

CALANGO: A phylogeny-aware comparative genomics tool for discovering quantitative genotype-phenotype associations across species

Living species vary significantly in phenotype and genomic content. Sophisticated statistical methods linking genes with phenotypes within a species have led to breakthroughs in complex genetic diseases and genetic breeding. Despite the abundance of genomic and phenotypic data available for thousands of species, finding genotype-phenotype associations across species is challenging due to the non-independence of species data resulting from common ancestry. To address this, we present CALANGO (comparative analysis with annotation-based genomic components), a phylogeny-aware comparative genomics tool to find homologous regions and biological roles associated with quantitative phenotypes across species. In two case studies, CALANGO identified both known and previously unidentified genotype-phenotype associations. The first study revealed unknown aspects of the ecological interaction between Escherichia coli, its integrated bacteriophages, and the pathogenicity phenotype. The second identified an association between maximum height in angiosperms and the expansion of a reproductive mechanism that prevents inbreeding and increases genetic diversity, with implications for conservation biology and agriculture.

An Emerging Animal Model for Querying the Role of Whole Genome Duplication in Development, Evolution, and Disease

Whole genome duplication (WGD) or polyploidization can occur at the cellular, tissue, and organismal levels. At the cellular level, tetraploidization has been proposed as a driver of aneuploidy and genome instability and correlates strongly with cancer progression, metastasis, and the development of drug resistance. WGD is also a key developmental strategy for regulating cell size, metabolism, and cellular function. In specific tissues, WGD is involved in normal development (e.g., organogenesis), tissue homeostasis, wound healing, and regeneration. At the organismal level, WGD propels evolutionary processes such as adaptation, speciation, and crop domestication. An essential strategy to further our understanding of the mechanisms promoting WGD and its effects is to compare isogenic strains that differ only in their ploidy. Caenorhabditis elegans (C. elegans) is emerging as an animal model for these comparisons, in part because relatively stable and fertile tetraploid strains can be produced rapidly from nearly any diploid strain. Here, we review the use of Caenorhabditis polyploids as tools to understand important developmental processes (e.g., sex determination, dosage compensation, and allometric relationships) and cellular processes (e.g., cell cycle regulation and chromosome dynamics during meiosis). We also discuss how the unique characteristics of the C. elegans WGD model will enable significant advances in our understanding of the mechanisms of polyploidization and its role in development and disease.

Plant NADPH-dependent thioredoxin reductases are crucial for the metabolism of sink leaves and plant acclimation to elevated CO2

Plants contain three NADPH-thioredoxin reductases (NTR) located in the cytosol/mitochondria (NTRA/B) and the plastid (NTRC) with important metabolic functions. However, mutants deficient in all NTRs remained to be investigated. Here, we generated and characterised the triple Arabidopsis ntrabc mutant alongside with ntrc single and ntrab double mutants under different environmental conditions. Both ntrc and ntrabc mutants showed reduced growth and substantial metabolic alterations, especially in sink leaves and under high CO₂ (HC), as compared to the wild type. However, ntrabc showed higher effective quantum yield of PSII under both constant and fluctuating light conditions, altered redox states of NADH/NAD⁺ and glutathione (GSH/GSSG) and lower potential quantum yield of PSII in sink leaves in ambient but not high CO₂ concentrations, as compared to ntrc, suggesting a functional interaction between chloroplastic and extra-chloroplastic NTRs in photosynthesis regulation depending on leaf development and environmental conditions. Our results unveil a previously unknown role of the NTR system in regulating sink leaf metabolism and plant acclimation to HC, while it is not affecting full plant development, indicating that the lack of the NTR system can be compensated, at least to some extent, by other redox mechanisms.

Two-step model of paleohexaploidy, ancestral genome reshuffling and plasticity of heat shock response in Asteraceae

An ancient hexaploidization event in the most but not all Asteraceae plants, may have been responsible for shaping the genomes of many horticultural, ornamental, and medicinal plants that promoting the prosperity of the largest angiosperm family on the earth. However, the duplication process of this hexaploidy, as well as the genomic and phenotypic diversity of extant Asteraceae plants caused by paleogenome reorganization, are still poorly understood. We analyzed 11 genomes from 10 genera in Asteraceae, and redated the Asteraceae common hexaploidization (ACH) event ~70.7-78.6 million years ago (Mya) and the Asteroideae specific tetraploidization (AST) event ~41.6-46.2 Mya. Moreover, we identified the genomic homologies generated from the ACH, AST and speciation events, and constructed a multiple genome alignment framework for Asteraceae. Subsequently, we revealed biased fractionations between the paleopolyploidization produced subgenomes, suggesting the ACH and AST both are allopolyplodization events. Interestingly, the paleochromosome reshuffling traces provided clear evidence for the two-step duplications of ACH event in Asteraceae. Furthermore, we reconstructed ancestral Asteraceae karyotype (AAK) that has 9 paleochromosomes, and revealed a highly flexible reshuffling of Asteraceae paleogenome. Of specific significance, we explored the genetic diversity of Heat Shock Transcription Factors (Hsfs) associated with recursive whole-genome polyploidizations, gene duplications, and paleogenome reshuffling, and revealed that the expansion of Hsfs gene families enable heat shock plasticity during the genome evolution of Asteraceae. Our study provides insights on polyploidy and paleogenome remodeling for the successful establishment of Asteraceae, and is helpful for further communication and exploration of the diversification of plant families and phenotypes.

Genome-Wide Identification and Analysis of OsSPXs Revealed Its Genetic Influence on Cold Tolerance of Dongxiang Wild Rice (DXWR)

SPX-domain proteins (small proteins with only the SPX domain) have been proven to be involved in phosphate-related signal transduction and regulation pathways. Except for OsSPX1 research showing that it plays a role in the process of rice adaptation to cold stress, the potential functions of other SPX genes in cold stress are unknown. Therefore, in this study, we identified six OsSPXs from the whole genome of DXWR. The phylogeny of OsSPXs has a strong correlation with its motif. Transcriptome data analysis showed that OsSPXs were highly sensitive to cold stress, and real-time PCR verified that the levels of OsSPX1, OsSPX2, OsSPX4, and OsSPX6 in cold-tolerant materials (DXWR) during cold treatment were higher than that of cold-sensitive rice (GZX49). The promoter region of DXWR OsSPXs contains a large number of cis-acting elements related to abiotic stress tolerance and plant hormone response. At the same time, these genes have expression patterns that are highly similar to cold-tolerance genes. This study provides useful information about OsSPXs, which is helpful for the gene-function research of DXWR and genetic improvements during breeding.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Global transcriptome analysis of allopolyploidization reveals large-scale repression of the D-subgenome in synthetic hexaploid wheat

Synthetic hexaploid wheat (SHW) lines are created as pre-breeding germplasm to diversify the D subgenome of hexaploid wheat and capitalize upon the untapped genetic diversity of the Aegilops tauschii gene pool. However, the phenotypes observed in the Ae. tauschii parents are not always recovered in the SHW lines, possibly due to inter-subgenome interactions. To elucidate this post-polyploidization genome reprogramming phenomenon, we performed RNA-seq of four SHW lines and their corresponding tetraploid and diploid parents, across ten tissues and three biological replicates. Homoeologue expression bias (HEB) analysis using more than 18,000 triads suggests massive suppression of homoeoalleles of the D subgenome in SHWs. Comparative transcriptome analysis of the whole-genome gene set further corroborated this finding. Alternative splicing analysis of the high-confidence genes indicates an additional layer of complexity where all five splice events are identified, and retained intron is predominant. Homoeologue expression upon resynthesis of hexaploid wheat has implications to the usage and handling of this germplasm in breeding as it relates to capturing the effects of epistatic interaction across subgenomes upon polyploidization. Special considerations must be given to this germplasm in pre-breeding activities to consider the extent of the inter-subgenome interactions on gene expression and their impact on traits for crop improvement.

Maize β-amylase7 encodes two proteins using alternative transcriptional start sites: nuclear BAM7 and plastidic BAM2

An unusual β-amylase7 (BAM7) gene in some angiosperms, including grasses such as maize (Zea mays), appears to encode two functionally distinct proteins; a nuclear-localized transcription factor (BAM7) and a plastid-localized starch hydrolase (BAM2). In Arabidopsis (Arabidopsis thaliana), these two proteins are encoded by separate genes on different chromosomes but their physiological functions are not well established. Using the maize BAM7 gene as a model, we detected two populations of transcripts by 5'-RACE that encode the predicted proteins. The two transcripts are apparently synthesized independently using separate core promoters about 1 kb apart, the second of which is located in the first intron of the full-length gene. The N-terminus of the shorter protein, ZmBAM7-S, begins near the 3' end of the first intron of ZmBAM7-L and starts with a predicted chloroplast transit peptide. We previously showed that ZmBAM7-S is catalytically active with properties like those of AtBAM2. Here, we report that ZmBAM7-S targets green fluorescent protein to plastids. The transcript encoding the longer protein, ZmBAM7-L, encodes an additional DNA-binding domain containing a functional nuclear localization signal. This putative dual-function gene originated at least 400 Mya, prior to the emergence of ferns, and has persisted in some angiosperms that lack a separate BAM2 gene. It appears to have been duplicated and subfunctionalized in at least four lineages of land plants, resulting in two genes resembling Arabidopsis BAM2 and BAM7. Targeting of two products from a single gene to different subcellular locations is not uncommon in plants, but it is unusual when they are predicted to serve completely different functions in the two locations.

Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis

BACKGROUND

Homoeologs are defined as homologous genes resulting from allopolyploidy. Bread wheat, Triticum aestivum, is an allohexaploid species with many homoeologs. Homoeolog expression bias, referring to the relative contribution of homoeologs to the transcriptome, is critical for determining the traits that influence wheat growth and development. Asymmetric transcription of homoeologs has been so far investigated in a tissue or organ-specific manner, which could be misleading due to a mixture of cell types.

RESULTS

Here, we perform single nuclei RNA sequencing and ATAC sequencing of wheat root to study the asymmetric gene transcription, reconstruct cell differentiation trajectories and cell-type-specific gene regulatory networks. We identify 22 cell types. We then reconstruct cell differentiation trajectories that suggest different origins between epidermis/cortex and endodermis, distinguishing bread wheat from Arabidopsis. We show that the ratio of asymmetrically transcribed triads varies greatly when analyzing at the single-cell level. Hub transcription factors determining cell type identity are also identified. In particular, we demonstrate that TaSPL14 participates in vasculature development by regulating the expression of BAM1. Combining single-cell transcription and chromatin accessibility data, we construct the pseudo-time regulatory network driving root hair differentiation. We find MYB3R4, REF6, HDG1, and GATAs as key regulators in this process.

CONCLUSIONS

Our findings reveal the transcriptional landscape of root organization and asymmetric gene transcription at single-cell resolution in polyploid wheat.

The use of Anchored Hybrid Enrichment data to resolve higher-level phylogenetic relationships: A proof-of-concept applied to Asterales (Eudicotyledoneae; Angiosperms)

Anchored Hybrid Enrichment (AHE) is a tool for capturing orthologous regions of the nuclear genome shared in low or single copy across lineages. Despite the increasing number of studies using this method, its usefulness to estimate relationships at deeper taxonomic levels in plants has not been fully explored. Here we present a proof of concept about the performance of nuclear loci obtained with AHE to infer phylogenetic relationships and explore the use of gene sampling schemes to estimate divergence times in Asterales. We recovered low-copy nuclear loci using the AHE method from herbarium material and silica-preserved samples. Maximum likelihood, Bayesian inference, and coalescence approaches were used to reconstruct phylogenomic relationships. Dating analyses were conducted under a multispecies coalescent approach by jointly inferring species tree and divergence times with random gene sampling schemes and multiple calibrations. We recovered 403 low-copy nuclear loci for 63 species representing nine out of eleven families of Asterales. Phylogenetic hypotheses were congruent among the applied methods and previously published results. Analyses with concatenated datasets were strongly supported, but coalescence-based analyses showed low support for the phylogenetic position of families Argophyllaceae and Alseuosmiaceae. Estimated family ages were congruent among gene sampling schemes, with the mean age for Asterales around 130 Myr. Our study documents the usefulness of AHE for resolving phylogenetic relationships at deep phylogenetic levels in Asterales. Observed phylogenetic inconsistencies were possibly due to the non-inclusion of families Phellinceae and Pentaphragmataceae. Random gene sampling schemes produced consistent age estimates with coalescence and species tree relaxed clock approaches.

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.

Insights into the differentiation and adaptation within Circaeasteraceae from Circaeaster agrestis genome sequencing and resequencing

Circaeaster agrestis and Kingdonia uniflora are sister species that reproduce sexually and mainly asexually respectively, providing a good system for comparative genome evolution between taxa with different reproductive models. Comparative genome analyses revealed the two species have similar genome size, but C. agrestis encodes many more genes. The gene families specific to C. agrestis show significant enrichment of genes associated with defense response, while those gene families specific to K. uniflora are enriched in genes regulating root system development. Collinearity analyses revealed C. agrestis experienced two rounds of whole-genome duplication. Fst outlier test across 25 C. agrestis populations uncovered a close inter-relationship between abiotic stress and genetic variability. Genetic feature comparisons showed K. uniflora presents much higher genome heterozygosity, transposable element load, linkage disequilibrium degree, and π_N/π_S ratio. This study provides new insights into understanding the genetic differentiation and adaptation within ancient lineages characterized by multiple reproductive models.

Evolutionary and Genomic Diversity of True Polyploidy in Tetrapods

True polyploid organisms have more than two chromosome sets in their somatic and germline cells. Polyploidy is a major evolutionary force and has played a significant role in the early genomic evolution of plants, different invertebrate taxa, chordates, and teleosts. However, the contribution of polyploidy to the generation of new genomic, ecological, and species diversity in tetrapods has traditionally been underestimated. Indeed, polyploidy represents an important pathway of genomic evolution, occurring in most higher-taxa tetrapods and displaying a variety of different forms, genomic configurations, and biological implications. Herein, we report and discuss the available information on the different origins and evolutionary and ecological significance of true polyploidy in tetrapods. Among the main tetrapod lineages, modern amphibians have an unparalleled diversity of polyploids and, until recently, they were considered to be the only vertebrates with closely related diploid and polyploid bisexual species or populations. In reptiles, polyploidy was thought to be restricted to squamates and associated with parthenogenesis. In birds and mammals, true polyploidy has generally been considered absent (non-tolerated). These views are being changed due to an accumulation of new data, and the impact as well as the different evolutionary and ecological implications of polyploidy in tetrapods, deserve a broader evaluation.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Phylotranscriptomic analyses reveal multiple whole-genome duplication events, the history of diversification and adaptations in the Araceae

BACKGROUND AND AIMS

The Araceae are one of the most diverse monocot families with numerous morphological and ecological novelties. Plastid and mitochondrial genes have been used to investigate the phylogeny and to interpret shifts in the pollination biology and biogeography of the Araceae. In contrast, the role of whole-genome duplication (WGD) in the evolution of eight subfamilies remains unclear.

METHODS

New transcriptomes or low-depth whole-genome sequences of 65 species were generated through Illumina sequencing. We reconstructed the phylogenetic relationships of Araceae using concatenated and species tree methods, and then estimated the age of major clades using TreePL. We inferred the WGD events by Ks and gene tree methods. We investigated the diversification patterns applying time-dependent and trait-dependent models. The expansions of gene families and functional enrichments were analysed using CAFE and InterProScan.

KEY RESULTS

Gymnostachydoideae was the earliest diverging lineage followed successively by Orontioideae, Lemnoideae and Lasioideae. In turn, they were followed by the clade of 'bisexual climbers' comprised of Pothoideae and Monsteroideae, which was resolved as the sister to the unisexual flowers clade of Zamioculcadoideae and Aroideae. A special WGD event ψ (psi) shared by the True-Araceae clade occurred in the Early Cretaceous. Net diversification rates first declined and then increased through time in the Araceae. The best diversification rate shift along the stem lineage of the True-Araceae clade was detected, and net diversification rates were enhanced following the ψ-WGD. Functional enrichment analyses revealed that some genes, such as those encoding heat shock proteins, glycosyl hydrolase and cytochrome P450, expanded within the True-Araceae clade.

CONCLUSIONS

Our results improve our understanding of aroid phylogeny using the large number of single-/low-copy nuclear genes. In contrast to the Proto-Araceae group and the lemnoid clade adaption to aquatic environments, our analyses of WGD, diversification and functional enrichment indicated that WGD may play a more important role in the evolution of adaptations to tropical, terrestrial environments in the True-Araceae clade. These insights provide us with new resources to interpret the evolution of the Araceae.

Inference of Ancient Polyploidy Using Transcriptome Data

Polyploidizations, or whole-genome duplications (WGDs), in plants have increased biological complexity, facilitated evolutionary innovation, and likely enabled adaptation under harsh conditions. Besides genomic data, transcriptome data have been widely employed to detect WGDs, due to their efficient accessibility to the gene space of a species. Age distributions based on synonymous substitutions (so-called K_S age distributions) for paralogs assembled from transcriptome data have identified numerous WGDs in plants, paving the way for further studies on the importance of WGDs for the evolution of seed and flowering plants. However, it is still unclear how transcriptome-based age distributions compare to those based on genomic data. In this chapter, we implemented three different de novo transcriptome assembly pipelines with two popular assemblers, i.e., Trinity and SOAPdenovo-Trans. We selected six plant species with published genomes and transcriptomes to evaluate how assembled transcripts from different pipelines perform when using K_S distributions to detect previously documented WGDs in the six species. Further, using genes predicted in each genome as references, we evaluated the effects of missing genes, gene family clustering, and de novo assembled transcripts on the transcriptome-based K_S distributions. Our results show that, although the transcriptome-based K_S distributions differ from the genome-based ones with respect to their shapes and scales, they are still reasonably reliable for unveiling WGDs, except in species where most duplicates originated from a recent WGD. We also discuss how to overcome some possible pitfalls when using transcriptome data to identify WGDs.

Revisiting ancient polyploidy in leptosporangiate ferns

Ferns, and particularly homosporous ferns, have long been assumed to have experienced recurrent whole-genome duplication (WGD) events because of their substantially large genome sizes, surprisingly high chromosome numbers, and high degrees of polyploidy among many extant members. As the number of sequenced fern genomes is limited, recent studies have employed transcriptome data to find evidence for WGDs in ferns. However, they have reached conflicting results concerning the occurrence of ancient polyploidy, for instance, in the lineage of leptosporangiate ferns. Because identifying WGDs in a phylogenetic context is the foremost step in studying the contribution of ancient polyploidy to evolution, we here revisited earlier identified WGDs in leptosporangiate ferns, mainly the core leptosporangiate ferns, by building K_S -age distributions and applying substitution rate corrections and by conducting statistical gene tree-species tree reconciliation analyses. Our integrative analyses not only identified four ancient WGDs in the sampled core leptosporangiate ferns but also identified false positives and false negatives for WGDs that recent studies have reported earlier. In conclusion, we underscore the significance of substitution rate corrections and uncertainties in gene tree-species tree reconciliations in calling WGD events and advance an exemplar workflow to overcome such often-overlooked issues.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.