The pineapple genome and the evolution of CAM photosynthesis

Unravelling the novel sex determination genotype with 'ZY' and a distinctive 2.15-2.95 Mb inversion among poplar species through haplotype-resolved genome assembly and comparative genomics analysis

Populus tomentosa, an indigenous tree species, is widely distributed and cultivated over 1,000,000 km2 in China, contributing significantly to forest production, ecological conservation and urban-rural greening. Although a reference genome is available for P. tomentosa, the intricate interspecific hybrid origins, chromosome structural variations (SVs) and sex determination mechanisms remain confusion and unclear due to its broad and even overlapping geographical distribution, extensive morphological variations and cross infiltration among white poplar species. We conducted a haplotype-resolved de novo assembly of P. tomentosa elite individual GM107, which comprises subgenomes a and b with a total genome size of 714.9 Mb. We then analysed the formation of hybrid species and the phylogenetic evolution and sex differentiation across the entire genus. Phylogenomic analyses suggested that GM107 likely originated from a hybridisation event between P. alba (♀) and P. davidiana (♂) approximately 3.8 Mya. A total of 1551 chromosome SVs were identified between the two subgenomes. More noteworthily, a distinctive inversion structure spanning 2.15-2.95 Mb was unveiled among Populus, Tacamahaca, Turaga, Aigeiros poplar species and Salix, highlighting a unique evolutionary feature. Intriguingly, a novel sex genotype of the ZY type, which represents a crossover between XY and ZW systems, was identified and confirmed through both natural and artificial hybrids populations. These novel insights offer significant theoretical value for the study of the species' evolutionary origins and serve as a valuable resource for ecological genetics and forest biotechnology.

Pathogen resistance was negatively regulated by the NAC transcription factor ScATAF1 in sugarcane

The NAC (NAM, ATAF, and CUC) is one of the largest transcription factor gene families in plants. In this study, 180, 141, and 131 NAC family members were identified from Saccharum complex, including S. officinarum, S. spontaneum, and Erianthus rufipilus. The Ka/Ks ratio of ATAF subfamily was all less than 1. Besides, 52 ATAF members from 12 representative plants were divided into three clades and there was only a significant expansion in maize. Surprisingly, ABA and JA cis-elements were abundant in hormonal response factor, followed by transcriptional regulator and abiotic stressor. The ATAF subfamily was differentially expressed in various tissues, under low temperature and smut pathogen treatments. Further, the ScATAF1 gene, with high expression in leaves, stem epidermis, and buds, was isolated. The encoded protein, lack of self-activation activity, was situated in the cell nucleus. Moreover, SA and JA stresses down-regulated the expression of this gene, while ABA, NaCl, and 4°C treatments led to its up-regulation. Interestingly, its expression in the smut susceptible sugarcane cultivars was much higher than the smut resistant ones. Notably, the colors presented slight brown in tobacco transiently overexpressing ScATAF1 at 1 d after DAB staining, while the symptoms were more obvious at 3 d after inoculation with Ralstonia solanacearum, with ROS, JA, and SA signaling pathway genes significantly up-regulated. We thus speculated ScATAF1 gene could negatively mediate hypersensitive reactions and produce ROS by JA and SA signaling pathways. These findings lay the groundwork for in-depth investigation on the biological roles of ATAF subfamily in sugarcane.

AcMYB266, a key regulator of the red coloration in pineapple peel: a case of subfunctionalization in tandem duplicated genes

Red fruit peel is an attractive target for pineapple breeding. Various pineapple accessions with distinct red coloration patterns exist; however, the precise molecular mechanism accounting for these differences remains unknown, which hinders the pineapple breeding process from combining high fruit quality with red peel. In this study, we characterized a transcription factor, AcMYB266, which is preferentially expressed in pineapple peel and positively regulates anthocyanin accumulation. Transgenic pineapple, Arabidopsis, and tobacco plants overexpressing AcMYB266 exhibited significant anthocyanin accumulation. Conversely, transient silencing of this gene led to decreased anthocyanin accumulation in pineapple red bracts. In-depth analysis indicated that variations of AcMYB266 sequences in the promoter instead of the protein-coding region seem to contribute to different red coloration patterns in peels of three representative pineapple varieties. In addition, we found that AcMYB266 was located in a cluster of four MYB genes exclusive to and conserved in Ananas species. Of this cluster, each was proved to regulate anthocyanin synthesis in different pineapple tissues, illustrating an interesting case of gene subfunctionalization after tandem duplication. In summary, we have characterized AcMYB266 as a key regulator of pineapple red fruit peel and identified an MYB cluster whose members were subfunctionalized to specifically regulate the red coloration of different pineapple tissues. The present study will assist in establishing a theoretical mechanism for pineapple breeding for red fruit peel and provide an interesting case for the investigation of gene subfunctionalization in plants.

Shifts in evolutionary lability underlie independent gains and losses of root-nodule symbiosis in a single clade of plants

Root nodule symbiosis (RNS) is a complex trait that enables plants to access atmospheric nitrogen converted into usable forms through a mutualistic relationship with soil bacteria. Pinpointing the evolutionary origins of RNS is critical for understanding its genetic basis, but building this evolutionary context is complicated by data limitations and the intermittent presence of RNS in a single clade of ca. 30,000 species of flowering plants, i.e., the nitrogen-fixing clade (NFC). We developed the most extensive de novo phylogeny for the NFC and an RNS trait database to reconstruct the evolution of RNS. Our analysis identifies evolutionary rate heterogeneity associated with a two-step process: An ancestral precursor state transitioned to a more labile state from which RNS was rapidly gained at multiple points in the NFC. We illustrate how a two-step process could explain multiple independent gains and losses of RNS, contrary to recent hypotheses suggesting one gain and numerous losses, and suggest a broader phylogenetic and genetic scope may be required for genome-phenome mapping.

Proteogenomic Gene Structure Validation in the Pineapple Genome

MD2 pineapple (Ananas comosus) is the second most important tropical crop that preserves crassulacean acid metabolism (CAM), which has high water-use efficiency and is fast becoming the most consumed fresh fruit worldwide. Despite the significance of environmental efficiency and popularity, until very recently, its genome sequence has not been determined and a high-quality annotated proteome has not been available. Here, we have undertaken a pilot proteogenomic study, analyzing the proteome of MD2 pineapple leaves using liquid chromatography-mass spectrometry (LC-MS/MS), which validates 1781 predicted proteins in the annotated F153 (V3) genome. In addition, a further 603 peptide identifications are found that map exclusively to an independent MD2 transcriptome-derived database but are not found in the standard F153 (V3) annotated proteome. Peptide identifications derived from these MD2 transcripts are also cross-referenced to a more recent and complete MD2 genome annotation, resulting in 402 nonoverlapping peptides, which in turn support 30 high-quality gene candidates novel to both pineapple genomes. Many of the validated F153 (V3) genes are also supported by an independent proteomics data set collected for an ornamental pineapple variety. The contigs and peptides have been mapped to the current F153 genome build and are available as bed files to display a custom gene track on the Ensembl Plants region viewer. These analyses add to the knowledge of experimentally validated pineapple genes and demonstrate the utility of transcript-derived proteomics to discover both novel genes and genetic structure in a plant genome, adding value to its annotation.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

TIP aquaporins in Cyperus esculentus: genome-wide identification, expression profiles, subcellular localizations, and interaction patterns

BACKGROUND

Tonoplast intrinsic proteins (TIPs), which typically mediate water transport across vacuolar membranes, play an essential role in plant growth, development, and stress responses. However, their characterization in tigernut (Cyperus esculentus L.), an oil-bearing tuber plant of the Cyperaceae family, is still in the infancy.

RESULTS

In this study, a first genome-wide characterization of the TIP subfamily was conducted in tigernut, resulting in ten members representing five previously defined phylogenetic groups, i.e., TIP1-5. Although the gene amounts are equal to that present in two model plants Arabidopsis and rice, the group composition and/or evolution pattern were shown to be different. Except for CeTIP1;3 that has no counterpart in both Arabidopsis and rice, complex orthologous relationships of 1:1, 1:2, 1:3, 2:1, and 2:2 were observed. Expansion of the CeTIP subfamily was contributed by whole-genome duplication (WGD), transposed, and dispersed duplications. In contrast to the recent WGD-derivation of CeTIP3;1/-3;2, synteny analyses indicated that TIP4 and - 5 are old WGD repeats of TIP2, appearing sometime before monocot-eudicot divergence. Expression analysis revealed that CeTIP genes exhibit diverse expression profiles and are subjected to developmental and diurnal fluctuation regulation. Moreover, when transiently overexpressed in tobacco leaves, CeTIP1;1 was shown to locate in the vacuolar membrane and function in homo/heteromultimer, whereas CeTIP2;1 is located in the cell membrane and only function in heteromultimer. Interestingly, CeTIP1;1 could mediate the tonoplast-localization of CeTIP2;1 via protein interaction, implying complex regulatory patterns.

CONCLUSIONS

Our findings provide a global view of CeTIP genes, which provide valuable information for further functional analysis and genetic improvement through manipulating key members in tigernut.

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

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

The Cissus quadrangularis genome reveals its adaptive features in an arid habitat

Cissus quadrangularis is a tetraploid species belonging to the Vitaceae family and is known for the Crassulacean acid metabolism (CAM) pathway in the succulent stem, while the leaves perform C3 photosynthesis. Here, we report a high-quality genome of C. quadrangularis comprising a total size of 679.2 Mb which was phased into two subgenomes. Genome annotation identified 51 857 protein-coding genes, while approximately 47.75% of the genome was composed of repetitive sequences. Gene expression ratios of two subgenomes demonstrated that the sub-A genome as the dominant subgenome played a vital role during the drought tolerance. Genome divergence analysis suggests that the tetraploidization event occurred around 8.9 million years ago. Transcriptome data revealed that pathways related to cutin, suberine, and wax metabolism were enriched in the stem during drought treatment, suggesting that these genes contributed to the drought adaption. Additionally, a subset of CAM-related genes displayed diurnal expression patterns in the succulent stems but not in leaves, indicating that stem-biased expression of existing genes contributed to the CAM evolution. Our findings provide insights into the mechanisms of drought adaptation and photosynthesis transition in plants.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Bringing CAM photosynthesis to the table: Paving the way for resilient and productive agricultural systems in a changing climate

Modern agricultural systems are directly threatened by global climate change and the resulting freshwater crisis. A considerable challenge in the coming years will be to develop crops that can cope with the consequences of declining freshwater resources and changing temperatures. One approach to meeting this challenge may lie in our understanding of plant photosynthetic adaptations and water use efficiency. Plants from various taxa have evolved crassulacean acid metabolism (CAM), a water-conserving adaptation of photosynthetic carbon dioxide fixation that enables plants to thrive under semi-arid or seasonally drought-prone conditions. Although past research on CAM has led to a better understanding of the inner workings of plant resilience and adaptation to stress, successful introduction of this pathway into C3 or C4 plants has not been reported. The recent revolution in molecular, systems, and synthetic biology, as well as innovations in high-throughput data generation and mining, creates new opportunities to uncover the minimum genetic tool kit required to introduce CAM traits into drought-sensitive crops. Here, we propose four complementary research avenues to uncover this tool kit. First, genomes and computational methods should be used to improve understanding of the nature of variations that drive CAM evolution. Second, single-cell 'omics technologies offer the possibility for in-depth characterization of the mechanisms that trigger environmentally controlled CAM induction. Third, the rapid increase in new 'omics data enables a comprehensive, multimodal exploration of CAM. Finally, the expansion of functional genomics methods is paving the way for integration of CAM into farming systems.

The Calvin-Benson-Bassham cycle in C4 and Crassulacean acid metabolism species

The Calvin-Benson-Bassham (CBB) cycle is the ancestral CO2 assimilation pathway and is found in all photosynthetic organisms. Biochemical extensions to the CBB cycle have evolved that allow the resulting pathways to act as CO2 concentrating mechanisms, either spatially in the case of C4 photosynthesis or temporally in the case of Crassulacean acid metabolism (CAM). While the biochemical steps in the C4 and CAM pathways are known, questions remain on their integration and regulation with CBB cycle activity. The application of omic and transgenic technologies is providing a more complete understanding of the biochemistry of C4 and CAM species and will also provide insight into the CBB cycle in these plants. As the global population increases, new solutions are required to increase crop yields and meet demands for food and other bioproducts. Previous work in C3 species has shown that increasing carbon assimilation through genetic manipulation of the CBB cycle can increase biomass and yield. There may also be options to improve photosynthesis in species using C4 photosynthesis and CAM through manipulation of the CBB cycle in these plants. This is an underexplored strategy and requires more basic knowledge of CBB cycle operation in these species to enable approaches for increased productivity.

Identification and Analysis of PEPC Gene Family Reveals Functional Diversification in Orchidaceae and the Regulation of Bacterial-Type PEPC

Phosphoenolpyruvate carboxylase (PEPC) gene family plays a crucial role in both plant growth and response to abiotic stress. Approximately half of the Orchidaceae species are estimated to perform CAM pathway, and the availability of sequenced orchid genomes makes them ideal subjects for investigating the PEPC gene family in CAM plants. In this study, a total of 33 PEPC genes were identified across 15 orchids. Specifically, one PEPC gene was found in Cymbidium goeringii and Platanthera guangdongensis; two in Apostasia shenzhenica, Dendrobium chrysotoxum, D. huoshanense, Gastrodia elata, G. menghaiensis, Phalaenopsis aphrodite, Ph. equestris, and Pl. zijinensis; three in C. ensifolium, C. sinense, D. catenatum, D. nobile, and Vanilla planifolia. These PEPC genes were categorized into four subgroups, namely PEPC-i, PEPC-ii, and PEPC-iii (PTPC), and PEPC-iv (BTPC), supported by the comprehensive analyses of their physicochemical properties, motif, and gene structures. Remarkably, PEPC-iv contained a heretofore unreported orchid PEPC gene, identified as VpPEPC4. Differences in the number of PEPC homolog genes among these species were attributed to segmental duplication, whole-genome duplication (WGD), or gene loss events. Cis-elements identified in promoter regions were predominantly associated with light responsiveness, and circadian-related elements were observed in each PEPC-i and PEPC-ii gene. The expression levels of recruited BTPC, VpPEPC4, exhibited a lower expression level than other VpPEPCs in the tested tissues. The expression analyses and RT-qPCR results revealed diverse expression patterns in orchid PEPC genes. Duplicated genes exhibited distinct expression patterns, suggesting functional divergence. This study offered a comprehensive analysis to unveil the evolution and function of PEPC genes in Orchidaceae.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Review: Structural-functional relationship of WRKY transcription factors: Unfolding the role of WRKY in plants

WRKY as the name suggests, are the transcription factors (TFs) that contain the signature WRKY domains, hence named after it. Since their discovery in 1994, they have been well studied in plants with exploration of approximately 74 WRKY genes in the model plant, Arabidopsis alone. However, the study of these transcription factors (TFs) is not just limited to model plant now. They have been studied widely in crop plants as well, because of their tremendous contribution in stress as well as in growth and development. Here, in this review, we describe the story of WRKY TFs from their identification to their origin, the binding mechanisms, structure and their contribution in regulating plant development and stress physiology. High throughput transcriptomics-based data also opened a doorway to understand the comprehensive and detailed functioning of WRKY TFs in plants. Indeed, the detailed functional role of each and every WRKY member in regulating the gene expression is required to pave the path to develop holistic understanding of their role in stress physiology and developmental processes in plants.

Genome-Wide Identification and Characterization of Gibberellic Acid-Stimulated Arabidopsis Gene Family in Pineapple (Ananas comosus)

The gibberellic acid-stimulated Arabidopsis (GASA) gene family plays a crucial role in growth, development, and stress response, and it is specific to plants. This gene family has been extensively studied in various plant species, and its functional role in pineapple has yet to be characterized. In this study, 15 AcGASA genes were identified in pineapple through a genome-wide scan and categorized into three major branches based on a phylogenetic tree. All AcGASA proteins share a common structural domain with 12 cysteine residues, but they exhibit slight variations in their physicochemical properties and motif composition. Predictions regarding subcellular localization suggest that AcGASA proteins are present in the cell membrane, Golgi apparatus, nucleus, and cell wall. An analysis of gene synteny indicated that both tandem and segmental repeats have a significant impact on the expansion of the AcGASA gene family. Our findings demonstrate the differing regulatory effects of these hormones (GA, NAA, IAA, MeJA, and ABA) on the AcGASA genes. We analyzed the expression profiles of GASA genes in different pineapple tissue parts, and the results indicated that AcGASA genes exhibit diverse expression patterns during the development of different plant tissues, particularly in the regulation of floral organ development. This study provides a comprehensive understanding of GASA family genes in pineapple. It serves as a valuable reference for future studies on the functional characterization of GASA genes in other perennial herbaceous plants.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

The CAM lineages of planet Earth

BACKGROUND AND SCOPE

The growth of experimental studies of crassulacean acid metabolism (CAM) in diverse plant clades, coupled with recent advances in molecular systematics, presents an opportunity to re-assess the phylogenetic distribution and diversity of species capable of CAM. It has been more than two decades since the last comprehensive lists of CAM taxa were published, and an updated survey of the occurrence and distribution of CAM taxa is needed to facilitate and guide future CAM research. We aimed to survey the phylogenetic distribution of these taxa, their diverse morphology, physiology and ecology, and the likely number of evolutionary origins of CAM based on currently known lineages.

RESULTS AND CONCLUSIONS

We found direct evidence (in the form of experimental or field observations of gas exchange, day-night fluctuations in organic acids, carbon isotope ratios and enzymatic activity) for CAM in 370 genera of vascular plants, representing 38 families. Further assumptions about the frequency of CAM species in CAM clades and the distribution of CAM in the Cactaceae and Crassulaceae bring the currently estimated number of CAM-capable species to nearly 7 % of all vascular plants. The phylogenetic distribution of these taxa suggests a minimum of 66 independent origins of CAM in vascular plants, possibly with dozens more. To achieve further insight into CAM origins, there is a need for more extensive and systematic surveys of previously unstudied lineages, particularly in living material to identify low-level CAM activity, and for denser sampling to increase phylogenetic resolution in CAM-evolving clades. This should allow further progress in understanding the functional significance of this pathway by integration with studies on the evolution and genomics of CAM in its many forms.

Lessons from the history of Agave: ecological and cultural context for valuation of CAM

BACKGROUND AND SCOPE

Crassulacean acid metabolism (CAM) is an intriguing physiological adaptation in plants that are widespread throughout many ecosystems. Despite the relatively recent mechanistic understanding of CAM in plant physiology, evidence from historical records suggests that ancient cultures in the Americas also recognized the value of CAM plants. Agave species, in particular, have a rich cultural legacy that provides a foundation for commercially valued products. Here, we review that legacy and potential relationships between ancient values and the needs of modern-day climate adaptation strategies.

CONCLUSIONS

There are many products that can be produced from Agave species, including food, sugar, fibre and medicines. Traditional knowledge about agricultural management and preparation of plant products can be combined with new ecophysiological knowledge and agronomic techniques to develop these resources in the borderland region of the southwestern USA and Mexico. Historical records of pre-Columbian practices in the Sonoran desert and remnants of centuries-old agriculture in Baja California and Sonora demonstrate the climate resilience of Agave agriculture. Commercial growth of both tequila and bacanora indicates the potential for large-scale production today, but also underscores the importance of adopting regenerative agricultural practices to accomplish environmentally sustainable production. Recent international recognition of the Appellation of Origin for several Agave species produced for spirits in Mexico might provide opportunities for agricultural diversification. In contrast, fibre is currently produced from several Agave species on many continents. Projections of growth with future climate change suggest that Agave spp. will be viable alternatives for commodity crops that suffer declines during drought and increased temperatures. Historical cultivation of Agave affirms that these CAM plants can supply sugar, soft and hard fibres, medicines and food supplements.

Prospects and perspectives: inferring physiological and regulatory targets for CAM from molecular and modelling approaches

BACKGROUND AND SCOPE

This review summarizes recent advances in our understanding of Crassulacean Acid Metabolism (CAM) by integrating evolutionary, ecological, physiological, metabolic and molecular perspectives. A number of key control loops which moderate the expression of CAM phases, and their metabolic and molecular control, are explored. These include nocturnal stomatal opening, activation of phosphoenolpyruvate carboxylase by a specific protein kinase, interactions with circadian clock control, as well as daytime decarboxylation and activation of Rubisco. The vacuolar storage and release of malic acid and the interplay between the supply and demand for carbohydrate reserves are also key metabolic control points.

FUTURE OPPORTUNITIES

We identify open questions and opportunities, with experimentation informed by top-down molecular modelling approaches allied with bottom-up mechanistic modelling systems. For example, mining transcriptomic datasets using high-speed systems approaches will help to identify targets for future genetic manipulation experiments to define the regulation of CAM (whether circadian or metabolic control). We emphasize that inferences arising from computational approaches or advanced nuclear sequencing techniques can identify potential genes and transcription factors as regulatory targets. However, these outputs then require systematic evaluation, using genetic manipulation in key model organisms over a developmental progression, combining gene silencing and metabolic flux analysis and modelling to define functionality across the CAM day-night cycle. From an evolutionary perspective, the origins and function of CAM succulents and responses to water deficits are set against the mesophyll and hydraulic limitations imposed by cell and tissue succulence in contrasting morphological lineages. We highlight the interplay between traits across shoots (3D vein density, mesophyll conductance and cell shrinkage) and roots (xylem embolism and segmentation). Thus, molecular, biophysical and biochemical processes help to curtail water losses and exploit rapid rehydration during restorative rain events. In the face of a changing climate, we hope such approaches will stimulate opportunities for future research.

The Agavoideae: an emergent model clade for CAM evolutionary biology

Crassulacean acid metabolism - or CAM photosynthesis - was described in the early to mid-20th century, and our understanding of this metabolic pathway was later expanded upon through detailed biochemical analyses of carbon balance. Soon after, scientists began to study the ecophysiological implications of CAM, and a large part of this early work was conducted in the genus Agave, in the subfamily Agavoideae of the family Asparagaceae. Today, the Agavoideae continues to be important for the study of CAM photosynthesis, from the ecophysiology of CAM species, to the evolution of the CAM phenotype and to the genomics underlying CAM traits. Here we review past and current work on CAM in the Agavoideae, in particular highlighting the work of Park Nobel in Agave, and focusing on the powerful comparative system the Agavoideae has become for studying the origins of CAM. We also highlight new genomics research and the potential for studying intraspecific variation within species of the Agavoideae, particularly species in the genus Yucca. The Agavoideae has served as an important model clade for CAM research for decades, and undoubtedly will continue to help push our understanding of CAM biology and evolution in the future.

Reconciling continuous and discrete models of C4 and CAM evolution

BACKGROUND

A current argument in the CAM biology literature has focused on the nature of the CAM evolutionary trajectory: whether there is a smooth continuum of phenotypes between plants with C3 and CAM photosynthesis or whether there are discrete steps of phenotypic evolutionary change such as has been modelled for the evolution of C4 photosynthesis. A further implication is that a smooth continuum would increase the evolvability of CAM, whereas discrete changes would make the evolutionary transition from C3 to CAM more difficult.

SCOPE

In this essay, I attempt to reconcile these two viewpoints, because I think in many ways this is a false dichotomy that is constraining progress in understanding how both CAM and C4 evolved. In reality, the phenotypic space connecting C3 species and strong CAM/C4 species is both a continuum of variably expressed quantitative traits and yet also contains certain combinations of traits that we are able to identify as discrete, recognizable phenotypes. In this sense, the evolutionary mechanics of CAM origination are no different from those of C4 photosynthesis, nor from the evolution of any other complex trait assemblage.

CONCLUSIONS

To make progress, we must embrace the concept of discrete phenotypic phases of CAM evolution, because their delineation will force us to articulate what aspects of phenotypic variation we think are significant. There are some current phenotypic gaps that are limiting our ability to build a complete CAM evolutionary model: the first is how a rudimentary CAM biochemical cycle becomes established, and the second is how the 'accessory' CAM cycle in C3+CAM plants is recruited into a primary metabolism. The connections to the C3 phenotype we are looking for are potentially found in the behaviour of C3 plants when undergoing physiological stress - behaviour that, strangely enough, remains essentially unexplored in this context.

Developmental dynamics of crassulacean acid metabolism (CAM) in Opuntia ficus-indica

BACKGROUND AND AIMS

The relative contributions of C3 photosynthesis and crassulacean acid metabolism (CAM) during the earliest stages of development were investigated to assess how much each might contribute to cactus pear (Opuntia ficus-indica) productivity.

METHODS

The developmental progression of C3 photosynthesis and CAM was assessed in seedlings and daughter cladodes of mature plants by titratable acidity, δ13C isotopic values and diel gas exchange measurements.

KEY RESULTS

Nocturnal acidification was observed in seedling cladodes and cotyledons at the earliest stages of development and became highly significant by 75 days of development. Seedling cotyledons showed mean δ13C values of -21.4 and -17.1 ‰ at 30 and 100 days of age, respectively. Seedling cladodes showed mean δ13C values of -19.4 and -14.5 ‰ at 30 and 100 days of age, respectively. These values are typical of CAM plants. Net CO2 assimilation was negative, then occurred in both the day and the night, with nighttime fixation becoming predominant once the primary cladode reached 5 cm in size. Emergent daughter cladodes growing on mature plants showed nocturnal titratable acidity at the earliest stages of development, which became significant when daughter cladodes were >2.5-5 cm in height. Emergent daughter cladodes showed mean δ13C values of -14.5 to -15.6 ‰, typical of CAM plants. CO2 assimilation studies revealed that net CO2 uptake was negative in daughter cladodes <12 cm in length, but then exhibited net positive CO2 assimilation in both the day and the night, with net nocturnal CO2 assimilation predominating once the daughter cladode grew larger.

CONCLUSIONS

Developing O. ficus-indica primary and daughter cladodes begin as respiring sink tissues that transition directly to performing CAM once net positive CO2 fixation is observed. Overall, these results demonstrate that CAM is the primary form of photosynthetic carbon assimilation for O. ficus-indica even at the earliest stages of seedling or daughter cladode development.

Comparative transcriptome analysis reveals candidate genes for cold stress response and early flowering in pineapple

Pineapple originates from tropical regions in South America and is therefore significantly impacted by cold stress. Periodic cold events in the equatorial regions where pineapple is grown may induce early flowering, also known as precocious flowering, resulting in monetary losses due to small fruit size and the need to make multiple passes for harvesting a single field. Currently, pineapple is one of the most important tropical fruits in the world in terms of consumption, and production losses caused by weather can have major impacts on worldwide exportation potential and economics. To further our understanding of and identify mechanisms for low-temperature tolerance in pineapple, and to identify the relationship between low-temperature stress and flowering time, we report here a transcriptomic analysis of two pineapple genotypes in response to low-temperature stress. Using meristem tissue collected from precocious flowering-susceptible MD2 and precocious flowering-tolerant Dole-17, we performed pairwise comparisons and weighted gene co-expression network analysis (WGCNA) to identify cold stress, genotype, and floral organ development-specific modules. Dole-17 had a greater increase in expression of genes that confer cold tolerance. The results suggested that low temperature stress in Dole-17 plants induces transcriptional changes to adapt and maintain homeostasis. Comparative transcriptomic analysis revealed differences in cuticular wax biosynthesis, carbohydrate accumulation, and vernalization-related gene expression between genotypes. Cold stress induced changes in ethylene and abscisic acid-mediated pathways differentially between genotypes, suggesting that MD2 may be more susceptible to hormone-mediated early flowering. The differentially expressed genes and module hub genes identified in this study are potential candidates for engineering cold tolerance in pineapple to develop new varieties capable of maintaining normal reproduction cycles under cold stress. In addition, a total of 461 core genes involved in the development of reproductive tissues in pineapple were also identified in this study. This research provides an important genomic resource for understanding molecular networks underlying cold stress response and how cold stress affects flowering time in pineapple.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Differential gene expression during floral transition in pineapple

Pineapple (Ananas comosus var. comosus) and ornamental bromeliads are commercially induced to flower by treatment with ethylene or its analogs. The apex is transformed from a vegetative to a floral meristem and shows morphological changes in 8 to 10 days, with flowers developing 8 to 10 weeks later. During eight sampling stages ranging from 6 h to 8 days after treatment, 7961 genes were found to exhibit differential expression (DE) after the application of ethylene. In the first 3 days after treatment, there was little change in ethylene synthesis or in the early stages of the ethylene response. Subsequently, three ethylene response transcription factors (ERTF) were up-regulated and the potential gene targets were predicted to be the positive flowering regulator CONSTANS-like 3 (CO), a WUSCHEL gene, two APETALA1/FRUITFULL (AP1/FUL) genes, an epidermal patterning gene, and a jasmonic acid synthesis gene. We confirm that pineapple has lost the flowering repressor FLOWERING LOCUS C. At the initial stages, the SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) was not significantly involved in this transition. Another WUSCHEL gene and a PHD homeobox transcription factor, though not apparent direct targets of ERTF, were up-regulated within a day of treatment, their predicted targets being the up-regulated CO, auxin response factors, SQUAMOSA, and histone H3 genes with suppression of abscisic acid response genes. The FLOWERING LOCUS T (FT), TERMINAL FLOWER (TFL), AGAMOUS-like APETELAR (AP2), and SEPETALA (SEP) increased rapidly within 2 to 3 days after ethylene treatment. Two FT genes were up-regulated at the apex and not at the leaf bases after treatment, suggesting that transport did not occur. These results indicated that the ethylene response in pineapple and possibly most bromeliads act directly to promote the vegetative to flower transition via APETALA1/FRUITFULL (AP1/FUL) and its interaction with SPL, FT, TFL, SEP, and AP2. A model based on AP2/ERTF DE and predicted DE target genes was developed to give focus to future research. The identified candidate genes are potential targets for genetic manipulation to determine their molecular role in flower transition.

Genomic analysis of the emergent aquatic plant Sparganium stoloniferum provides insights into its clonality, local adaptation and demographic history

Clonal propagation and extensive dispersal of seeds and asexual propagules are two important features of aquatic plants that help them adapt to aquatic environments. Accurate measurements of clonality and effective clonal dispersal are essential for understanding the evolution of aquatic plants. Here, we first assembled a high-quality chromosome-level genome of a widespread emergent aquatic plant Sparganium stoloniferum to provide a reference for its population genomic study. We then performed high-depth resequencing of 173 individuals from 20 populations covering different basins across its range in China. Population genomic analyses revealed three genetic lineages reflecting the northeast (NE), southwest (SW) and northwest (NW) of its geographical distribution. The NE lineage diverged in the middle Pleistocene while the SW and NW lineages diverged until about 2400 years ago. Clonal relationship analyses identified nine populations as monoclonal population. Dispersal of vegetative propagules was identified between five populations covering three basins in the NE lineage, and dispersal distance was up to 1041 km, indicating high dispersibility in emergent aquatic plant species. We also identified lineage-specific positively selected genes that are likely to be involved in adaptations to saline wetlands and high-altitude environments. Our findings accurately measure the clonality, determine the dispersal range and frequency of vegetative propagules, and detect genetic signatures of local adaptation in a widespread emergent aquatic plant species, providing new perspectives on the evolution of aquatic plants.

Diel dynamics of multi-omics in elkhorn fern provide new insights into weak CAM photosynthesis

Crassulacean acid metabolism (CAM) has high water-use efficiency (WUE) and is widely recognized to have evolved from C3 photosynthesis. Different plant lineages have convergently evolved CAM, but the molecular mechanism that underlies C3-to-CAM evolution remains to be clarified. Platycerium bifurcatum (elkhorn fern) provides an opportunity to study the molecular changes underlying the transition from C3 to CAM photosynthesis because both modes of photosynthesis occur in this species, with sporotrophophyll leaves (SLs) and cover leaves (CLs) performing C3 and weak CAM photosynthesis, respectively. Here, we report that the physiological and biochemical attributes of CAM in weak CAM-performing CLs differed from those in strong CAM species. We investigated the diel dynamics of the metabolome, proteome, and transcriptome in these dimorphic leaves within the same genetic background and under identical environmental conditions. We found that multi-omic diel dynamics in P. bifurcatum exhibit both tissue and diel effects. Our analysis revealed temporal rewiring of biochemistry relevant to the energy-producing pathway (TCA cycle), CAM pathway, and stomatal movement in CLs compared with SLs. We also confirmed that PHOSPHOENOLPYRUVATE CARBOXYLASE KINASE (PPCK) exhibits convergence in gene expression among highly divergent CAM lineages. Gene regulatory network analysis identified candidate transcription factors regulating the CAM pathway and stomatal movement. Taken together, our results provide new insights into weak CAM photosynthesis and new avenues for CAM bioengineering.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially βCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

Rethinking the potential productivity of crassulacean acid metabolism by integrating metabolic dynamics with shoot architecture, using the example of Agave tequilana

Terrestrial CAM plants typically occur in hot semiarid regions, yet can show high crop productivity under favorable conditions. To achieve a more mechanistic understanding of CAM plant productivity, a biochemical model of diel metabolism was developed and integrated with 3-D shoot morphology to predict the energetics of light interception and photosynthetic carbon assimilation. Using Agave tequilana as an example, this biochemical model faithfully simulated the four diel phases of CO2 and metabolite dynamics during the CAM rhythm. After capturing the 3-D form over an 8-yr production cycle, a ray-tracing method allowed the prediction of the light microclimate across all photosynthetic surfaces. Integration with the biochemical model thereby enabled the simulation of plant and stand carbon uptake over daily and annual courses. The theoretical maximum energy conversion efficiency of Agave spp. is calculated at 0.045-0.049, up to 7% higher than for C3 photosynthesis. Actual light interception, and biochemical and anatomical limitations, reduced this to 0.0069, or 15.6 Mg ha-1  yr-1 dry mass annualized over an 8-yr cropping cycle, consistent with observation. This is comparable to the productivity of many C3 crops, demonstrating the potential of CAM plants in climates where little else may be grown while indicating strategies that could raise their productivity.

Chromosome-scale genome sequence of Suaeda glauca sheds light on salt stress tolerance in halophytes

Soil salinity is a growing concern for global crop production and the sustainable development of humanity. Therefore, it is crucial to comprehend salt tolerance mechanisms and identify salt-tolerance genes to enhance crop tolerance to salt stress. Suaeda glauca, a halophyte species well adapted to the seawater environment, possesses a unique ability to absorb and retain high salt concentrations within its cells, particularly in its leaves, suggesting the presence of a distinct mechanism for salt tolerance. In this study, we performed de novo sequencing of the S. glauca genome. The genome has a size of 1.02 Gb (consisting of two sets of haplotypes) and contains 54 761 annotated genes, including alleles and repeats. Comparative genomic analysis revealed a strong synteny between the genomes of S. glauca and Beta vulgaris. Of the S. glauca genome, 70.56% comprises repeat sequences, with retroelements being the most abundant. Leveraging the allele-aware assembly of the S. glauca genome, we investigated genome-wide allele-specific expression in the analyzed samples. The results indicated that the diversity in promoter sequences might contribute to consistent allele-specific expression. Moreover, a systematic analysis of the ABCE gene families shed light on the formation of S. glauca's flower morphology, suggesting that dysfunction of A-class genes is responsible for the absence of petals in S. glauca. Gene family expansion analysis demonstrated significant enrichment of Gene Ontology (GO) terms associated with DNA repair, chromosome stability, DNA demethylation, cation binding, and red/far-red light signaling pathways in the co-expanded gene families of S. glauca and S. aralocaspica, in comparison with glycophytic species within the chenopodium family. Time-course transcriptome analysis under salt treatments revealed detailed responses of S. glauca to salt tolerance, and the enrichment of the transition-upregulated genes in the leaves associated with DNA repair and chromosome stability, lipid biosynthetic process, and isoprenoid metabolic process. Additionally, genome-wide analysis of transcription factors indicated a significant expansion of FAR1 gene family. However, further investigation is needed to determine the exact role of the FAR1 gene family in salt tolerance in S. glauca.

Defining Mechanisms of C3 to CAM Photosynthesis Transition toward Enhancing Crop Stress Resilience

Global climate change and population growth are persistently posing threats to natural resources (e.g., freshwater) and agricultural production. Crassulacean acid metabolism (CAM) evolved from C3 photosynthesis as an adaptive form of photosynthesis in hot and arid regions. It features the nocturnal opening of stomata for CO2 assimilation, diurnal closure of stomata for water conservation, and high water-use efficiency. To cope with global climate challenges, the CAM mechanism has attracted renewed attention. Facultative CAM is a specialized form of CAM that normally employs C3 or C4 photosynthesis but can shift to CAM under stress conditions. It not only serves as a model for studying the molecular mechanisms underlying the CAM evolution, but also provides a plausible solution for creating stress-resilient crops with facultative CAM traits. This review mainly discusses the recent research effort in defining the C3 to CAM transition of facultative CAM plants, and highlights challenges and future directions in this important research area with great application potential.

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

Acclimation to a combination of water deficit and nutrient deprivation through simultaneous increases in abscisic acid and bioactive jasmonates in the succulent plant Sempervivum tectorum L

Activation of hormonal responses defines the drought acclimation ability of plants and may condition their survival. However, aside ABA, little is known about the possible contribution of other phytohormones, such as jasmonates and salicylates, in the response of CAM plants to water deficit. Here, we aimed to study the physiological mechanisms underlying the stress tolerance of house leek (Sempervivum tectorum L.), a CAM plant adapted to survive harsh environments, to a combination of water deficit and nutrient deprivation. We exposed plants to the combination of these two abiotic stresses by withholding nutrient solution for 10 weeks and monitored their physiological response every two weeks by measuring various stress makers together with the accumulation of stress-related phytohormones and photoprotective molecules, such as tocopherols (vitamin E). Results showed that ABA content increased by 4.2-fold after four weeks of water deficit to keep later constant up to 10 weeks of stress, variations that occurred concomitantly with reductions in the relative leaf water content, which decreased by up to 20% only. The bioactive jasmonate, jasmonoyl-isoleucine was the other stress-related phytohormone that simultaneously increased under stress together with ABA. While contents of salicylic acid and the jasmonoyl-isoleucine precursors, 12-oxo-phytodienoic acid and jasmonic acid decreased with water deficit, those of jasmonoyl-isoleucine increased 3.6-fold at four weeks of stress. The contents of ABA and jasmonoyl-isoleucine correlated positively between them and with the content of α-tocopherol per unit of chlorophyll, thus suggesting a photoprotective activation role. It is concluded that S. tectorum not only withstands a combination of water deficit and nutrient deprivation for 10 weeks without any symptom of damage but also activates effective defense strategies through the simultaneous accumulation of ABA and the bioactive jasmonate form, jasmonoyl-isoleucine.

Identification of bromelain subfamily proteases encoded in the pineapple genome

Papain (aka C1A) family proteases, including bromelain enzymes, are widespread across the plant kingdom and play critical regulatory functions in protein turnover during development. The proteolytic activity exhibited by papain family proteases has led to their increased usage for a wide range of cosmetic, therapeutic, and medicinal purposes. Bromelain enzymes, or bromelains in short, are members of the papain family that are specific to the bromeliad plant family. The only major commercial extraction source of bromelain is pineapple. The importance of C1A family and bromelain subfamily proteases in pineapple development and their increasing economic importance led several researchers to utilize available genomic resources to identify protease-encoding genes in the pineapple genome. To date, studies are lacking in screening bromelain genes for targeted use in applied science studies. In addition, the bromelain genes coding for the enzymes present in commercially available bromelain products have not been identified and their evolutionary origin has remained unclear. Here, using the newly developed MD2 v2 pineapple genome, we aimed to identify bromelain-encoding genes and elucidate their evolutionary origin. Orthologous and phylogenetic analyses of all papain-family proteases encoded in the pineapple genome revealed a single orthogroup (189) and phylogenetic clade (XIII) containing the bromelain subfamily. Duplication mode and synteny analyses provided insight into the origin and expansion of the bromelain subfamily in pineapple. Proteomic analysis identified four bromelain enzymes present in two commercially available bromelain products derived from pineapple stem, corresponding to products of four putative bromelain genes. Gene expression analysis using publicly available transcriptome data showed that 31 papain-family genes identified in this study were up-regulated in specific tissues, including stem, fruit, and floral tissues. Some of these genes had higher expression in earlier developmental stages of different tissues. Similar expression patterns were identified by RT-qPCR analysis with leaf, stem, and fruit. Our results provide a strong foundation for future applicable studies on bromelain, such as transgenic approaches to increase bromelain content in pineapple, development of bromelain-producing bioreactors, and studies that aim to determine the medicinal and/or therapeutic viability of individual bromelain enzymes.

In-silico, evolutionary, and functional analysis of CHUP1 and its related proteins in Bienertia sinuspersici-a comparative study across C3, C4, CAM, and SCC4 model plants

Single-cell C₄ (SCC₄) plants with bienertioid anatomy carry out photosynthesis in a single cell. Chloroplast movement is the underlying phenomenon, where chloroplast unusual positioning 1 (CHUP1) plays a key role. This study aimed to characterize CHUP1 and CHUP1-like proteins in an SCC₄ photosynthetic plant, Bienertia sinuspersici. Also, a comparative analysis of SCC4 CHUP1 was made with C₃, C₄, and CAM model plants including an extant basal angiosperm, Amborella. The CHUP1 gene exists as a single copy from the basal angiosperms to SCC₄ plants. Our analysis identified that Chenopodium quinoa, a recently duplicated allotetraploid, has two copies of CHUP1. In addition, the numbers of CHUP1-like and its associated proteins such as CHUP1-like_a, CHUP1-like_b, HPR, TPR, and ABP varied between the species. Hidden Markov Model analysis showed that the gene size of CHUP1-like_a and CHUP1-like_b of SCC₄ species, Bienertia, and Suaeda were enlarged than other plants. Also, we identified that CHUP1-like_a and CHUP1-like_b are absent in Arabidopsis and Amborella, respectively. Motif analysis identified several conserved and variable motifs based on the orders (monocot and dicot) as well as photosynthetic pathways. For instance, CAM plants such as pineapple and cactus shared certain motifs of CHUP1-like_a irrespective of their distant phylogenetic relationship. The free ratio model showed that CHUP1 maintained purifying selection, whereas CHUP1-like_a and CHUP1-like_b have adaptive functions between SCC₄ plants and quinoa. Similarly, rice and maize branches displayed functional diversification on CHUP1-like_b. Relative gene expression data showed that during the subcellular compartmentalization process of Bienertia, CHUP1 and actin-binding proteins (ABP) genes showed a similar pattern of expression. Altogether, the results of this study provide insight into the evolutionary and functional details of CHUP1 and its associated proteins in the development of the SCC₄ system in comparison with other C₃, C₄, and CAM model plants.

Genome-wide mapping of DNase I hypersensitive sites in pineapple leaves

Pineapple [Ananas comosus (L.) Merr.] is the most economically important crop possessing crassulacean acid metabolism (CAM) photosynthesis which has a higher water use efficiency by control of nocturnal opening and diurnal closure of stomata. To provide novel insights into the diel regulatory landscape in pineapple leaves, we performed genome-wide mapping of DNase I hypersensitive sites (DHSs) in pineapple leaves at day (2a.m.) and night (10a.m.) using a simplified DNase-seq method. As a result, totally 33340 and 28753 DHSs were found in green-tip tissue, and 29597 and 40068 were identified in white-base tissue at 2a.m. and 10a.m., respectively. We observed that majority of the pineapple genes occupied less than two DHSs with length shorter than 1 kb, and the promotor DHSs showed a proximal trend to the transcription start site (>77% promotor DHSs within 1 kb). In addition, more intergenic DHSs were identified around transcription factors or transcription co-regulators (TFs/TCs) than other functional genes, indicating complex regulatory contexts around TFs/TCs. Through combined analysis of tissue preferential DHSs and genes, we respectively found 839 and 888 coordinately changed genes in green-tip at 2a.m. and 10a.m. (AcG2 and AcG10). Furthermore, AcG2-specific, AcG10-specific and common accessible DHSs were dissected from the total photosynthetic preferential DHSs, and the regulatory networks indicated dynamic regulations with multiple cis-regulatory elements occurred to genes preferentially expressed in photosynthetic tissues. Interestingly, binding motifs of several cycling TFs were identified in the DHSs of key CAM genes, revealing a circadian regulation to CAM coordinately diurnal expression. Our results provide a chromatin regulatory landscape in pineapple leaves during the day and night. This will provide important information to assist with deciphering the circadian regulation of CAM photosynthesis.

Unveiling the molecular mechanism involving anthocyanins in pineapple peel discoloration during fruit maturation

Peel color is a key factor that affects the fruit's aesthetic and economic values. Limited knowledge is available on the regulation of pineapple peel discoloration. Here, we report that a decrease in anthocyanin biosynthesis, particularly cyanidin, is predominantly associated with the pineapple peel color change during maturation. The findings suggest that the changes in the expression of key structural genes (early and late biosynthetic genes) of the anthocyanin (cyanidin) biosynthesis pathway are responsible for peel discoloration. Based on a gene co-expression analysis and a transient expression, two transcription factors i.e., AcHOX21 and AcMYB12, were identified, whose' downregulation leads to reduced anthocyanin accumulation with fruit maturation. The endogenous levels of jasmonic acid, gibberellic acid, and auxins are also involved in anthocyanin-content-led peel discoloration. Overall, the discovery of genes regulating anthocyanin biosynthesis in pineapple peel provides a theoretical basis for improving the fruit's aesthetic value through genetic engineering.

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.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Genomic and transcriptomic analysis of sacred fig (Ficus religiosa)

BACKGROUND

Peepal/Bodhi tree (Ficus religiosa L.) is an important, long-lived keystone ecological species. Communities on the Indian subcontinent have extensively employed the plant in Ayurveda, traditional medicine, and spiritual practices. The Peepal tree is often thought to produce oxygen both during the day and at night by Indian folks. The goal of our research was to produce molecular resources using whole-genome and transcriptome sequencing techniques.

RESULTS

The complete genome of the Peepal tree was sequenced using two next-generation sequencers Illumina HiSeq1000 and MGISEQ-2000. We assembled the draft genome of 406 Mb, using a hybrid assembly workflow. The genome annotation resulted in 35,093 protein-coding genes; 53% of its genome consists of repetitive sequences. To understand the physiological pathways in leaf tissues, we analyzed photosynthetically distinct conditions: bright sunny days and nights. The RNA-seq analysis supported the expression of 26,479 unigenes. The leaf transcriptomic analysis of the diurnal and nocturnal periods revealed the expression of the significant number of genes involved in the carbon-fixation pathway.

CONCLUSIONS

This study presents a draft hybrid genome assembly for F. religiosa and its functional annotated genes. The genomic and transcriptomic data-derived pathways have been analyzed for future studies on the Peepal tree.

Small Auxin Up RNA (SAUR) gene family identification and functional genes exploration during the floral organ and fruit developmental stages in pineapple (Ananas comosus L.) and its response to salinity and drought stresses

In plants, sexual reproduction relies on the proper development of floral organs that facilitate the successful development of fruits and seeds. Auxin responsive small auxin-up RNA (SAUR) genes play essential roles in floral organ formation and fruit development. However, little is known about the role of SAUR genes in pineapple floral organ formation and fruit development as well as stress responses. In this study, based on genome information and transcriptome datasets, 52 AcoSAUR genes were identified and grouped into 12 groups. The gene structure analysis revealed that most AcoSAUR genes did not have introns, although auxin-acting elements were abundant in the promoter region of AcoSAUR members. The expression analysis across the multiple flower and fruit development stages revealed differential expression of AcoSAUR genes, indicating a tissue and stage-specific function of AcoSAURs. Correlation analysis and pairwise comparisons between gene expression and tissue specificity identified stamen-, petal-, ovule-, and fruit-specific AcoSAURs involved in pineapple floral organs (AcoSAUR4/5/15/17/19) and fruit development (AcoSAUR6/11/36/50). RT-qPCR analysis revealed that AcoSAUR12/24/50 played positive roles in response to the salinity and drought treatment. This work provides an abundant genomic resource for functional analysis of AcoSAUR genes during the pineapple floral organs and fruit development stages. It also highlights the role of auxin signaling involved in pineapple reproductive organ growth.

Virus-Induced Gene Silencing (VIGS): A Powerful Tool for Crop Improvement and Its Advancement towards Epigenetics

Virus-induced gene silencing (VIGS) is an RNA-mediated reverse genetics technology that has evolved into an indispensable approach for analyzing the function of genes. It downregulates endogenous genes by utilizing the posttranscriptional gene silencing (PTGS) machinery of plants to prevent systemic viral infections. Based on recent advances, VIGS can now be used as a high-throughput tool that induces heritable epigenetic modifications in plants through the viral genome by transiently knocking down targeted gene expression. As a result of the progression of DNA methylation induced by VIGS, new stable genotypes with desired traits are being developed in plants. In plants, RNA-directed DNA methylation (RdDM) is a mechanism where epigenetic modifiers are guided to target loci by small RNAs, which play a major role in the silencing of the target gene. In this review, we described the molecular mechanisms of DNA and RNA-based viral vectors and the knowledge obtained through altering the genes in the studied plants that are not usually accessible to transgenic techniques. We showed how VIGS-induced gene silencing can be used to characterize transgenerational gene function(s) and altered epigenetic marks, which can improve future plant breeding programs.

High nucleotide similarity of three Copia lineage LTR retrotransposons among plant genomes

Transposable elements (TEs) are mobile elements found in the majority of eukaryotic genomes. TEs deeply impact the structure and evolution of chromosomes and can induce mutations affecting coding genes. In plants, the major group of TEs is long terminal repeat retrotransposons (LTR-RTs). They are classified into superfamilies (Gypsy, Copia) and subclassified into lineages. Horizontal transfer (HT), defined as the nonsexual transmission of genetic material between species, is a process allowing LTR-RTs to invade a new genome. Although this phenomenon was considered rare, recent studies demonstrate numerous transfers of LTR-RTs. This study aims to determine which LTR-RT lineages are shared with high similarity among 69 plant genomes. We identified and classified 88 450 LTR-RTs and determined 143 cases of high similarities between pairs of genomes. Most of them involved three Copia lineages (Oryco/Ivana, Retrofit/Ale, and Tork/Tar/Ikeros). A detailed analysis of three cases of high similarities involving Tork/Tar/Ikeros group shows an uneven distribution in the phylogeny of the elements and incongruence with between phylogenetic trees topologies, indicating they could be originated from HTs. Overall, our results suggest that LTR-RT Copia lineages share outstanding similarity between distant species and may likely be involved in HT mechanisms more frequent than initially estimated.

Characterization and expression profiles of the B-box gene family during plant growth and under low-nitrogen stress in Saccharum

BACKGROUND

B-box (BBX) zinc-finger transcription factors play crucial roles in plant growth, development, and abiotic stress responses. Nevertheless, little information is available on sugarcane (Saccharum spp.) BBX genes and their expression profiles.

RESULTS

In the present study, we characterized 25 SsBBX genes in the Saccharum spontaneum genome database. The phylogenetic relationships, gene structures, and expression patterns of these genes during plant growth and under low-nitrogen conditions were systematically analyzed. The SsBBXs were divided into five groups based on phylogenetic analysis. The evolutionary analysis further revealed that whole-genome duplications or segmental duplications were the main driving force for the expansion of the SsBBX gene family. The expression data suggested that many BBX genes (e.g., SsBBX1 and SsBBX13) may be helpful in both plant growth and low-nitrogen stress tolerance.

CONCLUSIONS

The results of this study offer new evolutionary insight into the BBX family members in how sugarcane grows and responds to stress, which will facilitate their utilization in cultivated sugarcane breeding.

Inference of Ancient Polyploidy Using Transcriptome Data

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

Revisiting ancient polyploidy in leptosporangiate ferns

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

Biotechnological interventions in reducing losses of tropical fruits and vegetables

Tropical fruits and vegetables are predominantly cultivated in warm climate zones, resulting in cultivar diversity in terms of structure, features, and physiology. These constitute a variety of bioactive ingredients such as vitamins, minerals, phenolic acids, anthocyanins, flavonoids, fatty acids, fiber, and their distinctive appearances attract customers across the world. The global production of fruit and vegetables has been attained a tremendous increase for the last few decades. However, huge losses at pre- and postharvest levels are major constraints in their judicious use. Traditional breeding strategies were used to minimize these losses, but their functionality is limited due to their time and labor intensiveness. Recent biotechnological, computational, and multiomics approaches not only address the losses concern but also aid in boosting crop productivity and nutritional values. This article emphasizes molecular tools that have been used to reduce losses of tropical fruits and vegetables at pre- and postharvest levels.

Large-scale analyses of angiosperm Flowering Locus T genes reveal duplication and functional divergence in monocots

FLOWERING LOCUS T (FT) are well-known key genes for initiating flowering in plants. Delineating the evolutionary history and functional diversity of FT genes is important for understanding the diversification of flowering time and how plants adapt to the changing surroundings. We performed a comprehensive phylogenetic analysis of FT genes in 47 sequenced flowering plants and the 1,000 Plant Transcriptomes (1KP) database with a focus on monocots, especially cereals. We revealed the evolutionary history of FT genes. The FT genes in monocots can be divided into three clades (I, II, and III), whereas only one monophyletic group was detected in early angiosperms, magnoliids, and eudicots. Multiple rounds of whole-genome duplications (WGD) events followed by gene retention contributed to the expansion and variation of FT genes in monocots. Amino acid sites in the clade II and III genes were preferentially under high positive selection, and some sites located in vital domain regions are known to change functions when mutated. Clade II and clade III genes exhibited high variability in important regions and functional divergence compared with clade I genes; thus, clade I is more conserved than clade II and III. Genes in clade I displayed higher expression levels in studied organs and tissues than the clade II and III genes. The co-expression modules showed that some of the FT genes might have experienced neofunctionalization and subfunctionalization, such as the acquisition of environmental resistance. Overall, FT genes in monocots might form three clades by the ancient gene duplication, and each clade was subsequently subjected to different selection pressures and amino acid substitutions, which eventually led to different expression patterns and functional diversification. Our study provides a global picture of FT genes' evolution in monocots, paving a road for investigating FT genes' function in future.

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

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