Complete chloroplast genome sequence of holoparasite Cistanche deserticola (Orobanchaceae) reveals gene loss and horizontal gene transfer from its host Haloxylon ammodendron (Chenopodiaceae)

The complete chloroplast genome of Pyankovia brachiata (amaranthaceae), an annual desert plant in China

Pyankovia brachiata (Pall.) Akhani & Roalson 2007, is an annual plant belonging to the genus Pyankovia, family Amaranthaceae, which is widely distributed in the inland deserts of Northwest China. P. brachiata was previously categorized under the genus Salsola in Salsoleae and has been a long-standing topic of debate. Therefore, the complete chloroplast genome of P. brachiata must be studied to provide a theoretical reference for species classification. In this study, we sequenced P. brachiata samples and determined the species' complete chloroplast genome. The complete chloroplast genome was 149,922 bp in length, with one large single copy (LSC: 83,565 bp), one small single copy (SSC: 18,535 bp), and two inverted repeat regions (IRa and IRb, 23,911 bp each). It contains 132 genes, including 87 protein-coding, eight rRNA, and 37 tRNA genes. The phylogenetic position showed that P. brachiata has the closest relationship with Caroxylon passerinum (accession number: NC057191.1). This study will provide genetic information and be beneficial to understanding the systematic position of P. brachiata within the Amaranthaceae.

Gene Losses and Homology of the Chloroplast Genomes of Taxillus and Phacellaria Species

Research on the chloroplast genome of parasitic plants is limited. In particular, the homology between the chloroplast genomes of parasitic and hyperparasitic plants has not been reported yet. In this study, three chloroplast genomes of Taxillus (Taxillus chinensis, Taxillus delavayi, and Taxillus thibetensis) and one chloroplast genome of Phacellaria (Phacellaria rigidula) were sequenced and analyzed, among which T. chinensis is the host of P. rigidula. The chloroplast genomes of the four species were 119,941-138,492 bp in length. Compared with the chloroplast genome of the autotrophic plant Nicotiana tabacum, all of the ndh genes, three ribosomal protein genes, three tRNA genes and the infA gene were lost in the three Taxillus species. Meanwhile, in P. rigidula, the trnV-UAC gene and the ycf15 gene were lost, and only one ndh gene (ndhB) existed. The results of homology analysis showed that the homology between P. rigidula and its host T. chinensis was low, indicating that P. rigidula grows on its host T. chinensis but they do not share the chloroplast genome. In addition, horizontal gene transfer was not found between P. rigidula and its host T. chinensis. Several candidate highly variable regions in the chloroplast genomes of Taxillus and Phacellaria species were selected for species identification study. Phylogenetic analysis revealed that the species of Taxillus and Scurrula were closely related and supported that Scurrula and Taxillus should be treated as congeneric, while species in Phacellaria had a close relationship with that in Viscum.

Plastome variations reveal the distinct evolutionary scenarios of plastomes in the subfamily Cereoideae (Cactaceae)

BACKGROUND

The cactus family (Cactaceae) has been reported to have evolved a minimal photosynthetic plastome size, with the loss of inverted-repeat (IR) regions and NDH gene suites. However, there are very limited genomic data on the family, especially Cereoideae, the largest subfamily of cacti.

RESULTS

In the present study, we assembled and annotated 35 plastomes, 33 of which were representatives of Cereoideae, alongside 2 previously published plastomes. We analyzed the organelle genomes of 35 genera in the subfamily. These plastomes have variations rarely observed in those of other angiosperms, including size differences (with ~ 30 kb between the shortest and longest), dramatic dynamic changes in IR boundaries, frequent plastome inversions, and rearrangements. These results suggested that cacti have the most complex plastome evolution among angiosperms.

CONCLUSION

These results provide unique insight into the dynamic evolutionary history of Cereoideae plastomes and refine current knowledge of the relationships within the subfamily.

Comparative Analyses of Chloroplast Genomes for Parasitic Species of Santalales in the Light of Two Newly Sequenced Species, Taxillus nigrans and Scurrula parasitica

When a flowering plant species changes its life history from self-supply to parasite, its chloroplast genomes may have experienced functional physical reduction, and gene loss. Most species of Santalales are hemiparasitic and few studies focus on comparing the chloroplast genomes of the species from this order. In this study, we collected and compared chloroplast genomes of 12 species of Santalales and sequenced the chloroplast genomes of Taxillus nigrans and Scurrula parasitica for the first time. The chloroplast genomes for these species showed typical quadripartite structural organization. Phylogenetic analysis suggested that these 12 species of Santalales clustered into three clades: Viscum (4 spp.) and Osyris (1 sp.) in the Santalaceae and Champereia (1 sp.) in the Opiliaceae formed one clade, while Taxillus (3 spp.) and Scurrula (1 sp.) in the Loranthaceae and Schoepfia (1 sp.) in the Schoepfiaceae formed another clade. Erythropalum (1 sp.), in the Erythropalaceae, appeared as a third, most distant, clade within the Santalales. In addition, both Viscum and Taxillus are monophyletic, and Scurrula is sister to Taxillus. A comparative analysis of the chloroplast genome showed differences in genome size and the loss of genes, such as the ndh genes, infA genes, partial ribosomal genes, and tRNA genes. The 12 species were classified into six categories by the loss, order, and structure of genes in the chloroplast genome. Each of the five genera (Viscum, Osyris, Champereia, Schoepfia, and Erythropalum) represented an independent category, while the three Taxillus species and Scurrula were classified into a sixth category. Although we found that different genes were lost in various categories, most genes related to photosynthesis were retained in the 12 species. Hence, the genetic information accorded with observations that they are hemiparasitic species. Our comparative genomic analyses can provide a new case for the chloroplast genome evolution of parasitic species.

Analysis of Complete Chloroplast Genome: Structure, Phylogenetic Relationships of Galega orientalis and Evolutionary Inference of Galegeae

Galega orientalis, a leguminous herb in the Fabaceae family, is an ecologically and economically important species widely cultivated for its strong stress resistance and high protein content. However, genomic information of Galega orientalis has not been reported, which limiting its evolutionary analysis. The small genome size makes chloroplast relatively easy to obtain genomic sequence for phylogenetic studies and molecular marker development. Here, the chloroplast genome of Galega orientalis was sequenced and annotated. The results showed that the chloroplast genome of G. orientalis is 125,280 bp in length with GC content of 34.11%. A total of 107 genes were identified, including 74 protein-coding genes, 29 tRNAs and four rRNAs. One inverted repeat (IR) region was lost in the chloroplast genome of G. orientalis. In addition, five genes (rpl22, ycf2, rps16, trnE-UUC and pbf1) were lost compared with the chloroplast genome of its related species G. officinalis. A total of 84 long repeats and 68 simple sequence repeats were detected, which could be used as potential markers in the genetic studies of G. orientalis and related species. We found that the Ka/Ks values of three genes petL, rpl20, and ycf4 were higher than one in the pairwise comparation of G. officinalis and other three Galegeae species (Calophaca sinica, Caragana jubata, Caragana korshinskii), which indicated those three genes were under positive selection. A comparative genomic analysis of 15 Galegeae species showed that most conserved non-coding sequence regions and two genic regions (ycf1 and clpP) were highly divergent, which could be used as DNA barcodes for rapid and accurate species identification. Phylogenetic trees constructed based on the ycf1 and clpP genes confirmed the evolutionary relationships among Galegeae species. In addition, among the 15 Galegeae species analyzed, Galega orientalis had a unique 30-bp intron in the ycf1 gene and Tibetia liangshanensis lacked two introns in the clpP gene, which is contrary to existing conclusion that only Glycyrrhiza species in the IR lacking clade (IRLC) lack two introns. In conclusion, for the first time, the complete chloroplast genome of G. orientalis was determined and annotated, which could provide insights into the unsolved evolutionary relationships within the genus Galegeae.

Variation in Chloroplast Genome Size: Biological Phenomena and Technological Artifacts

The development of bioinformatic solutions is guided by biological knowledge of the subject. In some cases, we use unambiguous biological models, while in others we rely on assumptions. A commonly used assumption for genomes is that related species have similar genome sequences. This is even more obvious in the case of chloroplast genomes due to their slow evolution. We investigated whether the lengths of complete chloroplast sequences are closely related to the taxonomic proximity of the species. The study was performed using all available RefSeq sequences from the asterid and rosid clades. In general, chloroplast length distributions are narrow at both the family and genus levels. In addition, clear biological explanations have already been reported for families and genera that exhibit particularly wide distributions. The main factors responsible for the length variations are parasitic life forms, IR loss, IR expansions and contractions, and polyphyly. However, the presence of outliers in the distribution at the genus level is a strong indication of possible inaccuracies in sequence assembly.

The complete chloroplast genome of Striga asiatica (L.) Kuntze 1891 (Orobanchaceae), a hemiparasitic weed from Guangxi China

Striga asiatica (L.) Kuntze 1891 is a hemiparasitic plant native to Asia and Africa. It is invasive and causes yield losses in crops such as corn, rice and sorghum. Lack of chloroplast genomic data has limited research into its obligate parasitic lifestyle. In this study, the complete chloroplast genome of Striga asiatica was sequenced and characterized. It is a quadripartite structure with a total length of 191,085 bp and a GC content of 37.86%. It has a large single copy region (LSC) of 51,406 bp, a small single copy region (SSC) of 273 bp, and two copies of the reverse repeat sequence (IRA and IRB) of 69,703 bp. A total of 122 protein-coding genes, 8 rRNA genes, and 44 tRNA genes were annotated in the chloroplast genome. There were a lot of ndh gene deletions and pseudogenizations in this chloroplast genome. For example, ndhA, D, E, H, I, and K were all pseudogenes because they were missing the 5' end start codon. ndhB, C, and J had shorter gene lengths than their homologs, and ndhF and ndhG were missing genes. The phylogenetic tree reveals that all Striga species form a clade, and a bootstrap value of 100 indicates that S. asiatica is closely related to Striga hermonthica and Striga sepera. The comprehensive chloroplast genomic resource of S. asiatica would assist researchers in comprehending hemiparasitic mechanisms, molecular markers, and evolutionary patterns of the genus Striga.

ITS1 Barcode and Phytochemical Analysis by Gas Chromatography-Mass Spectrometry of Corynaea crassa Hook. f (Balanophoraceae) from Ecuador and Peru

The use of medicinal plants is the basis of traditional healthcare. Recently, the use of herbal medicine has been increasing among consumers due to availability, economy, and less side effect. For instance, the hemiparasite plant Corynaea crassa has medicinal properties and could be found in some regions of America, from Costa Rica to Bolivia. Phytochemical and genetic characterization of medicinal plants is needed for proper identification of metabolites responsible for medicinal properties and for genotyping, respectively. Moreover, characterization of medicinal plants through the use of DNA barcodes is an important tool for phylogenetic analysis and identification of species; furthermore, complemented with phytochemical analysis, both are useful for identification of plant species and quality control of medicinal products. The objective of this study was to analyze the species of C. crassa collected in Ecuador and Peru from the phylogenetic and phytochemical point of view. Polymerase chain reaction (PCR) was performed for amplification of the internal transcribed spacer 1 (ITS1) region after DNA extraction of samples of C. crassa. Blast analysis was performed in the GenBank database with the ITS1 sequences obtained from two accessions of C. crassa from Ecuador (GenBank accession numbers OM471920 and OM471919 for isolates CIBE-17 and CIBE-18, respectively) and three from Peru (GenBank accession numbers OM471921, OM471922, and OM471923 for isolates CIBE-13, CIBE-14, and CIBE-15, respectively). The accessions available in the GenBank were used for phylogenetic analysis. For the phytochemical analysis, hydroalcoholic extracts were obtained by maceration using 80% ethanol as solvent, followed by a derivatization process and analysis by gas chromatography-mass spectrometry. Based on the phylogenetic analysis of the C. crassa samples, the ITS1 sequence could be used to differentiate C. crassa of different locations. The samples of C. crassa from Ecuador and Peru are more similar between them than with other clades including Helosis spp. The phytochemical study revealed differences in the presence and relative abundance of some metabolites; mainly eugenol, 1,4-lactone arabinonic acid, dimethoxyrabelomycin and azelaic acid, which are reported for the first time for the species under study and the genus Corynaea. These results are the first findings on the combined analysis using genetic and phytochemical analysis for C. crassa, which could be used as a useful tool for quality control of the C. crassa species in medicinal products.

Large-scale mRNA transfer between Haloxylon ammodendron (Chenopodiaceae) and herbaceous root holoparasite Cistanche deserticola (Orobanchaceae)

Exchanges of mRNA were shown between host and stem parasites but not root parasites. Cistanche deserticola (Orobanchaceae) is a holoparasitic herb which parasitizes on the roots of woody plant Haloxylon ammodendron (Chenopodiaceae). We used transcriptome sequencing and bioinformatic analyses to identify nearly ten thousand mobile mRNAs. Transcript abundance appears to be a driving force for transfer event and mRNA exchanges occur through haustorial junction. Mobility of selected mRNAs was confirmed in situ and in sunflower-Orobanche cumana heterologous parasitic system. Four C. deserticola →H. ammodendron mobile mRNAs appear to facilitate haustorium development. Of interest, two mobile mRNAs of putative resistance genes CdNLR1 and CdNLR2 cause root-specific hypersensitive response and retard parasite development, which might contribute to parasitic equilibrium. The present study provides evidence for the large-scale mRNA transfer event between a woody host and a root parasite, and demonstrates the functional relevance of six C. deserticola genes in host-parasite interactions.

De novo full length transcriptome analysis and gene expression profiling to identify genes involved in phenylethanol glycosides biosynthesis in Cistanche tubulosa

Background

The dried stem of Cistanche, is a famous Chinese traditional medicine. The main active pharmacodynamic components are phenylethanol glycosides (PhGs). Cistanche tubulosa produces higher level of PhGs in its stems than that of Cistanche deserticola. However, the key genes in the PhGs biosynthesis pathway is not clear in C. tubulosa.

Results

In this study, we performed the full-length transcriptome sequencing and gene expression profiling of C. tubulosa using PacBio combined with BGISEQ-500 RNA-seq technology. Totally, 237,772 unique transcripts were obtained, ranging from 199 bp to 31,857 bp. Among the unique transcripts, 188,135 (79.12%) transcripts were annotated. Interestingly, 1080 transcripts were annotated as 22 enzymes related to PhGs biosynthesis. We measured the content of echinacoside, acteoside and total PhGs at two development stages, and found that the content of PhGs was 46.74% of dry matter in young fleshy stem (YS1) and then decreased to 31.22% at the harvest stage (HS2). To compare with YS1, 13,631 genes were up-regulated, and 15,521 genes were down regulated in HS2. Many differentially expressed genes (DEGs) were identified to be involved in phenylpropanoid biosynthesis pathway, phenylalanine metabolism pathway, and tyrosine metabolism pathway.

Conclusions

This is the first report of transcriptome study of C. tubulosa which provided the foundation for understanding of PhGs biosynthesis. Based on these results, we proposed a potential model for PhGs biosynthesis in C. tubulosa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08921-x.

Structural mutations of small single copy (SSC) region in the plastid genomes of five Cistanche species and inter-species identification

BACKGROUND

Cistanche is an important genus of Orobanchaceae, with critical medicinal, economic, and desertification control values. However, the phylogenetic relationships of Cistanche genus remained obscure. To date, no effective molecular markers have been reported to discriminate effectively the Cistanche closely related species reported here. In this study, we obtained and characterized the plastomes of four Cistanche species from China, to clarify the phylogenetic relationship within the genus, and to develop molecular markers for species discrimination.  RESULTS: Four Cistanche species (Cistanche deserticola, Cistanche salsa, Cistanche tubulosa and Cistanche sinensis), were deep-sequenced with Illumina. Their plastomes were assembled using SPAdes and annotated using CPGAVAS2. The plastic genomes were analyzed in detail, finding that all showed the conserved quadripartite structure (LSC-IR-SSC-IR) and with full sizes ranging from 75 to 111 Kbp. We observed a significant contraction of small single copy region (SSC, ranging from 0.4-29 Kbp) and expansion of inverted repeat region (IR, ranging from 6-30 Kbp), with C. deserticola and C. salsa showing the smallest SSCs with only one gene (rpl32). Compared with other Orobanchaceae species, Cistanche species showed extremely high rates of gene loss and pseudogenization, as reported for other parasitic Orobanchaceae species. Furthermore, analysis of sequence divergence on protein-coding genes showed the three genes (rpl22, clpP and ycf2) had undergone positive selection in the Cistanche species under study. In addition, by comparison of all available Cistanche plastomes we found 25 highly divergent intergenic spacer (IGS) regions that were used to predict two DNA barcode markers (Cis-mk01 and Cis-mk02 based on IGS region trnR-ACG-trnN-GUU) and eleven specific DNA barcode markers using Ecoprimer software. Experimental validation showed 100% species discrimination success rate with both type of markers.

CONCLUSION

Our findings have shown that Cistanche species are an ideal model to investigate the structure variation, gene loss and pseudogenization during the process of plastome evolution in parasitic species, providing new insights into the evolutionary relationships among the Cistanche species. In addition, the developed DNA barcodes markers allow the proper species identification, ensuring the effective and safe use of Cistanche species as medicinal products.

Analysis of microbial diversity and community structure of rhizosphere soil of Cistanche salsa from different host plants

Host plants influence rhizosphere microorganism composition through root secretions, and rhizosphere associated microorganisms influence Cistanche seeds germination. At present, little is known about effects of different host plants on soil bacteria and fungi in the rhizosphere of Cistanche salsa. High-throughput sequencing was used here to reveal the similarities and differences in the structural composition of the soil microbial community of C. salsa from six host plants (i.e., Halocnemum strobilaceum, Atriplex patens, Kalidium foliatum, Caroxylon passerinum, Anabasis aphylla, Krascheninnikovia ceratoides). We discovered that Krascheninnikovia ceratoides-parasitizing C. salsa (YRCR6) had the highest diversity of rhizosphere bacterial communities, and Anabasis aphylla -parasitizing C. salsa (YRCR5) had the highest diversity of rhizosphere fungal communities. Fungal communities were more influenced by the host plant than bacterial communities. In addition, we discovered certain rhizosphere microorganisms that may be associated with Cistanche seeds germination, including Mortierella, Aspergillus alliaceus, and Cladosporium, which are account for a relatively high proportion in Halocnemum strobilaceum, Atriplex patens and Anabasis aphylla -parasitizing C. salsa. Redundancy analysis results also revealed that AP, HCO₃^–, pH, Ca²⁺, SO₄^2–, and K⁺ had a highly significant impact on the bacterial community structure (P < 0.01), while pH and SO₄^2– had a significant impact on the fungal community structure (P < 0.05). Conclusively, differences were noted in the structure of rhizosphere bacterial and fungal communities of C. salsa parasitizing different plants in the same habit and the difference may be related to the host plant. This result can provide a new ideas for the selection of host plants and the cultivation of C. salsa.

Chromosome-level genome assembly of a xerophytic plant, Haloxylon ammodendron

Haloxylon ammodendron is a xerophytic perennial shrub or small tree that has a high ecological value in anti-desertification due to its high tolerance to drought and salt stress. Here, we report a high-quality, chromosome-level genome assembly of H. ammodendron by integrating PacBio's high-fidelity sequencing and Hi-C technology. The assembled genome size was 685.4 Mb, of which 99.6% was assigned to nine pseudochromosomes with a contig N50 value of 23.6 Mb. Evolutionary analysis showed that both the recent substantial amplification of long terminal repeat retrotransposons and tandem gene duplication may have contributed to its genome size expansion and arid adaptation. An ample amount of low-GC genes was closely related to functions that may contribute to the desert adaptation of H. ammodendron. Gene family clustering together with gene expression analysis identified differentially expressed genes that may play important roles in the direct response of H. ammodendron to water-deficit stress. We also identified several genes possibly related to the degraded scaly leaves and well-developed root system of H. ammodendron. The reference-level genome assembly presented here will provide a valuable genomic resource for studying the genome evolution of xerophytic plants, as well as for further genetic breeding studies of H. ammodendron.

Predicting the suitable habitats of parasitic desert species based on a niche model with Haloxylon ammodendron and Cistanche deserticola as examples

Haloxylon ammodendron, an excellent tree species for sand fixation and afforestation in the desert areas of western China, is threatened by climate change and anthropogenic activities. The suitable habitat of this species is shrinking at a remarkable rate, although conservation measures have been implemented. Cistanche deserticola is an entirely parasitic herb that occurs in deserts, is a source of "desert ginseng" worldwide, and has extremely high medicinal value. Little is known about using niche models to simulate habitat suitability and evaluate important environmental variables related to parasitic species. In this study, we modeled the current suitable habitat of H. ammodendron and C. deserticola by MaxEnt based on occurrence record data of the distributions of these two species in China. We grouped H. ammodendron and C. deserticola into three groups according to the characteristics of parasitic species and modeled them with environmental factors. The results showed that bioclimate was the most important environmental parameter affecting the H. ammodendron and C. deserticola distribution. Precipitations, such as annual precipitation, precipitation seasonality, and precipitation in the driest quarter, were identified as the most critical parameters. The slope, diurnal temperature range, water vapor pressure, ground-frost frequency, and solar radiation also substantially contributed to the distribution of the two species. The proportions of the most suitable areas for Groups 1, 2, and 3 were 1.2%, 1.3%, and 1.7%, respectively, in China. When combined with cultural geography, five hot spot conservation areas were determined within the distribution of H. ammodendron and C. deserticola. The comprehensive analysis indicated that by using MaxEnt to model the suitable habitat of parasitic species, we further improved the accuracy of the prediction and coupled the error of the distribution of a single species. This study provides a useful reference for the protection of H. ammodendron forests and the management of C. deserticola plantations.

Comparative and phylogenetic analyses of the chloroplast genomes of species of Paeoniaceae

Plants belonging to family Paeoniaceae are not only economically important ornamental plants but also medicinal plants used as an important source of traditional Chinese medicine. Owing to the complex network evolution and polyploidy evolution of this family, its systematics and taxonomy are controversial and require a detailed investigation. In this study, three complete chloroplast genomes of sect. Paeonia, one of the sections of Paeonia, were sequenced and then analysed together with 16 other published chloroplast genomes of Paeoniaceae species. The total lengths of the chloroplast genomes of these species were 152,153–154,405 bp. A total of 82–87 protein-coding genes, 31–40 tRNA genes and 8 rRNA genes were annotated. Bioinformatics analysis revealed 61–74 simple sequence repeats (SSRs) in the chloroplast genomes, most of which have A/T base preference. Codon usage analysis showed that A/U-ending codons were more positive than C/G-ending codons, and a slight bias in codon usage was observed in these species. A comparative analysis of these 19 species of Paeoniaceae was then conducted. Fourteen highly variable regions were selected for species relationship study. Phylogenetic analysis revealed that the species of sect. Paeonia gathered in one branch and then divided into different small branches. P. lactiflora, P. anomala, P. anomala subsp. veitchii and P. mairei clustered together. P. intermedia was related to P. obovata and P. obovata subsp. willmottiae. P. emodi was the sister to all other species in the sect. Paeonia.

Identification and phylogenetic analysis of five Crataegus species (Rosaceae) based on complete chloroplast genomes

The chloroplast genomes of the five Crataegus species were shown to have a conserved genome structure. Complete chloroplast genome sequences were more suitable than highly variable regions for the identification and phylogenetic analysis of Crataegus species. Hawthorn, which is commonly used as a traditional Chinese medicine, is one of the most popular sour fruits and has high economic value. Crataegus pinnatifida var. pinnatifida and C. pinnatifida var. major are frequently adulterated with other Crataegus species on the herbal medicine market. However, most Crataegus plants are difficult to identify using traditional morphological methods. Here, we compared five Crataegus chloroplast (CP) genomes comprising two newly sequenced (i.e., C. pinnatifida var. pinnatifida and C. pinnatifida var. major) and three previously published CP genomes. The CP genomes of the five Crataegus species had a conserved genome structure, gene content and codon usage. The total length of the CP genomes was 159,654-159,865 bp. A total of 129-130 genes, including 84-85 protein-coding genes, 37 tRNA genes and 8 rRNA genes, were annotated. Bioinformatics analysis revealed 96-103 simple sequence repeats (SSRs) and 48-70 long repeats in the five CP genomes. Combining the results of mVISTA and nucleotide diversity, five highly variable regions were screened for species identification and relationship studies. Maximum likelihood trees were constructed on the basis of complete CP genome sequences and highly variable regions. The results showed that the former had higher discriminatory power for Crataegus species, indicating that the complete CP genome could be used as a super-barcode to accurately authenticate the five Crataegus species.

Complete Plastid and Mitochondrial Genomes of Aeginetia indica Reveal Intracellular Gene Transfer (IGT), Horizontal Gene Transfer (HGT), and Cytoplasmic Male Sterility (CMS)

Orobanchaceae have become a model group for studies on the evolution of parasitic flowering plants, and Aeginetia indica, a holoparasitic plant, is a member of this family. In this study, we assembled the complete chloroplast and mitochondrial genomes of A. indica. The chloroplast and mitochondrial genomes were 56,381 bp and 401,628 bp long, respectively. The chloroplast genome of A. indica shows massive plastid genes and the loss of one IR (inverted repeat). A comparison of the A. indica chloroplast genome sequence with that of a previous study demonstrated that the two chloroplast genomes encode a similar number of proteins (except atpH) but differ greatly in length. The A. indica mitochondrial genome has 53 genes, including 35 protein-coding genes (34 native mitochondrial genes and one chloroplast gene), 15 tRNA (11 native mitochondrial genes and four chloroplast genes) genes, and three rRNA genes. Evidence for intracellular gene transfer (IGT) and horizontal gene transfer (HGT) was obtained for plastid and mitochondrial genomes. ψndhB and ψcemA in the A. indica mitogenome were transferred from the plastid genome of A. indica. The atpH gene in the plastid of A. indica was transferred from another plastid angiosperm plastid and the atpI gene in mitogenome A. indica was transferred from a host plant like Miscanthus siensis. Cox2 (orf43) encodes proteins containing a membrane domain, making ORF (Open Reading Frame) the most likely candidate gene for CMS development in A. indica.

The Complete Chloroplast Genome of Plukenetia volubilis Provides Insights Into the Organelle Inheritance

Plukenetia volubilis L. (Malpighiales: Euphorbiaceae), also known as Sacha inchi, is considered a promising crop due to its high seed content of unsaturated fatty acids (UFAs), all of them highly valuable for food and cosmetic industries, but the genetic basis of oil biosynthesis of this non-model plant is still insufficient. Here, we sequenced the total DNA of Sacha inchi by using Illumina and Nanopore technologies and approached a de novo reconstruction of the whole nucleotide sequence and the organization of its 164,111 bp length of the chloroplast genome, displaying two copies of an inverted repeat sequence [inverted repeat A (IRA) and inverted repeat B (IRB)] of 28,209 bp, each one separating a small single copy (SSC) region of 17,860 bp and a large single copy (LSC) region of 89,833 bp. We detected two large inversions on the chloroplast genome that were not presented in the previously reported sequence and studied a promising cpDNA marker, useful in phylogenetic approaches. This chloroplast DNA (cpDNA) marker was used on a set of five distinct Colombian cultivars of P. volubilis from different geographical locations to reveal their phylogenetic relationships. Thus, we evaluated if it has enough resolution to genotype cultivars, intending to crossbreed parents and following marker's trace down to the F1 generation. We finally elucidated, by using molecular and cytological methods on cut flower buds, that the inheritance mode of P. volubilis cpDNA is maternally transmitted and proposed that it occurs as long as it is physically excluded during pollen development. This de novo chloroplast genome will provide a valuable resource for studying this promising crop, allowing the determination of the organellar inheritance mechanism of some critical phenotypic traits and enabling the use of genetic engineering in breeding programs to develop new varieties.

Folk nomenclature of plants in Cistanche deserticola-associated community in South Gobi, Mongolia

Cistanche deserticola is an important medicinal plant in Mongolia. Despite its significant role in local healing systems, little traditional knowledge had been reported. The present study investigated folk names of C. deserticola and other species of the same community in Umnugobi Province, South Gobi region of Mongolia, based on ethnobotanical approaches. The high correspondence between folk names and scientific names of plant species occurring in Cistanche-associated community shows the scientific meaning of folk nomenclature and classification in Mongolia. The Mongolian and folk names of plants were formed on the basis of observations and understanding of wild plants including their morphology, phenology and traditional uses as well. Results from this study will support the conservation of C. deserticola itself, a rare and endangered plant species listed in the Monglian Red Data Book. Our documentation of folk nomenclature based on 96 plant species in the Cistanche community, as a part of traditional knowledge associated with biodiversity, will be very helpful for making strategy of plant biodiversity conservation in Mongolia.

Extensive gene loss in the plastome of holoparasitic plant Cistanche tubulosa (Orobanchaceae)

Extensive photosynthetic gene loss and rapid evolutionary rate occur in the plastomes of parasitic plants. The holoparasitic plant Cistanche tubulosa of Orobanchaceae is an important medicinal resources that are distributed in arid areas. In this study, the complete plastome of C. tubulosa has been sequenced, assembled and analyzed. The total plastome of C. tubulosa was 75,375 bp in length, consisting of a pair of inverted repeats (IRs, 6,593 bp), a large single-copy region (LSC, 32,470 bp) and a small single-copy region (SSC, 29,719 bp). It contained 24 intact protein coding genes, nine pseudogenes, and 44 missing genes. In addition, all the protein-coding genes, which were related to photosynthesis and energy production, were pseudogenised or lost. Four rRNA genes and 24 tRNA genes were intact meanwhile five tRNA genes were missing. Phylogenetic tree indicated that C. tubulosa was closely related to C. phelypaea. Our results may improve understanding of the plastome organization, classification, and evolution of parasitic plants.

Comparing and phylogenetic analysis chloroplast genome of three Achyranthes species

In this study, the chloroplast genome sequencing of the Achyranthes longifolia, Achyranthes bidentata and Achyranthes aspera were performed by Next-generation sequencing technology. The results revealed that there were a length of 151,520 bp (A. longifolia), 151,284 bp (A. bidentata), 151,486 bp (A. aspera), respectively. These chloroplast genome have a highly conserved structure with a pair of inverted repeat (IR) regions (25,150 bp; 25,145 bp; 25,150 bp), a large single copy (LSC) regions (83,732 bp; 83,933 bp; 83,966 bp) and a small single copy (SSC) regions (17,252 bp; 17,263 bp; 17,254 bp) in A. bidentate, A. aspera and A. longifolia. There were 127 genes were annotated, which including 8 rRNA genes, 37 tRNA genes and 82 functional genes. The phylogenetic analysis strongly revealed that Achyranthes is monophyletic, and A. bidentata was the closest relationship with A. aspera and A. longifolia. A. bidentata and A. longifolia were clustered together, the three Achyranthes species had the same origin, then the gunes of Achyranthes is the closest relative to Alternanthera, and that forms a group with Alternanthera philoxeroides. The research laid a foundation and provided relevant basis for the identification of germplasm resources in the future.

A multilocus phylogeny of the non-photosynthetic parasitic plant Cistanche (Orobanchaceae) refutes current taxonomy and identifies four major morphologically distinct clades

Phylogenetic relationships of and within non-photosynthetic parasitic lineages are notoriously poorly known, which negatively affects our understanding of parasitic plants. This is also the case for Cistanche (Orobanchaceae), an Old World genus with about two dozen species, whose relationships have not yet been addressed using molecular phylogenetic approaches. Here we infer phylogenetic relationships within the genus, employing a taxonomically and geographically broad sampling covering all previously distinguished infrageneric groups and most of the currently recognized species. A combined matrix of three plastid markers (trnL-trnF, including the trnL intron and the intergenic spacer (IGS), trnS-trnfM IGS and psbA-trnH IGS) and one nuclear marker (ITS) was analyzed using maximum parsimony, maximum likelihood and Bayesian inference. Cistanche falls into four well-supported and geographically differentiated clades: East Asian Clade, Northwest African Clade, Southwest Asian Clade and Widespread Clade. Of those, only the East Asian Clade corresponds to a previously recognized taxonomic section, whereas the others either contain members of two or three sections (Widespread Clade and Southwest Asian Clade, respectively) or have not been taxonomically recognized so far (Northwest African Clade). Whereas the Southwest Asian Clade exhibits strong phylogenetic structure among and partly within species (the East Asian Clade and the Northwest African Clade are monospecific), phylogenetic resolution within the Widespread Clade is often low and hampered by discrepancies between nuclear and plastid markers. Both molecular and morphological data indicate that species diversity in Cistanche is currently underestimated.

The Complete Plastomes of Five Hemiparasitic Plants (Osyris wightiana, Pyrularia edulis, Santalum album, Viscum liquidambaricolum, and V. ovalifolium): Comparative and Evolutionary Analyses Within Santalales

Most species of Santalales (the sandalwood order) are hemiparasites, including both facultative and obligate hemiparasites. Despite its rich diversity, only a small fraction of the species in the sandalwood order have sequenced plastomes. The evolution of parasitism-associated plastome reduction in Santalales remains under-studied. Here, we report the complete plastomes of three facultative hemiparasites (Pyrularia edulis, Cervantesiaceae; Osyris wightiana, and Santalum album, Santalaceae), and two obligate hemiparasites (Viscum liquidambaricolum and Viscum ovalifolium, Viscaceae). Coupled with publicly available data, we investigated the dynamics of plastome degradation in Santalales hemiparasites. Our results indicate that these hemiparasites can be characterized by various degrees of plastome downsizing, structural rearrangement, and gene loss. The loss or pseudogenization of ndh genes was commonly observed in Santalales hemiparasites, which may be correlated to the lifestyle shift from photoautotroph to hemiparasitism. However, the obligate hemiparasites did not exhibit a consistently higher level of gene loss or pseudogenization compared to facultative hemiparasites, which suggests that the degree of plastome reduction is not correlated with the trophic level facultative or obligate hemiparasitism. Instead, closely related taxa tend to possess highly similar plastome size, structure, and gene content. This implies the parasitism-associated plastome degradation in Santalales may evolve in a lineage-specific manner.

The Complete Chloroplast Genome of Two Important Annual Clover Species, Trifolium alexandrinum and T. resupinatum: Genome Structure, Comparative Analyses and Phylogenetic Relationships with Relatives in Leguminosae

Trifolium L., which belongs to the IR lacking clade (IRLC), is one of the largest genera in the Leguminosae and contains several economically important fodder species. Here, we present whole chloroplast (cp) genome sequencing and annotation of two important annual grasses, Trifolium alexandrinum (Egyptian clover) and T. resupinatum (Persian clover). Abundant single nucleotide polymorphisms (SNPs) and insertions/deletions (In/Dels) were discovered between those two species. Global alignment of T. alexandrinum and T. resupinatum to a further thirteen Trifolium species revealed a large amount of rearrangement and repetitive events in these fifteen species. As hypothetical cp open reading frame (ORF) and RNA polymerase subunits, ycf1 and rpoC2 in the cp genomes both contain vast repetitive sequences and observed high Pi values (0.7008, 0.3982) between T. alexandrinum and T. resupinatum. Thus they could be considered as the candidate genes for phylogenetic analysis of Trifolium species. In addition, the divergence time of those IR lacking Trifolium species ranged from 84.8505 Mya to 4.7720 Mya. This study will provide insight into the evolution of Trifolium species.

Diverse trajectories of plastome degradation in holoparasitic Cistanche and genomic location of the lost plastid genes

The plastid genomes (plastomes) of non-photosynthetic plants generally undergo gene loss and pseudogenization. Despite massive plastomes reported in different parasitism types of the broomrape family (Orobanchaceae), more plastomes representing different degradation patterns in a single genus are expected to be explored. Here, we sequence and assemble the complete plastomes of three holoparasitic Cistanche species (C. salsa, C. mongolica, and C. sinensis) and compare them with the available plastomes of Orobanchaceae. We identified that the diverse degradation trajectories under purifying selection existed among three Cistanche clades, showing obvious size differences in the entire plastome, long single copy region, and non-coding region, and different patterns of the retention/loss of functional genes. With few exceptions of putatively functional genes, massive plastid fragments, which have been lost and transferred into the mitochondrial or nuclear genomes, are non-functional. In contrast to the equivalents of the Orobanche species, some plastid-derived genes with diverse genomic locations are found in Cistanche. The early and initially diverged clades in different genera such as Cistanche and Aphyllon possess obvious patterns of plastome degradation, suggesting that such key lineages should be considered prior to comparative analysis of plastome evolution, especially in the same genus.

Gene Losses and Variations in Chloroplast Genome of Parasitic Plant Macrosolen and Phylogenetic Relationships within Santalales

Macrosolen plants are parasitic shrubs, several of which are important medicinal plants, that are used as folk medicine in some provinces of China. However, reports on Macrosolen are limited. In this study, the complete chloroplast genome sequences of Macrosolen cochinchinensis, Macrosolen tricolor and Macrosolen bibracteolatus are reported. The chloroplast genomes were sequenced by Illumina HiSeq X. The length of the chloroplast genomes ranged from 129,570 bp (M. cochinchinensis) to 126,621 bp (M. tricolor), with a total of 113 genes, including 35 tRNA, eight rRNA, 68 protein-coding genes, and two pseudogenes (ycf1 and rpl2). The simple sequence repeats are mainly comprised of A/T mononucleotide repeats. Comparative genome analyses of the three species detected the most divergent regions in the non-coding spacers. Phylogenetic analyses using maximum parsimony and maximum likelihood strongly supported the idea that Loranthaceae and Viscaceae are monophyletic clades. The data obtained in this study are beneficial for further investigations of Macrosolen in respect to evolution and molecular identification.

Comparative analysis of the complete chloroplast genomes of seven Populus species: Insights into alternative female parents of Populus tomentosa

Populus tomentosa, of section Populus, is distributed mainly in northern China. This species has high resistance to many diseases and insects, and it plays key roles in shelterbelts and urban afforestation in northern China. It has long been suspected to be a hybrid, but its parents remain unknown. In the present study, we report four newly sequenced complete cp genomes from section Populus and comparative genomic analyses of these new sequences and three published cp genome sequences. The seven cp genomes ranged from 155,853 bp (P. tremula var. davidiana) to 156,746 bp (P. adenopoda) in length, and their gene orders, gene numbers and GC contents were similar. We analyzed SNPs, indels, SSRs and repeats among the seven cp genomes, and eight small inversions were detected in the ndhC-trnV, rbcL-accD, petA-psbJ, trnW-trnP, rpl16-rps3, trnL-ycf15, ycf15-trnL, and ndhF-trnL intergenic regions. Furthermore, seven divergent regions (trnH-psbA, matK, psbM-psbD, ndhC-trnV, ycf1, ndhF-ccsA and ccsA-ndhD) were found in more highly variable regions. The phylogenetic tree reveals that P. tomentosa is closely related to P. alba and P. alba var. pyramidalis. Hence, P. alba was involved in the formation of P. tomentosa.

Chloroplast genome analysis of box-ironbark Eucalyptus

Eucalyptus L'Hérit. (Myrtaceae) is a taxonomically complex and highly speciose genus that dominates much of Australia's woody vegetation. However, very little information is available about the molecular biology and chloroplast diversity of certain groups, such as Eucalyptus section Adnataria, which is found in many woodland habitats of eastern Australia. We report four new complete chloroplast genomes of Eucalyptus, including three genomes from species previously lacking any chloroplast reference sequences. Plastomes of E. albens, E. conica, E. crebra and E. melliodora assembled using a de novo approach were shown to be largely identical to each other, and similar in size and structure to previously published chloroplast genomes from Eucalyptus. A total of 132 genes (114 single-copy genes and 18 duplicated genes in the IR regions) were identified, and shown to be highly conserved in terms of gene order, content and organization. Slightly higher divergence in the intergenic spacers was identified through comparative genomic analyses. Chloroplast sequences of 35 additional individuals representing 12 species were assembled using a reference guided approach. Rates of nucleotide substitution varied among the protein coding genes, with 17 genes under possible positive selection, and 29 invariant genes. Phylogenetic analysis of either the whole reconstructed plastome sequences or the individual genes revealed extreme discordance with expected species boundaries or higher-level relationships. Plastome relationships were better predicted by geography than by nuclear DNA or taxonomic relationships, suggesting a substantial influence of gene flow over and above the effects of incomplete lineage sorting. These results provide resources for future research and valuable insights into the prevalence of interspecific gene flow among Eucalyptus species.

Convergent Plastome Evolution and Gene Loss in Holoparasitic Lennoaceae

The Lennoaceae, a small monophyletic plant family of root parasites endemic to the Americas, are one of the last remaining independently evolved lineages of parasitic angiosperms lacking a published plastome. In this study, we present the assembled and annotated plastomes of two species spanning the crown node of Lennoaceae, Lennoa madreporoides and Pholisma arenarium, as well as their close autotrophic relative from the sister family Ehretiaceae, Tiquilia plicata. We find that the plastomes of L. madreporoides and P. arenarium are similar in size and gene content, and substantially reduced compared to T. plicata, consistent with trends seen in other holoparasitic lineages. In particular, most plastid genes involved in photosynthesis function have been lost, whereas housekeeping genes (ribosomal protein-coding genes, rRNAs, and tRNAs) are retained. One notable exception is the persistence of a rbcL open reading frame in P. arenarium but not L. madreporoides suggesting a nonphotosynthetic function for this gene. Of the retained coding genes, d_N/d_S ratios indicate that some remain under purifying selection, whereas others show relaxed selection. Overall, this study supports the mounting evidence for convergent plastome evolution in flowering plants following the shift to heterotrophy.

RNA-seq highlights parallel and contrasting patterns in the evolution of the nuclear genome of fully mycoheterotrophic plants

BACKGROUND

While photosynthesis is the most notable trait of plants, several lineages of plants (so-called full heterotrophs) have adapted to obtain organic compounds from other sources. The switch to heterotrophy leads to profound changes at the morphological, physiological and genomic levels.

RESULTS

Here, we characterize the transcriptomes of three species representing two lineages of mycoheterotrophic plants: orchids (Epipogium aphyllum and Epipogium roseum) and Ericaceae (Hypopitys monotropa). Comparative analysis is used to highlight the parallelism between distantly related fully heterotrophic plants. In both lineages, we observed genome-wide elimination of nuclear genes that encode proteins related to photosynthesis, while systems associated with protein import to plastids as well as plastid transcription and translation remain active. Genes encoding components of plastid ribosomes that have been lost from the plastid genomes have not been transferred to the nuclear genomes; instead, some of the encoded proteins have been substituted by homologs. The nuclear genes of both Epipogium species accumulated nucleotide substitutions twice as rapidly as their photosynthetic relatives; in contrast, no increase in the substitution rate was observed in H. monotropa.

CONCLUSIONS

Full heterotrophy leads to profound changes in nuclear gene content. The observed increase in the rate of nucleotide substitutions is lineage specific, rather than a universal phenomenon among non-photosynthetic plants.

Understanding plastome evolution in Hemiparasitic Santalales: Complete chloroplast genomes of three species, Dendrotrophe varians, Helixanthera parasitica, and Macrosolen cochinchinensis

Santalales is a large order, with over 2200 species, most of which are root or aerial (stem) hemiparasites. In this study, we report the newly assembled chloroplast genome of Dendrotrophe varians (140,666 bp) in the family Amphorogynaceae and the cp genomes of Helixanthera parasitica (124,881 bp) and Macrosolen cochinchinensis (122,986 bp), both in the family Loranthaceae. We compared the cp genomes of 11 Santalales including eight currently available cp genomes. Santalales cp genomes are slightly or not reduced in size (119-147 kb), similar to other hemiparasitic species, when compared with typical angiosperm cp genomes (120-170 kb). In a phylogeny examining gene content, the NADH dehydrogenase gene group is the only one among eight functional gene groups that lost complete functionally in all examined Santalales. This supports the idea that the functional loss of ndh genes is the initial stage in the evolution of the plastome of parasitic plants, but the loss has occurred independently multiple times in angiosperms, while they are not found in some parasites. This suggests that the functional loss of ndh genes is not essential for the transition from autotroph to parasite. We additionally examined the correlation between gene content and type of parasitism (obligate/facultative and stem/root parasites) of all hemiparasitic species in which cp genomes have been reported to date. Correlation was not found in any types of parasitism.

Dense infraspecific sampling reveals rapid and independent trajectories of plastome degradation in a heterotrophic orchid complex

Heterotrophic plants provide excellent opportunities to study the effects of altered selective regimes on genome evolution. Plastid genome (plastome) studies in heterotrophic plants are often based on one or a few highly divergent species or sequences as representatives of an entire lineage, thus missing important evolutionary-transitory events. Here, we present the first infraspecific analysis of plastome evolution in any heterotrophic plant. By combining genome skimming and targeted sequence capture, we address hypotheses on the degree and rate of plastome degradation in a complex of leafless orchids (Corallorhiza striata) across its geographic range. Plastomes provide strong support for relationships and evidence of reciprocal monophyly between C. involuta and the endangered C. bentleyi. Plastome degradation is extensive, occurring rapidly over a few million years, with evidence of differing rates of genomic change among the two principal clades of the complex. Genome skimming and targeted sequence capture differ widely in coverage depth overall, with depth in targeted sequence capture datasets varying immensely across the plastome as a function of GC content. These findings will help to fill a knowledge gap in models of heterotrophic plastid genome evolution, and have implications for future studies in heterotrophs.

Characterization and Comparative Analysis of the Complete Chloroplast Genome of the Critically Endangered Species Streptocarpus teitensis (Gesneriaceae)

Streptocarpus teitensis (Gesneriaceae) is an endemic species listed as critically endangered in the International Union for Conservation of Nature (IUCN) red list of threatened species. However, the sequence and genome information of this species remains to be limited. In this article, we present the complete chloroplast genome structure of Streptocarpus teitensis and its evolution inferred through comparative studies with other related species. S. teitensis displayed a chloroplast genome size of 153,207 bp, sheltering a pair of inverted repeats (IR) of 25,402 bp each split by small and large single-copy (SSC and LSC) regions of 18,300 and 84,103 bp, respectively. The chloroplast genome was observed to contain 116 unique genes, of which 80 are protein-coding, 32 are transfer RNAs, and four are ribosomal RNAs. In addition, a total of 196 SSR markers were detected in the chloroplast genome of Streptocarpus teitensis with mononucleotides (57.1%) being the majority, followed by trinucleotides (33.2%) and dinucleotides and tetranucleotides (both 4.1%), and pentanucleotides being the least (1.5%). Genome alignment indicated that this genome was comparable to other sequenced members of order Lamiales. The phylogenetic analysis suggested that Streptocarpus teitensis is closely related to Lysionotus pauciflorus and Dorcoceras hygrometricum.

Sequencing and Analysis of Strobilanthes cusia (Nees) Kuntze Chloroplast Genome Revealed the Rare Simultaneous Contraction and Expansion of the Inverted Repeat Region in Angiosperm

Ban-Lan-Gen, the root tissues derived from several morphologically indistinguishable plant species, have been used widely in traditional Chinese medicines for numerous years. The identification of reliable markers to distinguish various source plant species is critical for the effective and safe use of products containing Ban-Lan-Gen. Here, we analyzed and characterized the complete chloroplast (cp) genome sequence of Strobilanthes cusia (Nees) Kuntze to identify high-resolution markers for the species determination of Southern Ban-Lan-Gen. Total DNA was extracted and subjected to next-generation sequencing. The cp genome was then assembled, and the gaps were filled using PCR amplification and Sanger sequencing. Genome annotation was conducted using CpGAVAS web server. The genome was 144,133 bp in length, presenting a typical quadripartite structure of large (LSC; 91,666 bp) and small (SSC; 17,328 bp) single-copy regions separated by a pair of inverted repeats (IRs; 17,811 bp). The genome encodes 113 unique genes, including 79 protein-coding, 30 transfer RNA, and 4 ribosomal RNA genes. A total of 20 tandem, 2 forward, and 6 palindromic repeats were detected in the genome. A phylogenetic analysis based on 65 protein-coding genes showed that S. cusia was closely related to Andrographis paniculata and Ruellia breedlovei, which belong to the same family, Acanthaceae. One interesting feature is that the IR regions apparently undergo simultaneous contraction and expansion, resulting in the presence of single copies of rps19, rpl2, rpl23, and ycf2 in the LSC region and the duplication of psbA and trnH genes in the IRs. This study provides the first complete cp genome in the genus Strobilanthes, containing critical information for the classification of various Strobilanthes species in the future. This study also provides the foundation for precisely determining the plant sources of Ban-Lan-Gen.

Herba Cistanches: Anti-aging

The Cistanche species ("Rou Cong Rong" in Chinese) is an endangered wild species growing in arid or semi-arid areas. The dried fleshy stem of Cistanches has been used as a tonic in China for many years. Modern pharmacological studies have since demonstrated that Herba Cistanches possesses broad medicinal functions, especially for use in anti-senescence, anti-oxidation, neuroprotection, anti-inflammation, hepatoprotection, immunomodulation, anti-neoplastic, anti-osteoporosis and the promotion of bone formation. This review summarizes the up-to-date and comprehensive information on Herba Cistanches covering the aspects of the botany, traditional uses, phytochemistry and pharmacology, to lay ground for fully elucidating the potential mechanisms of Herba Cistanches' anti-aging effect and promote its clinical application as an anti-aging herbal medicine.

Independent degradation in genes of the plastid ndh gene family in species of the orchid genus Cymbidium (Orchidaceae; Epidendroideae)

In this paper, we compare ndh genes in the plastid genome of many Cymbidium species and three closely related taxa in Orchidaceae looking for evidence of ndh gene degradation. Among the 11 ndh genes, there were frequently large deletions in directly repeated or AT-rich regions. Variation in these degraded ndh genes occurs between individual plants, apparently at population levels in these Cymbidium species. It is likely that ndh gene transfers from the plastome to mitochondrial genome (chondriome) occurred independently in Orchidaceae and that ndh genes in the chondriome were also relatively recently transferred between distantly related species in Orchidaceae. Four variants of the ycf1-rpl32 region, which normally includes the ndhF genes in the plastome, were identified, and some Cymbidium species contained at least two copies of that region in their organellar genomes. The four ycf1-rpl32 variants seem to have a clear pattern of close relationships. Patterns of ndh degradation between closely related taxa and translocation of ndh genes to the chondriome in Cymbidium suggest that there have been multiple bidirectional intracellular gene transfers between two organellar genomes, which have produced different levels of ndh gene degradation among even closely related species.

Gene losses and partial deletion of small single-copy regions of the chloroplast genomes of two hemiparasitic Taxillus species

Numerous variations are known to occur in the chloroplast genomes of parasitic plants. We determined the complete chloroplast genome sequences of two hemiparasitic species, Taxillus chinensis and T. sutchuenensis, using Illumina and PacBio sequencing technologies. These species are the first members of the family Loranthaceae to be sequenced. The complete chloroplast genomes of T. chinensis and T. sutchuenensis comprise circular 121,363 and 122,562 bp-long molecules with quadripartite structures, respectively. Compared with the chloroplast genomes of Nicotiana tabacum and Osyris alba, all ndh genes as well as three ribosomal protein genes, seven tRNA genes, four ycf genes, and the infA gene of these two species have been lost. The results of the maximum likelihood and neighbor-joining phylogenetic trees strongly support the theory that Loranthaceae and Viscaceae are monophyletic clades. This research reveals the effect of a parasitic lifestyle on the chloroplast structure and genome content of T. chinensis and T. sutchuenensis, and enhances our understanding of the discrepancies in terms of assembly results between Illumina and PacBio.

Molecular Structure and Phylogenetic Analyses of Complete Chloroplast Genomes of Two Aristolochia Medicinal Species

The family Aristolochiaceae, comprising about 600 species of eight genera, is a unique plant family containing aristolochic acids (AAs). The complete chloroplast genome sequences of Aristolochia debilis and Aristolochia contorta are reported here. The results show that the complete chloroplast genomes of A. debilis and A. contorta comprise circular 159,793 and 160,576 bp-long molecules, respectively and have typical quadripartite structures. The GC contents of both species were 38.3% each. A total of 131 genes were identified in each genome including 85 protein-coding genes, 37 tRNA genes, eight rRNA genes and one pseudogene (ycf1). The simple-sequence repeat sequences mainly comprise A/T mononucletide repeats. Phylogenetic analyses using maximum parsimony (MP) revealed that A. debilis and A. contorta had a close phylogenetic relationship with species of the family Piperaceae, as well as Laurales and Magnoliales. The data obtained in this study will be beneficial for further investigations on A. debilis and A. contorta from the aspect of evolution, and chloroplast genetic engineering.

Degeneration of photosynthetic capacity in mixotrophic plants, Chimaphila japonica and Pyrola decorata (Ericaceae)

The evolution of photosynthesis is an important feature of mixotrophic plants. Previous inferences proposed that mixotrophic taxa tend to retain most genes relating to photosynthetic functions but vary in plastid gene content. However, no sequence data are available to test this hypothesis in Ericaceae. To investigate changes in plastid genomes that may result from a transition from autotrophy to mixotrophy, the plastomes of two mixotrophic plants, Pyrola decorata and Chimaphila japonica, were sequenced at Illumina's Genome Analyzer and compared to the published plastome of the autotrophic plant Rhododendron simsii, which also belongs to Ericaceae. The greatest discrepancy between mixotrophic and autotrophic plants was that ndh genes for both P. decorata and C. japonica plastomes have nearly all become pseudogenes. P. decorata and C. japonica also retained all genes directly involved in photosynthesis under strong selection. The calculated rate of nonsynonymous nucleotide substitutions and synonymous substitutions of protein-coding genes (dN/dS) showed that substitution rates in shade plants were apparently higher than those in sunlight plants. The two mixotrophic plastomes were generally very similar to that of non-parasitic plants, although ndh genes were largely pseudogenized. Photosynthesis genes under strong selection were retained in the two mixotrophs, however, with greatly increased substitution rates. Further research is needed to gain a clearer understanding of the evolution of autotrophy and mixotrophy in Ericaceae.

The complete chloroplast genome sequences of Lychnis wilfordii and Silene capitata and comparative analyses with other Caryophyllaceae genomes

The complete chloroplast genomes of Lychnis wilfordii and Silene capitata were determined and compared with ten previously reported Caryophyllaceae chloroplast genomes. The chloroplast genome sequences of L. wilfordii and S. capitata contain 152,320 bp and 150,224 bp, respectively. The gene contents and orders among 12 Caryophyllaceae species are consistent, but several microstructural changes have occurred. Expansion of the inverted repeat (IR) regions at the large single copy (LSC)/IRb and small single copy (SSC)/IR boundaries led to partial or entire gene duplications. Additionally, rearrangements of the LSC region were caused by gene inversions and/or transpositions. The 18 kb inversions, which occurred three times in different lineages of tribe Sileneae, were thought to be facilitated by the intermolecular duplicated sequences. Sequence analyses of the L. wilfordii and S. capitata genomes revealed 39 and 43 repeats, respectively, including forward, palindromic, and reverse repeats. In addition, a total of 67 and 56 simple sequence repeats were discovered in the L. wilfordii and S. capitata chloroplast genomes, respectively. Finally, we constructed phylogenetic trees of the 12 Caryophyllaceae species and two Amaranthaceae species based on 73 protein-coding genes using both maximum parsimony and likelihood methods.

The Chloroplast Genome of Passiflora edulis (Passifloraceae) Assembled from Long Sequence Reads: Structural Organization and Phylogenomic Studies in Malpighiales

The family Passifloraceae consists of some 700 species classified in around 16 genera. Almost all its members belong to the genus Passiflora. In Brazil, the yellow passion fruit (Passiflora edulis) is of considerable economic importance, both for juice production and consumption as fresh fruit. The availability of chloroplast genomes (cp genomes) and their sequence comparisons has led to a better understanding of the evolutionary relationships within plant taxa. In this study, we obtained the complete nucleotide sequence of the P. edulis chloroplast genome, the first entirely sequenced in the Passifloraceae family. We determined its structure and organization, and also performed phylogenomic studies on the order Malpighiales and the Fabids clade. The P. edulis chloroplast genome is characterized by the presence of two copies of an inverted repeat sequence (IRA and IRB) of 26,154 bp, each separating a small single copy region of 13,378 bp and a large single copy (LSC) region of 85,720 bp. The annotation resulted in the identification of 105 unique genes, including 30 tRNAs, 4 rRNAs, and 71 protein coding genes. Also, 36 repetitive elements and 85 SSRs (microsatellites) were identified. The structure of the complete cp genome of P. edulis differs from that of other species because of rearrangement events detected by means of a comparison based on 22 members of the Malpighiales. The rearrangements were three inversions of 46,151, 3,765 and 1,631 bp, located in the LSC region. Phylogenomic analysis resulted in strongly supported trees, but this could also be a consequence of the limited taxonomic sampling used. Our results have provided a better understanding of the evolutionary relationships in the Malpighiales and the Fabids, confirming the potential of complete chloroplast genome sequences in inferring evolutionary relationships and the utility of long sequence reads for generating very accurate biological information.

Limited mitogenomic degradation in response to a parasitic lifestyle in Orobanchaceae

In parasitic plants, the reduction in plastid genome (plastome) size and content is driven predominantly by the loss of photosynthetic genes. The first completed mitochondrial genomes (mitogenomes) from parasitic mistletoes also exhibit significant degradation, but the generality of this observation for other parasitic plants is unclear. We sequenced the complete mitogenome and plastome of the hemiparasite Castilleja paramensis (Orobanchaceae) and compared them with additional holoparasitic, hemiparasitic and nonparasitic species from Orobanchaceae. Comparative mitogenomic analysis revealed minimal gene loss among the seven Orobanchaceae species, indicating the retention of typical mitochondrial function among Orobanchaceae species. Phylogenetic analysis demonstrated that the mobile cox1 intron was acquired vertically from a nonparasitic ancestor, arguing against a role for Orobanchaceae parasites in the horizontal acquisition or distribution of this intron. The C. paramensis plastome has retained nearly all genes except for the recent pseudogenization of four subunits of the NAD(P)H dehydrogenase complex, indicating a very early stage of plastome degradation. These results lend support to the notion that loss of ndh gene function is the first step of plastome degradation in the transition to a parasitic lifestyle.

Chloroplast sequence of treegourd (Crescentia cujete, Bignoniaceae) to study phylogeography and domestication

PREMISE OF THE STUDY

Crescentia cujete (Bignoniaceae) fruit rinds are traditionally used for storage vessels and handicrafts. We assembled its chloroplast genome and identified single-nucleotide polymorphisms (SNPs).

METHODS AND RESULTS

Using a genome skimming approach, the whole chloroplast of C. cujete was assembled using 3,106,928 sequence reads of 150 bp. The chloroplast is 154,662 bp in length, structurally divided into a large single copy region (84,788 bp), a small single copy region (18,299 bp), and two inverted repeat regions (51,575 bp) with 88 genes annotated. By resequencing the whole chloroplast, we identified 66 SNPs in C. cujete (N = 30) and 68 SNPs in C. amazonica (N = 6). Nucleotide diversity was estimated at 1.1 × 10^-3 and 3.5 × 10^-3 for C. cujete and C. amazonica, respectively.

CONCLUSIONS

This broadened C. cujete genetic toolkit will be important to study the origin, domestication, diversity, and phylogeography of treegourds in the Neotropics.

Understanding the evolution of holoparasitic plants: the complete plastid genome of the holoparasite Cytinus hypocistis (Cytinaceae)

Background and Aims Plant plastid genomes are highly conserved in size, gene content and structure; however, parasitic plants are a noticeable exception to this evolutionary stability. Although the evolution of parasites could help to better understand plastome evolution in general, complete plastomes of parasites have been sequenced only for some lineages so far. Here we contribute to filling this gap by providing and analysing the complete plastome sequence of Cytinus hypocistis, the first parasite sequenced for Malvales and a species suspected to have an extremely small genome. Methods We sequenced and assembled de novo the plastid genome of Cytinus hypocistis using a shotgun approach on genomic DNA. Phylogenomic analyses based on coding regions were performed on Malvidae. For each coding region present in Cytinus, we tested for relaxation or intensification of selective pressures in the Cytinus lineage compared with autotrophic Malvales. Key Results Cytinus hypocistis has an extremely divergent genome that is among the smallest sequenced to date (19·4 kb), with only 23 genes and no inverted repeat regions. Phylogenomic analysis confirmed the position of Cytinus within Malvales. All coding regions of Cytinus plastome presented very high substitution rates compared with non-parasitic Malvales. Conclusions Some regions were inferred to be under relaxed negative selection in Cytinus, suggesting that further plastome reduction is occurring due to relaxed purifying selection associated with the loss of photosynthetic activity. On the other hand, increased selection intensity and strong positive selection were detected for rpl22 in the Cytinus lineage, which might indicate an evolutionary role in the host-parasite arms race, a point that needs further research.

Complete plastid genome sequence of Primula sinensis (Primulaceae): structure comparison, sequence variation and evidence for accD transfer to nucleus

Species-rich genus Primula L. is a typical plant group with which to understand genetic variance between species in different levels of relationships. Chloroplast genome sequences are used to be the information resource for quantifying this difference and reconstructing evolutionary history. In this study, we reported the complete chloroplast genome sequence of Primula sinensis and compared it with other related species. This genome of chloroplast showed a typical circular quadripartite structure with 150,859 bp in sequence length consisting of 37.2% GC base. Two inverted repeated regions (25,535 bp) were separated by a large single-copy region (82,064 bp) and a small single-copy region (17,725 bp). The genome consists of 112 genes, including 78 protein-coding genes, 30 tRNA genes and four rRNA genes. Among them, seven coding genes, seven tRNA genes and four rRNA genes have two copies due to their locations in the IR regions. The accD and infA genes lacking intact open reading frames (ORF) were identified as pseudogenes. SSR and sequence variation analyses were also performed on the plastome of Primula sinensis, comparing with another available plastome of P. poissonii. The four most variable regions, rpl36–rps8, rps16–trnQ, trnH–psbA and ndhC–trnV, were identified. Phylogenetic relationship estimates using three sub-datasets extracted from a matrix of 57 protein-coding gene sequences showed the identical result that was consistent with previous studies. A transcript found from P. sinensis transcriptome showed a high similarity to plastid accD functional region and was identified as a putative plastid transit peptide at the N-terminal region. The result strongly suggested that plastid accD has been functionally transferred to the nucleus in P. sinensis.

Drastic reduction of plastome size in the mycoheterotrophic Thismia tentaculata relative to that of its autotrophic relative Tacca chantrieri

PREMISE OF THE STUDY

Heterotrophic angiosperms tend to have reduced plastome sizes relative to those of their autotrophic relatives because genes that code for proteins involved in photosynthesis are lost. However, some plastid-encoded proteins may have vital nonphotosynthetic functions, and the plastome therefore may be retained after the loss of photosynthesis.

METHODS

We sequenced the plastome of the mycoheterotrophic species Thismia tentaculata and a representative of its sister genus, Tacca chantrieri, using next-generation technology, and we compared sequences and structures of genes and genomes of these species.

KEY RESULTS

The plastome of Tacca chantrieri is similar to those of other autotrophic taxa of Dioscoreaceae, except in a few local rearrangements and one gene loss. The plastome of Thismia tentaculata is ca. 16 kbp long with a quadripartite structure and is among the smallest known plastomes. Synteny is minimal between the plastomes of Tacca chantrieri and Thismia tentaculata. The latter includes only 12 candidate genes, with all except accD involved in protein synthesis. Of the 12 genes, trnE, trnfM, and accD are frequently among the few that remain in depauperate plastomes.

CONCLUSIONS

The plastome of Thismia tentaculata, like those of most other heterotrophic plants, includes a small number of genes previously suggested to be essential to plastome survival.

Massive intracellular gene transfer during plastid genome reduction in nongreen Orobanchaceae

Plastid genomes (plastomes) of nonphotosynthetic plants experience extensive gene losses and an acceleration of molecular evolutionary rates. Here, we inferred the mechanisms and timing of reductive genome evolution under relaxed selection in the broomrape family (Orobanchaceae). We analyzed the plastomes of several parasites with a major focus on the genus Orobanche using genome-descriptive and Bayesian phylogenetic-comparative methods. Besides this, we scanned the parasites' other cellular genomes to trace the fate of all genes that were purged from their plastomes. Our analyses indicate that the first functional gene losses occurred within 10 Myr of the transition to obligate parasitism in Orobanchaceae, and that the physical plastome reduction proceeds by small deletions that accumulate over time. Evolutionary rate shifts coincide with the genomic reduction process in broomrapes, suggesting that the shift of selectional constraints away from photosynthesis to other molecular processes alters the plastid rate equilibrium. Most of the photosynthesis-related genes or fragments of genes lost from the plastomes of broomrapes have survived in their nuclear or mitochondrial genomes as the results of multiple intracellular transfers and subsequent fragmentation. Our findings indicate that nonessential DNA is eliminated much faster in the plastomes of nonphotosynthetic parasites than in their other cellular genomes.

Complete Plastid Genome of the Recent Holoparasite Lathraea squamaria Reveals Earliest Stages of Plastome Reduction in Orobanchaceae

Plants from the family Orobanchaceae are widely used as a model to study different aspects of parasitic lifestyle including host–parasite interactions and physiological and genomic adaptations. Among the latter, the most prominent are those that occurred due to the loss of photosynthesis; they include the reduction of the photosynthesis-related gene set in both nuclear and plastid genomes. In Orobanchaceae, the transition to non-photosynthetic lifestyle occurred several times independently, but only one lineage has been in the focus of evolutionary studies. These studies included analysis of plastid genomes and transcriptomes and allowed the inference of patterns and mechanisms of genome reduction that are thought to be general for parasitic plants. Here we report the plastid genome of Lathraea squamaria, a holoparasitic plant from Orobanchaceae, clade Rhinantheae. We found that in this plant the degree of plastome reduction is the least among non-photosynthetic plants. Like other parasites, Lathraea possess a plastome with elevated absolute rate of nucleotide substitution. The only gene lost is petL, all other genes typical for the plastid genome are present, but some of them–those encoding photosystem components (22 genes), cytochrome b₆/f complex proteins (4 genes), plastid-encoded RNA polymerase subunits (2 genes), ribosomal proteins (2 genes), ccsA and cemA–are pseudogenized. Genes for cytochrome b₆/f complex and photosystems I and II that do not carry nonsense or frameshift mutations have an increased ratio of non-synonymous to synonymous substitution rates, indicating the relaxation of purifying selection. Our divergence time estimates showed that transition to holoparasitism in Lathraea lineage occurred relatively recently, whereas the holoparasitic lineage Orobancheae is about two times older.