Anita Bhandari

Expression analysis of octopamine and tyramine receptors in Drosophila

The monoamines octopamine and tyramine, which are the invertebrate counterparts of epinephrine and norepinephrine, transmit their action through sets of G protein-coupled receptors. Four different octopamine receptors (Oamb, Octß1R, Octß2R, Octß3R) and 3 different tyramine receptors (TyrR, TyrRII, TyrRIII) are present in the fruit fly Drosophila melanogaster. Utilizing the presumptive promoter regions of all 7 octopamine and tyramine receptors, the Gal4/UAS system is utilized to elucidate their complete expression pattern in larvae as well as in adult flies. All these receptors show strong expression in the nervous system but their exact expression patterns vary substantially. Common to all octopamine and tyramine receptors is their expression in mushroom bodies, centers for learning and memory in insects. Outside the central nervous system, the differences in the expression patterns are more conspicuous. However, four of them are present in the tracheal system, where they show different regional preferences within this organ. On the other hand, TyrR appears to be the only receptor present in the heart muscles and TyrRII the only one expressed in oenocytes. Skeletal muscles express octß2R, Oamb and TyrRIII, with octß2R being present in almost all larval muscles. Taken together, this study provides comprehensive information about the sites of expression of all octopamine and tyramine receptors in the fruit fly, thus facilitating future research in the field.

Sequence, phylogenetic and variant analyses of antithrombin III

Antithrombin III (ATIII) performs a critical anticoagulant function by precluding the activation of blood clotting proteinases, aided by its two cofactors, heparin and heparan sulfate. Though several studies have been carried out on physiological, biochemical and structural perspectives on ATIII, so far there are limited studies on the molecular evolution of ATIII. Herein, we carried out molecular phylogenetic analyses of ATIII genes, combining gene structures, synteny and sequence-structural features for ATIII spanning 50 vertebrate genomes. ATIII is maintained over 450 MY on same genomic loci in vertebrates with few changes in ray-finned fishes and lost one intron 262c in tetrapods and coelacanth. In ray-finned fishes, ATIII gene has additional intron at the position 262c, which shared by group V1 members, corroborating that it is lost in other vertebrates and also in lobed-finned fish coelacanth (Latimeria chalumnae). We found that heparin binding basic residues, hD helix, three pairs of Cys-Cys salt bridges, N-glycosylation sites, serpin motifs and inhibitory reactive center loop (RCL) of ATIII protein are highly conserved. Using 1092 human genomes available from 1000G project, we also compiled 1997 ATIII variants, which reveals 76.2% single nucleotide polymorphisms (SNPs), 11.8% deletions and 8.1% insertions as three major classes of gene variations. These understandings may have medical importance as well.

Conservation of tubulin-binding sequences in TRPV1 throughout evolution.

Background Transient Receptor Potential Vanilloid sub type 1 (TRPV1), commonly known as capsaicin receptor can detect multiple stimuli ranging from noxious compounds, low pH, temperature as well as electromagnetic wave at different ranges. In addition, this receptor is involved in multiple physiological and sensory processes. Therefore, functions of TRPV1 have direct influences on adaptation and further evolution also. Availability of various eukaryotic genomic sequences in public domain facilitates us in studying the molecular evolution of TRPV1 protein and the respective conservation of certain domains, motifs and interacting regions that are functionally important. Methodology and Principal Findings Using statistical and bioinformatics tools, our analysis reveals that TRPV1 has evolved about ∼420 million years ago (MYA). Our analysis reveals that specific regions, domains and motifs of TRPV1 has gone through different selection pressure and thus have different levels of conservation. We found that among all, TRP box is the most conserved and thus have functional significance. Our results also indicate that the tubulin binding sequences (TBS) have evolutionary significance as these stretch sequences are more conserved than many other essential regions of TRPV1. The overall distribution of positively charged residues within the TBS motifs is conserved throughout evolution. In silico analysis reveals that the TBS-1 and TBS-2 of TRPV1 can form helical structures and may play important role in TRPV1 function. Conclusions and Significance Our analysis identifies the regions of TRPV1, which are important for structure – function relationship. This analysis indicates that tubulin binding sequence-1 (TBS-1) near the TRP-box forms a potential helix and the tubulin interactions with TRPV1 via TBS-1 have evolutionary significance. This interaction may be required for the proper channel function and regulation and may also have significance in the context of Taxol®-induced neuropathy.

Revising angiotensinogen from phylogenetic and genetic variants perspectives

Angiotensinogen (AGT) belongs to the serpin superfamily. It acts as the unique substrate of all angiotensin peptides, which generates a spectrum of angiotensin peptides in the renin-angiotensin system and regulates hypertension. This serpin belongs to the multiple member group V2 of the intron encoded vertebrate serpin classification. Despite huge advancements in the understanding of angiotensinogen based on biochemical properties and its roles in the RAS, phylogenetic history of AGT remains forgotten. To date, there is no comprehensive study illustrating the phylogenetic history of AGT. Herein, we investigated phylogenetic traits of AGT gene across vertebrates. Gene structures of AGT gene from selected ray-finned fishes varied in exon I and II with insertions of two novel introns in the core domain for ray-finned fishes at the position 77c and 233c. We that found AGT loci is conserved from lampreys to human and estimated to be older than 500 MY. By comparing AGT protein in 57 vertebrate genomes, we illustrated that the reactive center loop (RCL) of AGT protein became from inhibitory (in lampreys, GTEAKAETVVGIMPI†SMPPT) to non-inhibitory (in human, EREPTESTQQLNKPE†VLEVT) during period of 500 MY. We identified 690 AGT variants by analysis of 1092 human genomes with top three variation classes belongs to SNPs (89.7%), somatic SNVs (5.2%) and deletion (2.9%). There are 32 key residues out of 121 missense variants, which are deleterious for AGT protein, computed by combination of SIFT and PolyPhen V2 methods. These results may have clinical implications for understanding hypertension.

Molecular Phylogeny of OVOL Genes Illustrates a Conserved C2H2 Zinc Finger Domain Coupled by Hypervariable Unstructured Regions

OVO-like proteins (OVOL) are members of the zinc finger protein family and serve as transcription factors to regulate gene expression in various differentiation processes. Recent studies have shown that OVOL genes are involved in epithelial development and differentiation in a wide variety of organisms; yet there is a lack of comprehensive studies that describe OVOL proteins from an evolutionary perspective. Using comparative genomic analysis, we traced three different OVOL genes (OVOL1-3) in vertebrates. One gene, OVOL3, was duplicated during a whole-genome-duplication event in fish, but only the copy (OVOL3b) was retained. From early-branching metazoa to humans, we found that a core domain, comprising a tetrad of C2H2 zinc fingers, is conserved. By domain comparison of the OVOL proteins, we found that they evolved in different metazoan lineages by attaching intrinsically-disordered (ID) segments of N/C-terminal extensions of 100 to 1000 amino acids to this conserved core. These ID regions originated independently across different animal lineages giving rise to different types of OVOL genes over the course of metazoan evolution. We illustrated the molecular evolution of metazoan OVOL genes over a period of 700 million years (MY). This study both extends our current understanding of the structure/function relationship of metazoan OVOL genes, and assembles a good platform for further characterization of OVOL genes from diverged organisms.

A model species for agricultural pest genomics: The genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

Genetic variants and evolutionary analyses of heparin cofactor II.

Heparin cofactor II (HCII) belongs to serpin superfamily and it acts as a thrombin inhibitor in the coagulation cascade, in a glycosaminoglycan-dependent pathway using the release of a sequestered hirudin-like N-terminal tail for interaction with thrombin. This serpin belongs to multiple member group V2 of vertebrate serpin classification. However, there is no comprehensive study illustrating the exact phylogenetic history of HCII, to date. Herein, we explored phylogenetic traits of HCII genes. Structures of HCII gene from selected ray-finned fishes and lamprey varied in exon I and II with insertions of novel introns of which one in core domain for ray-finned fishes in exon II at the position 241c. We found HCII remain nested in the largest intron of phosphatidylinositol (PI) 4-kinase (PIK4CA) gene (genetic variants of this gene cause schizophrenia) at the origin of vertebrates, dated about 500MY old. We found that sequence features such as two acidic repeats (AR1-II), GAG-binding helix-D, three serpin motifs and inhibitory reactive center loop (RCL) of HCII protein are highly conserved in 55 vertebrates analyzed. We identified 985 HCII variants by analysis of 1092 human genomes with top three variation classes belongs to SNPs (84.3%), insertion (7.1%) and deletion (5.0%). We identified 37 deleterious mutations in the human HCII protein and we have described these mutations in relation to HCII sequence-structure-function relationships. These understandings may have clinical and medical importance as well.

Molecular phylogeny of C1 inhibitor depicts two immunoglobulin-like domains fusion in fishes and ray-finned fishes specific intron insertion after separation from zebrafish

C1 inhibitor (C1IN) is a multi-facet serine protease inhibitor in the plasma cascades, inhibiting several proteases, notably, regulates both complement and contact system activation. Despite huge advancements in the understanding of C1IN based on biochemical properties and its roles in the plasma cascades, the phylogenetic history of C1IN remains uncharacterized. To date, there is no comprehensive study illustrating the phylogenetic history of C1IN. Herein, we explored phylogenetic history of C1IN gene in vertebrates. Fishes have C1IN with two immunoglobulin like domains attached in the N-terminal region. The RCL regions of CIIN from fishes and tetrapod genomes have variations at the positions P2 and P1. Gene structures of C1IN gene from selected ray-finned fishes varied in the Ig domain region with creation of novel intron splitting exon Im2 into Im2a and Im2b. This intron is limited to ray-finned fishes with genome size reduced below 1 Gb. Hence, we suggest that genome compaction and associated double-strand break repairs are behind this intron gain. This study reveals the evolutionary history of C1IN and confirmed that this gene remains the same locus for ∼450 MY in 52 vertebrates analysed, but it is not found in frogs and lampreys.

Understanding V(D)J recombination initiator RAG1 gene using molecular phylogenetic and genetic variant analyses and upgrading missense and non-coding variants of clinical importance.

The recombination-activating genes (RAGs) encode for V(D)J recombinases responsible for rearrangements of antigen-receptor genes during T and B cell development, and RAG expression is known to correlate strictly with the process of rearrangement. There have been several studies of RAG1 illustrating biochemical, physiological and immunological properties. Hitherto, there are limited studies on RAG1 focusing molecular phylogenetic analyses, evolutionary traits, and genetic variants in human populations. Hence, there is a need of a comprehensive study on this topic. In the current report, we have shed light into insights of evolutionary traits and genetic variants of human RAG1 gene using 1092 genomes from human populations. Syntenic analyses revealed that two RAG genes are physically linked and conserved on the same locus in head-to-head orientation from sea urchin to human for about 550 MY. Spliceosomal introns have been in invaded in fishes and sea urchin, whereas gene structures of RAG1 gene from tetrapods remained single exon architecture. We compiled 751 genetic variants in human RAG1 gene using 1092 human genomes; where major stockholders of variant classes are 79% single nucleotide polymorphisms (SNPs), 12.2% somatic single nucleotide variants (somatic SNVs) and 6.8% deletion. Out of 267 missense variants, 140 are deleterious mutations. We identified 284 non-coding variants with 94% regulatory in nature.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

Background The Colorado potato beetle, Leptinotarsa decemlineata Say, is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to not only rapidly adapt to a broad range of solanaceaeous plants and variable climates during its global invasion, but, most notably, to rapidly evolve resistance to insecticides (over 50 different compounds in all major classes, in some cases within the first year of use). To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species, using whole-genome sequencing, transcriptome sequencing, and a large community-driven annotation effort. Results We present a 140x coverage whole genome sequence from a single female L. decemlineata, with a reference gene set of 24,740 genes. Transposable elements comprise at least 17% of the genome, and are heavily represented in an analysis of rapidly evolving gene families compared to other Coleoptera. Population genetic analyses provide evidence of high levels of nucleotide diversity, local geographic structure, and recent population growth in pest populations, pointing to the availability of considerable standing genetic variation. These factors may play an important role in rapid evolutionary change. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes (e.g. cysteine peptidase genes) in gut tissues, as well as expansions of the gustatory receptors for bitter tasting plants in the nightshade family, Solanaceae. Despite its notoriety for adapting to insecticides, L. decemlineata has a similar suite of genes involved in resistance (metabolic detoxification and cuticle penetration) compared to other beetles, although expansions in specific cytochrome P450 subfamilies are known to be associated with insecticide resistance. Finally, this beetle has interesting duplications in RNAi genes that might be linked to its high sensitivity to RNAi and could be important in the future development of gene targeted pesticides. Conclusions As a representative of one of the most evolutionarily diverse lineages, the L. decemlineata genome will undoubtedly provide new opportunities for deeper understanding on the ecology, evolution, and management of this species, as well as new opportunities to leverage genomic technologies to understand the basis of a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.