GreeNC 2.0: a comprehensive database of plant long non-coding RNAs.

The Green Non-Coding Database (GreeNC) is one of the reference databases for the study of plant long non-coding RNAs (lncRNAs). Here we present our most recent update where 16 species have been updated, while 78 species have been added, resulting in the annotation of more than 495 000 lncRNAs. Moreover, sequence clustering was applied providing information about sequence conservation and gene families. The current version of the database is available at: http://greenc.sequentiabiotech.com/wiki2/Main_Page.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development.

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Azadirachta indica MicroRNAs: Genome-Wide Identification, Target Transcript Prediction, and Expression Analyses.

MicroRNAs are short, endogenous, non-coding RNAs, liable for essential regulatory function. Numerous miRNAs have been identified and studied in plants with known genomic or small RNA resources. Despite the availability of genomic and transcriptomic resources, the miRNAs have not been reported in the medicinal tree Azadirachta indica (Neem) till date. Here for the first time, we report extensive identification of miRNAs and their possible targets in A. indica which might help to unravel their therapeutic potential. A comprehensive search of miRNAs in the A. indica genome by C-mii tool was performed. Overall, 123 miRNAs classified into 63 families and their stem-loop hairpin structures were predicted. The size of the A. indica (ain)-miRNAs ranged between 19 and 23 nt in length, and their corresponding ain-miRNA precursor sequence MFEI value averaged as -1.147 kcal/mol. The targets of ain-miRNAs were predicted in A. indica as well as Arabidopsis thaliana plant. The gene ontology (GO) annotation revealed the involvement of ain-miRNA targets in developmental processes, transport, stress, and metabolic processes including secondary metabolism. Stem-loop qRT-PCR was carried out for 25 randomly selected ain-miRNAs and differential expression patterns were observed in different A. indica tissues. Expression of miRNAs and its targets shows negative correlation in a dependent manner.

A potential endogenous gibberellin-mediated signaling cascade regulated floral transition in Magnolia × soulangeana 'Changchun'.

The floral transition is a critical developmental switch in plants, and has profound effects on the flower production and yield. Magnolia × soulangeana 'Changchun' is known as a woody ornamental plant, which can bloom in spring and summer, respectively. In this study, anatomical observation, physiological measurement, transcriptome, and small RNA sequencing were performed to investigate potential endogenous regulatory mechanisms underlying floral transition in 'Changchun'. Transition of the shoot apical meristem from vegetative to reproductive growth occurred between late April and early May. During this specific developmental process, a total of 161,645 unigenes were identified, of which 73,257 were significantly differentially expressed, while a number of these two categories of miRNAs were 299 and 148, respectively. Further analysis of differentially expressed genes (DEGs) revealed that gibberellin signaling could regulate floral transition in 'Changchun' in a DELLA-dependent manner. In addition, prediction and analysis of miRNA targeted genes suggested that another potential molecular regulatory module was mediated by the miR172 family and other several novel miRNAs (Ms-novel_miR139, Ms-novel_miR229, and Ms-novel_miR232), with the participation of up- or down-regulating genes, including MsSVP, MsAP2, MsTOE3, MsAP1, MsGATA6, MsE2FA, and MsMDS6. Through the integrated analysis of mRNA and miRNA, our research results will facilitate the understanding of the potential molecular mechanism underlying floral transition in 'Changchun', and also provide basic experimental data for the plant germplasm resources innovation in Magnolia.

PLncDB V2.0: a comprehensive encyclopedia of plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with little or no protein coding potential. The expanding list of lncRNAs and accumulating evidence of their functions in plants have necessitated the creation of a comprehensive database for lncRNA research. However, currently available plant lncRNA databases have some deficiencies, including the lack of lncRNA data from some model plants, uneven annotation standards, a lack of visualization for expression patterns, and the absence of epigenetic information. To overcome these problems, we upgraded our Plant Long noncoding RNA Database (PLncDB, http://plncdb.tobaccodb.org/), which was based on a uniform annotation pipeline. PLncDB V2.0 currently contains 1 246 372 lncRNAs for 80 plant species based on 13 834 RNA-Seq datasets, integrating lncRNA information from four other resources including EVLncRNAs, RNAcentral and etc. Expression patterns and epigenetic signals can be visualized using multiple tools (JBrowse, eFP Browser and EPexplorer). Targets and regulatory networks for lncRNAs are also provided for function exploration. In addition, PLncDB V2.0 is hierarchical and user-friendly and has five built-in search engines. We believe PLncDB V2.0 is useful for the plant lncRNA community and data mining studies and provides a comprehensive resource for data-driven lncRNA research in plants.

Identification of miRNAs linked to peanut nodule functional processes.

microRNAs (miRNAs) are non-coding small RNAs that regulate gene expression at post-transcriptional level. Thousands of miRNAs have been identified in legumes, but studies about miRNAs linked to peanut nodule functionality are scarce. In this work we analyzed transcriptional changes in peanut nodules to identify miRNAs involved in functional processes of these organs. We found 32 miRNAs precursors differentially expressed in nodules compared with roots, and predicted the potential targets of their corresponding mature miRNAs. Among them, 20 belong to 14 conserved miRNAs families and 12 are Arachis hypogaea-specific miRNAs. Expression levels of 3 miRNAs (ahy-miR399, ahy-miR159 and ahy-miR3508) were confirmed experimentally by qPCR. We also demonstrated that the expression of these miRNAs was not affected by inoculation of a biocontrol bacterium or a fungal pathogen. The catalogue of differentially expressed miRNA precursors and the expression of the corresponding mature miRNA potential targets in the nodules of A. hypogaea obtained in this work is a database of strong candidates, including A. hypogaea-specific miRNAs, for the regulation of the nodule functionality. The analysis of their role in this process will certainly lead to the characterization of essential regulators in these particular aeschynomenoid nodules.

Grape-RNA: A Database for the Collection, Evaluation, Treatment, and Data Sharing of Grape RNA-Seq Datasets.

Since its inception, RNA sequencing (RNA-seq) has become the most effective way to study gene expression. After more than a decade of development, numerous RNA-seq datasets have been created, and the full utilization of these datasets has emerged as a major issue. In this study, we built a comprehensive database named Grape-RNA, which is focused on the collection, evaluation, treatment, and data sharing of grape RNA-seq datasets. This database contains 1529 RNA-seq samples, 112 microRNA samples from the public platform, and 485 RNA-seq in-house datasets sequenced by our lab. We classified these data into 25 conditions and provide the sample information, cleaned raw data, expression level, assembled unigenes, useful tools, and other relevant information to the users. Thus, this study provides data and tools that should be beneficial for researchers by allowing them to easily use the RNA-seq. The provided information can greatly contribute to grape breeding and genomic and biological research. This study may improve the usage of RNA-seq.

Phylogenomics, co-evolution of ecological niche and morphology, and historical biogeography of buckeyes, horsechestnuts, and their relatives (Hippocastaneae, Sapindaceae) and the value of RAD-Seq for deep evolutionary inferences back to the Late Cretaceous.

In this study, we used RAD-seq data to resolve the phylogeny of the tribe Hippocastaneae (Sapindaceae) and conducted comparative analyses to gain insights into the evolution and biogeography of the group that had fossils dating back to the late Cretaceous. Hippocastaneae, including the horsechestnuts and buckeyes, is a well-supported clade in Sapindaceae that comprises 12-14 species in Aesculus, two in Billia, and one in Handeliodendron. Most species in the tribe are distributed in Eurasia and North America and exhibit a classic pattern of intercontinental disjunction in the Northern Hemisphere, while Billia occurs from southern Mexico to northern South America. The earliest fossils of Aesculus date back to at least the earliest Paleocene of eastern Asia and western North America, where there are also putative occurrences from the latest Cretaceous. The group provides an excellent system for understanding floristic disjunction in the Northern Hemisphere extending to the Neotropics. However, a strongly supported and well resolved phylogeny is presently lacking for the tribe. Previous phylogenetic studies using several gene regions revealed five well-supported clades in Aesculus, largely corresponding to five recognized taxonomic sections, but relationships among these clades and among Aesculus, Billia, and Handeliodendron were not well supported. In this study, we used RAD-seq data from 68 samples representing all clades and species of Hippocastaneae except Billia, for which we used one of two species, to further resolve relationships within the tribe. Our phylogenomic analyses showed strong support for a sister relationship between Aesculus and Handeliodendron, in contrast to previous findings which supported Billia as sister to Aesculus. Within Aesculus, relationships among sections were strongly supported as (sect. Calothyrsus, (sect. Aesculus, (sect. Macrothyrsus, (sect. Parryana, sect. Pavia)))). We found that the traditionally recognized section Calothyrsus was monophyletic, with all eastern Asian species sister to the western North American species, A. californica. Analyses of divergence times combined with biogeographic analyses suggested a Late Cretaceous origin of Hippocastaneae, in eastern Asia, western North America, and Central America (including southern Mexico), followed by isolation of Billia in Central America, extinction of the tribe ancestor in western North America, and divergence of Aesculus from Handeliodendron in eastern Asia. A Late Cretaceous origin of the common ancestor of Aesculus in eastern Asia was followed by dispersals into western North America, Europe, and eastern North America during the Late Cretaceous and the Paleogene. Our results support Aesculus as a relic of the boreotropical flora and subsequent intercontinental spread of the genus through the Bering and North Atlantic land bridges. We performed character mapping analyses, which revealed that biogeographic isolation and niche divergence may have played important roles in driving morphological evolution and lineage divergence in Aesculus. Our study demonstrates the value of RAD-seq data for reconstructing phylogeny back to the Late Cretaceous.

Maize transposable elements contribute to long non-coding RNAs that are regulatory hubs for abiotic stress response.

BACKGROUND: Several studies have mined short-read RNA sequencing datasets to identify long non-coding RNAs (lncRNAs), and others have focused on the function of individual lncRNAs in abiotic stress response. However, our understanding of the complement, function and origin of lncRNAs - and especially transposon derived lncRNAs (TE-lncRNAs) - in response to abiotic stress is still in its infancy. RESULTS: We utilized a dataset of 127 RNA sequencing samples that included total RNA datasets and PacBio fl-cDNA data to discover lncRNAs in maize. Overall, we identified 23,309 candidate lncRNAs from polyA+ and total RNA samples, with a strong discovery bias within total RNA. The majority (65%) of the 23,309 lncRNAs had sequence similarity to transposable elements (TEs). Most had similarity to long-terminal-repeat retrotransposons from the Copia and Gypsy superfamilies, reflecting a high proportion of these elements in the genome. However, DNA transposons were enriched for lncRNAs relative to their genomic representation by ~ 2-fold. By assessing the fraction of lncRNAs that respond to abiotic stresses like heat, cold, salt and drought, we identified 1077 differentially expressed lncRNA transcripts, including 509 TE-lncRNAs. In general, the expression of these lncRNAs was significantly correlated with their nearest gene. By inferring co-expression networks across our large dataset, we found that 39 lncRNAs are as major hubs in co-expression networks that respond to abiotic stress, and 18 appear to be derived from TEs. CONCLUSIONS: Our results show that lncRNAs are enriched in total RNA samples, that most (65%) are derived from TEs, that at least 1077 are differentially expressed during abiotic stress, and that 39 are hubs in co-expression networks, including a small number that are evolutionary conserved. These results suggest that lncRNAs, including TE-lncRNAs, may play key regulatory roles in moderating abiotic responses.

Integrated transcriptome and miRNA analysis uncovers molecular regulators of aerial stem-to-rhizome transition in the medical herb Gynostemma pentaphyllum.

BACKGROUND: Gynostemma pentaphyllum is an important perennial medicinal herb belonging to the family Cucurbitaceae. Aerial stem-to-rhizome transition before entering the winter is an adaptive regenerative strategy in G. pentaphyllum that enables it to survive during winter. However, the molecular regulation of aerial stem-to-rhizome transition is unknown in plants. Here, integrated transcriptome and miRNA analysis was conducted to investigate the regulatory network of stem-to-rhizome transition. RESULTS: Nine transcriptome libraries prepared from stem/rhizome samples collected at three stages of developmental stem-to-rhizome transition were sequenced and a total of 5428 differentially expressed genes (DEGs) were identified. DEGs associated with gravitropism, cell wall biosynthesis, photoperiod, hormone signaling, and carbohydrate metabolism were found to regulate stem-to-rhizome transition. Nine small RNA libraries were parallelly sequenced, and seven significantly differentially expressed miRNAs (DEMs) were identified, including four known and three novel miRNAs. The seven DEMs targeted 123 mRNAs, and six pairs of miRNA-target showed significantly opposite expression trends. The GpmiR166b-GpECH2 module involved in stem-to-rhizome transition probably promotes cell expansion by IBA-to-IAA conversion, and the GpmiR166e-GpSGT-like module probably protects IAA from degradation, thereby promoting rhizome formation. GpmiR156a was found to be involved in stem-to-rhizome transition by inhibiting the expression of GpSPL13A/GpSPL6, which are believed to negatively regulate vegetative phase transition. GpmiR156a and a novel miRNA Co.47071 co-repressed the expression of growth inhibitor GpRAV-like during stem-to-rhizome transition. These miRNAs and their targets were first reported to be involved in the formation of rhizomes. In this study, the expression patterns of DEGs, DEMs and their targets were further validated by quantitative real-time PCR, supporting the reliability of sequencing data. CONCLUSIONS: Our study revealed a comprehensive molecular network regulating the transition of aerial stem to rhizome in G. pentaphyllum. These results broaden our understanding of developmental phase transitions in plants.

Computational tools for plant small RNA detection and categorization.

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: insights on nitrogen starvation responses.

Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.

Grass phasiRNAs and male fertility.

Recent studies have indicated that a special type of small noncoding RNAs, phased small-interfering RNAs (phasiRNAs) play crucial roles in many cellular processes of plant development. PhasiRNAs are generated from long RNA precursors at intervals of 21 or 24 nt in plants, and they are produced from both protein-coding gene and long noncoding RNA genes. Different from those in eudicots, grass phasiRNAs include a special class of small RNAs that are specifically expressed in reproductive organs. These grass phasiRNAs are associated with gametogenesis, especially with anther development and male fertility. In this review, we summarized current knowledge on these small noncoding RNAs in male germ cells and their possible biological functions and mechanisms in grass species.

Survey of 800+ data sets from human tissue and body fluid reveals xenomiRs are likely artifacts.

miRNAs are small 22-nucleotide RNAs that can post-transcriptionally regulate gene expression. It has been proposed that dietary plant miRNAs can enter the human bloodstream and regulate host transcripts; however, these findings have been widely disputed. We here conduct the first comprehensive meta-study in the field, surveying the presence and abundances of cross-species miRNAs (xenomiRs) in 824 sequencing data sets from various human tissues and body fluids. We find that xenomiRs are commonly present in tissues (17%) and body fluids (69%); however, the abundances are low, comprising 0.001% of host human miRNA counts. Further, we do not detect a significant enrichment of xenomiRs in sequencing data originating from tissues and body fluids that are exposed to dietary intake (such as liver). Likewise, there is no significant depletion of xenomiRs in tissues and body fluids that are relatively separated from the main bloodstream (such as brain and cerebro-spinal fluids). Interestingly, the majority (81%) of body fluid xenomiRs stem from rodents, which are a rare human dietary contribution but common laboratory animals. Body fluid samples from the same studies tend to group together when clustered by xenomiR compositions, suggesting technical batch effects. Last, we performed carefully designed and controlled animal feeding studies, in which we detected no transfer of plant miRNAs into rat blood, or bovine milk sequences into piglet blood. In summary, our comprehensive computational and experimental results indicate that xenomiRs originate from technical artifacts rather than dietary intake.

Genome-wide identification of potato long intergenic noncoding RNAs responsive to Pectobacterium carotovorum subspecies brasiliense infection.

BACKGROUND: Long noncoding RNAs (lncRNAs) represent a class of RNA molecules that are implicated in regulation of gene expression in both mammals and plants. While much progress has been made in determining the biological functions of lncRNAs in mammals, the functional roles of lncRNAs in plants are still poorly understood. Specifically, the roles of long intergenic nocoding RNAs (lincRNAs) in plant defence responses are yet to be fully explored. RESULTS: In this study, we used strand-specific RNA sequencing to identify 1113 lincRNAs in potato (Solanum tuberosum) from stem tissues. The lincRNAs are expressed from all 12 potato chromosomes and generally smaller in size compared to protein-coding genes. Like in other plants, most potato lincRNAs possess single exons. A time-course RNA-seq analysis between a tolerant and a susceptible potato cultivar showed that 559 lincRNAs are responsive to Pectobacterium carotovorum subsp. brasiliense challenge compared to mock-inoculated controls. Moreover, coexpression analysis revealed that 17 of these lincRNAs are highly associated with 12 potato defence-related genes. CONCLUSIONS: Together, these results suggest that lincRNAs have potential functional roles in potato defence responses. Furthermore, this work provides the first library of potato lincRNAs and a set of novel lincRNAs implicated in potato defences against P. carotovorum subsp. brasiliense, a member of the soft rot Enterobacteriaceae phytopathogens.

A naïve Bayesian classifier for identifying plant microRNAs.

MicroRNAs (miRNAs) are important regulatory molecules in eukaryotic organisms. Existing methods for the identification of mature miRNA sequences in plants rely extensively on the search for stem-loop structures, leading to high false negative rates. Here, we describe a probabilistic method for ranking putative plant miRNAs using a naïve Bayes classifier and its publicly available implementation. We use a number of properties to construct the classifier, including sequence length, number of observations, existence of detectable predicted miRNA* sequences, the distribution of nearby reads and mapping multiplicity. We apply the method to small RNA sequence data from soybean, peach, Arabidopsis and rice and provide experimental validation of several predictions in soybean. The approach performs well overall and strongly enriches for known miRNAs over other types of sequences. By utilizing a Bayesian approach to rank putative miRNAs, our method is able to score miRNAs that would be eliminated by other methods, such as those that have low counts or lack detectable miRNA* sequences. As a result, we are able to detect several soybean miRNA candidates, including some that are 24 nucleotides long, a class that is almost universally eliminated by other methods.

Long noncoding RNA transcriptome of plants.

Since their discovery more than two decades ago, animal long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes. Recently, a large number of lncRNAs have also been identified in higher plants, and here, we review their identification, classification and known regulatory functions in various developmental events and stress responses. Knowledge gained from a deeper understanding of this special group of noncoding RNAs may lead to biotechnological improvement of crops. Some possible examples in this direction are discussed.

Genome-wide discovery and analysis of phased small interfering RNAs in Chinese sacred lotus.

Phased small interfering RNA (phasiRNA) generating loci (briefly as PHAS) in plants are a novel class of genes that are normally regulated by microRNAs (miRNAs). Similar to miRNAs, phasiRNAs encoded by PHAS play important regulatory roles by targeting protein coding transcripts in plant species. We performed a genome-wide discovery of PHAS loci in Chinese sacred lotus and identified a total of 106 PHAS loci. Of these, 47 loci generate 21 nucleotide (nt) phasiRNAs and 59 loci generate 24 nt phasiRNAs, respectively. We have also identified a new putative TAS3 and a putative TAS4 loci in the lotus genome. Our results show that some of the nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance proteins and MYB transcription factors potentially generate phasiRNAs. Furthermore, our results suggest that some large subunit (LSU) rRNAs can derive putative phasiRNAs, which is potentially resulted from crosstalk between small RNA biogenesis pathways that are employed to process rRNAs and PHAS loci, respectively. Some of the identified phasiRNAs have putative trans-targets with less than 4 mismatches, suggesting that the identified PHAS are involved in many different pathways. Finally, the discovery of 24 nt PHAS in lotus suggests that there are 24 nt PHAS in dicots.

A genome-wide identification and characterization of mircoRNAs and their targets in 'Suli' pear (Pyrus pyrifolia white pear group).

MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that are endogenous regulators of gene expression. miRNAs play a crucial role in cells via degradation of target mRNAs or by inhibition of target protein translation. In the present study, 186 new potentially conserved pear miRNAs belonging to 37 families were identified. The length of mature miRNAs ranged from 19 to 24 nt, and most of the miRNAs (154 out of 186) were 21 nt in length. The length of pre-miRNAs in pear was also found to vary from 62 to 282 nt with an average of 105 ± 43 nt. The potential miRNAs belonged to 29 clusters involving 20 different miRNA families. Using these potential miRNAs, we further scoured of the pear genome and found 326 potential target genes, which included transcription factors, stress responsive genes, and the genes involved in transmembrane transport and signal transduction. Gene ontology analysis of these potential targets suggested that 47 biological processes were potentially regulated by miRNAs, including oxidation-reduction, stress response, transport, etc. KEGG pathway analysis showed that the identified miRNAs were found in 15 metabolism networks which were related to starch and sucrose metabolism, and ascorbate and aldarate metabolism, among others. Our study will help in the further understanding of the essential role of miRNAs in growth and development and stress response of pear.

Comparative genomics of four Liliales families inferred from the complete chloroplast genome sequence of Veratrum patulum O. Loes. (Melanthiaceae).

The sequence of the chloroplast genome, which is inherited maternally, contains useful information for many scientific fields such as plant systematics, biogeography and biotechnology because its characteristics are highly conserved among species. There is an increase in chloroplast genomes of angiosperms that have been sequenced in recent years. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Veratrum patulum Loes. (Melanthiaceae, Liliales) was analyzed completely. The circular double-stranded DNA of 153,699 bp consists of two inverted repeat (IR) regions of 26,360 bp each, a large single copy of 83,372 bp, and a small single copy of 17,607 bp. This plastome contains 81 protein-coding genes, 30 distinct tRNA and four genes of rRNA. In addition, there are six hypothetical coding regions (ycf1, ycf2, ycf3, ycf4, ycf15 and ycf68) and two open reading frames (ORF42 and ORF56), which are also found in the chloroplast genomes of the other species. The gene orders and gene contents of the V. patulum plastid genome are similar to that of Smilax china, Lilium longiflorum and Alstroemeria aurea, members of the Smilacaceae, Liliaceae and Alstroemeriaceae (Liliales), respectively. However, the loss rps16 exon 2 in V. patulum results in the difference in the large single copy regions in comparison with other species. The base substitution rate is quite similar among genes of these species. Additionally, the base substitution rate of inverted repeat region was smaller than that of single copy regions in all observed species of Liliales. The IR regions were expanded to trnH_GUG in V. patulum, a part of rps19 in L. longiflorum and A. aurea, and whole sequence of rps19 in S. china. Furthermore, the IGS lengths of rbcL-accD-psaI region were variable among Liliales species, suggesting that this region might be a hotspot of indel events and the informative site for phylogenetic studies in Liliales. In general, the whole chloroplast genome of V. patulum, a potential medicinal plant, will contribute to research on the genetic applications of this genus.