The pineapple genome and the evolution of CAM photosynthesis

The highly continuous reference genome of a leaf-chimeric red pineapple (Ananas comosus var. bracteatus f. tricolor) provides insights into elaboration of leaf color

Ananas comosus var. bracteatus f. tricolor (GL1) is a red pineapple accession whose mostly green leaves with chimeric white leaf margins turn red in spring and autumn and during flowering. It is an important ornamental plant and ideal plant research model for anthocyanin metabolism, chimeric leaf development, and photosynthesis. Here, we generated a highly contiguous chromosome-scale genome assembly for GL1 and compared it with other 3 published pineapple assemblies (var. comosus accessions MD2 and F153, and var. bracteatus accession CB5). The GL1 assembly has a total size of ∼461 Mb, with a contig N50 of ∼2.97 Mb and Benchmarking Universal Single-Copy Ortholog score of 97.3%. More than 99% of the contigs are anchored to 25 pseudochromosomes. Compared with the other 3 published pineapple assemblies, the GL1 assembly was confirmed to be more continuous. Our evolutionary analysis showed that the Bromeliaceae and Poaceae diverged from their nearest common ancestor ∼82.36 million years ago (MYA). Population structure analysis showed that while GL1 has not undergone admixture, bracteatus accession CB5 has resulted from admixture of 3 species of Ananas. Through classification of orthogroups, analysis of genes under positive selection, and analysis of presence/absence variants, we identified a series of genes related to anthocyanin metabolism and development of chimeric leaves. The structure and evolution of these genes were compared among the published pineapple assemblies with reveal candidate genes for these traits. The GL1 genome assembly and its comparisons with other 3 pineapple genome assemblies provide a valuable resource for the genetic improvement of pineapple and serve as a model for understanding the genomic basis of important traits in different pineapple varieties and other pan-cereal crops.

Genome assembly, annotation, and comparative analysis of the cattail Typha latifolia

Cattails (Typha species) comprise a genus of emergent wetland plants with a global distribution. Typha latifolia and Typha angustifolia are two of the most widespread species, and in areas of sympatry can interbreed to produce the hybrid Typha × glauca. In some regions, the relatively high fitness of Typha × glauca allows it to outcompete and displace both parent species, while simultaneously reducing plant and invertebrate biodiversity, and modifying nutrient and water cycling. We generated a high-quality whole-genome assembly of T. latifolia using PacBio long-read and high coverage Illumina sequences that will facilitate evolutionary and ecological studies in this hybrid zone. Genome size was 287 Mb and consisted of 1158 scaffolds, with an N50 of 8.71 Mb; 43.84% of the genome were identified as repetitive elements. The assembly has a BUSCO score of 96.03%, and 27,432 genes and 2700 RNA sequences were putatively identified. Comparative analysis detected over 9000 shared orthologs with related taxa and phylogenomic analysis supporting T. latifolia as a divergent lineage within Poales. This high-quality scaffold-level reference genome will provide a useful resource for future population genomic analyses and improve our understanding of Typha hybrid dynamics.

The genome and preliminary single-nuclei transcriptome of Lemna minuta reveals mechanisms of invasiveness

The ability to trace every cell in some model organisms has led to the fundamental understanding of development and cellular function. However, in plants the complexity of cell number, organ size, and developmental time makes this a challenge even in the diminutive model plant Arabidopsis (Arabidopsis thaliana). Duckweed, basal nongrass aquatic monocots, provide an opportunity to follow every cell of an entire plant due to their small size, reduced body plan, and fast clonal growth habit. Here we present a chromosome-resolved genome for the highly invasive Lesser Duckweed (Lemna minuta) and generate a preliminary cell atlas leveraging low cell coverage single-nuclei sequencing. We resolved the 360 megabase genome into 21 chromosomes, revealing a core nonredundant gene set with only the ancient tau whole-genome duplication shared with all monocots, and paralog expansion as a result of tandem duplications related to phytoremediation. Leveraging SMARTseq2 single-nuclei sequencing, which provided higher gene coverage yet lower cell count, we profiled 269 nuclei covering 36.9% (8,457) of the L. minuta transcriptome. Since molecular validation was not possible in this nonmodel plant, we leveraged gene orthology with model organism single-cell expression datasets, gene ontology, and cell trajectory analysis to define putative cell types. We found that the tissue that we computationally defined as mesophyll expressed high levels of elemental transport genes consistent with this tissue playing a role in L. minuta wastewater detoxification. The L. minuta genome and preliminary cell map provide a paradigm to decipher developmental genes and pathways for an entire plant.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Current status and future perspectives on the evolution of cis-regulatory elements in plants

Cis-regulatory elements (CREs) are short stretches (∼5-15 base pairs) of DNA capable of being bound by a transcription factor and influencing the expression of nearby genes. These regions are of great interest to anyone studying the relationship between phenotype and genotype as these sequences often dictate genes' spatio-temporal expression. Indeed, several associative signals between genotype and phenotype are known to lie outside of protein-coding regions. Therefore, a key to understand evolutionary biology requires their characterization in current and future genome assemblies. In this review, we cover some recent examples of how CRE variation contributes to phenotypic evolution, discuss evidence for the selective pressures experienced by non-coding regions of the genome, and consider several studies on accessible chromatin regions in plants and what they can tell us about CREs. Finally, we discuss how current advances in sequencing technologies will improve our knowledge of CRE variation.

Horticultural innovation by viral-induced gene regulation of carotenogenesis

Multipartite viral vectors provide a simple, inexpensive and effective biotechnological tool to transiently manipulate (i.e. reduce or increase) gene expression in planta and characterise the function of genetic traits. The development of virus-induced gene regulation (VIGR) systems usually involve the targeted silencing or overexpression of genes involved in pigment biosynthesis or degradation in plastids, thereby providing rapid visual assessment of success in establishing RNA- or DNA-based VIGR systems in planta. Carotenoids pigments provide plant tissues with an array of yellow, orange, and pinkish-red colours. VIGR-induced transient manipulation of carotenoid-related gene expression has advanced our understanding of carotenoid biosynthesis, regulation, accumulation and degradation, as well as plastid signalling processes. In this review, we describe mechanisms of VIGR, the importance of carotenoids as visual markers of technology development, and knowledge gained through manipulating carotenogenesis in model plants as well as horticultural crops not always amenable to transgenic approaches. We outline how VIGR can be utilised in plants to fast-track the characterisation of gene function(s), accelerate fruit tree breeding programs, edit genomes, and biofortify plant products enriched in carotenoid micronutrients for horticultural innovation.

Chromosome-scale assembly and population diversity analyses provide insights into the evolution of Sapindus mukorossi

Sapindus mukorossi is an environmentally friendly plant and renewable energy source whose fruit has been widely used for biomedicine, biodiesel, and biological chemicals due to its richness in saponin and oil contents. Here, we report the first chromosome-scale genome assembly of S. mukorossi (covering ~391 Mb with a scaffold N50 of 24.66 Mb) and characterize its genetic architecture and evolution by resequencing 104 S. mukorossi accessions. Population genetic analyses showed that genetic diversity in the southwestern distribution area was relatively higher than that in the northeastern distribution area. Gene flow events indicated that southwest species may be the donor population for the distribution areas in China. Genome-wide selective sweep analysis showed that a large number of genes are involved in defense responses, growth and development, including SmRPS2, SmRPS4, SmRPS7, SmNAC2, SmNAC23, SmNAC102, SmWRKY6, SmWRKY26, and SmWRKY33. We also identified several candidate genes controlling six agronomic traits by genome-wide association studies, including SmPCBP2, SmbHLH1, SmCSLD1, SmPP2C, SmLRR-RKs, and SmAHP. Our study not only provides a rich genomic resource for further basic research on Sapindaceae woody trees but also identifies several economically significant genes for genomics-enabled improvements in molecular breeding.

ksrates: positioning whole-genome duplications relative to speciation events in KS distributions

SUMMARY

We present ksrates, a user-friendly command-line tool to position ancient whole-genome duplication events with respect to speciation events in a phylogeny by comparing paralog and ortholog KS distributions derived from genomic or transcriptomic sequences, while adjusting for substitution rate differences among the lineages involved.

AVAILABILITY AND IMPLEMENTATION

ksrates is implemented in Python 3 and as a Nextflow pipeline. The source code, Singularity and Docker containers, documentation and tutorial are available via https://github.com/VIB-PSB/ksrates.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Improved High-Quality Genome Assembly and Annotation of Pineapple (Ananas comosus) Cultivar MD2 Revealed Extensive Haplotype Diversity and Diversified FRS/FRF Gene Family

Pineapple (Ananas comosus (L.) Merr.) is the second most important tropical fruit crop globally, and ‘MD2’ is the most important cultivated variety. A high-quality genome is important for molecular-based breeding, but available pineapple genomes still have some quality limitations. Here, PacBio and Hi-C data were used to develop a new high-quality MD2 assembly and gene prediction. Compared to the previous MD2 assembly, major improvements included a 26.6-fold increase in contig N50 length, phased chromosomes, and >6000 new genes. The new MD2 assembly also included 161.6 Mb additional sequences and >3000 extra genes compared to the F153 genome. Over 48% of the predicted genes harbored potential deleterious mutations, indicating that the high level of heterozygosity in this species contributes to maintaining functional alleles. The genome was used to characterize the FAR1-RELATED SEQUENCE (FRS) genes that were expanded in pineapple and rice. Transposed and dispersed duplications contributed to expanding the numbers of these genes in the pineapple lineage. Several AcFRS genes were differentially expressed among tissue-types and stages of flower development, suggesting that their expansion contributed to evolving specialized functions in reproductive tissues. The new MD2 assembly will serve as a new reference for genetic and genomic studies in pineapple.

The Cymbidium genome reveals the evolution of unique morphological traits

The marvelously diverse Orchidaceae constitutes the largest family of angiosperms. The genus Cymbidium in Orchidaceae is well known for its unique vegetation, floral morphology, and flower scent traits. Here, a chromosome-scale assembly of the genome of Cymbidium ensifolium (Jianlan) is presented. Comparative genomic analysis showed that C. ensifolium has experienced two whole-genome duplication (WGD) events, the most recent of which was shared by all orchids, while the older event was the τ event shared by most monocots. The results of MADS-box genes analysis provided support for establishing a unique gene model of orchid flower development regulation, and flower shape mutations in C. ensifolium were shown to be associated with the abnormal expression of MADS-box genes. The most abundant floral scent components identified included methyl jasmonate, acacia alcohol and linalool, and the genes involved in the floral scent component network of C. ensifolium were determined. Furthermore, the decreased expression of photosynthesis-antennae and photosynthesis metabolic pathway genes in leaves was shown to result in colorful striped leaves, while the increased expression of MADS-box genes in leaves led to perianth-like leaves. Our results provide fundamental insights into orchid evolution and diversification.

Small parsimony for natural genomes in the DCJ-indel model

The Small Parsimony Problem (SPP) aims at finding the gene orders at internal nodes of a given phylogenetic tree such that the overall genome rearrangement distance along the tree branches is minimized. This problem is intractable in most genome rearrangement models, especially when gene duplication and loss are considered. In this work, we describe an Integer Linear Program algorithm to solve the SPP for natural genomes, i.e. genomes that contain conserved, unique, and duplicated markers. The evolutionary model that we consider is the DCJ-indel model that includes the Double-Cut and Join rearrangement operation and the insertion and deletion of genome segments. We evaluate our algorithm on simulated data and show that it is able to reconstruct very efficiently and accurately ancestral gene orders in a very comprehensive evolutionary model.

Integrative transcriptomic and metabolomic analysis of D-leaf of seven pineapple varieties differing in N-P-K% contents

BACKGROUND

Pineapple (Ananas comosus L. Merr.) is the third most important tropical fruit in China. In other crops, farmers can easily judge the nutritional requirements from leaf color. However, concerning pineapple, it is difficult due to the variation in leaf color of the cultivated pineapple varieties. A detailed understanding of the mechanisms of nutrient transport, accumulation, and assimilation was targeted in this study. We explored the D-leaf nitrogen (N), phosphorus (P), and potassium (K) contents, transcriptome, and metabolome of seven pineapple varieties.

RESULTS

Significantly higher N, P, and K% contents were observed in Bali, Caine, and Golden pineapple. The transcriptome sequencing of 21 libraries resulted in the identification of 14,310 differentially expressed genes in the D-leaves of seven pineapple varieties. Genes associated with N transport and assimilation in D-leaves of pineapple was possibly regulated by nitrate and ammonium transporters, and glutamate dehydrogenases play roles in N assimilation in arginine biosynthesis pathways. Photosynthesis and photosynthesis-antenna proteins pathways were also significantly regulated between the studied genotypes. Phosphate transporters and mitochondrial phosphate transporters were differentially regulated regarding inorganic P transport. WRKY, MYB, and bHLH transcription factors were possibly regulating the phosphate transporters. The observed varying contents of K% in the D-leaves was associated to the regulation of K⁺ transporters and channels under the influence of Ca²⁺ signaling. The UPLC-MS/MS analysis detected 873 metabolites which were mainly classified as flavonoids, lipids, and phenolic acids.

CONCLUSIONS

These findings provide a detailed insight into the N, P, K% contents in pineapple D-leaf and their transcriptomic and metabolomic signatures.

High-quality reference genome sequences of two coconut cultivars provide insights into evolution of monocot chromosomes and differentiation of fiber content and plant height

BACKGROUND

Coconut is an important tropical oil and fruit crop whose evolutionary position renders it a fantastic species for the investigation of the evolution of monocot chromosomes and the subsequent differentiation of ancient plants.

RESULTS

Here, we report the assembly and annotation of reference-grade genomes of Cn. tall and Cn. dwarf, whose genome sizes are 2.40 Gb and 2.39 Gb, respectively. The comparative analysis reveals that the two coconut subspecies diverge about 2-8 Mya while the conserved Arecaceae-specific whole-genome duplication (ω WGD) occurs approximately 47-53 Mya. It additionally allows us to reconstruct the ancestral karyotypes of the ten ancient monocot chromosomes and the evolutionary trajectories of the 16 modern coconut chromosomes. Fiber synthesis genes in Cn. tall, related to lignin and cellulose synthesis, are found at a higher copy number and expression level than dwarf coconuts. Integrated multi-omics analysis reveals that the difference in coconut plant height is the result of altered gibberellin metabolism, with both the GA20ox copy number and a single-nucleotide change in the promoter together leading to the difference in plant height between Cn. tall and Cn. dwarf.

CONCLUSION

We provide high-quality coconut genomes and reveal the genetic basis of trait differences between two coconuts through multi-omics analysis. We also reveal that the selection of plant height has been targeted for the same gene for millions of years, not only in natural selection of ancient plant as illustrated in coconut, but also for artificial selection in cultivated crops such as rice and maize.

Underwater CAM photosynthesis elucidated by Isoetes genome

To conserve water in arid environments, numerous plant lineages have independently evolved Crassulacean Acid Metabolism (CAM). Interestingly, Isoetes, an aquatic lycophyte, can also perform CAM as an adaptation to low CO2 availability underwater. However, little is known about the evolution of CAM in aquatic plants and the lack of genomic data has hindered comparison between aquatic and terrestrial CAM. Here, we investigate underwater CAM in Isoetes taiwanensis by generating a high-quality genome assembly and RNA-seq time course. Despite broad similarities between CAM in Isoetes and terrestrial angiosperms, we identify several key differences. Notably, Isoetes may have recruited the lesser-known 'bacterial-type' PEPC, along with the 'plant-type' exclusively used in other CAM and C4 plants for carboxylation of PEP. Furthermore, we find that circadian control of key CAM pathway genes has diverged considerably in Isoetes relative to flowering plants. This suggests the existence of more evolutionary paths to CAM than previously recognized.

Characterization of germline development and identification of genes associated with germline specification in pineapple

Understanding germline specification in plants could be advantageous for agricultural applications. In recent decades, substantial efforts have been made to understand germline specification in several plant species, including Arabidopsis, rice, and maize. However, our knowledge of germline specification in many agronomically important plant species remains obscure. Here, we characterized the female germline specification and subsequent female gametophyte development in pineapple using callose staining, cytological, and whole-mount immunolocalization analyses. We also determined the male germline specification and gametophyte developmental timeline and observed male meiotic behavior using chromosome spreading assays. Furthermore, we identified 229 genes that are preferentially expressed at the megaspore mother cell (MMC) stage during ovule development and 478 genes that are preferentially expressed at the pollen mother cell (PMC) stage of anther development using comparative transcriptomic analysis. The biological functions, associated regulatory pathways and expression patterns of these genes were also analyzed. Our study provides a convenient cytological reference for exploring pineapple germline development and a molecular basis for the future functional analysis of germline specification in related plant species.

Starch branching enzymes as putative determinants of postharvest quality in horticultural crops

Starch branching enzymes (SBEs) are key determinants of the structure and amount of the starch in plant organs, and as such, they have the capacity to influence plant growth, developmental, and fitness processes, and in addition, the industrial end-use of starch. However, little is known about the role of SBEs in determining starch structure-function relations in economically important horticultural crops such as fruit and leafy greens, many of which accumulate starch transiently. Further, a full understanding of the biological function of these types of starches is lacking. Because of this gap in knowledge, this minireview aims to provide an overview of SBEs in horticultural crops, to investigate the potential role of starch in determining postharvest quality. A systematic examination of SBE sequences in 43 diverse horticultural species, identified SBE1, 2 and 3 isoforms in all species examined except apple, olive, and Brassicaceae, which lacked SBE1, but had a duplicated SBE2. Among our findings after a comprehensive and critical review of published data, was that as apple, banana, and tomato fruits ripens, the ratio of the highly digestible amylopectin component of starch increases relative to the more digestion-resistant amylose fraction, with parallel increases in SBE2 transcription, fruit sugar content, and decreases in starch. It is tempting to speculate that during the ripening of these fruit when starch degradation occurs, there are rearrangements made to the structure of starch possibly via branching enzymes to increase starch digestibility to sugars. We propose that based on the known action of SBEs, and these observations, SBEs may affect produce quality, and shelf-life directly through starch accumulation, and indirectly, by altering sugar availability. Further studies where SBE activity is fine-tuned in these crops, can enrich our understanding of the role of starch across species and may improve horticulture postharvest quality.

Genome-Wide Identification and Expression Analysis of MADS-Box Family Genes in Litchi (Litchi chinensis Sonn.) and Their Involvement in Floral Sex Determination

Litchi possesses unique flower morphology and adaptive reproduction strategies. Although previous attention has been intensively devoted to the mechanisms underlying its floral induction, the molecular basis of flower sex determination remains largely unknown. MADS-box genes are promising candidates for this due to their significant roles in various aspects of inflorescence and flower organogenesis. Here, we present a detailed overview of phylogeny and expression profiles of 101 MADS-box genes that were identified in litchi. These LcMADSs are unevenly located across the 15 chromosomes and can be divided into type I and type II genes. Fifty type I MADS-box genes are subdivided into Mα, Mβ and Mγ subgroups, while fifty-one type II LcMADSs consist of 37 MIKC^C -type and 14 MIKC *-type genes. Promoters of both types of LcMADS genes contain mainly ABA and MeJA response elements. Tissue-specific and development-related expression analysis reveal that LcMADS51 could be positively involved in litchi carpel formation, while six MADS-box genes, including LcMADS42/46/47/75/93/100, play a possible role in stamen development. GA is positively involved in the sex determination of litchi flowers by regulating the expression of LcMADS51 (LcSTK). However, JA down-regulates the expression of floral organ identity genes, suggesting a negative role in litchi flower development.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.

Genomic insights into the fast growth of paulownias and the formation of Paulownia witches' broom

Paulownias are among the fastest growing trees in the world, but they often suffer tremendous loss of wood production due to infection by Paulownia witches' broom (PaWB) phytoplasmas. In this study, we have sequenced and assembled a high-quality nuclear genome of Paulownia fortunei, a commonly cultivated paulownia species. The assembled genome of P. fortunei is 511.6 Mb in size, with 93.2% of its sequences anchored to 20 pseudo-chromosomes, and it contains 31 985 protein-coding genes. Phylogenomic analyses show that the family Paulowniaceae is sister to a clade composed of Phrymaceae and Orobanchaceae. Higher photosynthetic efficiency is achieved by integrating C3 photosynthesis and the crassulacean acid metabolism pathway, which may contribute to the extremely fast growth habit of paulownia trees. Comparative transcriptome analyses reveal modules related to cambial growth and development, photosynthesis, and defense responses. Additional genome sequencing of PaWB phytoplasma, combined with functional analyses, indicates that the effector PaWB-SAP54 interacts directly with Paulownia PfSPLa, which in turn causes the degradation of PfSPLa by the ubiquitin-mediated pathway and leads to the formation of witches' broom. Taken together, these results provide significant insights into the biology of paulownias and the regulatory mechanism for the formation of PaWB.

Convergent evolution of gene regulatory networks underlying plant adaptations to dry environments

Plants transitioned from an aquatic to a terrestrial lifestyle during their evolution. On land, fluctuations on water availability in the environment became one of the major problems they encountered. The appearance of morpho-physiological adaptations to cope with and tolerate water loss from the cells was undeniably useful to survive on dry land. Some of these adaptations, such as carbon concentrating mechanisms (CCMs), desiccation tolerance (DT) and root impermeabilization, appeared in multiple plant lineages. Despite being crucial for evolution on land, it has been unclear how these adaptations convergently evolved in the various plant lineages. Recent advances on whole genome and transcriptome sequencing are revealing that co-option of genes and gene regulatory networks (GRNs) is a common feature underlying the convergent evolution of these adaptations. In this review, we address how the study of CCMs and DT has provided insight into convergent evolution of GRNs underlying plant adaptation to dry environments, and how these insights could be applied to currently emerging understanding of evolution of root impermeabilization through different barrier cell types. We discuss examples of co-option, conservation and innovation of genes and GRNs at the cell, tissue and organ levels revealed by recent phylogenomic (comparative genomic) and comparative transcriptomic studies.

Genome-Wide Investigation of SBT Family Genes in Pineapple and Functional Analysis of AcoSBT1.12 in Floral Transition

SBT (Subtilisin-like serine protease), a clan of serine proteolytic enzymes, plays a versatile role in plant growth and defense. Although SBT family genes have been obtained from studies of dicots such as Arabidopsis, little is known about the potential functions of SBT in the monocots. In this study, 54 pineapple SBT genes (AcoSBTs) were divided into six subfamilies and then identified to be experienced strong purifying selective pressure and distributed on 25 chromosomes unevenly. Cis-acting element analysis indicated that almost all AcoSBTs promoters contain light-responsive elements. Further, the expression pattern via RNA-seq data showed that different AcoSBTs were preferentially expressed in different above-ground tissues. Transient expression in tobacco showed that AcoSBT1.12 was located in the plasma membrane. Moreover, Transgenic Arabidopsis ectopically overexpressing AcoSBT1.12 exhibited delayed flowering time. In addition, under the guidance of bioinformatic prediction, we found that AcoSBT1.12 could interact with AcoCWF19L, AcoPUF2, AcoCwfJL, Aco012905, and AcoSZF1 by yeast-two hybrid (Y2H). In summary, this study provided valuable information on pineapple SBT genes and illuminated the biological function of AcoSBT1.12 in floral transition.

Recent amplification of microsatellite-associated miniature inverted-repeat transposable elements in the pineapple genome

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous DNA transposable elements that play important roles in genome organization and evolution. Genome-wide identification and characterization of MITEs provide essential information for understanding genome structure and evolution.

RESULTS

We performed genome-wide identification and characterization of MITEs in the pineapple genome. The top two MITE families, accounting for 29.39% of the total MITEs and 3.86% of the pineapple genome, have insertion preference in (TA) n dinucleotide microsatellite regions. We therefore named these MITEs A. comosus microsatellite-associated MITEs (Ac-mMITEs). The two Ac-mMITE families, Ac-mMITE-1 and Ac-mMITE-2, shared sequence similarity in the terminal inverted repeat (TIR) regions, suggesting that these two Ac-mMITE families might be derived from a common or closely related autonomous elements. The Ac-mMITEs are frequently clustered via adjacent insertions. Among the 21,994 full-length Ac-mMITEs, 46.1% of them were present in clusters. By analyzing the Ac-mMITEs without (TA) n microsatellite flanking sequences, we found that Ac-mMITEs were likely derived from Mutator-like DNA transposon. Ac-MITEs showed highly polymorphic insertion sites between cultivated pineapples and their wild relatives. To better understand the evolutionary history of Ac-mMITEs, we filtered and performed comparative analysis on the two distinct groups of Ac-mMITEs, microsatellite-targeting MITEs (mt-MITEs) that are flanked by dinucleotide microsatellites on both sides and mutator-like MITEs (ml-MITEs) that contain 9/10 bp TSDs. Epigenetic analysis revealed a lower level of host-induced silencing on the mt-MITEs in comparison to the ml-MITEs, which partially explained the significantly higher abundance of mt-MITEs in pineapple genome. The mt-MITEs and ml-MITEs exhibited differential insertion preference to gene-related regions and RNA-seq analysis revealed their differential influences on expression regulation of nearby genes.

CONCLUSIONS

Ac-mMITEs are the most abundant MITEs in the pineapple genome and they were likely derived from Mutator-like DNA transposon. Preferential insertion in (TA) n microsatellite regions of Ac-mMITEs occurred recently and is likely the result of damage-limiting strategy adapted by Ac-mMITEs during co-evolution with their host. Insertion in (TA) n microsatellite regions might also have promoted the amplification of mt-MITEs. In addition, mt-MITEs showed no or negligible impact on nearby gene expression, which may help them escape genome control and lead to their amplification.

Inference of Gene Regulatory Network Uncovers the Linkage between Circadian Clock and Crassulacean Acid Metabolism in Kalanchoë fedtschenkoi

The circadian clock drives time-specific gene expression, enabling biological processes to be temporally controlled. Plants that conduct crassulacean acid metabolism (CAM) photosynthesis represent an interesting case of circadian regulation of gene expression as stomatal movement is temporally inverted relative to stomatal movement in C3 plants. The mechanisms behind how the circadian clock enabled physiological differences at the molecular level is not well understood. Recently, the rescheduling of gene expression was reported as a mechanism to explain how CAM evolved from C3. Therefore, we investigated whether core circadian clock genes in CAM plants were re-phased during evolution, or whether networks of phase-specific genes were simply re-wired to different core clock genes. We identified candidate core clock genes based on gene expression features and then applied the Local Edge Machine (LEM) algorithm to infer regulatory relationships between this new set of core candidates and known core clock genes in Kalanchoë fedtschenkoi. We further inferred stomata-related gene targets for known and candidate core clock genes and constructed a gene regulatory network for core clock and stomata-related genes. Our results provide new insight into the mechanism of circadian control of CAM-related genes in K. fedtschenkoi, facilitating the engineering of CAM machinery into non-CAM plants for sustainable crop production in water-limited environments.

Comparative transcriptome analysis of a fan-shaped inflorescence in pineapple using RNA-seq

Pineapple plant usually has a capitulum. However, a fan-shaped inflorescence was exceptionally evolved in pineapple, having multiple crown buds. In order to reveal the molecular mechanisms of the formation of the fan-shaped inflorescence, fruit traits and the transcriptional differences between the fan-shaped inflorescence and the wild-shaped inflorescence pineapples were analyzed in three tissues, i.e., the flower stem apex, the base of the inflorescence, and the inflorescence axis. The weight (i.e., individual yield) of fan-shaped fruit is 4.5 times that of wild-shaped fruit；and non-significant difference in soluble solids, soluble sugar, titratable acid, and Vitamin C was found. Between the fan-shaped inflorescence and wild-shaped inflorescence, a total of 5370 differentially expressed genes were identified across the three tissues. Of these genes, there were 489 overlapping differentially expressed genes in all three tissue comparisons. Between the two pineapples, functional analysis indicated that 444 transcription factors and 206 inflorescence development-related genes were differentially expressed in at least one tissue comparison, while 45 transcription factors and 21 inflorescence development-related genes were overlapped across three tissues. Among the 489 overlapping differentially expressed genes in the three tissue comparisons, excluding the inflorescence development-related genes and transcription factors, 80 of them revealed a higher percentage of involvement in the biological processes relating to response to auxin, and reproductive processes. RNA-seq value and real-time quantitative PCR analysis exhibited the similar gene expression patterns in the three tissues. Our result provided novel cues for understanding the molecular mechanisms of the formation of the fan-shaped inflorescence in pineapple, making a valuable resource for the study of plant breeding and the speciation of pineapple.

Random mutagenesis in vegetatively propagated crops: opportunities, challenges and genome editing prospects

In order to meet the growing human food and nutrition demand a perpetual process of crop improvement is idealized. It has seen changing trends and varying concepts throughout human history; from simple selection to complex gene-editing. Among these techniques, random mutagenesis has been shown to be a promising technology to achieve desirable genetic gain with less time and minimal efforts. Over the decade, several hundred varieties have been released through random mutagenesis, but the production is falling behind the demand. Several food crops like banana, potato, cassava, sweet potato, apple, citrus, and others are vegetatively propagated. Since such crops are not propagated through seed, genetic improvement through classical breeding is impractical for them. Besides, in the case of polyploids, accomplishment of allelic homozygosity requires a considerable land area, extensive fieldwork with huge manpower, and hefty funding for an extended period of time. Apart from induction, mapping of induced genes to facilitate the knowledge of biological processes has been performed only in a few selected facultative vegetative crops like banana and cassava which can form a segregating population. During the last few decades, there has been a shift in the techniques used for crop improvement. With the introduction of the robust technologies like meganucleases, zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats (CRISPR) more and more crops are being subjected to gene editing. However, more work needs to be done in case of vegetatively propagated crops.

Comparative phylogenetic analysis of CBL reveals the gene family evolution and functional divergence in Saccharum spontaneum

BACKGROUND

The identification and functional analysis of genes that improve tolerance to low potassium stress in S. spontaneum is crucial for breeding sugarcane cultivars with efficient potassium utilization. Calcineurin B-like (CBL) protein is a calcium sensor that interacts with specific CBL-interacting protein kinases (CIPKs) upon plants' exposure to various abiotic stresses.

RESULTS

In this study, nine CBL genes were identified from S. spontaneum. Phylogenetic analysis of 113 CBLs from 13 representative plants showed gene expansion and strong purifying selection in the CBL family. Analysis of CBL expression patterns revealed that SsCBL01 was the most commonly expressed gene in various tissues at different developmental stages. Expression analysis of SsCBLs under low K⁺ stress indicated that potassium deficiency moderately altered the transcription of SsCBLs. Subcellular localization showed that SsCBL01 is a plasma membrane protein and heterologous expression in yeast suggested that, while SsCBL01 alone could not absorb K⁺, it positively regulated K⁺ absorption mediated by the potassium transporter SsHAK1.

CONCLUSIONS

This study provided insights into the evolution of the CBL gene family and preliminarily demonstrated that the plasma membrane protein SsCBL01 was involved in the response to low K⁺ stress in S. spontaneum.

Haplotype-resolved genome assembly and allele-specific gene expression in cultivated ginger

Ginger (Zingiber officinale) is one of the most valued spice plants worldwide; it is prized for its culinary and folk medicinal applications and is therefore of high economic and cultural importance. Here, we present a haplotype-resolved, chromosome-scale assembly for diploid ginger anchored to 11 pseudochromosome pairs with a total length of 3.1 Gb. Remarkable structural variation was identified between haplotypes, and two inversions larger than 15 Mb on chromosome 4 may be associated with ginger infertility. We performed a comprehensive, spatiotemporal, genome-wide analysis of allelic expression patterns, revealing that most alleles are coordinately expressed. The alleles that exhibited the largest differences in expression showed closer proximity to transposable elements, greater coding sequence divergence, more relaxed selection pressure, and more transcription factor binding site differences. We also predicted the transcription factors potentially regulating 6-gingerol biosynthesis. Our allele-aware assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in ginger.

Somatic variations led to the selection of acidic and acidless orange cultivars

Somatic variations are a major source of genetic diversification in asexual plants, and underpin clonal evolution and the breeding of asexual crops. Sweet orange is a model species for studying somatic variation because it reproduces asexually through apomixis and is propagated asexually through grafting. To dissect the genomic basis of somatic variation, we de novo assembled a reference genome of sweet orange with an average of three gaps per chromosome and a N50 contig of 24.2 Mb, as well as six diploid genomes of somatic mutants of sweet oranges. We then sequenced 114 somatic mutants with an average genome coverage of 41×. Categorization of the somatic variations yielded insights into the single-nucleotide somatic mutations, structural variations and transposable element (TE) transpositions. We detected 877 TE insertions, and found TE insertions in the transporter or its regulatory genes associated with variation in fruit acidity. Comparative genomic analysis of sweet oranges from three diversity centres supported a dispersal from South China to the Mediterranean region and to the Americas. This study provides a global view on the somatic variations, the diversification and dispersal history of sweet orange and a set of candidate genes that will be useful for improving fruit taste and flavour.

Characterization of auxin transporter AUX, PIN and PILS gene families in pineapple and evaluation of expression profiles during reproductive development and under abiotic stresses

Polar auxin transport in plant is mediated by influx and efflux transporters, which are encoded by AUX/LAX, PIN and PILS genes, respectively. The auxin transporter gene families have been characterized in several species from monocots and eudicots. However, a genome-wide overview of auxin transporter gene families in pineapple is not yet available. In this study, we identified a total of threeAcAUX genes, 12 AcPIN genes, and seven AcPILS genes in the pineapple genome, which were variably located on 15 chromosomes. The exon-intron structure of these genes and properties of deduced proteins were relatively conserved within the same family. Most protein motifs were widespread in the AUX, PIN or PILS proteins, whereas a few motifs were absent in only one or two proteins. Analysis of the expression profiles of these genes elucidated that several genes exhibited either preferential or tissue-specific expression patterns in vegetative and/or reproductive tissues. AcAUX2 was specifically expressed in the early developmental ovules, while AcPIN1b and AcPILS2 were strongly expressed in stamens and ovules. AcPIN9b, AcPILS1, AcPILS6a, 6b and 6c were abundantly expressed in stamens. Furthermore, qRT-PCR results showed that several genes in these families were responsive to various abiotic stresses. Comparative analysis indicated that the genes with close evolutionary relationships among pineapple, rice and Arabidopsis exhibited similar expression patterns. Overexpression of the AcAUX1 in Arabidopsis rescued the phenotype in aux1-T, and resulted in increased lateral roots in WT. These results will provide new insights into auxin transporter genes of pineapple and facilitate our understanding of their roles in pineapple growth and development.

Conversion between 100-million-year-old duplicated genes contributes to rice subspecies divergence

BACKGROUND

Duplicated gene pairs produced by ancient polyploidy maintain high sequence similarity over a long period of time and may result from illegitimate recombination between homeologous chromosomes. The genomes of Asian cultivated rice Oryza sativa ssp. indica (XI) and Oryza sativa ssp. japonica (GJ) have recently been updated, providing new opportunities for investigating ongoing gene conversion events and their impact on genome evolution.

RESULTS

Using comparative genomics and phylogenetic analyses, we evaluated gene conversion rates between duplicated genes produced by polyploidization 100 million years ago (mya) in GJ and XI. At least 5.19-5.77% of genes duplicated across the three rice genomes were affected by whole-gene conversion after the divergence of GJ and XI at ~ 0.4 mya, with more (7.77-9.53%) showing conversion of only portions of genes. Independently converted duplicates surviving in the genomes of different subspecies often use the same donor genes. The ongoing gene conversion frequency was higher near chromosome termini, with a single pair of homoeologous chromosomes, 11 and 12, in each rice genome being most affected. Notably, ongoing gene conversion has maintained similarity between very ancient duplicates, provided opportunities for further gene conversion, and accelerated rice divergence. Chromosome rearrangements after polyploidization are associated with ongoing gene conversion events, and they directly restrict recombination and inhibit duplicated gene conversion between homeologous regions. Furthermore, we found that the converted genes tended to have more similar expression patterns than nonconverted duplicates. Gene conversion affects biological functions associated with multiple genes, such as catalytic activity, implying opportunities for interaction among members of large gene families, such as NBS-LRR disease-resistance genes, contributing to the occurrence of the gene conversion.

CONCLUSION

Duplicated genes in rice subspecies generated by grass polyploidization ~ 100 mya remain affected by gene conversion at high frequency, with important implications for the divergence of rice subspecies.

Engineering of Crassulacean Acid Metabolism

Crassulacean acid metabolism (CAM) has evolved from a C₃ ground state to increase water use efficiency of photosynthesis. During CAM evolution, selective pressures altered the abundance and expression patterns of C₃ genes and their regulators to enable the trait. The circadian pattern of CO₂ fixation and the stomatal opening pattern observed in CAM can be explained largely with a regulatory architecture already present in C₃ plants. The metabolic CAM cycle relies on enzymes and transporters that exist in C₃ plants and requires tight regulatory control to avoid futile cycles between carboxylation and decarboxylation. Ecological observations and modeling point to mesophyll conductance as a major factor during CAM evolution. The present state of knowledge enables suggestions for genes for a minimal CAM cycle for proof-of-concept engineering, assuming altered regulation of starch synthesis and degradation are not critical elements of CAM photosynthesis and sufficient malic acid export from the vacuole is possible.

Genome-wide identification and expression profiling of DREB genes in Saccharum spontaneum

Background

The dehydration-responsive element-binding proteins (DREBs) are important transcription factors that interact with a DRE/CRT (C-repeat) sequence and involve in response to multiple abiotic stresses in plants. Modern sugarcane are hybrids from the cross between Saccharum spontaneum and Saccharum officinarum, and the high sugar content is considered to the attribution of S. officinaurm, while the stress tolerance is attributed to S. spontaneum. To understand the molecular and evolutionary characterization and gene functions of the DREBs in sugarcane, based on the recent availability of the whole genome information, the present study performed a genome-wide in silico analysis of DREB genes and transcriptome analysis in the polyploidy S. spontaneum.

Results

Twelve DREB1 genes and six DREB2 genes were identified in S. spontaneum genome and all proteins contained a conserved AP2/ERF domain. Eleven SsDREB1 allele genes were assumed to be originated from tandem duplications, and two of them may be derived after the split of S. spontaneum and the proximal diploid species sorghum, suggesting tandem duplication contributed to the expansion of DREB1-type genes in sugarcane. Phylogenetic analysis revealed that one DREB2 gene was lost during the evolution of sugarcane. Expression profiling showed different SsDREB genes with variable expression levels in the different tissues, indicating seven SsDREB genes were likely involved in the development and photosynthesis of S. spontaneum. Furthermore, SsDREB1F, SsDREB1L, SsDREB2D, and SsDREB2F were up-regulated under drought and cold condition, suggesting that these four genes may be involved in both dehydration and cold response in sugarcane.

Conclusions

These findings demonstrated the important role of DREBs not only in the stress response, but also in the development and photosynthesis of S. spontaneum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07799-5.

Transition From Proto-Kranz-Type Photosynthesis to HCO3 - Use Photosynthesis in the Amphibious Plant Hygrophila polysperma

Hygrophila polysperma is a heterophyllous amphibious plant. The growth of H. polysperma in submerged conditions is challenging due to the low CO2 environment, increased resistance to gas diffusion, and bicarbonate ion (HCO3 -) being the dominant dissolved inorganic carbon source. The submerged leaves of H. polysperma have significantly higher rates of underwater photosynthesis compared with the terrestrial leaves. 4,4'-Diisothiocyanatostilbene-2,2'-disulfonate (DIDS), an anion exchanger protein inhibitor, and ethoxyzolamide (EZ), an inhibitor of internal carbonic anhydrase, repressed underwater photosynthesis by the submerged leaves. These results suggested that H. polysperma acclimates to the submerged condition by using HCO3 - for photosynthesis. H. polysperma transports HCO3 - into the leaf by a DIDS-sensitive HCO3 - transporter and converted to CO2 by carbonic anhydrase. Additionally, proteome analysis revealed that submerged leaves accumulated fewer proteins associated with C4 photosynthesis compared with terrestrial leaves. This finding suggested that H. polysperma is capable of C4 and C3 photosynthesis in the terrestrial and submerged leaves, respectively. The ratio of phosphoenol pyruvate carboxylase to ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) in the submerged leaves was less than that in the terrestrial leaves. Upon anatomical observation, the terrestrial leaves exhibited a phenotype similar to the Kranz anatomy found among C4 plants; however, chloroplasts in the bundle sheath cells were not located adjacent to the vascular bundles, and the typical Kranz anatomy was absent in submerged leaves. These results suggest that H. polysperma performs proto-Kranz type photosynthesis in a terrestrial environment and shifts from a proto-Kranz type in terrestrial leaves to a HCO3 - use photosynthesis in the submerged environments.

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Draft genome of Puya raimondii (Bromeliaceae), the Queen of the Andes

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species.

The chromosome-level reference genome of Coptis chinensis provides insights into genomic evolution and berberine biosynthesis

Coptis chinensis Franch, a perennial herb, is mainly distributed in southeastern China. The rhizome of C. chinensis has been used as a traditional medicine for more than 2000 years in China and many other Asian countries. The pharmacological activities of C. chinensis have been validated by research. Here, we present a de novo high-quality genome of C. chinensis with a chromosome-level genome of ~958.20 Mb, a contig N50 of 1.58 Mb, and a scaffold N50 of 4.53 Mb. We found that the relatively large genome size of C. chinensis was caused by the amplification of long terminal repeat (LTR) retrotransposons. In addition, a whole-genome duplication event in ancestral Ranunculales was discovered. Comparative genomic analysis revealed that the tyrosine decarboxylase (TYDC) and (S)-norcoclaurine synthase (NCS) genes were expanded and that the aspartate aminotransferase gene (ASP5) was positively selected in the berberine metabolic pathway. Expression level and HPLC analyses showed that the berberine content was highest in the roots of C. chinensis in the third and fourth years. The chromosome-level reference genome of C. chinensis provides important genomic data for molecular-assisted breeding and active ingredient biosynthesis.

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

Genome Skimming Reveals Widespread Hybridization in a Neotropical Flowering Plant Radiation

The tropics hold at least an order of magnitude greater plant diversity than the temperate zone, yet the reasons for this difference are still subject to debate. Much of tropical plant diversity is in highly speciose genera and understanding the drivers of such high species richness will help solve the tropical diversity enigma. Hybridization has recently been shown to underlie many adaptive radiations, but its role in the evolution of speciose tropical plant genera has received little attention. Here, we address this topic in the hyperdiverse Bromeliaceae genus Vriesea using genome skimming data covering the three genomic compartments. We find evidence for hybridization in ca. 11% of the species in our dataset, both within the genus and between Vriesea and other genera, which is commensurate with hybridization underlying the hyperdiversity of Vriesea, and potentially other genera in Tillandsioideae. While additional genomic research will be needed to further clarify the contribution of hybridization to the rapid diversification of Vriesea, our study provides an important first data point suggesting its importance to the evolution of tropical plant diversity.

Genome-wide identification, classification, and expression analysis of the HSF gene family in pineapple (Ananas comosus)

Transcription factors (TFs), such as heat shock transcription factors (HSFs), usually play critical regulatory functions in plant development, growth, and response to environmental cues. However, no HSFs have been characterized in pineapple thus far. Here, we identified 22 AcHSF genes from the pineapple genome. Gene structure, motifs, and phylogenetic analysis showed that AcHSF families were distinctly grouped into three subfamilies (12 in Group A, seven in Group B, and four in Group C). The AcHSF promoters contained various cis-elements associated with stress, hormones, and plant development processes, for instance, STRE, WRKY, and ABRE binding sites. The majority of HSFs were expressed in diverse pineapple tissues and developmental stages. The expression of AcHSF-B4b/AcHSF-B4c and AcHSF-A7b/AcHSF-A1c were enriched in the ovules and fruits, respectively. Six genes (AcHSF-A1a , AcHSF-A2, AcHSF-A9a, AcHSF-B1a, AcHSF-B2a, and AcHSF-C1a) were transcriptionally modified by cold, heat, and ABA. Our results provide an overview and lay the foundation for future functional characterization of the pineapple HSF gene family.

Genome Sequencing and Assembly Strategies and a Comparative Analysis of the Genomic Characteristics in Penaeid Shrimp Species

Penaeid shrimp (family Penaeidae) represents one of the most economically and ecologically important groups of crustaceans. However, their genome sequencing and assembly have encountered extreme difficulties during the last 20 years. In this study, based on our previous genomic data, we investigated the genomic characteristics of four penaeid shrimp species and identified potential factors that result in their poor genome assembly, including heterozygosity, polyploidization, and repeats. Genome sequencing and comparison of somatic cells (diploid) of the four shrimp species and a single sperm cell (haploid) of Litopenaeus vannamei identified a common bimodal distribution of K-mer depths, suggesting either high heterozygosity or abundant homo-duplicated sequences present in their genomes. However, penaeids have not undergone whole-genome duplication as indicated by a series of approaches. Besides, the remarkable expansion of simple sequence repeats was another outstanding character of penaeid genomes, which also made the genome assembly highly fragmented. Due to this situation, we tried to assemble the genome of penaeid shrimp using various genome sequencing and assembly strategies and compared the quality. Therefore, this study provides new insights about the genomic characteristics of penaeid shrimps while improving their genome assemblies.

Layers of crosstalk between circadian regulation and environmental signalling in plants

Circadian regulation has a pervasive influence upon plant development, physiology and metabolism, impacting upon components of fitness and traits of agricultural importance. Circadian regulation is inextricably connected to the responses of plants to their abiotic environments, from the cellular to whole plant scales. Here, we review the crosstalk that occurs between circadian regulation and responses to the abiotic environment from the intracellular scale through to naturally fluctuating environments. We examine the spatial crosstalk that forms part of plant circadian regulation, at the subcellular, tissue, organ and whole-plant scales. This includes a focus on chloroplast and mitochondrial signalling, alternative splicing, long-distance circadian signalling and circadian regulation within natural environments. We also consider mathematical models for plant circadian regulation, to suggest future areas for advancing understanding of roles for circadian regulation in plant responses to environmental cues.

The absence of the caffeine synthase gene is involved in the naturally decaffeinated status of Coffea humblotiana, a wild species from Comoro archipelago

Caffeine is the most consumed alkaloid stimulant in the world. It is synthesized through the activity of three known N-methyltransferase proteins. Here we are reporting on the 422-Mb chromosome-level assembly of the Coffea humblotiana genome, a wild and endangered, naturally caffeine-free, species from the Comoro archipelago. We predicted 32,874 genes and anchored 88.7% of the sequence onto the 11 chromosomes. Comparative analyses with the African Robusta coffee genome (C. canephora) revealed an extensive genome conservation, despite an estimated 11 million years of divergence and a broad diversity of genome sizes within the Coffea genus. In this genome, the absence of caffeine is likely due to the absence of the caffeine synthase gene which converts theobromine into caffeine through an illegitimate recombination mechanism. These findings pave the way for further characterization of caffeine-free species in the Coffea genus and will guide research towards naturally-decaffeinated coffee drinks for consumers.

Genome-wide identification and transcriptional analysis of ammonium transporters in Saccharum

Ammonium transporters (AMTs) are plasma membrane proteins that exclusively transport ammonium/ammonia. It is essential for the nitrogen demand of plantsby AMT-mediated acquisition of ammonium from soils. The molecular characteristics and evolutionary history of AMTs in Saccharum spp. remain unclear. We comprehensively evaluated the AMT gene family in the latest release of the S. spontaneum genome and identified 6 novel AMT genes. These genes belong to 3 clusters: AMT2 (2 genes), AMT3 (3 genes), and AMT4 (one gene). Evolutionary analyses suggested that the S. spontaneum AMT gene family may have expanded via whole-genome duplication events. All of the 6 AMT genes are located on 5 chromosomes of S. spontaneum. Expression analyses revealed that AMT3;2 was highly expressed in leaves and in the daytime, and AMT2;1/3;2/4 were dynamic expressed in different leaf segments, as well as AMT2;1/3;2 demonstrated a high transcript accumulation level in leaves and roots and were significantly dynamic expressed under low-nitrogen conditions. The results suggest the functional roles of AMT genes on tissue expression and ammonium absorption in Saccharum. This study will provide some reference information for further elucidation of the functional mechanism and regulation of expression of the AMT gene family in Saccharum.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

Genome-wide characterization and analysis of bHLH transcription factors related to anthocyanin biosynthesis in spine grapes (Vitis davidii)

As one of the largest transcription factor family, basic helix-loop-helix (bHLH) transcription factor family plays an important role in plant metabolism, physiology and growth. Berry color is one of the important factors that determine grape quality. However, the bHLH transcription factor family's function in anthocyanin synthesis of grape berry has not been studied systematically. We identified 115 bHLH transcription factors in grape genome and phylogenetic analysis indicated that bHLH family could be classified into 25 subfamilies. First, we screened six candidate genes by bioinformatics analysis and expression analysis. We found one of the candidate genes VdbHLH037 belonged to III (f) subfamily and interacted with genes related to anthocyanin synthesis through phylogenetic analysis and interaction network prediction. Therefore, we speculated that VdbHLH037 participated in the anthocyanin synthesis process. To confirm this, we transiently expressed VdbHLH037 in grape and Arabidopsis transformation. Compared with the control, transgenic materials can accumulate more anthocyanins. These results provide a good base to study the function of the VdbHLH family in anthocyanin synthesis of grape berry.

The chromosome-level genome of dragon fruit reveals whole-genome duplication and chromosomal co-localization of betacyanin biosynthetic genes

Dragon fruits are tropical fruits economically important for agricultural industries. As members of the family of Cactaceae, they have evolved to adapt to the arid environment. Here we report the draft genome of Hylocereus undatus, commercially known as the white-fleshed dragon fruit. The chromosomal level genome assembly contains 11 longest scaffolds corresponding to the 11 chromosomes of H. undatus. Genome annotation of H. undatus found ~29,000 protein-coding genes, similar to Carnegiea gigantea (saguaro). Whole-genome duplication (WGD) analysis revealed a WGD event in the last common ancestor of Cactaceae followed by extensive genome rearrangements. The divergence time between H. undatus and C. gigantea was estimated to be 9.18 MYA. Functional enrichment analysis of orthologous gene clusters (OGCs) in six Cactaceae plants found significantly enriched OGCs in drought resistance. Fruit flavor-related functions were overrepresented in OGCs that are significantly expanded in H. undatus. The H. undatus draft genome also enabled the discovery of carbohydrate and plant cell wall-related functional enrichment in dragon fruits treated with trypsin for a longer storage time. Lastly, genes of the betacyanin (a red-violet pigment and antioxidant with a very high concentration in dragon fruits) biosynthetic pathway were found to be co-localized on a 12 Mb region of one chromosome. The consequence may be a higher efficiency of betacyanin biosynthesis, which will need experimental validation in the future. The H. undatus draft genome will be a great resource to study various cactus plants.

How to resolve the enigma of diurnal malate remobilisation from the vacuole in plants with crassulacean acid metabolism?

Opening of stomata in plants with crassulacean acid metabolism (CAM) is mainly shifted to the night period when atmospheric CO₂ is fixed by phosphoenolpyruvate carboxylase and stored as malic acid in the vacuole. As such, CAM plants ameliorate transpirational water losses and display substantially higher water-use efficiency compared with C₃ and C₄ plants. In the past decade significant technical advances have allowed an unprecedented exploration of genomes, transcriptomes, proteomes and metabolomes of CAM plants and efforts are ongoing to engineer the CAM pathway in C₃ plants. Whilst research efforts have traditionally focused on nocturnal carboxylation, less information is known regarding the drivers behind diurnal malate remobilisation from the vacuole that liberates CO₂ to be fixed by RuBisCo behind closed stomata. To shed more light on this process, we provide a stoichiometric analysis to identify potentially rate-limiting steps underpinning diurnal malate mobilisation and help direct future research efforts. Within this remit we address three key questions: Q1 Does light-dependent assimilation of CO₂ via RuBisCo dictate the rate of malate mobilisation? Q2: Do the enzymes responsible for malate decarboxylation limit daytime mobilisation from the vacuole? Q3: Does malate efflux from the vacuole set the pace of decarboxylation?

The genome sequence of Aloe vera reveals adaptive evolution of drought tolerance mechanisms

Aloe vera is a species from Asphodelaceae family having characteristics like drought resistance and numerous medicinal properties. However, the genetic basis of these phenotypes is yet unknown primarily due to unavailability of its genome sequence. Thus, we report the first Aloe vera genome sequence comprising of 12.93 Gbp and harboring 86,177 protein-coding genes. It is the first genome from Asphodelaceae family and the largest angiosperm genome sequenced and assembled till date. We also report the first genome-wide phylogeny of monocots including Aloe vera to resolve its phylogenetic position. The comprehensive comparative analysis of Aloe vera with other available high-quality monocot genomes revealed adaptive evolution in several genes of drought stress response, CAM pathway, and circadian rhythm and positive selection in DNA damage response genes in Aloe vera. This study provides clues on the genetic basis of evolution of drought stress tolerance capabilities of Aloe vera.