Distinctive plastome evolution in carnivorous angiosperms.

BACKGROUND: Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive. RESULTS: We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection. CONCLUSION: Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

Haplotype-resolved genomes of wild octoploid progenitors illuminate genomic diversifications from wild relatives to cultivated strawberry.

Strawberry is an emerging model for studying polyploid genome evolution and rapid domestication of fruit crops. Here we report haplotype-resolved genomes of two wild octoploids (Fragaria chiloensis and Fragaria virginiana), the progenitor species of cultivated strawberry. Substantial variation is identified between species and between haplotypes. We redefine the four subgenomes and track the genetic contributions of diploid species by additional sequencing of the diploid F. nipponica genome. We provide multiple lines of evidence that F. vesca and F. iinumae, rather than other described extant species, are the closest living relatives of these wild and cultivated octoploids. In response to coexistence with quadruplicate gene copies, the octoploid strawberries have experienced subgenome dominance, homoeologous exchanges and coordinated expression of homoeologous genes. However, some homoeologues have substantially altered expression bias after speciation and during domestication. These findings enhance our understanding of the origin, genome evolution and domestication of strawberries.

The Evolution of Mitochondrial Genomes between Two Cymbidium Sister Species: Dozens of Circular Chromosomes and the Maintenance and Deterioration of Genome Synteny.

Plant mitochondrial genomes (mitogenomes) exhibit fluid genome architectures, which could lead to the rapid erosion of genome synteny over a short evolutionary time scale. Among the species-rich orchid family, the leafy Cymbidium lancifolium and leafless Cymbidium macrorhizon are sister species with remarkable differences in morphology and nutritional physiology. Although our understanding of the evolution of mitochondria is incomplete, these sister taxa are ideal for examining this subject. In this study, the complete mitogenomes of C. lancifolium and C. macrorhizon, totaling 704,244 bp and 650,751 bp, respectively, were assembled. In the 2 mitogenomes, 38 protein-coding genes, 18 cis- and 6 trans-spliced introns, and approximately 611 Kb of homologous sequences are identical; overall, they have 99.4% genome-wide similarity. Slight variations in the mitogenomes of C. lancifolium and C. macrorhizon in repeat content (21.0 Kb and 21.6 Kb, respectively) and mitochondrial DNA of plastid origin (MIPT; 38.2 Kb and 37.5 Kb, respectively) were observed. The mitogenome architectures of C. lancifolium and C. macrorhizon are complex and comprise 23 and 22 mini-circular chromosomes, respectively. Pairwise comparisons indicate that the two mitogenomes are largely syntenic, and the disparity in chromosome numbers is likely due to repeat-mediated rearrangements among different chromosomes. Notably, approximately 93.2 Kb C. lancifolium mitochondrial sequences lack any homology in the C. macrorhizon mitogenome, indicating frequent DNA gains and losses, which accounts mainly for the size variation. Our findings provide unique insights into mitogenome evolution in leafy and leafless plants of sister species and shed light on mitogenome dynamics during the transition from mixotrophy to mycoheterotrophy.

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility.

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

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.

Pan-mitogenomics reveals the genetic basis of cytonuclear conflicts in citrus hybridization, domestication, and diversification.

Although interactions between the cytoplasmic and nuclear genomes occurred during diversification of many plants, the evolutionary conflicts due to cytonuclear interactions are poorly understood in crop breeding. Here, we constructed a pan-mitogenome and identified chimeric open reading frames (ORFs) generated by extensive structural variations (SVs). Meanwhile, short reads from 184 accessions of citrus species were combined to construct three variation maps for the nuclear, mitochondrial, and chloroplast genomes. The population genomic data showed discordant topologies between the cytoplasmic and nuclear genomes because of differences in mutation rates and levels of heteroplasmy from paternal leakage. An analysis of species-specific SVs indicated that mitochondrial heteroplasmy was common and that chloroplast heteroplasmy was undetectable. Interestingly, we found a prominent divergence in the mitogenomes and the highest genetic load in the, which may provide the basis for cytoplasmic male sterility (CMS) and thus influence the reshuffling of the cytoplasmic and nuclear genomes during hybridization. Using cytoplasmic replacement experiments, we identified a type of species-specific CMS in mandarin related to two chimeric mitochondrial genes. Our analyses indicate that cytoplasmic genomes from mandarin have rarely been maintained in hybrids and that paternal leakage produced very low levels of mitochondrial heteroplasmy in mandarin. A genome-wide association study (GWAS) provided evidence for three nuclear genes that encode pentatricopeptide repeat (PPR) proteins contributing to the cytonuclear interactions in the Citrus genus. Our study demonstrates the occurrence of evolutionary conflicts between cytoplasmic and nuclear genomes in citrus and has important implications for genetics and breeding.

Phylogenetic Analysis of Wild Species and the Maternal Origin of Cultivars in the Genus Lilium Using 114 Plastid Genomes.

Lilies are one of the most important ornamental flowers worldwide with approximately 100 wild species and numerous cultivars, but the phylogenetic relationships among wild species and their contributions to these cultivars are poorly resolved. We collected the major Lilium species and cultivars and assembled their plastome sequences. Our phylogenetic reconstruction using 114 plastid genomes, including 70 wild species representing all sections and 42 cultivars representing six hybrid divisions and two outgroups, uncovered well-supported genetic relationships within Lilium. The wild species were separated into two distinct groups (groups A and B) associated with geographical distribution, which further diversified into eight different clades that were phylogenetically well supported. Additional support was provided by the distributions of indels and single-nucleotide variants, which were consistent with the topology. The species of sections Archelirion, Sinomartagon III, and Leucolirion 6a and 6b were the maternal donors for Oriental hybrids, Asiatic hybrids, Trumpet hybrids, and Longiflorum hybrids, respectively. The maternal donors of the OT hybrids originated from the two sections Archelirion and Leucolirion 6a, and LA hybrids were derived from the two sections Leucolirion 6b and Sinomartagon. Our study provides an important basis for clarifying the infrageneric classification and the maternal origin of cultivars in Lilium.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions.

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

Plastid RNA editing reduction accompanied with genetic variations in Cymbidium, a genus with diverse lifestyle modes.

Recent sequencing efforts have broadly uncovered the evolutionary trajectory of plastid genomes (plastomes) of flowering plants in diverse habitats, yet our knowledge of the evolution of plastid posttranscriptional modifications is limited. In this study, we generated 11 complete plastomes and performed ultra-deep transcriptome sequencing to investigate the co-evolution of plastid RNA editing and genetic variation in Cymbidium, a genus with diverse trophic lifestyles. Genome size and gene content is reduced in terrestrial and green mycoheterotrophic orchids relative to their epiphytic relatives. This could be partly due to extensive losses and pseudogenization of ndh genes for the plastid NADH dehydrogenase-like complex, but independent pseudogenization of ndh genes has also occurred in the epiphyte C. mannii, which was reported to use strong crassulacean acid metabolism photosynthesis. RNA editing sites are abundant but variable in number among Cymbidium plastomes. The nearly twofold variation in editing abundance is mainly due to extensive reduction of ancestral editing sites in ndh transcripts of terrestrial, mycoheterotrophic, and C. mannii plastomes. The co-occurrence of editing reduction and pseudogenization in ndh genes suggests functional constraints on editing machinery may be relaxed, leading to nonrandom loss of ancestral edited sites via reduced editing efficiency. This study represents the first systematic examination of RNA editing evolution linked to plastid genome variation in a single genus. We also propose an explanation for how genomic and posttranscriptional variations might be affected by lifestyle-associated ecological adaptation strategies in Cymbidium.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing.

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Differential expressions of anthocyanin synthesis genes underlie flower color divergence in a sympatric Rhododendron sanguineum complex.

BACKGROUND: The Rhododendron sanguineum complex is endemic to alpine mountains of northwest Yunnan and southeast Tibet of China. Varieties in this complex exhibit distinct flower colors even at the bud stage. However, the underlying molecular regulations for the flower color variation have not been well characterized. Here, we investigated this via measuring flower reflectance profiles and comparative transcriptome analyses on three coexisting varieties of the R. sanguineum complex, with yellow flush pink, bright crimson, and deep blackish crimson flowers respectively. We compared the expression levels of differentially-expressed-genes (DEGs) of the anthocyanin / flavonoid biosynthesis pathway using RNA-seq and qRT-PCR data. We performed clustering analysis based on transcriptome-derived Single Nucleotide Polymorphisms (SNPs) data, and finally analyzed the promoter architecture of DEGs. RESULTS: Reflectance spectra of the three color morphs varied distinctively in the range between 400 and 700 nm, with distinct differences in saturation, brightness, hue, and saturation/hue ratio, an indirect measurement of anthocyanin content. We identified 15,164 orthogroups that were shared among the three varieties. The SNP clustering analysis indicated that the varieties were not monophyletic. A total of 40 paralogous genes encoding 12 enzymes contributed to the flower color polymorphism. These anthocyanin biosynthesis-related genes were associated with synthesis, modification and transportation properties (RsCHS, RsCHI, RsF3H, RsF3'H, RsFLS, RsANS, RsAT, RsOMT, RsGST), as well as genes involved in catabolism and degradation (RsBGLU, RsPER, RsCAD). Variations in sequence and cis-acting elements of these genes might correlate with the anthocyanin accumulation, thus may contribute to the divergence of flower color in the R. sanguineum complex. CONCLUSIONS: Our results suggested that the varieties are very closely related and flower color variations in the R. sanguineum complex correlate tightly with the differential expression levels of genes involved in the anabolic and catabolic synthesis network of anthocyanin. Our study provides a scenario involving intricate relationships between genetic mechanisms for floral coloration accompanied by gene flow among the varieties that may represent an early case of pollinator-mediated incipient sympatric speciation.

Temporal Regulation of the Metabolome and Proteome in Photosynthetic and Photorespiratory Pathways Contributes to Maize Heterosis.

Heterosis or hybrid vigor is widespread in plants and animals. Although the molecular basis for heterosis has been extensively studied, metabolic and proteomic contributions to heterosis remain elusive. Here we report an integrative analysis of time-series metabolome and proteome data in maize (Zea mays) hybrids and their inbred parents. Many maize metabolites and proteins are diurnally regulated, and many of these show nonadditive abundance in the hybrids, including key enzymes and metabolites involved in carbon assimilation. Compared with robust trait heterosis, metabolic heterosis is relatively mild. Interestingly, most amino acids display negative mid-parent heterosis (MPH), i.e., having lower values than the average of the parents, while sugars, alcohols, and nucleoside metabolites show positive MPH. From the network perspective, metabolites in the photosynthetic pathway show positive MPH, whereas metabolites in the photorespiratory pathway show negative MPH, which corresponds to nonadditive protein abundance and enzyme activities of key enzymes in the respective pathways in the hybrids. Moreover, diurnally expressed proteins that are upregulated in the hybrids are enriched in photosynthesis-related gene-ontology terms. Hybrids may more effectively remove toxic metabolites generated during photorespiration, and thus maintain higher photosynthetic efficiency. These metabolic and proteomic resources provide unique insight into heterosis and its utilization for high yielding maize and other crop plants.

Episodic and guanine-cytosine-biased bursts of intragenomic and interspecific synonymous divergence in Ajugoideae (Lamiaceae) mitogenomes.

Synonymous substitution rates in plant mitochondrial genomes vary by orders of magnitude among species, whereas synonymous rates among genes within a genome are generally consistent. Exceptionally, genes within the Ajuga reptans (Lamiaceae) mitochondrial genome exhibit unprecedented intragenomic heterogeneity in synonymous sequence divergence, but the biological mechanisms underlying this rate variation remain unclear. We tracked the origin and evolutionary trajectory of mitochondrial rate variations by dense sampling in Ajugoideae and found differences in the timing and magnitude of rate acceleration for particular genes. The most divergent genes accelerated earlier, retained a high rate across Ajugoideae, and are generally devoid of RNA editing, whereas moderately diverged genes accelerated later and retained relatively higher RNA editing frequency. The acceleration of mutation rates correlates with increased guanine-cytosine (GC) content, suggesting a key role for GC-biased gene conversion and/or repair after the breakage of ancestral gene clusters.

Evolution of self-compatibility by a mutant Sm-RNase in citrus.

Self-incompatibility (SI) is an important mechanism that prevents self-fertilization and inbreeding in flowering plants. The most widespread SI system utilizes S ribonucleases (S-RNases) and S-locus F-boxes (SLFs) as S determinants. In citrus, SI is ancestral, and Citrus maxima (pummelo) is self-incompatible, while Citrus reticulata (mandarin) and its hybrids are self-compatible (SC). Here, we identify nine highly polymorphic pistil-specific, developmentally expressed S-RNases from pummelo that segregate with S haplotypes in a gametophytic manner and cluster with authentic S-RNases. We provide evidence that these S-RNases function as the female S determinants in citrus. Moreover, we show that each S-RNase is linked to approximately nine SLFs. In an analysis of 117 citrus SLF and SFL-like (SLFL) genes, we reveal that they cluster into 12 types and that the S-RNases and intra-haplotypic SLF and SLFL genes co-evolved. Our data support the notion that citrus have a S locus comprising a S-RNase and several SLFs that fit the non-self-recognition model. We identify a predominant single nucleotide mutation, Sm-RNase, in SC citrus, which provides a 'natural' loss of function. We show that SI-SC transitions due to the Sm-RNase initially arose in mandarin, spreading to its hybrids and became fixed. Identification of an evolutionarily distant new genus utilizing the S-RNase-based SI system, >100 million years separated from the nearest S-RNase family, is a milestone for evolutionary comparative studies.

Extensive Shifts from Cis- to Trans-splicing of Gymnosperm Mitochondrial Introns.

Hundreds of plant mitogenomes have been sequenced from angiosperms, but relatively few mitogenomes are available from its sister lineage, gymnosperms. To examine mitogenomic diversity among extant gymnosperms, we generated draft mitogenomes from 11 diverse species and compared them with four previously published mitogenomes. Examined mitogenomes from Pinaceae and cycads retained all 41 protein genes and 26 introns present in the common ancestor of seed plants, whereas gnetophyte and cupressophyte mitogenomes experienced extensive gene and intron loss. In Pinaceae and cupressophyte mitogenomes, an unprecedented number of exons are distantly dispersed, requiring trans-splicing of 50-70% of mitochondrial introns to generate mature transcripts. RNAseq data confirm trans-splicing of these dispersed exons in Pinus. The prevalence of trans-splicing in vascular plant lineages with recombinogenic mitogenomes suggests that genomic rearrangement is the primary cause of shifts from cis- to trans-splicing in plant mitochondria.

Complete loss of RNA editing from the plastid genome and most highly expressed mitochondrial genes of Welwitschia mirabilis.

Comparative genomics among gymnosperms suggested extensive loss of mitochondrial RNA editing sites from Welwitschia mirabilis based on predictive analysis. However, empirical or transcriptome data to confirm this massive loss event are lacking, and the potential mechanisms of RNA site loss are unclear. By comparing genomic sequences with transcriptomic and reverse-transcription PCR sequencing data, we performed a comprehensive analysis of the pattern of RNA editing in the mitochondrial and plastid genomes (mitogenome and plastome, respectively) of W. mirabilis and a second gymnosperm, Ginkgo biloba. For W. mirabilis, we found only 99 editing sites located in 13 protein-coding genes in the mitogenome and a complete loss of RNA editing from the plastome. The few genes having high editing frequency in the Welwitschia mitogenome showed a strong negative correlation with gene expression level. Comparative analyses with G. biloba, containing 1,405 mitochondrial and 345 plastid editing sites, revealed that the editing loss from W. mirabilis is mainly due to the substitution of editable cytidines to thymidines at the genomic level, which could be caused by retroprocessing. Our result is the first study to uncover massive editing loss from both the mitogenome and plastome in a single genus. Furthermore, our results suggest that gene expression level and retroprocessing both contributed to the evolution of RNA editing in plant organellar genomes.

Physiological and Transcriptome Analyses Reveal Short-Term Responses and Formation of Memory Under Drought Stress in Rice.

In some plants, exposure to stress can induce a memory response, which appears to play an important role in adaptation to recurrent stress environments. However, whether rice exhibits drought stress memory and the molecular mechanisms that might underlie this process have remained unclear. Here, we ensured that rice drought memory was established after cycles of mild drought and re-watering treatment, and studied gene expression by whole-transcriptome strand-specific RNA sequencing (ssRNA-seq). We detected 6,885 transcripts and 238 lncRNAs involved in the drought memory response, grouped into 16 distinct patterns. Notably, the identified genes of dosage memory generally did not respond to the initial drought treatment. Our results demonstrate that stress memory can be developed in rice under appropriate water deficient stress, and lncRNA, DNA methylation and endogenous phytohormones (especially abscisic acid) participate in rice short-term drought memory, possibly acting as memory factors to activate drought-related memory transcripts in pathways such as photosynthesis and proline biosynthesis, to respond to the subsequent stresses.

Alternative splicing coupled to nonsense-mediated mRNA decay contributes to the high-altitude adaptation of maca (Lepidium meyenii).

Alpine plants remain the least studied plant communities in terrestrial ecosystems. However, how they adapt to high-altitude environments is far from clear. Here, we used RNA-seq to investigate a typical alpine plant maca (Lepidium meyenii) to understand its high-altitude adaptation at transcriptional and post-transcriptional level. At transcriptional level, we found that maca root significantly up-regulated plant immunity genes in day-time comparing to night-time, and up-regulated abiotic (cold/osmotic) stress response genes in Nov and Dec comparing to Oct. In addition, 17 positively selected genes were identified, which could be involved in mitochondrion. At post-transcriptional level, we found that maca had species-specific characterized alternative splicing (AS) profile which could be influenced by stress environments. For example, the alternative 3' splice site events (A3SS, 39.62%) were predominate AS events in maca, rather than intron retention (IR, 23.17%). Interestingly, besides serine/arginine-rich (SR) proteins and long non-coding RNAs (lncRNAs), a lot of components in nonsense-mediated mRNA decay (NMD) were identified under differential alternative splicing (DAS), supporting AS coupled to NMD as essential mechanisms for maca's stress responses and high-altitude adaptation. Taken together, we first attempted to unveil maca's high-altitude adaptation mechanisms based on transcriptome and post-transcriptome evidence. Our data provided valuable insights to understand the high-altitude adaptation of alpine plants.

The terpene limonene induced the green mold of citrus fruit through regulation of reactive oxygen species (ROS) homeostasis in Penicillium digitatum spores.

Herein, the dosage effect of limonene on the P. digitatum spore germination and its regulatory mechanisms were investigated. Results showed that limonene only at low concentrations displayed a stimulatory role, with the optimal concentration being 0.25 µL/mL. GC-MS and GC analysis revealed that limonene contents remained relative stable and no evidence of transformation was observed at stimulatory concentrations. Metabolomics analysis showed that 61 metabolites including organic acids, amino acids, sugars, nucleosides, fatty acids, and their derivatives, were significantly changed (P < 0.05), suggesting the transitions between soluble sugars and energy-related metabolisms. Proteomic analysis demonstrated that proteins in energy-related pathways and ROS homeostasis were also influenced. These were further confirmed by the activities of catalase (CAT), superoxide dismutase (SOD), and glutathione-S-transferase (GST), the contents of reactive oxygen species (ROS), hydrogen peroxide (H2O2), and glutathione (GSH). Our present research indicates that ROS homeostasis is involved in the limonene induced spore germination of P. digitatum.

Phylogenomic evidence for ancient recombination between plastid genomes of the Cupressus-Juniperus-Xanthocyparis complex (Cupressaceae).

BACKGROUND: Phylogenetic relationships among Eastern Hemisphere cypresses, Western Hemisphere cypresses, junipers, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 11 new plastid genomes (plastomes) from the five variously described genera in this complex (Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis) and compared them with additional plastomes from diverse members of Cupressaceae. RESULTS: Phylogenetic analysis of protein-coding genes recovered a topology in which Juniperus is sister to Cupressus, whereas a tree based on whole plastomes indicated that the Callitropsis-Hesperocyparis-Xanthocyparis (CaHX) clade is sister to Cupressus. A sliding window analysis of site-specific phylogenetic support identified a ~ 15 kb region, spanning the genes ycf1 and ycf2, which harbored an anomalous signal relative to the rest of the genome. After excluding these genes, trees based on the remainder of the genes and genome consistently recovered a topology grouping the CaHX clade and Cupressus with strong bootstrap support. In contrast, trees based on the ycf1 and ycf2 region strongly supported a sister relationship between Cupressus and Juniperus. CONCLUSIONS: These results demonstrate that standard phylogenomic analyses can result in strongly supported but conflicting trees. We suggest that the conflicting plastomic signals result from an ancient introgression event involving ycf1 and ycf2 that occurred in an ancestor of this species complex. The introgression event was facilitated by plastomic recombination in an ancestral heteroplasmic individual carrying distinct plastid haplotypes, offering further evidence that recombination occurs between plastomes. Finally, we provide strong support for previous proposals to recognize five genera in this species complex: Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis.