Genome sequencing and functional analysis of a multipurpose medicinal herb Tinospora cordifolia (Giloy)

Tinospora cordifolia (Willd.) Hook.f. & Thomson, also known as Giloy, is among the most important medicinal plants that have numerous therapeutic applications in human health due to the production of a diverse array of secondary metabolites. To gain genomic insights into the medicinal properties of T. cordifolia, the genome sequencing was carried out using 10× Genomics linked read and Nanopore long-read technologies. The draft genome assembly of T. cordifolia was comprised of 1.01 Gbp, which is the genome sequenced from the plant family Menispermaceae. We also performed the genome size estimation for T. cordifolia, which was found to be 1.13 Gbp. The deep sequencing of transcriptome from the leaf tissue was also performed. The genome and transcriptome assemblies were used to construct the gene set, resulting in 17,245 coding gene sequences. Further, the phylogenetic position of T. cordifolia was also positioned as basal eudicot by constructing a genome-wide phylogenetic tree using multiple species. Further, a comprehensive comparative evolutionary analysis of gene families contraction/expansion and multiple signatures of adaptive evolution was performed. The genes involved in benzyl iso-quinoline alkaloid, terpenoid, lignin and flavonoid biosynthesis pathways were found with signatures of adaptive evolution. These evolutionary adaptations in genes provide genomic insights into the presence of diverse medicinal properties of this plant. The genes involved in the common symbiosis signalling pathway associated with endosymbiosis (Arbuscular Mycorrhiza) were found to be adaptively evolved. The genes involved in adventitious root formation, peroxisome biogenesis, biosynthesis of phytohormones, and tolerance against abiotic and biotic stresses were also found to be adaptively evolved in T. cordifolia.

Genome sequencing of Syzygium cumini (jamun) reveals adaptive evolution in secondary metabolism pathways associated with its medicinal properties

Syzygium cumini, also known as jambolan or jamun, is an evergreen tree widely known for its medicinal properties, fruits, and ornamental value. To understand the genomic and evolutionary basis of its medicinal properties, we sequenced S. cumini genome for the first time from the world's largest tree genus Syzygium using Oxford Nanopore and 10x Genomics sequencing technologies. We also sequenced and assembled the transcriptome of S. cumini in this study. The tetraploid and highly heterozygous draft genome of S. cumini had a total size of 709.9 Mbp with 61,195 coding genes. The phylogenetic position of S. cumini was established using a comprehensive genome-wide analysis including species from 18 Eudicot plant orders. The existence of neopolyploidy in S. cumini was evident from the higher number of coding genes and expanded gene families resulting from gene duplication events compared to the other two sequenced species from this genus. Comparative evolutionary analyses showed the adaptive evolution of genes involved in the phenylpropanoid-flavonoid (PF) biosynthesis pathway and other secondary metabolites biosynthesis such as terpenoid and alkaloid in S. cumini, along with genes involved in stress tolerance mechanisms, which was also supported by leaf transcriptome data generated in this study. The adaptive evolution of secondary metabolism pathways is associated with the wide range of pharmacological properties, specifically the anti-diabetic property, of this species conferred by the bioactive compounds that act as nutraceutical agents in modern medicine.

Genome sequencing and de novo and reference-based genome assemblies of Bos indicus breeds

BACKGROUND

Indian cattle breeds (Bos indicus) are known for their remarkable adaptability to hot and humid climates, higher nutritious quality of milk, better disease tolerance, and greater ability to perform in poor feed compared to taurine cattle (Bos taurus). Distinct phenotypic differences are observed among the B. indicus breeds; however, the whole genome sequences were unavailable for these indigenous breeds.

OBJECTIVE

We aimed to perform whole genome sequencing to construct the draft genome assemblies of four B. indicus breeds; Ongole, Kasargod Dwarf, Kasargod Kapila, and Vechur (the smallest cattle of the world).

METHODS

We sequenced the whole genomes using Illumina short-read technology, and constructed de novo and reference-based genome assemblies of these native B. indicus breeds for the first time.

RESULTS

The draft de novo genome assemblies of B. indicus breeds ranged from 1.98 to 3.42 Gbp. We also constructed the mitochondrial genome assemblies (~ 16.3 Kbp), and yet unavailable 18S rRNA marker gene sequences of these B. indicus breeds. The genome assemblies helped to identify the bovine genes related to distinct phenotypic characteristics and other biological processes for this species compared to B. taurus, which are plausibly responsible for providing better adaptive traits. We also identified the genes that showed sequence variation in dwarf and non-dwarf breeds of B. indicus compared to B. taurus.

CONCLUSIONS

The genome assemblies of these Indian cattle breeds, the 18S rRNA marker genes, and identification of the distinct genes in B. indicus breeds compared to B. taurus will help in future studies on these cattle species.

Genome of Phyllanthus emblica: the medicinal plant Amla with super antioxidant properties

Phyllanthus emblica or Indian gooseberry, commonly known as amla, is an important medicinal horticultural plant used in traditional and modern medicines. It bears stone fruits with immense antioxidant properties due to being one of the richest natural sources of vitamin C and numerous flavonoids. This study presents the first genome sequencing of this species performed using 10x Genomics and Oxford Nanopore Technology. The draft genome assembly was 519 Mbp in size and consisted of 4,384 contigs, N50 of 597 Kbp, 98.4% BUSCO score, and 37,858 coding sequences. This study also reports the genome-wide phylogeny of this species with 26 other plant species that resolved the phylogenetic position of P. emblica. The presence of three ascorbate biosynthesis pathways including L-galactose, galacturonate, and myo-inositol pathways was confirmed in this genome. A comprehensive comparative evolutionary genomic analysis including gene family expansion/contraction and identification of multiple signatures of adaptive evolution provided evolutionary insights into ascorbate and flavonoid biosynthesis pathways and stone fruit formation through lignin biosynthesis. The availability of this genome will be beneficial for its horticultural, medicinal, dietary, and cosmetic applications and will also help in comparative genomics analysis studies.

Mapping of rDNA on the chromosomes of Eleusine species by fluorescence in situ hybridization

Mapping of rDNA sites on the chromosomes of four diploid and two tetraploid species of Eleusine has provided valuable information on genome relationship between the species. Presence of 18S-5.8S-26S rDNA on the largest pair of the chromosomes, location of 5S rDNA at four sites on two pairs of chromosomes and presence of 18S-5.8S-26S and 5S rDNA at same location on one pair of chromosomes have clearly differentiated E. multiflora from rest of the species of Eleusine. The two tetraploid species, E. coracana and E. africana have the same number of 18S-5.8S-26S and 5S rDNA sites and located at similar position on the chromosomes. Diploid species, E. indica, E. floccifolia and E. tristachya have the same 18S-5.8S-26S sites and location on the chromosomes which also resembled with the two pairs of 18S-5.8S-26S rDNA locations in tetraploid species, E. coracana and E. africana. The 5S rDNA sites on chromosomes of E. indica and E. floccifolia were also comparable to the 5S rDNA sites of E. africana and E. coracana. The similarity of the rDNA sites and their location on chromosomes in the three diploid and two polyploid species also supports the view that genome donors to tetraploid species may be from these diploid species.

One year outcome of babies with severe neonatal hyperbilirubinemia and reversible abnormality in brainstem auditory evoked responses

Brainstem auditory evoked responses (BAER) were longitudinally recorded prospectively in 18 term infants with neonatal hyperbilirubinemia (NHB) (total serum bilirubin > 15 mg/dl). Seven neonates had abnormal BAER. Wave complex IV-V was absent in eight recordings in NHB group while they were normal in the control group (p < 0.001). Prolongation of latency of waves I and V and interwave conduction time (wave I-V) occurred in jaundiced infants especially when unconjugated serum bilirubin level rose above 22 mg/dl. The abnormalities in BAER reversed to normal in all seven neonates after exchange blood transfusion indicating transient nature of bilirubin toxicity to the brain. All seven neonates in the study and control group had normal hearing, development quotient and were free of neurological sequelae on follow up for one year.

Auditory brain stem evoked response in normal neonates and infants

Auditory brainstem evoked responses (ABER) were recorded to provide normative data regarding the latency values, brainstem conduction times in term neonates and infants up to 24 months of age. The mean latency for wave I - 1.86 +/- 0.11 m sec, wave III - 5.11 +/- 0.31 m sec, wave V - 7.10 +/- 0.30 m sec progressively decreased at birth to wave I - 1.64 +/- 0.17 m sec, wave III - 4.40 +/- 0.27 m sec, wave V - 6.15 +/- 0.21 m sec at 12 months of age by stimulus strength of 75 decibels at a rate 20 clicks/sec. Interwave conduction times improved to I-III - 2.76 +/- 0.29 m sec, I-V - 4.51 +/- 0.29 m sec, III-V - 1.75 +/- 0.17 m sec at 12 months from I-III - 3.25 +/- 0.26 m sec, I-V - 5.23 +/- 0.23 m sec, III-V - 1.99 +/- 0.35 m sec at birth.

Involvement of central nervous system in diabetes mellitus

Brainstem auditory evoked responses were recorded in 22 diabetic patients with a variable duration of illness (mean 5.8 years) and 14 normal healthy controls of comparable age. The initial 10 millisecond components, found to be most consistent and reproducible, were analysed. Variations in the form of individual wave latency, interpeak latencies and V wave amplitude were compared in both the groups. No difference was found in any of the parameters. It was concluded that central neural pathways are not involved at least initially in diabetes mellitus.