The RNA-Seq approach to studying the expression of mosquito mitochondrial genes

MITGARD: an automated pipeline for mitochondrial genome assembly in eukaryotic species using RNA-seq data

MOTIVATION

Over the past decade, the field of next-generation sequencing (NGS) has seen dramatic advances in methods and a decrease in costs. Consequently, a large expansion of data has been generated by NGS, most of which have originated from RNA-sequencing (RNA-seq) experiments. Because mitochondrial genes are expressed in most eukaryotic cells, mitochondrial mRNA sequences are usually co-sequenced within the target transcriptome, generating data that are commonly underused or discarded. Here, we present MITGARD, an automated pipeline that reliably recovers the mitochondrial genome from RNA-seq data from various sources. The pipeline identifies mitochondrial sequence reads based on a phylogenetically related reference, assembles them into contigs, and extracts a complete mtDNA for the target species.

RESULTS

We demonstrate that MITGARD can reconstruct the mitochondrial genomes of several species throughout the tree of life. We noticed that MITGARD can recover the mitogenomes in different sequencing schemes and even in a scenario of low-sequencing depth. Moreover, we showed that the use of references from congeneric species diverging up to 30 million years ago (MYA) from the target species is sufficient to recover the entire mitogenome, whereas the use of species diverging between 30 and 60 MYA allows the recovery of most mitochondrial genes. Additionally, we provide a case study with original data in which we estimate a phylogenetic tree of snakes from the genus Bothrops, further demonstrating that MITGARD is suitable for use on biodiversity projects. MITGARD is then a valuable tool to obtain high-quality information for studies focusing on the phylogenetic and evolutionary aspects of eukaryotes and provides data for easily identifying a sample using barcoding, and to check for cross-contamination using third-party tools.

Genomes, expression profiles, and diversity of mitochondria of the White-footed Deermouse Peromyscus leucopus, reservoir of Lyme disease and other zoonoses

The cricetine rodents Peromyscus leucopus and P. maniculatus are key reservoirs for several zoonotic diseases in North America. We determined the complete circular mitochondrial genome sequences of representatives of 3 different stock colonies of P. leucopus, one stock colony of P. maniculatus and two wild populations of P. leucopus. The genomes were syntenic with that of the murids Mus musculus and Rattus norvegicus. Phylogenetic analysis confirmed that these two Peromyscus species are sister taxa in a clade with P. polionotus and also uncovered a distinction between P. leucopus populations in the eastern and the central United States. In one P. leucopus lineage four extended regions of mitochondrial pseudogenes were identified in the nuclear genome. RNA-seq analysis revealed transcription of the entire genome and differences from controls in the expression profiles of mitochondrial genes in the blood, but not in liver or brain, of animals infected with the zoonotic pathogen Borrelia hermsii. PCR and sequencing of the D-loop of the mitochondrion identified 32 different haplotypes among 118 wild P. leucopus at a Connecticut field site. These findings help to further establish P. leucopus as a model organism for studies of emerging infectious diseases, ecology, and in other disciplines.

Mitochondrial genome characterization of Melipona bicolor: Insights from the control region and gene expression data

The stingless bee Melipona bicolor is the only bee in which true polygyny occurs. Its mitochondrial genome was first sequenced in 2008, but it was incomplete and no information about its transcription was known. We combined short and long reads of M. bicolor DNA with RNASeq data to obtain insights about mitochondrial evolution and gene expression in bees. The complete genome has 15,001 bp, including a control region of 255 bp that contains all conserved structures described in honeybees with the highest AT content reported so far for bees (98.1%), displaying a compact but functional region. Gene expression control is similar to other insects however unusual patterns of expression may suggest the existence of different isoforms for the mitochondrially encoded 12S rRNA. Results reveal unique and shared features of the mitochondrial genome in terms of sequence evolution and gene expression making M. bicolor an interesting model to study mitochondrial genomic evolution.

Variation and evolution of polyadenylation profiles in sauropsid mitochondrial mRNAs as deduced from the high-throughput RNA sequencing

BACKGROUND

Genes encoded in vertebrate mitochondrial DNAs are transcribed as a polycistronic transcript for both strands, which is later processed into individual mRNAs, rRNAs and tRNAs, followed by modifications, such as polyadenylation at the 3' end of mRNAs. Although mechanisms of the mitochondrial transcription and RNA processing have been extensively studied using some model organisms, structural variability of mitochondrial mRNAs across different groups of vertebrates is poorly understood. We conducted the high-throughput RNA sequencing to identify major polyadenylation sites for mitochondrial mRNAs in the Japanese grass lizard, Takydromus tachydromoides and compared the polyadenylation profiles with those identified similarly for 23 tetrapod species, featuring sauropsid taxa (reptiles and birds).

RESULTS

As compared to the human, a major polyadenylation site for the NADH dehydrogenase subunit 5 mRNA of the grass lizard was located much closer to its stop codon, resulting in considerable truncation of the 3' untranslated region for the mRNA. Among the other sauropsid taxa, several distinct polyadenylation profiles from the human counterpart were found for different mRNAs. They included various truncations of the 3' untranslated region for NADH dehydrogenase subunit 5 mRNA in four taxa, bird-specific polyadenylation of the light-strand-transcribed NADH dehydrogenase subunit 6 mRNA, and the combination of the ATP synthase subunit 8/6 mRNA with a neighboring mRNA into a tricistronic mRNA in the side-necked turtle Pelusios castaneus. In the last case of P. castaneus, as well as another example for NADH dehydrogenase subunit 1 mRNAs of some birds, the association between the polyadenylation site change and the gene overlap was highlighted. The variations in the polyadenylation profile were suggested to have arisen repeatedly in diverse sauropsid lineages. Some of them likely occurred in response to gene rearrangements in the mitochondrial DNA but the others not.

CONCLUSIONS

These results demonstrate structural variability of mitochondrial mRNAs in sauropsids. The efficient and comprehensive characterization of the mitochondrial mRNAs will contribute to broaden our understanding of their structural and functional evolution.

De novo sequencing and transcriptome analysis of venom glands of endoparasitoid Aenasius arizonensis (Girault) (=Aenasius bambawalei Hayat) (Hymenoptera, Encyrtidae)

Aenasius bambawalei Hayat (Encyrtidae: Hymenoptera) has been synonymized with Aenasius arizonensis (Girault) is a small, newly discovered endoparasitoid of the cotton mealybug Phenacoccuss solenopsis Tinsley (Pseudococcidae: Hemiptera), which completes its life cycle inside the body of its host and it is a potential insect control tool. Despite the acquired knowledge regarding host-parasitoid interaction, little information is available on the factors of parasitoid origin able to modulate mealybug physiology. The components of A. arizonensis venom have not been well studied but venom from other parasitoids and wasps contain biologically active proteins that have potential applications in pest management or may be of medicinal importance. To provide an insight into the transcripts expressed in the venom gland of A. arizonensis, a transcriptomic database was developed utilizing high throughput RNA sequencing approaches to analyze the genes expressed in venom glands of this endoparasitic wasp. The resulting A. arizonensis RNA sequences were assembled de-novo with contigs then blasted against the NCBI non-redundant sequence database. Contigs which matched database sequences were mostly homologous to genes from hymenopteran parasitoids such as Nasonia vitripennis, Copidosoma floridanum, Fopius arsenus and Pteromalas puparium. Further analysis of the A. arizonensis database was then performed which focused on selected genes encoding proteins potentially involved in host developmental arrest, disrupting the host immune system, host paralysis, and transcripts that support these functions. Sequenced mRNAS predicted to encode full length ORFs of Calreticulin, Serine Protease Precursor and Arginine kinase proteins were identified and the tissue specific expression of these putative venom genes was analyzed by RT-PCR. In addition, results also demonstrate that de novo transcriptome assembly allows useful venom gene expression analysis in a species lacking a genome sequence database and may provide useful information for devising control tools for insect pests and other applications.

RNA-Seq analysis using de novo transcriptome assembly as a reference for the salmon louse Caligus rogercresseyi

Despite the economic and environmental impacts that sea lice infestations have on salmon farming worldwide, genomic data generated by high-throughput transcriptome sequencing for different developmental stages, sexes, and strains of sea lice is still limited or unknown. In this study, RNA-seq analysis was performed using de novo transcriptome assembly as a reference for evidenced transcriptional changes from six developmental stages of the salmon louse Caligus rogercresseyi. EST-datasets were generated from the nauplius I, nauplius II, copepodid and chalimus stages and from female and male adults using MiSeq Illumina sequencing. A total of 151,788,682 transcripts were yielded, which were assembled into 83,444 high quality contigs and subsequently annotated into roughly 24,000 genes based on known proteins. To identify differential transcription patterns among salmon louse stages, cluster analyses were performed using normalized gene expression values. Herein, four clusters were differentially expressed between nauplius I–II and copepodid stages (604 transcripts), five clusters between copepodid and chalimus stages (2,426 transcripts), and six clusters between female and male adults (2,478 transcripts). Gene ontology analysis revealed that the nauplius I–II, copepodid and chalimus stages are mainly annotated to aminoacid transfer/repair/breakdown, metabolism, molting cycle, and nervous system development. Additionally, genes showing differential transcription in female and male adults were highly related to cytoskeletal and contractile elements, reproduction, cell development, morphogenesis, and transcription-translation processes. The data presented in this study provides the most comprehensive transcriptome resource available for C. rogercresseyi, which should be used for future genomic studies linked to host-parasite interactions.

The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector

Mosquitoes are vectors of a number of important human and animal diseases. The development of novel vector control strategies requires a thorough understanding of mosquito biology. To facilitate this, we used RNA-seq to identify novel genes and provide the first high-resolution view of the transcriptome throughout development and in response to blood feeding in a mosquito vector of human disease, Aedes aegypti, the primary vector for Dengue and yellow fever. We characterized mRNA expression at 34 distinct time points throughout Aedes development, including adult somatic and germline tissues, by using polyA+ RNA-seq. We identify a total of 14,238 novel new transcribed regions corresponding to 12,597 new loci, as well as many novel transcript isoforms of previously annotated genes. Altogether these results increase the annotated fraction of the transcribed genome into long polyA+ RNAs by more than twofold. We also identified a number of patterns of shared gene expression, as well as genes and/or exons expressed sex-specifically or sex-differentially. Expression profiles of small RNAs in ovaries, early embryos, testes, and adult male and female somatic tissues also were determined, resulting in the identification of 38 new Aedes-specific miRNAs, and ~291,000 small RNA new transcribed regions, many of which are likely to be endogenous small-interfering RNAs and Piwi-interacting RNAs. Genes of potential interest for transgene-based vector control strategies also are highlighted. Our data have been incorporated into a user-friendly genome browser located at www.Aedes.caltech.edu, with relevant links to Vectorbase (www.vectorbase.org)

Global analyses of Ceratocystis cacaofunesta mitochondria: from genome to proteome

Background

The ascomycete fungus Ceratocystis cacaofunesta is the causal agent of wilt disease in cacao, which results in significant economic losses in the affected producing areas. Despite the economic importance of the Ceratocystis complex of species, no genomic data are available for any of its members. Given that mitochondria play important roles in fungal virulence and the susceptibility/resistance of fungi to fungicides, we performed the first functional analysis of this organelle in Ceratocystis using integrated “omics” approaches.

Results

The C. cacaofunesta mitochondrial genome (mtDNA) consists of a single, 103,147-bp circular molecule, making this the second largest mtDNA among the Sordariomycetes. Bioinformatics analysis revealed the presence of 15 conserved genes and 37 intronic open reading frames in C. cacaofunesta mtDNA. Here, we predicted the mitochondrial proteome (mtProt) of C. cacaofunesta, which is comprised of 1,124 polypeptides - 52 proteins that are mitochondrially encoded and 1,072 that are nuclearly encoded. Transcriptome analysis revealed 33 probable novel genes. Comparisons among the Gene Ontology results of the predicted mtProt of C. cacaofunesta, Neurospora crassa and Saccharomyces cerevisiae revealed no significant differences. Moreover, C. cacaofunesta mitochondria were isolated, and the mtProt was subjected to mass spectrometric analysis. The experimental proteome validated 27% of the predicted mtProt. Our results confirmed the existence of 110 hypothetical proteins and 7 novel proteins of which 83 and 1, respectively, had putative mitochondrial localization.

Conclusions

The present study provides the first partial genomic analysis of a species of the Ceratocystis genus and the first predicted mitochondrial protein inventory of a phytopathogenic fungus. In addition to the known mitochondrial role in pathogenicity, our results demonstrated that the global function analysis of this organelle is similar in pathogenic and non-pathogenic fungi, suggesting that its relevance in the lifestyle of these organisms should be based on a small number of specific proteins and/or with respect to differential gene regulation. In this regard, particular interest should be directed towards mitochondrial proteins with unknown function and the novel protein that might be specific to this species. Further functional characterization of these proteins could enhance our understanding of the role of mitochondria in phytopathogenicity.

Complete mitochondrial genome of Concholepas concholepas inferred by 454 pyrosequencing and mtDNA expression in two mollusc populations

Despite the great relevance of mitochondrial genome analysis in evolutionary studies, there is scarce information on how the transcripts associated with the mitogenome are expressed and their role in the genetic structuring of populations. This work reports the complete mitochondrial genome of the marine gastropod Concholepas concholepas, obtained by 454 pryosequencing, and an analysis of mitochondrial transcripts of two populations 1000 km apart along the Chilean coast. The mitochondrion of C. concholepas is 15,495 base pairs (bp) in size and contains the 37 subunits characteristic of metazoans, as well as a non-coding region of 330 bp. In silico analysis of mitochondrial gene variability showed significant differences among populations. In terms of levels of relative abundance of transcripts associated with mitochondrion in the two populations (assessed by qPCR), the genes associated with complexes III and IV of the mitochondrial genome had the highest levels of expression in the northern population while transcripts associated with the ATP synthase complex had the highest levels of expression in the southern population. Moreover, fifteen polymorphic SNPs were identified in silico between the mitogenomes of the two populations. Four of these markers implied different amino acid substitutions (non-synonymous SNPs). This work contributes novel information regarding the mitochondrial genome structure and mRNA expression levels of C. concholepas.

Comparative transcriptome analyses of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from Kenya by RNA-Seq

Malaria causes more than 300 million clinical cases and 665,000 deaths each year, and the majority of the mortality and morbidity occurs in sub-Saharan Africa. Due to the lack of effective vaccines and wide-spread resistance to antimalarial drugs, mosquito control is the primary method of malaria prevention and control. Currently, malaria vector control relies on the use of insecticides, primarily pyrethroids. The extensive use of insecticides has imposed strong selection pressures for resistance in the mosquito populations. Consequently, resistance to pyrethroids in Anopheles gambiae, the main malaria vector in sub-Saharan Africa, has become a major obstacle for malaria control. A key element of resistance management is the identification of resistance mechanisms and subsequent development of reliable resistance monitoring tools. Field-derived An. gambiae from Western Kenya were phenotyped as deltamethrin-resistant or -susceptible by the standard WHO tube test, and their expression profile compared by RNA-seq. Based on the current annotation of the An. gambiae genome, a total of 1,093 transcripts were detected as significantly differentially accumulated between deltamethrin-resistant and -susceptible mosquitoes. These transcripts are distributed over the entire genome, with a large number mapping in QTLs previously linked to pyrethorid resistance, and correspond to heat-shock proteins, metabolic and transport functions, signal transduction activities, cytoskeleton and others. The detected differences in transcript accumulation levels between resistant and susceptible mosquitoes reflect transcripts directly or indirectly correlated with pyrethroid resistance. RNA-seq data also were used to perform a de-novo Cufflinks assembly of the An. gambiae genome.

Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration

Background

The naid annelids contain a number of species that vary in their ability to regenerate lost body parts, making them excellent candidates for evolution of regeneration studies. However, scant sequence data exists to facilitate such studies. We constructed a cDNA library from the naid Pristina leidyi, a species that is highly regenerative and also reproduces asexually by fission, using material from a range of regeneration and fission stages for our library. We then sequenced the transcriptome of P. leidyi using 454 technology.

Results

454 sequencing produced 1,550,174 reads with an average read length of 376 nucleotides. Assembly of 454 sequence reads resulted in 64,522 isogroups and 46,679 singletons for a total of 111,201 unigenes in this transcriptome. We estimate that over 95% of the transcripts in our library are present in our transcriptome. 17.7% of isogroups had significant BLAST hits to the UniProt database and these include putative homologs of a number of genes relevant to regeneration research. Although many sequences are incomplete, the mean sequence length of transcripts (isotigs) is 707 nucleotides. Thus, many sequences are large enough to be immediately useful for downstream applications such as gene expression analyses. Using in situ hybridization, we show that two Wnt/β-catenin pathway genes (homologs of frizzled and β-catenin) present in our transcriptome are expressed in the regeneration blastema of P. leidyi, demonstrating the usefulness of this resource for regeneration research.

Conclusions

454 sequencing is a rapid and efficient approach for identifying large numbers of genes in an organism that lacks a sequenced genome. This transcriptome dataset will be a valuable resource for molecular analyses of regeneration in P. leidyi and will serve as a starting point for comparisons to non-regenerating naids. It also contributes significantly to the still limited genomic resources available for annelids and lophotrochozoans more generally.

RNA sequencing reveals differential expression of mitochondrial and oxidation reduction genes in the long-lived naked mole-rat when compared to mice

The naked mole-rat (Heterocephalus glaber) is a long-lived, cancer resistant rodent and there is a great interest in identifying the adaptations responsible for these and other of its unique traits. We employed RNA sequencing to compare liver gene expression profiles between naked mole-rats and wild-derived mice. Our results indicate that genes associated with oxidoreduction and mitochondria were expressed at higher relative levels in naked mole-rats. The largest effect is nearly 300-fold higher expression of epithelial cell adhesion molecule (Epcam), a tumour-associated protein. Also of interest are the protease inhibitor, alpha2-macroglobulin (A2m), and the mitochondrial complex II subunit Sdhc, both ageing-related genes found strongly over-expressed in the naked mole-rat. These results hint at possible candidates for specifying species differences in ageing and cancer, and in particular suggest complex alterations in mitochondrial and oxidation reduction pathways in the naked mole-rat. Our differential gene expression analysis obviated the need for a reference naked mole-rat genome by employing a combination of Illumina/Solexa and 454 platforms for transcriptome sequencing and assembling transcriptome contigs of the non-sequenced species. Overall, our work provides new research foci and methods for studying the naked mole-rat's fascinating characteristics.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

Targeted gene expression in the transgenic Aedes aegypti using the binary Gal4-UAS system

In this study, we report the establishment of the binary Gal4/UAS system for the yellow fever mosquito Aedes aegypti. We utilized the 1.8-kb 5' upstream region of the vitellogenin gene (Vg) to genetically engineer mosquito lines with the Vg-Gal4 activator and established UAS-EGFP responder transgenic mosquito lines to evaluate the binary Gal4/UAS system. The results show that the Vg-Gal4 driver leads to a high level of tissue-, stage- and sex-specific expression of the EGFP reporter in the fat body of Vg-Gal4/UAS-EGFP hybrids after blood-meal activation. In addition, the applicability of this system to study hormonal regulation of gene expression was demonstrated in in vitro organ culture experiments in which the EGFP reporter was highly activated in isolated fat bodies of previtellogenic Vg-Gal4/UAS-EGFP females incubated in the presence of 20-hydroxyecdysone (20E). Hence, this study has opened the door for further refinement of genetic tools in mosquitoes.

Transcriptome profiling of chemosensory appendages in the malaria vector Anopheles gambiae reveals tissue- and sex-specific signatures of odor coding

Background

Chemosensory signal transduction guides the behavior of many insects, including Anopheles gambiae, the major vector for human malaria in sub-Saharan Africa. To better understand the molecular basis of mosquito chemosensation we have used whole transcriptome RNA sequencing (RNA-seq) to compare transcript expression profiles between the two major chemosensory tissues, the antennae and maxillary palps, of adult female and male An. gambiae.

Results

We compared chemosensory tissue transcriptomes to whole body transcriptomes of each sex to identify chemosensory enhanced genes. In the six data sets analyzed, we detected expression of nearly all known chemosensory genes and found them to be highly enriched in both olfactory tissues of males and females. While the maxillary palps of both sexes demonstrated strict chemosensory gene expression overlap, we observed acute differences in sensory specialization between male and female antennae. The relatively high expression levels of chemosensory genes in the female antennae reveal its role as an organ predominately assigned to chemosensation. Remarkably, the expression of these genes was highly conserved in the male antennae, but at much lower relative levels. Alternatively, consistent with a role in mating, the male antennae displayed significant enhancement of genes involved in audition, while the female enhancement of these genes was observed, but to a lesser degree.

Conclusions

These findings suggest that the chemoreceptive spectrum, as defined by gene expression profiles, is largely similar in female and male An. gambiae. However, assuming sensory receptor expression levels are correlated with sensitivity in each case, we posit that male and female antennae are perceptive to the same stimuli, but possess inverse receptive prioritizations and sensitivities. Here we have demonstrated the use of RNA-seq to characterize the sensory specializations of an important disease vector and grounded future studies investigating chemosensory processes.