Polyploidy-associated genome modifications during land plant evolution

Phylogeny and Functional Differentiation of the Terpene Synthase Gene Family in Angiosperms with Emphasis on Rosa chinensis

Terpenes are pivotal for plant growth, development, and adaptation to environmental stresses. With the advent of extensive genomic data and sophisticated bioinformatics tools, new insights into the evolutionary dynamics and functional diversification of terpene synthases (TPSs) have emerged. Despite genome-wide identifications of the TPS family in certain species, comprehensive cross-species analyses remain scarce. In this study, we conducted a genome-wide identification and subgroup classification of TPS families across 115 angiosperms with available genomic sequences. Our phylogenomic synteny network analysis elucidated the complex evolutionary history of TPS genes, revealing notable expansions and contractions among subgroups. Specifically, TPS-a showed significant expansion, while TPS-b was variably lost in some Poaceae, indicating adaptive responses. TPS-c maintained considerable conservation across species, whereas TPS-e/f diverged into distinct evolutionary trajectories despite functional overlap, with TPS-e further splitting into two angiosperm-specific clades. The TPS-g subgroup displayed lineage-restricted distribution, primarily in super-rosids and monocots. Notably, TPS-d and TPS-h subgroups were absent in angiosperms. Employing Rosa chinensis as a case study, we identified RcTPS23, a conserved bifunctional terpene synthase, highlighting the utility of cross-species synteny data in functional prediction. This comprehensive analysis elucidates the phylogenetic and functional landscape of TPS subgroups in angiosperms, providing a robust framework for predicting TPS function and guiding further functional investigations.

Identification of the BBX gene family in blueberry at different chromosome ploidy levels and fruit development and response under stress

BACKGROUND

Blueberry (Vaccinium spp.) fruits are rich in flavonoids such as anthocyanins and have a high nutritional value. The zinc finger protein transcription factor B-box (BBX) plays important roles in plant growth and development, hormone response, abiotic stress, and anthocyanin accumulation. However, studies on the BBX family in blueberry are lacking.

RESULTS

In total, 83 VcBBX and 24 VdBBX genes were identified in tetraploid and diploid blueberry, respectively. A correlation was observed between the number of BBX genes in blueberry and chromosome ploidy. Gene loss and specific replication during blueberry evolution may lead to an imbalance of quantitative relationship between VcBBX and VdBBX genes. The analysis of transcriptome and quantitative reverse transcription-polymerase chain reaction data revealed that the expression pattern of BBX genes depended on the developmental stage of blueberry fruit. Gibberellin inhibited the expression of most VcBBX genes. Abscisic acid promoted the expression of some members of the BBX family. The expression levels of VcBBX15b4, VcBBX21a1, and VcBBX30a in blueberry leaves were significantly downregulated under blue light treatment, whereas that of VcBBX15c3 was significantly upregulated under red light treatment.

CONCLUSION

In total, 83 VcBBX and 24 VdBBX genes were identified in 2 types of blueberries. Fruit development and transcription profiles under different stresses were analyzed. These findings will support further investigation of how BBX genes are involved in regulating hormone treatment and light stress during the growth and development of blueberry.

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

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Phenanthrene metabolism in Panicum miliaceum: anatomical adaptations, degradation pathway, and computational analysis of a dioxygenase enzyme

Polycyclic aromatic compounds (PAHs) are persistent organic pollutants of environmental concern due to their potential impacts on food chain, with plants being particularly vulnerable. While plants can uptake, transport, and transform PAHs, the precise mechanisms underlying their localization and degradation are not fully understood. Here, a cultivation experiment conducted with Panicum miliaceum exposed different concentrations of phenanthrene (PHE). Intermediate PHE degradation compounds were identified via GC-MS analysis, leading to the proposal of a phytodegradation pathway featuring three significant benzene ring cleavage steps. Our results showed that P. miliaceum exhibited the ability to effectively degrade high levels of PHE, resulting in the production of various intermediate products through several chemical changes. Examination of the localization and anatomical characteristics revealed structural alterations linked to PHE stress, with an observed enhancement in PHE accumulation density in both roots and shoots as treatment levels increased. Following a 2-week aging period, a decrease in the amount of PHE accumulation was observed, along with a change in its localization. Bioinformatics analysis of the P. miliaceum 2-oxoglutarate-dependent dioxygenase (2-ODD) DAO-like protein revealed a 299 amino acid structure with two highly conserved domains, namely 2OG-FeII_Oxy and DIOX_N. Molecular docking analysis aligned with experimental results, strongly affirming the potential link and direct action of 2-ODD DAO-like protein with PHE. Our study highlights P. miliaceum capacity for PAHs degradation and elucidates the mechanisms behind enhanced degradation efficiency. By integrating experimental evidence with bioinformatics analysis, we offer valuable insights into the potential applications of plant-based remediation strategies for PAHs-contaminated environments.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Trends in the evolution of intronless genes in Poaceae

Intronless genes (IGs), which are a feature of prokaryotes, are a fascinating group of genes that are also present in eukaryotes. In the current study, a comparison of Poaceae genomes revealed that the origin of IGs may have involved ancient intronic splicing, reverse transcription, and retrotranspositions. Additionally, IGs exhibit the typical features of rapid evolution, including recent duplications, variable copy numbers, low divergence between paralogs, and high non-synonymous to synonymous substitution ratios. By tracing IG families along the phylogenetic tree, we determined that the evolutionary dynamics of IGs differed among Poaceae subfamilies. IG families developed rapidly before the divergence of Pooideae and Oryzoideae and expanded slowly after the divergence. In contrast, they emerged gradually and consistently in the Chloridoideae and Panicoideae clades during evolution. Furthermore, IGs are expressed at low levels. Under relaxed selection pressure, retrotranspositions, intron loss, and gene duplications and conversions may promote the evolution of IGs. The comprehensive characterization of IGs is critical for in-depth studies on intron functions and evolution as well as for assessing the importance of introns in eukaryotes.

Genome-Wide Characterization and Evolutionary Expansion of Poplar NAC Transcription Factors and Their Tissue-Specific Expression Profiles under Drought

The NAC (NAM, ATAF1/2 and CUC2) is a large gene family of plant-specific transcription factors that play a pivotal role in various physiological processes and abiotic stresses. Due to the lack of genome-wide characterization, intraspecific and interspecific synteny, and drought-responsive expression pattern of NAC genes in poplar, the functional characterization of drought-related NAC genes have been scarcely reported in Populus species. Here, we identified a total of 170 NAC domain-containing genes in the P. trichocarpa genome, 169 of which were unevenly distributed on its nineteen chromosomes. These NAC genes were phylogenetically divided into twenty subgroups, some of which exhibited a similar pattern of exon-intron architecture. The synteny and Ka/Ks analysis indicated that the expansion of NAC genes in poplar was mainly due to gene duplication events occurring before and after the divergence of Populus and Salix. Ten PdNAC (P. deltoids × P. euramericana cv.'Nanlin895') genes were randomly selected and cloned. Their drought-responsive expression profiles showed a tissue-specific pattern. The transcription factor PdNAC013 was verified to be localized in the nucleus. Our research results provide genomic information for the expansion of NAC genes in the poplar genome, and for further characterizing putative poplar NAC genes associated with water-deficit.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Paleopolyploidies and Genomic Fractionation in Major Eudicot Clades

Eudicots account for ~75% of living angiosperms, containing important food and energy crops. Recently, high-quality genome sequences of several eudicots including Aquilegia coerulea and Nelumbo nucifera have become available, providing an opportunity to investigate the early evolutionary characteristics of eudicots. We performed genomic hierarchical and event-related alignments to infer homology within and between representative species of eudicots. The results provide strong evidence for multiple independent polyploidization events during the early diversification of eudicots, three of which are likely to be allopolyploids: The core eudicot-common hexaploidy (ECH), Nelumbo-specific tetraploidy (NST), and Ranunculales-common tetraploidy (RCT). Using different genomes as references, we constructed genomic alignment to list the orthologous and paralogous genes produced by polyploidization and speciation. This could provide a fundamental framework for studying other eudicot genomes and gene(s) evolution. Further, we revealed significantly divergent evolutionary rates among these species. By performing evolutionary rate correction, we dated RCT to be ~118-134 million years ago (Mya), after Ranunculales diverged with core eudicots at ~123-139 Mya. Moreover, we characterized genomic fractionation resulting from gene loss and retention after polyploidizations. Notably, we revealed a high degree of divergence between subgenomes. In particular, synonymous nucleotide substitutions at synonymous sites (Ks) and phylogenomic analyses implied that A. coerulea might provide the subgenome(s) for the gamma-hexaploid hybridization.

The mosaic oat genome gives insights into a uniquely healthy cereal crop

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia^1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

Assembly of the hexaploid oat genome and its diploid and tetraploid relatives clarifies the evolutionary history of oat and allows mapping of genes for agronomic traits.

In-silico evolutionary analysis of plant-OBERON proteins during compatible MYMV infection in respect of improving host resistance

Yellow mosaic disease (YMD) of pulses caused by mungbean yellow mosaic virus is a major threat to crop production. An infection that is compatible with regulating and interacting host proteins and the virus causes YMD. Oberon families of proteins OBE1-4 and VIN1-4 are imperative for plants, functions in meristem and vascular development, and were also regulated during compatible disease infection. Furthermore, in-silico expression results suggested the involvement of OBE1 and OBE2 proteins during virus infection of Vigna, Arabidopsis and soybean. Moreover, a common ancestor for the meristem and virus movement related Oberons was inferred through phylogenetic analysis. Protein interaction studies showed three amino acids (Aspartate, glutamate and lysine) in the plant homeodomain (PHD), involved in interaction with the N-terminal region of the virus movement protein and were also conserved in both monocot and dicots. Additionally, major differences in the nuclear localization signal (NLS) showing clade specific conservation and significant variation between dicots and monocots were ascertained in meristem and virus movement related Oberons. Consequently, a combination of PHD, CCD and their interactions with the VPg viral domain increases the susceptibility to YMD. Further, modification in the NLS regions of the viral movement clade Oberons, to knock out allele generation in the OBE1 and OBE2 homologs through genome-editing approaches could be established as alternate strategies for the improvement of host resistance and control yellow mosaic disease in plants, especially in pulse crops.

Role of histone deacetylase inhibitors in androgenic callus induction of Oryza sativa sub indica, in sight into evolution and mode of action of histone deacetylase genes

BACKGROUND

The potential of paddy breeding has reached its pinnacle, and hybrids have been the principal research outcome. Hence, our hypothesis was based on improvising the callus induction efficiency of recalcitrant Oryza sativa sub. indica hybrids by intervening into their cellular functions like cell division and histone regulation for the production of doubled haploids, a better output compared to hybrids.

METHODOLOGY AND RESULTS

Insight into the mechanism of cell division is the foremost concern in altering the same and hence studies on evolution, expression and action of histone deacetylase and its 12 genes (9 HDA and 3 HD-tunin genes) were chosen in the hypothesis. Expression of HDA genes at three stages (anther dehiscence, 1st callusing and second callusing stages) with inhibitor (trichostatin-A) interventions indicated 1st callusing stage as the most important in influencing callus induction and also the genes HDA19, 6, 15 and 5 were the most important. TSA alone had a significant impact on the regulation of the genes HDT 702, HDA19, HDA9, and HDA6. Higher expression of HDA19 and HDA6 was involved in maximizing callus induction; HDA15 had an antagonistic expression compared to HDA19/6 and might be involved in chlorophyll regulation during regeneration. Results of evolutionary analysis on histone deacetylases indicated a long and single lineage of origin denoting its importance in the basic cellular functions. The tubulin deacetylation gene HDA5, which was exclusively found in dicotyledons, had a recent evolutionary history only from terrestrial plants, and also had significant conservation in its motifs and NLS region.

CONCLUSION

By combating the recalcitrant nature of Indica cultivars, molecular editing on a combination of HDA genes will enhance the callus induction and regeneration efficiency of the next generation of doubled haploids, therby improving the total yield.

Comprehensive evolutionary analysis of the TCP gene family: Further insights for its origin, expansion, and diversification

TCP proteins are plant-specific transcription factors, which are involved in a broad range of physiological processes of plant growth and development. However, the origin and evolutionary history of this gene family is not fully resolved. Here, we present a genome-wide survey of TCP genes in 59 species (including 42 genomes and 17 transcriptomes) covering all main lineages of green plants, and reconstruct the evolutionary history of this gene family. Our results suggested that the origin of TCP genes predated the emergence of land plants, possibly in the common ancestor of Phragmoplastophyta. The TCP gene family gradually experienced a continuous expansion and grew from a few members in algae, moss and lycophytes to dozens, and sometimes over 50 members in angiosperms. Phylogenetic analysis indicated that at least four subclades (Class I and three subclades of Class II) have been occurred in the ancestor of spermatophyte (seed plant). Both dispersed duplication and segmental duplication or whole-genome duplication (WGD) contributed significantly to the expansion of the TCP gene family over the course of evolution. Our findings provide a comprehensive evolutionary analysis of the TCP gene family and highlight the importance of gene duplications in the evolution of this plant-specific transcription factors.

The Chloranthus sessilifolius genome provides insight into early diversification of angiosperms

Most extant angiosperms belong to Mesangiospermae, which comprises eudicots, monocots, magnoliids, Chloranthales and Ceratophyllales. However, phylogenetic relationships between these five lineages remain unclear. Here, we report the high-quality genome of a member of the Chloranthales lineage (Chloranthus sessilifolius). We detect only one whole genome duplication within this species and find that polyploidization events in different Mesangiospermae lineage are mutually independent. We also find that the members of all floral development-related gene lineages are present in C. sessilifolius despite its extremely simplified flower. The AP1 and PI genes, however, show a weak floral tissue-specialized expression. Our phylogenomic analyses suggest that Chloranthales and magnoliids are sister groups, and both are together sister to the clade comprising Ceratophyllales and eudicots, while the monocot lineage is sister to all other Mesangiospermae. Our findings suggest that in addition to hybridization, incomplete lineage sorting may largely account for phylogenetic inconsistencies between the observed gene trees.

An Evolutionary Analysis of B-Box Transcription Factors in Strawberry Reveals the Role of FaBBx28c1 in the Regulation of Flowering Time

Flowering connects vegetative and generative developmental phases and plays a significant role in strawberry production. The mechanisms that regulate strawberry flowering time are unclear. B-box transcription factors (BBXs) play important roles in the flowering time regulation of plants. Nevertheless, BBXs in octoploid cultivated strawberry (Fragaria ananassa) and their functions in flowering time regulation have not been identified. Here, we identified 51 FaBBXs from cultivated strawberry and 16 FvBBXs from diploid wild strawberry (Fragaria vesca), which were classified into five groups according to phylogenetic analysis. Further evolutionary analysis showed that whole-genome duplication or segmental duplication is a crucial factor that leads to the expansion of the BBX gene family in two strawberry species. Moreover, some loss and acquisition events of FaBBX genes were identified in the genome of cultivated strawberry that could have affected traits of agronomic interest, such as fruit quality. The promoters of FaBBX genes showed an enrichment in light-responsive, cis-regulatory elements, with 16 of these genes showing changes in their transcriptional activity in response to blue light treatment. On the other hand, FaBBX28c1, whose transcriptional activity is reduced in response to blue light, showed a delay in flowering time in Arabidopsis transgenic lines, suggesting its role in the regulation of flowering time in cultivated strawberry. Our results provide new evolutionary insight into the BBX gene family in cultivated strawberry and clues regarding their function in flowering time regulation in plants.

Insights into angiosperm evolution, floral development and chemical biosynthesis from the Aristolochia fimbriata genome

Aristolochia, a genus in the magnoliid order Piperales, has been famous for centuries for its highly specialized flowers and wide medicinal applications. Here, we present a new, high-quality genome sequence of Aristolochia fimbriata, a species that, similar to Amborella trichopoda, lacks further whole-genome duplications since the origin of extant angiosperms. As such, the A. fimbriata genome is an excellent reference for inferences of angiosperm genome evolution, enabling detection of two novel whole-genome duplications in Piperales and dating of previously reported whole-genome duplications in other magnoliids. Genomic comparisons between A. fimbriata and other angiosperms facilitated the identification of ancient genomic rearrangements suggesting the placement of magnoliids as sister to monocots, whereas phylogenetic inferences based on sequence data we compiled yielded ambiguous relationships. By identifying associated homologues and investigating their evolutionary histories and expression patterns, we revealed highly conserved floral developmental genes and their distinct downstream regulatory network that may contribute to the complex flower morphology in A. fimbriata. Finally, we elucidated the genetic basis underlying the biosynthesis of terpenoids and aristolochic acids in A. fimbriata.

Evolution of the MLO gene families in octoploid strawberry (Fragaria ×ananassa) and progenitor diploid species identified potential genes for strawberry powdery mildew resistance

Powdery mildew (PM) caused by Podosphaera aphanis is a major fungal disease of cultivated strawberry. Mildew Resistance Locus O (MLO) is a gene family described for having conserved seven-transmembrane domains. Induced loss-of-function in specific MLO genes can confer durable and broad resistance against PM pathogens. However, the genomic structure and potential role of MLO genes for PM resistance have not been characterized yet in the octoploid cultivated strawberry. In the present study, MLO gene families were characterized in four diploid progenitor species (Fragaria vesca, F. iinumae, F. viridis, and F. nipponica) and octoploid cultivated (Fragaria ×ananassa) strawberry, and potential sources of MLO-mediated susceptibility were identified. Twenty MLO sequences were identified in F. vesca and 68 identified in F. ×ananassa. Phylogenetic analysis divided diploid and octoploid strawberry MLO genes into eight different clades, in which three FveMLO (MLO10, MLO17, and MLO20) and their twelve orthologs of FaMLO were grouped together with functionally characterized MLO genes conferring PM susceptibility. Copy number variations revealed differences in MLO composition among homoeologous chromosomes, supporting the distinct origin of each subgenome during the evolution of octoploid strawberry. Dissecting genomic sequence and structural variations in candidate FaMLO genes revealed their potential role associated with genetic controls and functionality in strawberry against PM pathogen. Furthermore, the gene expression profiling and RNAi silencing of putative FaMLO genes in response to the pathogen indicate the function in PM resistance. These results are a critical first step in understanding the function of strawberry MLO genes and will facilitate further genetic studies of PM resistance in cultivated strawberry.

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

From time immemorial, humans have exploited plants as a source of food and medicines. The World Health Organization (WHO) has recorded 21,000 plants with medicinal value out of 300,000 species available worldwide. The promising modern "multi-omics" platforms and tools have been proven as functional platforms able to endow us with comprehensive knowledge of the proteome, genome, transcriptome, and metabolome of medicinal plant systems so as to reveal the novel connected genetic (gene) pathways, proteins, regulator sequences and secondary metabolite (molecule) biosynthetic pathways of various drug and protein molecules from a variety of plants with therapeutic significance. This review paper endeavors to abridge the contemporary advancements in research areas of multi-omics and the information involved in decoding its prospective relevance to the utilization of plants with medicinal value in the present global scenario. The crosstalk of medicinal plants with genomics, transcriptomics, proteomics, and metabolomics approaches will be discussed.

Genome-wide analysis of the AREB/ABF gene lineage in land plants and functional analysis of TaABF3 in Arabidopsis

BACKGROUND

Previous studies have shown that ABFs (abscisic acid-responsive transcription factors) are important ABA-signaling components that participate in abiotic stress response. However, little is known about the function of ABFs in Triticum aestivum. In addition, although various ABFs have been identified in other species, the phylogenetic relationship between ABF transcription factors has not been systemically investigated in land plants.

RESULTS

In this study, we systemically collected ABFs from land plants and analyzed the phylogenetic relationship of these ABF genes. The ABF genes are present in all the land plants we investigated, including moss, lycophyte, monocots, and eudicots. Furthermore, these ABF genes are phylogenetically divided into seven subgroups, differentiations that are supported by variation in the gene structure, protein properties, and motif patterns. We further demonstrated that the expression of ABF genes varies among different tissues and developmental stages, and are induced by one or more environmental stresses. Furthermore, we found that three wheat ABFs (TaABF1, TaABF2, and TaABF3) were significantly induced by drought stress. Compared with wild-type (WT) plants, transgenic Arabidopsis plants overexpressing TaABF3 displayed enhanced drought tolerance.

CONCLUSIONS

These results provide important ground work for understanding the phylogenetic relationships between plant ABF genes. Our results also indicate that TaABFs may participate in regulating plant response to abiotic stresses.

Distinct nucleotide patterns among three subgenomes of bread wheat and their potential origins during domestication after allopolyploidization

BACKGROUND

The speciation and fast global domestication of bread wheat have made a great impact on three subgenomes of bread wheat. DNA base composition is an essential genome feature, which follows the individual-strand base equality rule and [AT]-increase pattern at the genome, chromosome, and polymorphic site levels among thousands of species. Systematic analyses on base compositions of bread wheat and its wild progenitors could facilitate further understanding of the evolutionary pattern of genome/subgenome-wide base composition of allopolyploid species and its potential causes.

RESULTS

Genome/subgenome-wide base-composition patterns were investigated by using the data of polymorphic site in 93 accessions from worldwide populations of bread wheat, its diploid and tetraploid progenitors, and their corresponding reference genome sequences. Individual-strand base equality rule and [AT]-increase pattern remain in recently formed hexaploid species bread wheat at the genome, subgenome, chromosome, and polymorphic site levels. However, D subgenome showed the fastest [AT]-increase across polymorphic site from Aegilops tauschii to bread wheat than that on A and B subgenomes from wild emmer to bread wheat. The fastest [AT]-increase could be detected almost all chromosome windows on D subgenome, suggesting different mechanisms between D and other two subgenomes. Interestingly, the [AT]-increase is mainly contributed by intergenic regions at non-selective sweeps, especially the fastest [AT]-increase of D subgenome. Further transition frequency and sequence context analysis indicated that three subgenomes shared same mutation type, but D subgenome owns the highest mutation rate on high-frequency mutation type. The highest mutation rate on D subgenome was further confirmed by using a bread-wheat-private SNP set. The exploration of loci/genes related to the [AT] value of D subgenome suggests the fastest [AT]-increase of D subgenome could be involved in DNA repair systems distributed on three subgenomes of bread wheat.

CONCLUSIONS

The highest mutation rate is detected on D subgenome of bread wheat during domestication after allopolyploidization, leading to the fastest [AT]-increase pattern of D subgenome. The phenomenon may come from the joint action of multiple repair systems inherited from its wild progenitors.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

Phylogenomic analysis of the APETALA2 transcription factor subfamily across angiosperms reveals both deep conservation and lineage-specific patterns

The APETALA2 (AP2) subfamily of transcription factors are key regulators of angiosperm root, shoot, flower and embryo development. The broad diversity of anatomical and morphological structures is potentially associated with the genomic dynamics of the AP2 subfamily. However, a comprehensive phylogenomic analysis of the AP2 subfamily across angiosperms is lacking. We combined phylogenetic and synteny analysis of distinct AP2 subclades in the completed genomes of 107 angiosperm species. We identified major changes in copy number variation and genomic context within subclades across lineages, and discuss how these changes may have contributed to the evolution of lineage-specific traits. Multiple AP2 subclades show highly conserved patterns of copy number and synteny across angiosperms, while others are more dynamic and show distinct lineage-specific patterns. As examples of lineage-specific morphological divergence due to AP2 subclade dynamics, we hypothesize that loss of PLETHORA1/2 in monocots correlates with the absence of taproots, whereas independent lineage-specific changes of PLETHORA4/BABY BOOM and WRINKLED1 genes in Brassicaceae and monocots point towards regulatory divergence of embryogenesis between these lineages. Additionally, copy number expansion of TOE1 and TOE3/AP2 in asterids is implicated with differential regulation of flower development. Moreover, we show that the genomic context of AP2s is in general highly specialized per angiosperm lineage. To our knowledge, this study is the first to shed light on the evolutionary divergence of the AP2 subfamily subclades across major angiosperm lineages and emphasizes the need for lineage-specific characterization of developmental networks to understand trait variability further.

Global Role of Crop Genomics in the Face of Climate Change

The development of climate change resilient crops is necessary if we are to meet the challenge of feeding the growing world's population. We must be able to increase food production despite the projected decrease in arable land and unpredictable environmental conditions. This review summarizes the technological and conceptual advances that have the potential to transform plant breeding, help overcome the challenges of climate change, and initiate the next plant breeding revolution. Recent developments in genomics in combination with high-throughput and precision phenotyping facilitate the identification of genes controlling critical agronomic traits. The discovery of these genes can now be paired with genome editing techniques to rapidly develop climate change resilient crops, including plants with better biotic and abiotic stress tolerance and enhanced nutritional value. Utilizing the genetic potential of crop wild relatives (CWRs) enables the domestication of new species and the generation of synthetic polyploids. The high-quality crop plant genome assemblies and annotations provide new, exciting research targets, including long non-coding RNAs (lncRNAs) and cis-regulatory regions. Metagenomic studies give insights into plant-microbiome interactions and guide selection of optimal soils for plant cultivation. Together, all these advances will allow breeders to produce improved, resilient crops in relatively short timeframes meeting the demands of the growing population and changing climate.

Balanced Genome Triplication in Wheat Causes Premature Growth Arrest and an Upheaval of Genome-Wide Gene Regulation

Polyploidy, or whole genome duplication (WGD), is a driving evolutionary force across the tree of life and has played a pervasive role in the evolution of the plant kingdom. It is generally believed that a major genetic attribute contributing to the success of polyploidy is increased gene and genome dosage. The evolution of polyploid wheat has lent support to this scenario. Wheat has evolved at three ploidal levels: diploidy, tetraploidy, and hexaploidy. Ample evidence testifies that the evolutionary success, be it with respect to evolvability, natural adaptability, or domestication has dramatically increased with each elevation of the ploidal levels. A long-standing question is what would be the outcome if a further elevation of ploidy is superimposed on hexaploid wheat? Here, we characterized a spontaneously occurring nonaploid wheat individual in selfed progenies of synthetic hexaploid wheat and compared it with its isogenic hexaploid siblings at the phenotypic, cytological, and genome-wide gene-expression levels. The nonaploid manifested severe defects in growth and development, albeit with a balanced triplication of the three wheat subgenomes. Transcriptomic profiling of the second leaf of nonaploid, taken at a stage when phenotypic abnormality was not yet discernible, already revealed significant dysregulation in global-scale gene expression with ca. 25.2% of the 49,436 expressed genes being differentially expressed genes (DEGs) at a twofold change cutoff relative to the hexaploid counterpart. Both up- and downregulated DEGs were identified in the nonaploid vs. hexaploid, including 457 genes showing qualitative alteration, i.e., silencing or activation. Impaired functionality at both cellular and organismal levels was inferred from gene ontology analysis of the DEGs. Homoeologous expression analysis of 9,574 sets of syntenic triads indicated that, compared with hexaploid, the proportions showing various homeologous expression patterns were highly conserved in the nonaploid although gene identity showed moderate reshuffling among some of the patterns in the nonaploid. Together, our results suggest hexaploidy is likely the upper limit of ploidy level in wheat; crossing this threshold incurs severe ploidy syndrome that is preceded by disruptive dysregulation of global gene expression.

A cotton NAC transcription factor GhirNAC2 plays positive roles in drought tolerance via regulating ABA biosynthesis

NAC protein is a large plant specific transcription factor family, which plays important roles in the response to abiotic stresses. However, the regulation mechanism of most NAC proteins in drought stress remains to be further uncovered. In this study, we elucidated the molecular functions of a NAC protein, GhirNAC2, in response to drought stress in cotton. GhirNAC2 was greatly induced by drought and phytohormone abscisic acid (ABA). Subcellular localization demonstrated that GhirNAC2 was located in the nucleus. Co-suppression of GhirNAC2 in cotton led to larger stomata aperture, elevated water loss and finally reduced transgenic plants tolerance to drought stress. Furthermore, the endogenous ABA content was significantly lower in GhirNAC2-suppressed transgenic plant leaves compared to wild type. in vivo and in vitro studies showed that GhirNAC2 directly binds to the promoter of GhNCED3a/3c, key genes in ABA biosynthesis, which were both down-regulated in GhirNAC2-suppressed transgenic lines. Transient silencing of GhNCED3a/3c also significantly reduced the resistance to drought stress in cotton plants. However, ectopic expression of GhirNAC2 in tobacco significantly enhanced seed germination, root growth and plant survival under drought stress. Taken together, GhirNAC2 plays a positive role in cotton drought tolerance, which functions by modulating ABA biosynthesis and stomata closure via regulating GhNCED3a/3c expression.

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

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

Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization

Cold stress profoundly affects plant growth and development and is a key factor affecting the geographic distribution and evolution of plants. Plants have evolved adaptive mechanisms to cope with cold stress. Here, through the genomic analysis of Arabidopsis, three Brassica species and 17 other representative species, we found that both cold-related genes (CRGs) and their collinearity were preferentially retained after polyploidization followed by genome instability, while genome-wide gene sets exhibited a variety of other expansion mechanisms. The cold-related regulatory network was increased in Brassicaceae genomes, which were recursively affected by polyploidization. By combining our findings regarding the selective retention of CRGs from this ecological genomics study with the available knowledge of cold-induced chromosome doubling, we hypothesize that cold stress may have contributed to the success of polyploid plants through both increasing polyploidization and selectively maintaining CRGs during evolution. This hypothesis requires further biological and ecological exploration to obtain solid supporting evidence, which will potentially contribute to understanding the generation of polyploids and to the field of ecological genomics.

Pangenomics Comes of Age: From Bacteria to Plant and Animal Applications

The pangenome refers to a collection of genomic sequence found in the entire species or population rather than in a single individual; the sequence can be core, present in all individuals, or accessory (variable or dispensable), found in a subset of individuals only. While pangenomic studies were first undertaken in bacterial species, developments in genome sequencing and assembly approaches have allowed construction of pangenomes for eukaryotic organisms, fungi, plants, and animals, including two large-scale human pangenome projects. Analysis of the these pangenomes revealed key differences, most likely stemming from divergent evolutionary histories, but also surprising similarities.

The Aquilegia genome reveals a hybrid origin of core eudicots

Background

Whole-genome duplications (WGDs) have dominated the evolutionary history of plants. One consequence of WGD is a dramatic restructuring of the genome as it undergoes diploidization, a process under which deletions and rearrangements of various sizes scramble the genetic material, leading to a repacking of the genome and eventual return to diploidy. Here, we investigate the history of WGD in the columbine genus Aquilegia, a basal eudicot, and use it to illuminate the origins of the core eudicots.

Results

Within-genome synteny confirms that columbines are ancient tetraploids, and comparison with the grape genome reveals that this tetraploidy appears to be shared with the core eudicots. Thus, the ancient gamma hexaploidy found in all core eudicots must have involved a two-step process: first, tetraploidy in the ancestry of all eudicots, then hexaploidy in the ancestry of core eudicots. Furthermore, the precise pattern of synteny sharing suggests that the latter involved allopolyploidization and that core eudicots thus have a hybrid origin.

Conclusions

Novel analyses of synteny sharing together with the well-preserved structure of the columbine genome reveal that the gamma hexaploidy at the root of core eudicots is likely a result of hybridization between a tetraploid and a diploid species.

The allotetraploid origin and asymmetrical genome evolution of the common carp Cyprinus carpio

Common carp (Cyprinus carpio) is an allotetraploid species derived from recent whole genome duplication and provides a model to study polyploid genome evolution in vertebrates. Here, we generate three chromosome-level reference genomes of C. carpio and compare to related diploid Cyprinid genomes. We identify a Barbinae lineage as potential diploid progenitor of C. carpio and then divide the allotetraploid genome into two subgenomes marked by a distinct genome similarity to the diploid progenitor. We estimate that the two diploid progenitors diverged around 23 Mya and merged around 12.4 Mya based on the divergence rates of homoeologous genes and transposable elements in two subgenomes. No extensive gene losses are observed in either subgenome. Instead, we find gene expression bias across surveyed tissues such that subgenome B is more dominant in homoeologous expression. CG methylation in promoter regions may play an important role in altering gene expression in allotetraploid C. carpio.

The common carp is derived from recent whole genome duplication and represents a model for polyploid genome evolution, rare in vertebrates. Here, the authors generate and analyse chromosome-level reference genomes for common carp, and describe subgenome gene expression changes.

A draft genome and transcriptome of common milkweed (Asclepias syriaca) as resources for evolutionary, ecological, and molecular studies in milkweeds and Apocynaceae

Milkweeds (Asclepias) are used in wide-ranging studies including floral development, pollination biology, plant-insect interactions and co-evolution, secondary metabolite chemistry, and rapid diversification. We present a transcriptome and draft nuclear genome assembly of the common milkweed, Asclepias syriaca. This reconstruction of the nuclear genome is augmented by linkage group information, adding to existing chloroplast and mitochondrial genomic resources for this member of the Apocynaceae subfamily Asclepiadoideae. The genome was sequenced to 80.4× depth and the draft assembly contains 54,266 scaffolds ≥1 kbp, with N50 = 3,415 bp, representing 37% (156.6 Mbp) of the estimated 420 Mbp genome. A total of 14,474 protein-coding genes were identified based on transcript evidence, closely related proteins, and ab initio models, and 95% of genes were annotated. A large proportion of gene space is represented in the assembly, with 96.7% of Asclepias transcripts, 88.4% of transcripts from the related genus Calotropis, and 90.6% of proteins from Coffea mapping to the assembly. Scaffolds covering 75 Mbp of the Asclepias assembly formed 11 linkage groups. Comparisons of these groups with pseudochromosomes in Coffea found that six chromosomes show consistent stability in gene content, while one may have a long history of fragmentation and rearrangement. The progesterone 5β-reductase gene family, a key component of cardenolide production, is likely reduced in Asclepias relative to other Apocynaceae. The genome and transcriptome of common milkweed provide a rich resource for future studies of the ecology and evolution of a charismatic plant family.

Next-generation transcriptome assembly and analysis: Impact of ploidy

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

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

Myc-like transcriptional factors in wheat: structural and functional organization of the subfamily I members

BACKGROUND

Myc-like regulatory factors carrying the basic helix-loop-helix (bHLH) domain belong to a large superfamily of transcriptional factors (TFs) present in all eukaryotic kingdoms. In plants, the representatives of this superfamily regulate diverse biological processes including growth and development as well as response to various stresses. As members of the regulatory MBW complexes, they participate in biosynthesis of flavonoids. In wheat, only one member (TaMyc1) of the Myc-like TFs family has been studied, while structural and functional organization of further members remained uncharacterized. From two Myc-subfamilies described recently in the genomes of Triticeae tribe species, we investigated thoroughly the members of the subfamily I which includes the TaMyc1 gene.

RESULTS

Comparison of the promoter regions of the Myc subfamily I members in wheat suggested their division into two groups (likely homoeologous sets): TaMyc-1 (TaMyc-A1/TaMyc1, TaMyc-B1, TaMyc-D1) and TaMyc-2 (TaMyc-A2 and TaMyc-D2). It was demonstrated that the TaMyc-D1 copy has lost its functionality due to the frame shift mutation. The study of functional features of the other four copies suggested some of them to be involved in the biosynthesis of anthocyanins. In particular, TaMyc-B1 is assumed to be a co-regulator of the gene TaC1-A1 (encoding R2R3-Myb factor) in the MBW regulatory complex activating anthocyanin synthesis in wheat coleoptile. The mRNA levels of the TaMyc-A1, TaMyc-B1, TaMyc-A2 and TaMyc-D2 genes increased significantly in wheat seedlings exposed to osmotic stress. Salinity stress induced expression of TaMyc-B1 and TaMyc-A2, while TaMyc-A1 was repressed.

CONCLUSIONS

The features of the structural and functional organization of the members of subfamily I of Myc-like TFs in wheat were determined. Myc-like co-regulator (TaMyc-B1) of anthocyanin synthesis in wheat coleoptile was described for the first time. The Myc-encoding genes presumably involved in response to drought and salinity were determined in wheat. The results obtained are important for further manipulations with Myc genes, aimed on increasing wheat adaptability.

Floral MADS-box protein interactions in the early diverging angiosperm Aristolochia fimbriata Cham. (Aristolochiaceae: Piperales)

Floral identity MADS-box A, B, C, D, E, and AGL6 class genes are predominantly single copy in Magnoliids, and predate the whole genome duplication (WGD) events in monocots and eudicots. By comparison with the model species Arabidopsis thaliana, the expression patterns of B-, C-, and D-class genes in stamen, carpel, and ovules are conserved in Aristolochia fimbriata, whereas A-, E-class, and AGL6 genes have different expression patterns. Nevertheless, the interactions of these proteins that act through multimeric complexes remain poorly known in early divergent angiosperms. This study evaluates protein interactions among all floral MADS-box A. fimbriata proteins using the Yeast Two Hybrid System (Y2H). We found no homodimers and less heterodimers formed by AfimFUL when compared to AfimAGL6, which allowed us to suggest AGL6 homodimers in combination with AfimSEP2 as the most likely tetramer in sepal identity. We found AfimAP3-AfimPI obligate heterodimers and AfimAG-AfimSEP2 protein interactions intact suggesting conserved stamen and carpel tetrameric complexes in A. fimbriata. We observed a broader interaction partner set for AfimSEP2 than for its paralog AfimSEP1. We show conserved and exclusive MADS-box protein interactions in A. fimbriata in comparison with other eudicot and monocot model species in order to establish plesiomorphic MADS-box protein floral networks in angiosperms.

Divergent gene expression levels between diploid and autotetraploid Tolmiea relative to the total transcriptome, the cell, and biomass

PREMISE OF THE STUDY

Studies of gene expression and polyploidy are typically restricted to characterizing differences in transcript concentration. Using diploid and autotetraploid Tolmiea, we present an integrated approach for cross-ploidy comparisons that account for differences in transcriptome size and cell density and make multiple comparisons of transcript abundance.

METHODS

We use RNA spike-in standards in concert with cell size and density to identify and correct for differences in transcriptome size and compare levels of gene expression across multiple scales: per transcriptome, per cell, and per biomass.

KEY RESULTS

In total, ~17% of all loci were identified as differentially expressed (DEGs) between the diploid and autopolyploid species. The per-transcriptome normalization, the method researchers typically use, captured the fewest DEGs (58% of total DEGs) and failed to detect any DEGs not found by the alternative normalizations. When transcript abundance was normalized per biomass and per cell, ~66% and ~82% of the total DEGs were recovered, respectively. The discrepancy between per-transcriptome and per-cell recovery of DEGs occurs because per-transcriptome normalizations are concentration-based and therefore blind to differences in transcriptome size.

CONCLUSIONS

While each normalization enables valid comparisons at biologically relevant scales, a holistic comparison of multiple normalizations provides additional explanatory power not available from any single approach. Notably, autotetraploid loci tend to conserve diploid-like transcript abundance per biomass through increased gene expression per cell, and these loci are enriched for photosynthesis-related functions.

GenoDup Pipeline: a tool to detect genome duplication using the dS-based method

Understanding whole genome duplication (WGD), or polyploidy, is fundamental to investigating the origin and diversification of organisms in evolutionary biology. The wealth of genomic data generated by next generation sequencing (NGS) has resulted in an urgent need for handy and accurate tools to detect WGD. Here, I present a useful and user-friendly pipeline called GenoDup for inferring WGD using the dS-based method. I have successfully applied GenoDup to identify WGD in empirical data from both plants and animals. The GenoDup Pipeline provides a reliable and useful tool to infer WGD from NGS data.

Unbiased subgenome evolution following a recent whole-genome duplication in pear (Pyrus bretschneideri Rehd.)

Genome fractionation (also known as diploidization) frequently occurs following paleopolyploidization events. Biased fractionation between subgenomes has been found in some paleo-allopolyploids, while this phenomenon is absent in paleo-autopolyploids. Pear (Pyrus bretschneideri Rehd.) experienced a recent whole-genome duplication (WGD, ~30 million years ago); however, the evolutionary fate of the two subgenomes derived from this WGD event is not clear. In this study, we identified the two paleo-subgenomes in pear using peach (Prunus persica) as an outgroup and investigated differences in the gene loss rate, evolutionary rate, gene expression level, and DNA methylation level between these two subgenomes. Fractionation bias was not found between the two pear subgenomes, which evolved at similar evolutionary rates. The DNA methylation level of the two subgenomes showed little bias, and we found no expression dominance between the subgenomes. However, we found that singleton genes and homeologous genes within each subgenome showed divergent evolutionary patterns of selective constraints, expression and epigenetic modification. These results provide insights into subgenome evolution following paleopolyploidization in pear.

Comparative Analysis of Homologous Sequences of Saccharum officinarum and Saccharum spontaneum Reveals Independent Polyploidization Events

Sugarcane (Saccharum spp. hybrids) is an economically important crop widely grown in tropical and subtropical regions for sugar and ethanol production. However, the large genome size, high ploidy level, interspecific hybridization and aneuploidy make sugarcane one of the most complex genomes and have long hampered genome research in sugarcane. Modern sugarcane cultivars are derived from interspecific hybridization between S. officinarum and S. spontaneum with 80-90% of the genome from S. officinarum and 10-20% of the genome from S. spontaneum. We constructed bacterial artificial chromosome (BAC) libraries of S. officinarum variety LA Purple (2n = 8x = 80) and S. spontaneum haploid clone AP85-441 (2n = 4x = 32), and selected and sequenced 97 BAC clones from the two Saccharum BAC libraries. A total of 5,847,280 bp sequence from S. officinarum and 5,011,570 bp from S. spontaneum were assembled and 749 gene models were annotated in these BACs. A relatively higher gene density and lower repeat content were observed in S. spontaneum BACs than in S. officinarum BACs. Comparative analysis of syntenic regions revealed a high degree of collinearity in genic regions between Saccharum and Sorghum bicolor and between S. officinarum and S. spontaneum. In the syntenic regions, S. spontaneum showed expansion relative to S. officinarum, and both S. officinarum and S. spontaneum showed expansion relative to sorghum. Among the 75 full-length LTR retrotransposons identified in the Saccharum BACs, none of them are older than 2.6 mys and no full-length LTR elements are shared between S. officinarum and S. spontaneum. In addition, divergence time estimated using a LTR junction marker and a syntenic gene shared by 3 S. officinarum and 1 S. spontaneum BACs revealed that the S. spontaneum intergenic region was distant to those from the 3 homologous regions in S. officinarum. Our results suggested that S. officinarum and S. spontaneum experienced at least two rounds of independent polyploidization in each lineage after their divergence from a common ancestor.

Functional Diversity of Diterpene Synthases in the Biofuel Crop Switchgrass

Diterpenoids constitute a diverse class of metabolites with critical functions in plant development, defense, and ecological adaptation. Major monocot crops, such as maize (Zea mays) and rice (Oryza sativa), deploy diverse blends of specialized diterpenoids as core components of biotic and abiotic stress resilience. Here, we describe the genome-wide identification and functional characterization of stress-related diterpene synthases (diTPSs) in the dedicated bioenergy crop switchgrass (Panicum virgatum). Mining of the allotetraploid switchgrass genome identified an expansive diTPS family of 31 members, and biochemical analysis of 11 diTPSs revealed a modular metabolic network producing a diverse array of diterpenoid metabolites. In addition to ent-copalyl diphosphate (CPP) and ent-kaurene synthases predictably involved in gibberellin biosynthesis, we identified syn-CPP and ent-labda-13-en-8-ol diphosphate (LPP) synthases as well as two diTPSs forming (+)-labda-8,13E-dienyl diphosphate (8,13-CPP) and ent-neo-cis-trans-clerodienyl diphosphate (CT-CLPP) scaffolds not observed previously in plants. Structure-guided mutagenesis of the (+)-8,13-CPP and ent-neo-CT-CLPP synthases revealed residue substitutions in the active sites that altered product outcome, representing potential neofunctionalization events that occurred during diversification of the switchgrass diTPS family. The conversion of ent-CPP, ent-LPP, syn-CPP, and ent-neo-CT-CLPP by promiscuous diTPSs further yielded distinct labdane-type diterpene olefins and alcohols. Of these metabolites, the formation of 9β-hydroxy-syn-pimar-15-ene and the expression of the corresponding genes were induced in roots and leaves in response to oxidative stress and ultraviolet irradiation, indicating their possible roles in abiotic stress adaptation. Together, these findings expand the known chemical space of diterpenoid metabolism in monocot crops toward systematically investigating and ultimately improving stress resilience traits in crop species.

Evolution of floral diversity: genomics, genes and gamma

A salient feature of flowering plant diversification is the emergence of a novel suite of floral features coinciding with the origin of the most species-rich lineage, Pentapetalae. Advances in phylogenetics, developmental genetics and genomics, including new analyses presented here, are helping to reconstruct the specific evolutionary steps involved in the evolution of this clade. The enormous floral diversity among Pentapetalae appears to be built on a highly conserved ground plan of five-parted (pentamerous) flowers with whorled phyllotaxis. By contrast, lability in the number and arrangement of component parts of the flower characterize the early-diverging eudicot lineages subtending Pentapetalae. The diversification of Pentapetalae also coincides closely with ancient hexaploidy, referred to as the gamma whole-genome triplication, for which the phylogenetic timing, mechanistic details and molecular evolutionary consequences are as yet not fully resolved. Transcription factors regulating floral development often persist in duplicate or triplicate in gamma-derived genomes, and both individual genes and whole transcriptional programmes exhibit a shift from broadly overlapping to tightly defined expression domains in Pentapetalae flowers. Investigations of these changes associated with the origin of Pentapetalae can lead to a more comprehensive understanding of what is arguably one of the most important evolutionary diversification events within terrestrial plants.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation

Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.

Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch

Silver birch (Betula pendula) is a pioneer boreal tree that can be induced to flower within 1 year. Its rapid life cycle, small (440-Mb) genome, and advanced germplasm resources make birch an attractive model for forest biotechnology. We assembled and chromosomally anchored the nuclear genome of an inbred B. pendula individual. Gene duplicates from the paleohexaploid event were enriched for transcriptional regulation, whereas tandem duplicates were overrepresented by environmental responses. Population resequencing of 80 individuals showed effective population size crashes at major points of climatic upheaval. Selective sweeps were enriched among polyploid duplicates encoding key developmental and physiological triggering functions, suggesting that local adaptation has tuned the timing of and cross-talk between fundamental plant processes. Variation around the tightly-linked light response genes PHYC and FRS10 correlated with latitude and longitude and temperature, and with precipitation for PHYC. Similar associations characterized the growth-promoting cytokinin response regulator ARR1, and the wood development genes KAK and MED5A.

Revisiting the phosphatidylethanolamine-binding protein (PEBP) gene family reveals cryptic FLOWERING LOCUS T gene homologs in gymnosperms and sheds new light on functional evolution

Angiosperms and gymnosperms are two major groups of extant seed plants. It has been suggested that gymnosperms lack FLOWERING LOCUS T (FT), a key integrator at the core of flowering pathways in angiosperms. Taking advantage of newly released gymnosperm genomes, we revisited the evolutionary history of the plant phosphatidylethanolamine-binding protein (PEBP) gene family through phylogenetic reconstruction. Expression patterns in three gymnosperm taxa and heterologous expression in Arabidopsis were studied to investigate the functions of gymnosperm FT-like and TERMINAL FLOWER 1 (TFL1)-like genes. Phylogenetic reconstruction suggests that an ancient gene duplication predating the divergence of seed plants gave rise to the FT and TFL1 genes. Expression patterns indicate that gymnosperm TFL1-like genes play a role in the reproductive development process, while GymFT1 and GymFT2, the FT-like genes resulting from a duplication event in the common ancestor of gymnosperms, function in both growth rhythm and sexual development pathways. When expressed in Arabidopsis, both spruce FT-like and TFL1-like genes repressed flowering. Our study demonstrates that gymnosperms do have FT-like and TFL1-like genes. Frequent gene and genome duplications contributed significantly to the expansion of the plant PEBP gene family. The expression patterns of gymnosperm PEBP genes provide novel insight into the functional evolution of this gene family.

Deep into the Aristolochia Flower: Expression of C, D, and E-Class Genes in Aristolochia fimbriata (Aristolochiaceae)

Aristolochia fimbriata (Aristolochiaceae) is a member of an early diverging lineage of flowering plants and a promising candidate for evo-devo studies. Aristolochia flowers exhibit a unique floral synorganization that consists of a monosymmetric and petaloid calyx formed by three congenitally fused sepals, and a gynostemium formed by the congenital fusion between stamens and the stigmatic region of the carpels. This floral ground plan atypical in the magnoliids can be used to evaluate the role of floral organ identity MADS-box genes during early flower evolution. In this study, we present in situ hybridization experiments for the homologs of the canonical C-, D-, and E-class genes. Spatiotemporal expression of the C-class gene AfimAG is restricted to stamens, ovary, and ovules, suggesting a conserved stamen and carpel identity function, consistent with that reported in core-eudicots and monocots. The D-class gene AfimSTK is detected in the anthers, the stigmas, the ovary, the ovules, the fruit, and the seeds, suggesting conserved roles in ovule and seed identity and unique roles in stamens, ovary, and fruit development. In addition, AfimSTK expression patterns in areas of organ abscission and dehiscence zones suggest putative roles linked to senescence processes. We found that both E-class genes are expressed in the anthers and the ovary; however, AfimSEP2 exhibits higher expression compared to AfimSEP1. These findings provide a comprehensive picture of the ancestral expression patterns of the canonical MADS-box floral organ identity genes and the foundations for further comparative analyses in other magnoliids.