Elucidation of the 1-phenethylisoquinoline pathway from an endemic conifer Cephalotaxus hainanensis.

Cephalotaxines harbor great medical potential, but their natural source, the endemic conifer Cephalotaxus is highly endangered, creating a conflict between biotechnological valorization and preservation of biodiversity. Here, we construct the whole biosynthetic pathway to the 1-phenethylisoquinoline scaffold, as first committed compound for phenylethylisoquinoline alkaloids (PIAs), combining metabolic modeling, and transcriptome mining of Cephalotaxus hainanensis to infer the biosynthesis for PIA precursor. We identify a novel protein, ChPSS, driving the Pictet-Spengler condensation and show that this enzyme represents the branching point where PIA biosynthesis diverges from the concurrent benzylisoquinoline-alkaloids pathway. We also pinpoint ChDBR as crucial step to form 4-hydroxydihydrocinnamaldehyde diverging from lignin biosynthesis. The elucidation of the early PIA pathway represents an important step toward microbe-based production of these pharmaceutically important alkaloids resolving the conflict between biotechnology and preservation of biodiversity.

Full-Length Transcriptome Analysis of Four Different Tissues of Cephalotaxus oliveri.

Cephalotaxus oliveri is a tertiary relict conifer endemic to China, regarded as a national second-level protected plant in China. This species has experienced severe changes in temperature and precipitation in the past millions of years, adapting well to harsh environments. In view of global climate change and its endangered conditions, it is crucial to study how it responds to changes in temperature and precipitation for its conservation work. In this study, single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing were combined to generate the complete transcriptome of C. oliveri. Using the RNA-seq data to correct the SMRT sequencing data, the four tissues obtained 63,831 (root), 58,108 (stem), 33,013 (leaf) and 62,436 (male cone) full-length unigenes, with a N50 length of 2523, 3480, 3181, and 3267 bp, respectively. Additionally, 35,887, 11,306, 36,422, and 25,439 SSRs were detected for the male cone, leaf, root, and stem, respectively. The number of long non-coding RNAs predicted from the root was the largest (11,113), and the other tissues were 3408 (stem), 3193 (leaf), and 3107 (male cone), respectively. Functional annotation and enrichment analysis of tissue-specific expressed genes revealed the special roles in response to environmental stress and adaptability in the different four tissues. We also characterized the gene families and pathways related to abiotic factors. This work provides a comprehensive transcriptome resource for C. oliveri, and this resource will facilitate further studies on the functional genomics and adaptive evolution of C. oliveri.

Plastome phylogenomics of Cephalotaxus (Cephalotaxaceae) and allied genera.

BACKGROUND AND AIMS: Cephalotaxus is a paleo-endemic genus in East Asia that consists of about 7-9 conifer species. Despite its great economic and ecological importance, the relationships between Cephalotaxus and related genera, as well as the interspecific relationships within Cephalotaxus, have long been controversial, resulting in contrasting taxonomic proposals in delimitation of Cephalotaxaceae and Taxaceae. Based on plastome data, this study aims to reconstruct a robust phylogeny to infer the systematic placement and the evolutionary history of Cephalotaxus. METHODS: A total of 11 plastomes, representing all species currently recognized in Cephalotaxus and two Torreya species, were sequenced and assembled. Combining these with previously published plastomes, we reconstructed a phylogeny of Cephalotaxaceae and Taxaceae with nearly full taxonomic sampling. Under a phylogenetic framework and molecular dating, the diversification history of Cephalotaxus and allied genera was explored. KEY RESULTS: Phylogenetic analyses of 81 plastid protein-coding genes recovered robust relationships between Cephalotaxus and related genera, as well as providing a well-supported resolution of interspecific relationships within Cephalotaxus, Taxus, Torreya and Amentotaxus. Divergence time estimation indicated that most extant species of these genera are relatively young, although fossil and other molecular evidence consistently show that these genera are ancient plant lineages. CONCLUSIONS: Our results justify the taxonomic proposal that recognizes Cephalotaxaceae as a monotypic family, and contribute to a clear-cut delineation between Cephalotaxaceae and Taxaceae. Given that extant species of Cephalotaxus are derived from recent divergence events associated with the establishment of monsoonal climates in East Asia and Pleistocene climatic fluctuations, they are not evolutionary relics.

Characterisation, expression and functional analysis of PAL gene family in Cephalotaxus hainanensis.

Phenylalanine ammonia lyase (PAL) is the first committed step in the formation of phenylpropanoids, and catalyses the deamination of L-phenylalanine (L-Phe) to yield cinnamic acid. While PALs are common in plants, PAL genes involved in alkaloid biosynthesis in Cephalotaxus hainanensis have never been described. To obtain better knowledge of PAL genes and their number and function involved in Cephalotaxus alkaloid biosynthesis four PAL genes were screened and cloned. In vitro enzymatic analysis showed that all four PAL recombinant proteins could convert L-Phe to product trans-cinnamic acid, and showed strict substrate specificity. Moreover, the expression profiles of four ChPALs were analysed using qRT-PCR, which showed that they had higher transcript levels in roots and stems, and that different ChPALs displayed different response sensitivities and change patterns in response to stimuli. Several metabolic compounds were measured in stimulated leaves using UPLC-MS, and indicating the concentration of Cephalotaxus alkaloids and cinnamic acid in leaves subjected to different conditions. These concentrations increased significantly after treatment with 100 mM NaCl, 100 mM mannitol, 100 µM SA and 10 µM ABA. The expression levels of four PAL genes showed indications of upregulation after treatment. These results supply an important foundation for further research on candidate genes involved in the biosynthesis of Cephalotaxus alkaloids.

Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China.

Cephalotaxus oliveri is an endangered tertiary relict conifer endemic to China. The species survives in a wide range from west to east with heterogeneous climatic conditions. Precipitation and temperature are main restrictive factors for distribution of C. oliveri. In order to comprehend the mechanism of adaptive evolution to climate variation, we employed ISSR markers to detect adaptive evolution loci, to identify the association between variation in temperature and precipitation and adaptive loci, and to investigate the genetic structure for 22 C. oliveri natural populations. In total, 14 outlier loci were identified, of which five were associated with temperature and precipitation. Among outlier loci, linkage disequilibrium (LD) was high (42.86%), which also provided strong evidence for selection. In addition, C. oliveri possessed high genetic variation (93.31%) and population differentiation, which may provide raw material to evolution and accelerate local adaptation, respectively. Ecological niche modeling showed that global warming will cause a shift for populations of C. oliveri from south to north with a shrinkage of southern areas. Our results contribute to understand the potential response of conifers to climatic changes, and provide new insights for conifer resource management and conservation strategies.

De Novo characterization of a Cephalotaxus hainanensis transcriptome and genes related to paclitaxel biosynthesis.

Cephalotaxus hainanensis, an endangered plant, is known to contain several metabolites with anti-cancer activity. Despite its clinical impact, the alkaloid metabolism of this species has remained largely uncharacterized. The potential of Cephalotaxus for metabolic engineering of medically interesting compounds has, so far, not been exploited, due to the almost complete lack of molecular information. We have therefore performed a high throughput RNA-seq analysis and assembled the transcriptome de novo. Raw reads comprising 4.3 Gbp were assembled de novo into 39,416 unique sequences (unigenes) with a mean length of 1,089.8 bp and a total assembly size of 45.8 Mbp, which equals to more than 50 times the number of Cephalotaxaceae sequences currently deposited in the GenBank (as of August 2013). As proof of principle for medically interesting pathways, gene fragments related to paclitaxel biosynthesis were searched and detected. To verify their functionality, the metabolic product paclitaxel, and its precursor baccatin III, were identified in the leaves of C. hainanensis by HPLC, and shown to be induced by MeJA. This finding demonstrates exemplarily the potential of the annotated transcriptome as information resource for the biotechnological exploitation of plant secondary metabolism.

Phylogeography of Cephalotaxus oliveri (Cephalotaxaceae) in relation to habitat heterogeneity, physical barriers and the uplift of the Yungui Plateau.

Habitat heterogeneity, physical barriers, and the uplift of the Yungui Plateau were found to deeply affect the phylogeographic pattern and evolutionary history of Cephalotaxus oliveri, a perennial conifer endemic to China. In this study, we explored the phylogeography using three chloroplast sequences (trnL-trnF, trnT-trnD and atpB-rbcL) in 22 natural populations of C. oliveri distributed throughout its range. The Yungui Plateau populations of C. oliveri were revealed to origin ca. 9.15Ma by molecular clock estimation, which is consistent with rapid uplift of the Qinghai-Tibetan Plateau (QTP) ca. 8-10Ma. Additionally, geological effects of the Yungui Plateau were suggested to promote the rapid intra-specific differentiation of C. oliveri in the Pliocene and Early Pleistocene. The relatively low level of genetic diversity (h=0.719, θ=1.17×10(-3)) and high population differentiation (NST=0.771 and GST=0.642) implied restricted gene flow among populations, which was confirmed by the Nested Clade Analysis (NCA). Mismatch distribution and haplotypes network provided evidences of recent demographic population expansion. Furthermore, the statistical dispersal-vicariance analysis indicated that the center of origin was in Central China. The comparison of haplotype distribution patterns indicated that the regions of HNHPS and HBLD were the potential refugia during the Pleistocene ice ages. Our results highlighted that habitat heterogeneity and physical barriers presenting in a species range can predict genetic patterns.

The complete chloroplast genome sequence of Cephalotaxus oliveri (Cephalotaxaceae): evolutionary comparison of cephalotaxus chloroplast DNAs and insights into the loss of inverted repeat copies in gymnosperms.

We have determined the complete chloroplast (cp) genome sequence of Cephalotaxus oliveri. The genome is 134,337 bp in length, encodes 113 genes, and lacks inverted repeat (IR) regions. Genome-wide mutational dynamics have been investigated through comparative analysis of the cp genomes of C. oliveri and C. wilsoniana. Gene order transformation analyses indicate that when distinct isomers are considered as alternative structures for the ancestral cp genome of cupressophyte and Pinaceae lineages, it is not possible to distinguish between hypotheses favoring retention of the same IR region in cupressophyte and Pinaceae cp genomes from a hypothesis proposing independent loss of IRA and IRB. Furthermore, in cupressophyte cp genomes, the highly reduced IRs are replaced by short repeats that have the potential to mediate homologous recombination, analogous to the situation in Pinaceae. The importance of repeats in the mutational dynamics of cupressophyte cp genomes is also illustrated by the accD reading frame, which has undergone extreme length expansion in cupressophytes. This has been caused by a large insertion comprising multiple repeat sequences. Overall, we find that the distribution of repeats, indels, and substitutions is significantly correlated in Cephalotaxus cp genomes, consistent with a hypothesis that repeats play a role in inducing substitutions and indels in conifer cp genomes.

Characterization of 15 polymorphic microsatellite loci for Cephalotaxus oliveri (Cephalotaxaceae), a conifer of medicinal importance.

Cephalotaxus oliveri is a scarce medicinal conifer endemic to the south central region of China and Vietnam. A small fragmented population presently exists due to anthropogenic disturbance. C. oliveri has been used for its alkaloids harringtonine and homoharringtonine, which are effective against leucocythemia and lymphadenosarcoma. Monoecious plants have been detected in nature, although they were understood to be dioecious. In order to study the mating system, population genetics and the genetic effects of habitat fragmentation on C. oliveri, 15 polymorphic and 12 monomorphic microsatellite loci were developed for C. oliveri by using the Fast Isolation by AFLP of Sequences Containing repeats (FIASCO) protocol. The polymorphisms were assessed in 96 individuals from three natural populations (32 individuals per population). The number of alleles per locus ranged from two to 33, the observed and expected heterozygosity per locus ranged from 0.000 to 1.000 and from 0.000 to 0.923, respectively. These loci would facilitate a comprehensive understanding of the genetic dynamics on C. oliveri, which will be useful for establishing effective conservation strategies for this species.

Development of microsatellite loci for Cephalotaxus oliveri (Cephalotaxaceae) and cross-amplification in Cephalotaxus.

PREMISE OF THE STUDY: Microsatellite loci were developed for Cephalotaxus oliveri, an endemic and endangered conifer in China, which will allow assessment of the levels of genetic diversity and a means to understand the genetic consequences of habitat fragmentation. METHODS AND RESULTS: Using the Fast Isolation by AFLP of Sequences COntaining (FIASCO) Repeats protocol, 19 microsatellite loci were identified in C. oliveri, 13 of which were polymorphic within a sample of 52 individuals representing five natural populations. The actual number of alleles per locus ranged from one to five. Twelve polymorphic loci were also successfully amplified in C. fortunei. CONCLUSIONS: These microsatellite loci will provide a useful tool for further investigating genetic variation in natural populations of C. oliveri, which will inform future conservation and management strategies. Additionally, cross-amplification in C. fortunei suggested the potential utility of these loci in this and other congeneric species.

Molecular cloning and characterization of a 2C-methyl-D: -erythritol 2,4-cyclodiphosphate synthase gene from Cephalotaxus harringtonia.

The full-length MECPS cDNA sequence (designated as Chmecps, GenBank Accession No.: DQ415658) was isolated by rapid amplification of cDNA ends (RACE) for the first time from Cephalotaxus harringtonia. The full-length cDNA of Chmecps was 1,146 bp containing a 753 bp open reading frame (ORF) encoding a polypeptide of 250 amino acids with a calculated mass of 26.67 kDa and an isoelectric point of 9.35. Comparative and bioinformatics analyses revealed that ChMECPS showed extensive homology with MECPSs from other plant species. Phylogenetic analysis indicated ChMECPS was more ancient than other plant MECPSs. Southern hybridization analysis of the genomic DNA showed that Chmecps was a single copy gene. Tissue expression pattern analysis revealed that ChMECPS expressed strongly in root and leaf, weakly in stem.