Chloroplast Genome Sequence of Clusterbean (Cyamopsis tetragonoloba L.): Genome Structure and Comparative Analysis

The chromosome-scale genome assembly of cluster bean provides molecular insight into edible gum (galactomannan) biosynthesis family genes

Cluster bean (Cyamopsis tetragonoloba (L.) Taub 2n = 14, is commonly known as Guar. Apart from being a vegetable crop, it is an abundant source of a natural hetero-polysaccharide called guar gum or galactomannan. Here, we are reporting a chromosome-scale reference genome assembly of a popular cluster bean cultivar RGC-936, by combining sequencing data from Illumina, 10X Genomics, Oxford Nanopore technologies. An initial assembly of 1580 scaffolds with an N50 value of 7.12 Mb was generated and these scaffolds were anchored to a high density SNP linkage map. Finally, a genome assembly of 550.31 Mb (94% of the estimated genome size of ~ 580 Mb (through flow cytometry) with 58 scaffolds was obtained, including 7 super scaffolds with a very high N50 value of 78.27 Mb. Phylogenetic analysis using single copy orthologs among 12 angiosperms showed that cluster bean shared a common ancestor with other legumes 80.6 MYA. No evidence of recent whole genome duplication event in cluster bean was found in our analysis. Further comparative transcriptomics analyses revealed pod-specific up-regulation of genes encoding enzymes involved in galactomannan biosynthesis. The high-quality chromosome-scale cluster bean genome assembly will facilitate understanding of the molecular basis of galactomannan biosynthesis and aid in genomics-assisted improvement of cluster bean.

Comparative genomics analysis of endangered wild Egyptian Moringa peregrina (Forssk.) Fiori plastome, with implications for the evolution of Brassicales order

Moringa is a mono-genus belonging to the Moringaceae family, which includes 13 species. Among them, Moringa peregrina is plant species native to the Arabian Peninsula, Southern Sinai in Egypt, and the Horn of Africa, and comprehensive studies on its nutritional, industrial, and medicinal values have been performed. Herein, we sequenced and analyzed the initial complete chloroplast genome of Moringa peregrina. Concurrently, we analyzed the new chloroplast genome along with 25 chloroplast genomes related to species representing eight families in the Brassicales order. The results indicate that the plastome sequence of M. peregrina consists of 131 genes, with an average GC content of 39.23%. There is a disparity in the IR regions of the 26 species ranging from 25,804 to 31,477 bp. Plastome structural variations generated 20 hotspot regions that could be considered prospective DNA barcode locations in the Brassicales order. Tandem repeats and SSR structures are reported as significant evidence of structural variations among the 26 tested specimens. Furthermore, selective pressure analysis was performed to estimate the substitution rate within the Moringaceae family, which revealing that the ndhA and accD genes are under positive selective pressure. The phylogenetic analysis of the Brassicales order produced an accurate monophyletic annotation cluster of the Moringaceae and Capparaceae species, offering unambiguous identification without overlapping groups between M. oleifera and M. peregrina, which are genetically strongly associated. Divergence time estimation suggests that the two Moringa species recently diversified, 0.467 Ma. Our findings highlight the first complete plastome of the Egyptian wild-type of M. peregrina, which can be used for determining plastome phylogenetic relationships and systematic evolution history within studies on the Moringaceae family.

Guar gum-based nanoformulations: Implications for improving drug delivery

Poorly soluble drugs are reported to easily degrade in the gastrointestinal tract and contribute in limiting the effect of drug to its targeted site. Oral administration of drug is one of the prominent ways to deliver a drug, although, it experiences barriers like acidic pH, presence of microflora and enzymes in the gastrointestinal tract. Collectively all of these participate in the degradation of drug before it reaches its target site and thus, they impede the sustained effect of drug. A quest of choosing a polymer with good stability profile and releasing the drug to its targeted site is always been a challenge for the scientists worldwide. Many polymers have been reported to prevent the degradation of drug and one such naturally occurring biocompatible polymer is guar gum. Guar gum-based nanoformulations have been extensively used in past decades to achieve controlled drug release which defines its importance. The coating of guar gum over the drug improves the bioavailability of the drug and thus helps in minimizing the risk of drug degradation. This review intends to highlight the beneficial role of guar gum-based nanoformulations to improve drug delivery by ameliorating the bioavailibility.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization

BACKGROUND

Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences.

RESULTS

The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens.

CONCLUSION

G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

Plastome evolution and phylogeny of the tribe Ruteae (Rutaceae)

Rutaceae is a large family, and the genus-level classification in the subfamilies or tribes of this family is not unified based on different taxonomic treatments. Until now, phylogenetic relationships of some genera in traditional tribe Ruteae have not been clearly resolved. In this study, seven new complete plastomes of this tribe were sequenced, and a comparative analysis was performed to investigate their plastome characteristics and evolution. In addition, we inferred the phylogenetic relationships of Ruteae based on complete plastome and nuclear ITS data. All plastomes exhibited a typical quadripartite structure and were relatively conserved in their structure and gene arrangement. Their genome sizes ranged from 154,656 bp to 160,677 bp, and the size variation was found to be associated with differences in IR expansion and gene loss. A total of 112 to 114 genes were identified in the genomes, including 78 to 79 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 2 pseudogenes. Sequence divergence analysis indicated that non-coding regions exhibited a higher percentage of variable characters, and nine non-coding and six coding regions were identified as divergent hotspots. Phylogenetic results based on different datasets showed that this tribe was divided into three reciprocally exclusive groups. The phylogenetic analyses between plastome and nuclear ITS data were partly incongruent with each other. This study provides new insights into plastome evolution of Ruteae as well as Rutaceae. The availability of these plastomes provides useful genomic resources for molecular DNA barcodes and phylogenetically informative markers and deepens our understanding of the phylogeny in Ruteae.

The chloroplast genome of black pepper (Piper nigrum L.) and its comparative analysis with related Piper species

Piper nigrum, also known as black pepper, is an economically and ecologically important crop of the genus Piper. It has been titled as the king of spices due to its wide consumption throughout the world. In the present investigation, the chloroplast genome of P. nigrum has been assembled from a whole genome sequence by integrating the short and long reads generated through Illumina and PacBio platforms, respectively. The chloroplast genome was observed to be 161,522 bp in size, having a quadripartite structure with a large single copy (LSC) region of 89,153 bp and a small single copy (SSC) region of 18,255 bp separated by a copy of inverted repeats (IRs), each 27,057 bp in length. Taking into consideration all the duplicated genes, a total of 131 genes were observed, which included 81 protein-coding genes, 37 tRNAs, 4 rRNAs, and 1 pseudogene. Individually, the LSC region consisted of 83 genes, the SSC region had 13 genes, and 18 genes were present in each IR region. Additionally, 216 SSRs were detected and 11 of these were validated through amplification in 12 species of Piper. The features of the chloroplast genome have been compared with those of the genus Piper. Our results provide useful insights into evolutionary and molecular studies of black pepper which will contribute to its further genetic improvement and breeding.

In silico Identification of miRNAs and Their Targets in Cluster Bean for Their Role in Development and Physiological Responses

Cluster bean popularly known as “guar” is a drought-tolerant, annual legume that has recently emerged as an economically important crop, owing to its high protein and gum content. The guar gum has wide range of applications in food, pharma, and mining industries. India is the leading exporter of various cluster bean-based products all across the globe. Non-coding RNAs (miRNAs) are involved in regulating the expression of the target genes leading to variations in the associated pathways or final protein concentrations. The understanding of miRNAs and their associated targets in cluster bean is yet to be used to its full potential. In the present study, cluster bean EST (Expressed Sequence Tags) database was exploited to identify the miRNA and their predicted targets associated with metabolic and biological processes especially response to diverse biotic and abiotic stimuli using in silico approach. Computational analysis based on cluster bean ESTs led to the identification of 57 miRNAs along with their targets. To the best of our knowledge, this is the first report on identification of miRNAs and their targets using ESTs in cluster bean. The miRNA related to gum metabolism was also identified. Most abundant miRNA families predicted in our study were miR156, miR172, and miR2606. The length of most of the mature miRNAs was found to be 21nt long and the range of minimal folding energy (MFE) was 5.8–177.3 (−kcal/mol) with an average value of 25.4 (−kcal/mol). The identification of cluster bean miRNAs and their targets is predicted to hasten the miRNA discovery, resulting in better knowledge of the role of miRNAs in cluster bean development, physiology, and stress responses.

Comparative transcriptomic analysis and antioxidant defense mechanisms in clusterbean (Cyamopsis tetragonoloba (L.) Taub.) genotypes with contrasting drought tolerance

To understand drought tolerance mechanism(s) in clusterbean (Cyamopsis tetragonoloba), we conducted physiological, biochemical, and de novo comparative transcriptome analysis of drought-tolerant (RGC-1002) and drought-sensitive (RGC-1066) genotypes subjected to 30 days of drought stress. Relative water content (RWC) was maintained in tolerant genotype but was reduced in sensitive genotype. Leaf pigment concentrations were higher in tolerant genotype. Net photosynthesis was significantly decreased in sensitive genotype but insignificant reduction was found in tolerant genotype. Enzymatic antioxidant (GR, APX, DHAR) activities were enhanced in tolerant genotype, while there were insignificant changes in these enzymes in sensitive genotype. The ratios of antioxidant molecules (ASC/DHA and GSH/GSSG) were higher in tolerant genotype as compared to sensitive genotype. In sensitive genotype, 6625 differentially expressed genes (DEGs) were upregulated and 5365 genes were downregulated. In tolerant genotype, 5206 genes were upregulated and 2793 genes were downregulated. In tolerant genotype, transketolase family protein, phosphoenolpyruvate carboxylase 3, temperature-induced lipocalin, and cytochrome oxidase were highly upregulated. Moreover, according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the drought tolerance may be attributed to upregulated starch and sucrose metabolism-related genes in tolerant genotype. Finally, quantitative real-time PCR confirmed the reproducibility of the RNA-seq data.

The complete chloroplast genome sequence of Rubus hirsutus Thunb. and a comparative analysis within Rubus species

Rubus hirsutus is a type of tonifying kidney-essence herb that belongs to the Rosaceae family, and has been commonly used to treat multiple diseases, such as polyuria, impotence, and infertility. In this study, we determined the complete chloroplast sequence of R. hirsutus and conduced a comparative analysis within the genus Rubus. The assembled chloroplast (cp.) genome is 156,380 bp in length with a GC content of 37.0% and shares a conserved quadripartite structure within the other cp. genomes in this genus. A total of 132 unique genes were annotated in the cp. genome of R. hirsutus, which contained 87 protein-coding genes, 37 tRNAs, and eight rRNAs. Seventeen duplicated genes were identified in the inverted repeats region. Furthermore, 70 simple sequence repeats and 35 long repeats were detected in total in the R. hirsutus chloroplast genome. Eight mutational hotspots were identified in the cp. genome of this species with higher nucleotide variations in non-coding regions than those of coding regions. Furthermore, the gene order, codon usage, and repeat sequence distribution were highly consistent in Rubus according to the results of a comparative analysis. A phylogenetic analysis indicated that there was a sister relationship between R. hirsutus and R. chingii. Overall, the complete chloroplast genome of R. hirsutus and the comparative analysis will help to further the evolutionary study, conservation, phylogenetic reconstruction, and development of molecular barcodes for the genus Rubus.

The complete chloroplast genome sequence of a hybrid blackberry (Rubus spp.) cultivar

Blackberry (Rubus spp.) is an important hybrid fruit crop popular in the US Pacific Northwest and the European region with complex origins. In this study, we report the complete chloroplast genome sequence of a hybrid blackberry cultivar 'Arapohol' using next-generation sequencing technology. The complete chloroplast genome size is 156,621 bp. The genome contains 134 genes, including 40 tRNA genes, 86 protein-coding genes, and 8 rRNA genes. Phylogenetic analysis based on 11 complete chloroplast genomes revealed that taxa is closely related to Rubus niveus. The complete chloroplast genome of this Rubus sp. provides valuable information for understanding the origination of this crop species.

Comparative Analysis of Chloroplast Genomes of Four Medicinal Capparaceae Species: Genome Structures, Phylogenetic Relationships and Adaptive Evolution

This study presents for the first time the complete chloroplast genomes of four medicinal species in the Capparaceae family belonging to two different genera, Cadaba and Maerua (i.e., C. farinosa, C. glandulosa, M. crassifolia and M. oblongifolia), to investigate their evolutionary process and to infer their phylogenetic positions. The four species are considered important medicinal plants, and are used in the treatment of many diseases. In the genus Cadaba, the chloroplast genome ranges from 156,481 bp to 156,560 bp, while that of Maerua ranges from 155,685 bp to 155,436 bp. The chloroplast genome of C. farinosa, M. crassifolia and M. oblongifolia contains 138 genes, while that of C. glandulosa contains 137 genes, comprising 81 protein-coding genes, 31 tRNA genes and 4 rRNA genes. Out of the total genes, 116–117 are unique, while the remaining 19 are replicated in inverted repeat regions. The psbG gene, which encodes for subunit K of NADH dehydrogenase, is absent in C. glandulosa. A total of 249 microsatellites were found in the chloroplast genome of C. farinosa, 251 in C. glandulosa, 227 in M. crassifolia and 233 in M. oblongifolia, the majority of which are mononucleotides A/T found in the intergenic spacer. Comparative analysis revealed variable hotspot regions (atpF, rpoC2, rps19 and ycf1), which can be used as molecular markers for species authentication and as regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. The monophyly of Capparaceae and other families under Brassicales, as well as the phylogenetic positions of the studied species, are highly supported by all the relationships in the phylogenetic tree. The cp genomes reported in this study will provide resources for studying the genetic diversity of Capparaceae, as well as resolving phylogenetic relationships within the family.

Complete chloroplast genome sequencing of Job's tears (Coix L.): genome structure, comparative analysis, and phylogenetic relationships

Job's tears, also known as adlay, is a valuable plant that has commonly been used in traditional Chinese medicine, as well as an edible food. Due to the lack of knowledge of its genetics and gaps in its evolutionary analysis, breeding of adlay has been hindered. Here, we report five complete chloroplast genomes of various species and varieties in the genus by Illumina sequencing, while their genome structure, comparative analysis, and phylogenetic relationships were conducted. Genome sizes ranged from 140,860 to 140,864 bp in length, GC contents were 38.43%, and genome architecture was of a typical quadripartite structure. We annotated 82～83 protein-coding genes and 46～47 non-coding RNA genes in each genome and they functionally associated with self-replication, photosynthesis, cytochrome synthesis and other unknown functions. Three codons that encoded tryptophan, arginine and leucine were used frequently at rates of 41.42, 37.98, and 32.28% respectively. The preferred codons consistently ended with A or T. A total of 146 simple sequence repeats (SSR), 9 insertions and deletions (InDels) and 143 single nucleotide polymorphisms (SNPs) were observed among genomes. The InDel and SNP variations were mostly distributed in intergenic regions. It confirmed that Coix, Sorghum, Saccharum, Zea, Tripsacum and Saccharum were closely genera and the genetic distance of Sorghum to Coix was closer than Zea to Coix. These results give us more insight into the evolution of Coix in a wide range of evolutionary studies.

Complete chloroplast genome sequences of Dipterygium glaucum and Cleome chrysantha and other Cleomaceae Species, comparative analysis and phylogenetic relationships

This current study presents, for the first time, the complete chloroplast genome of two Cleomaceae species: Dipterygium glaucum and Cleome chrysantha in order to evaluate the evolutionary relationship. The cp genome is 158,576 bp in length with 35.74% GC content in D. glaucum and 158,111 bp with 35.96% GC in C. chrysantha. Inverted repeats IR 26,209 bp, 26,251 bp each, LSC of 87,738 bp, 87,184 bp and SSC of 18,420 bp, 18,425 bp respectively. There are 136 genes in the genome, which includes 80 protein coding genes, 31 tRNA genes and four rRNA genes were observed in both chloroplast genomes. 117 genes are unique while the remaining 19 genes are duplicated in IR regions. The analysis of repeats shows that the cp genome includes all types of repeats with more frequent occurrences of palindromic; Also, this analysis indicates that the total number of simple sequence repeats (SSR) were 323 in D. glaucum, and 313 in C. chrysantha, of which the majority of the SSRs in these plastid genomes were mononucleotide repeats A/T which are located in the intergenic spacer. Moreover, the comparative analysis of the four cp sequences revealed four hotspot genes (atpF, rpoC2, rps19, and ycf1), these variable regions could be used as molecular makers for the species authentication as well as resources for inferring phylogenetic relationships of the species. All the relationships in the phylogenetic tree are with high support, this indicate that the complete chloroplast genome is a useful data for inferring phylogenetic relationship within the Cleomaceae and other families. The simple sequence repeats identified will be useful for identification, genetic diversity, and other evolutionary studies of the species. This study reported the first cp genome of the genus Dipterygium and Cleome. The finding of this study will be beneficial for biological disciplines such as evolutionary and genetic diversity studies of the species within the core Cleomaceae.

Chloroplast genome sequence of Chongming lima bean (Phaseolus lunatus L.) and comparative analyses with other legume chloroplast genomes

Background

Lima bean (Phaseolus lunatus L.) is a member of subfamily Phaseolinae belonging to the family Leguminosae and an important source of plant proteins for the human diet. As we all know, lima beans have important economic value and great diversity. However, our knowledge of the chloroplast genome level of lima beans is limited.

Results

The chloroplast genome of lima bean was obtained by Illumina sequencing technology for the first time. The Cp genome with a length of 150,902 bp, including a pair of inverted repeats (IRA and IRB 26543 bp each), a large single-copy (LSC 80218 bp) and a small single-copy region (SSC 17598 bp). In total, 124 unique genes including 82 protein-coding genes, 34 tRNA genes, and 8 rRNA genes were identified in the P. lunatus Cp genome. A total of 61 long repeats and 290 SSRs were detected in the lima bean Cp genome. It has a typical 50 kb inversion of the Leguminosae family and an 70 kb inversion to subtribe Phaseolinae. rpl16, accD, petB, rsp16, clpP, ndhA, ndhF and ycf1 genes in coding regions was found significant variation, the intergenic regions of trnk-rbcL, rbcL-atpB, ndhJ-rps4, psbD-rpoB, atpI-atpA, atpA-accD, accD-psbJ, psbE-psbB, rsp11-rsp19, ndhF-ccsA was found in a high degree of divergence. A phylogenetic analysis showed that P. lunatus appears to be more closely related to P. vulgaris, V.unguiculata and V. radiata.

Conclusions

The characteristics of the lima bean Cp genome was identified for the first time, these results will provide useful insights for species identification, evolutionary studies and molecular biology research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07467-8.

Comparative analysis and phylogenetic investigation of Hong Kong Ilex chloroplast genomes

Ilex is a monogeneric plant group (containing approximately 600 species) in the Aquifoliaceae family and one of the most commonly used medicinal herbs. However, its taxonomy and phylogenetic relationships at the species level are debatable. Herein, we obtained the complete chloroplast genomes of all 19 Ilex types that are native to Hong Kong. The genomes are conserved in structure, gene content and arrangement. The chloroplast genomes range in size from 157,119 bp in Ilex graciliflora to 158,020 bp in Ilex kwangtungensis. All these genomes contain 125 genes, of which 88 are protein-coding and 37 are tRNA genes. Four highly varied sequences (rps16-trnQ, rpl32-trnL, ndhD-psaC and ycf1) were found. The number of repeats in the Ilex genomes is mostly conserved, but the number of repeating motifs varies. The phylogenetic relationship among the 19 Ilex genomes, together with eight other available genomes in other studies, was investigated. Most of the species could be correctly assigned to the section or even series level, consistent with previous taxonomy, except Ilex rotunda var. microcarpa, Ilex asprella var. tapuensis and Ilex chapaensis. These species were reclassified; I. rotunda was placed in the section Micrococca, while the other two were grouped with the section Pseudoaquifolium. These studies provide a better understanding of Ilex phylogeny and refine its classification.

Exploring the edible gum (galactomannan) biosynthesis and its regulation during pod developmental stages in clusterbean using comparative transcriptomic approach

Galactomannan is a polymer of high economic importance and is extracted from the seed endosperm of clusterbean (C. tetragonoloba). In the present study, we worked to reveal the stage-specific galactomannan biosynthesis and its regulation in clusterbean. Combined electron microscopy and biochemical analysis revealed high protein and gum content in RGC-936, while high oil bodies and low gum content in M-83. A comparative transcriptome study was performed between RGC-936 (high gum) and M-83 (low gum) varieties at three developmental stages viz. 25, 39, and 50 days after flowering (DAF). Total 209,525, 375,595 and 255,401 unigenes were found at 25, 39 and 50 DAF respectively. Differentially expressed genes (DEGs) analysis indicated a total of 5147 shared unigenes between the two genotypes. Overall expression levels of transcripts at 39DAF were higher than 50DAF and 25DAF. Besides, 691 (RGC-936) and 188 (M-83) candidate unigenes that encode for enzymes involved in the biosynthesis of galactomannan were identified and analyzed, and 15 key enzyme genes were experimentally validated by quantitative Real-Time PCR. Transcription factor (TF) WRKY was observed to be co-expressed with key genes of galactomannan biosynthesis at 39DAF. We conclude that WRKY might be a potential biotechnological target (subject to functional validation) for developing high gum content varieties.

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.

Amazing symmetrical clustering in chloroplast genomes

Background

Previously, a seven-cluster pattern claiming to be a universal one in bacterial genomes has been reported. Keeping in mind the most popular theory of chloroplast origin, we checked whether a similar pattern is observed in chloroplast genomes.

Results

Surprisingly, eight cluster structure has been found, for chloroplasts. The pattern observed for chloroplasts differs rather significantly, from bacterial one, and from that latter observed for cyanobacteria. The structure is provided by clustering of the fragments of equal length isolated within a genome so that each fragment is converted in triplet frequency dictionary with non-overlapping triplets with no gaps in frame tiling. The points in 63-dimensional space were clustered due to elastic map technique. The eight cluster found in chloroplasts comprises the fragments of a genome bearing tRNA genes and exhibiting excessively high GC-content, in comparison to the entire genome.

Conclusion

Chloroplasts exhibit very specific symmetry type in distribution of coding and non-coding fragments of a genome in the space of triplet frequencies: this is mirror symmetry. Cyanobacteria may have both mirror symmetry, and the rotational symmetry typical for other bacteria.

Methods and Tools for Plant Organelle Genome Sequencing, Assembly, and Downstream Analysis

Organelles play an important role in a eukaryotic cell. Among them, the two organelles, chloroplast and mitochondria, are responsible for the critical function of photosynthesis and aerobic respiration. Organellar genomes are also very important for plant systematic studies. Here we have described the methods for isolation of the mitochondrial and plastid DNA and its subsequent sequencing with the help of NGS technology. We have also discussed in detail the various tools available for assembly, annotation, and visualization of the organelle genome sequence.

New Insights Into the Plastome Evolution of the Millettioid/Phaseoloid Clade (Papilionoideae, Leguminosae)

The Millettioid/Phaseoloid (MP) clade from the subfamily Papilionoideae (Leguminosae) consists of six tribes and ca. 3,000 species. Previous studies have revealed some plastome structural variations (PSVs) within this clade. However, many deep evolutionary relationships within the clade remain unresolved. Due to limited taxon sampling and few genetic markers in previous studies, our understanding of the evolutionary history of this clade is limited. To address this issue, we sampled 43 plastomes (35 newly sequenced) representing all the six tribes of the MP clade to examine genomic structural variations and phylogenetic relationships. Plastomes of the species from the MP clade were typically quadripartite (size ranged from 140,029 to 160,040 bp) and contained 109-111 unique genes. We revealed four independent gene losses (ndhF, psbI, rps16, and trnS-GCU), multiple IR-SC boundary shifts, and six inversions in the tribes Desmodieae, Millettieae, and Phaseoleae. Plastomes of the species from the MP clade have experienced significant variations which provide valuable information on the evolution of the clade. Plastid phylogenomic analyses using Maximum Likelihood and Bayesian methods yielded a well-resolved phylogeny at the tribal and generic levels within the MP clade. This result indicates that plastome data is useful and reliable data for resolving the evolutionary relationships of the MP clade. This study provides new insights into the phylogenetic relationships and PSVs within this clade.

Elucidation of Galactomannan Biosynthesis Pathway Genes through Transcriptome Sequencing of Seeds Collected at Different Developmental Stages of Commercially Important Indian Varieties of Cluster Bean (Cyamopsis tetragonoloba L.)

Cyamopsis tetragonoloba (L) endosperm predominantly contains guar gum a polysaccharide, which has tremendous industrial applications in food, textile, paper, oil drilling and water treatment. In order to understand the genes controlling galactomannan biosynthesis, mRNA was isolated from seeds collected at different developmental stages; young pods, mature pods and young leaf from two guar varieties, HG365 and HG870 and subjected to Illumina sequencing. De novo assembly of fourteen individual read files from two varieties of guar representing seven developmental stages gave a total of 1,13,607 contigs with an N50 of 1,244 bases. Annotation of assemblies with GO mapping revealed three levels of distribution, namely, Biological Processes, Molecular Functions and Cellular Components. GO studies identified major genes involved in galactomannan biosynthesis: Cellulose synthase D1 (CS D1) and GAUT-like gene families. Among the polysaccharide biosynthetic process (GO:0000271) genes the transcript abundance for CS was found to be predominantly more in leaf samples, whereas, the transcript abundance for GAUT-like steadily increased from 65% to 90% and above from stage1 to stage5 indicating accumulation of galactomannan in developing seeds; and validated by qRT-PCR analysis. Galactomannan quantification by HPLC showed HG365 (12.98–20.66%) and HG870 (7.035–41.2%) gradually increasing from stage1 to stage 5 (10–50 DAA) and highest accumulation occurred in mature and dry seeds with 3.8 to 7.1 fold increase, respectively. This is the first report of transcriptome sequencing and complete profiling of guar seeds at different developmental stages, young pods, mature pods and young leaf material from two commercially important Indian varieties and elucidation of galactomannan biosynthesis pathway. It is envisaged that the data presented herein will be very useful for improvement of guar through biotechnological interventions in future.

The genome size of clusterbean (Cyamopsis tetragonoloba) is significantly smaller compared to its wild relatives as estimated by flow cytometry

Clusterbean (C. tetragonoloba) is an important, leguminous vegetable and industrial crop with vast genetic diversity but meager genetic, cytological and genomic information. In the present study, an optimized procedure of flow cytometry was used to estimate the genome size of three clusterbean species, represented by C. tetragonoloba (cv. RGC-936) and two wild relatives (C. serreta and C. senegalensis). For accurate estimation of genomic content, singlet G₀/G₁ populations of multiple tissues such as leaves, hypocotyl, and matured seeds were determined and used along with three different plant species viz. Pisum sativum (as primary), Oryza sativa, and Glycine max (secondary), as external and internal reference standards. Seed tissue of the test sample and G. max provided the best estimate of nuclear DNA content in comparison to other sample tissues and reference standards. The genome size of C. tetragonoloba was detemined at 580.9±0.02Mbp (1C), while that of C. serreta and C. senegalensis was estimated at 979.6±0.02Mbp (1C) and 943.4±0.03Mbp (1C), respectively. Thus, the wild relatives harbor, nearly double the genome content of the cultivated cluster bean. Findings of this study will enrich genomic database of the legume family and can serve as the starting point for clusterbean evolutionary and genomics studies.

Comparison of Four Complete Chloroplast Genomes of Medicinal and Ornamental Meconopsis Species: Genome Organization and Species Discrimination

High-throughput sequencing of chloroplast genomes has been used to gain insight into the evolutionary relationships of plant species. In this study, we sequenced the complete chloroplast genomes of four species in the Meconopsis genus: M. racemosa, M. integrifolia (Maxim.) Franch, M. horridula and M. punicea. These plants grow in the wild and are recognized as having important medicinal and ornamental applications. The sequencing results showed that the size of the Meconopsis chloroplast genome ranges from 151864 to 153816 bp. A total of 127 genes comprising 90 protein-coding genes, 37 tRNA genes and 8 rRNA genes were observed in all four chloroplast genomes. Comparative analysis of the four chloroplast genomes revealed five hotspot regions (matK, rpoC2, petA, ndhF, and ycf1), which could potentially be used as unique molecular markers for species identification. In addition, the ycf1 gene may also be used as an effective molecular marker to distinguish Papaveraceae and determine the evolutionary relationships among plant species in the Papaveraceae family. Futhermore, these four genomes can provide valuable genetic information for other related studies.

Comprehensive Analysis of Rhodomyrtus tomentosa Chloroplast Genome

In the last decade, several studies have relied on a small number of plastid genomes to deduce deep phylogenetic relationships in the species-rich Myrtaceae. Nevertheless, the plastome of Rhodomyrtus tomentosa, an important representative plant of the Rhodomyrtus (DC.) genera, has not yet been reported yet. Here, we sequenced and analyzed the complete chloroplast (CP) genome of R. tomentosa, which is a 156,129-bp-long circular molecule with 37.1% GC content. This CP genome displays a typical quadripartite structure with two inverted repeats (IRa and IRb), of 25,824 bp each, that are separated by a small single copy region (SSC, 18,183 bp) and one large single copy region (LSC, 86,298 bp). The CP genome encodes 129 genes, including 84 protein-coding genes, 37 tRNA genes, eight rRNA genes and three pseudogenes (ycf1, rps19, ndhF). A considerable number of protein-coding genes have a universal ATG start codon, except for psbL and ndhD. Premature termination codons (PTCs) were found in one protein-coding gene, namely atpE, which is rarely reported in the CP genome of plants. Phylogenetic analysis revealed that R. tomentosa has a sister relationship with Eugenia uniflora and Psidium guajava. In conclusion, this study identified unique characteristics of the R. tomentosa CP genome providing valuable information for further investigations on species identification and the phylogenetic evolution between R. tomentosa and related species.

Green Biotechnology: A Brief Update on Plastid Genome Engineering

Plant genetic engineering has become an inevitable tool in the molecular breeding of crops. Significant progress has been made in the generation of novel plastid transformation vectors and optimized transformation protocols. There are several advantages of plastid genome engineering over conventional nuclear transformation. Some of the advantages include multigene engineering by expression of biosynthetic pathway genes as operons, extremely high-level expression of protein accumulation, lack of transgene silencing, etc. Transgene containment owing to maternal inheritance is another important advantage of plastid genome engineering. Chloroplast genome modification usually results in alteration of several thousand plastid genome copies in a cell. Several therapeutic proteins, edible vaccines, antimicrobial peptides, and industrially important enzymes have been successfully expressed in chloroplasts so far. Here, we critically recapitulate the latest developments in plastid genome engineering. Latest advancements in plastid genome sequencing are briefed. In addition, advancement of extending the toolbox for plastid engineering for selected applications in the area of molecular farming and production of industrially important enzyme is briefed.

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.

High Quality Unigenes and Microsatellite Markers from Tissue Specific Transcriptome and Development of a Database in Clusterbean (Cyamopsis tetragonoloba, L. Taub)

Clusterbean (Cyamopsis tetragonoloba L. Taub), is an important industrial, vegetable and forage crop. This crop owes its commercial importance to the presence of guar gum (galactomannans) in its endosperm which is used as a lubricant in a range of industries. Despite its relevance to agriculture and industry, genomic resources available in this crop are limited. Therefore, the present study was undertaken to generate RNA-Seq based transcriptome from leaf, shoot, and flower tissues. A total of 145 million high quality Illumina reads were assembled using Trinity into 127,706 transcripts and 48,007 non-redundant high quality (HQ) unigenes. We annotated 79% unigenes against Plant Genes from the National Center for Biotechnology Information (NCBI), Swiss-Prot, Pfam, gene ontology (GO) and KEGG databases. Among the annotated unigenes, 30,020 were assigned with 116,964 GO terms, 9984 with EC and 6111 with 137 KEGG pathways. At different fragments per kilobase of transcript per millions fragments sequenced (FPKM) levels, genes were found expressed higher in flower tissue followed by shoot and leaf. Additionally, we identified 8687 potential simple sequence repeats (SSRs) with an average frequency of one SSR per 8.75 kb. A total of 28 amplified SSRs in 21 clusterbean genotypes resulted in polymorphism in 13 markers with average polymorphic information content (PIC) of 0.21. We also constructed a database named ‘ClustergeneDB’ for easy retrieval of unigenes and the microsatellite markers. The tissue specific genes identified and the molecular marker resources developed in this study is expected to aid in genetic improvement of clusterbean for its end use.