Dataset Information

Mapping 3' transcript ends in the bank vole (Clethrionomys glareolus) mitochondrial genome with RNA-Seq.

ABSTRACT: Although posttranscriptional modification of mitochondrial (mt) transcripts plays key roles in completion of the coding information and in the expression of mtDNA-encoded genes, there is little experimental evidence on the polyadenylation status and the location of mt gene poly(A) sites for non-human mammals.Poly(A)-enriched RNA-Seq reads collected for two wild-caught bank voles (Clethrionomys glareolus) were mapped to the complete mitochondrial genome of that species. Transcript polyadenylation was detected as unmapped adenine residues at the ends of the mapped reads. Where the tRNA punctuation model applied, there was the expected polyadenylation, except for the nad5 transcript, whose polyadenylated 3' end is at an intergenic sequence/cytochrome b boundary. As in human, two pairs of bank vole genes, nad4l/nad4 and atp8/atp6, are expressed from bicistronic transcripts. TAA stop codons of four bank vole protein-coding genes (nad1, atp6, cox3 and nad4) are incompletely encoded in the DNA and are completed by polyadenylation. This is three genes (nad2, nad3 and cob) less than in human. The bank vole nad2 gene encodes a full stop codon (TAA in one vole and TAG in the other), which is followed by a 2 bp UTR and the gene conforms to the tRNA punctuation model. In contrast, the annotations of the reference mouse and some other rodent mt genomes in GenBank include complete TAG stop codons in both nad1 and nad2, which overlap downstream trnI and trnW, respectively. Thus the RNA-Seq data of bank voles provides a model for stop codons of mt-encoded genes in mammals comparable to humans, but at odds with some of the interpretation based purely on genomic data in mouse and other rodents.This work demonstrates how RNA-Seq data were useful to recover mtDNA transcriptome data in a non-model rodent and to shed more light on mammalian mtDNA transcriptome and post-transcriptional modification. Even though gene content and organisation of mtDNA are strongly conserved among mammals, annotations that neglect the transcriptome may be prone to errors in relation to the stop codons.

SUBMITTER: Markova S

PROVIDER: S-EPMC4624183 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Mapping 3' transcript ends in the bank vole (Clethrionomys glareolus) mitochondrial genome with RNA-Seq.

Marková Silvia S Filipi Karolína K Searle Jeremy B JB Kotlík Petr P

BMC genomics 20151026

<h4>Background</h4>Although posttranscriptional modification of mitochondrial (mt) transcripts plays key roles in completion of the coding information and in the expression of mtDNA-encoded genes, there is little experimental evidence on the polyadenylation status and the location of mt gene poly(A) sites for non-human mammals.<h4>Results</h4>Poly(A)-enriched RNA-Seq reads collected for two wild-caught bank voles (Clethrionomys glareolus) were mapped to the complete mitochondrial genome of that ...[more]

PMID: 26503603

Similar Datasets

Project description:BACKGROUND: Chemical communication in mammals involves globular lipocalins that protect and transport pheromones during their passage out of the body. Efficient communication via this protein - pheromone complex is essential for triggering multiple responses including aggression, mate choice, copulatory behaviour, and onset and synchronization of oestrus. The roles of lipocalins in communication were studied in many organisms and especially in mice (i.e. Mus musculus domesticus) which excrete Major Urinary Proteins (Mup) in excessive amounts in saliva and urine. Other mammals, however, often lack the genes for Mups or their expression is very low. Therefore, we aimed at characterization of candidate lipocalins in Myodes glareolus which are potentially linked to chemical communication. One of them is Aphrodisin which is a unique lipocalin that was previously described from hamster vaginal discharge and is known to carry pheromones stimulating copulatory behaviour in males. RESULTS: Here we show that Aphrodisin-like proteins exist in other species, belong to a group of Odorant Binding Proteins (Obp), and contrary to the expression of Aphrodisin only in hamster genital tract and parotid glands of females, we have detected these transcripts in both sexes of M. glareolus with the expression confirmed in various tissues including prostate, prepucial and salivary glands, liver and uterus. On the level of mRNA, we have detected three different gene variants. To assess their relevance for chemical communication we investigated the occurrence of particular proteins in saliva, urine and vaginal discharge. On the protein level we confirmed the presence of Obp2 and Obp3 in both saliva and urine. Appropriate bands in the range of 17-20 kDa from vaginal discharge were, however, beyond the MS detection limits. CONCLUSION: Our results demonstrate that three novel Obps (Obp1, Obp2, and Obp3) are predominant lipocalins in Myodes urine and saliva. On the protein level we have detected further variants and thus we assume that similarly as Major Urinary Proteins in mice, these proteins may be important in chemical communication in this Cricetid rodent.

Project description:Species-level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species-wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid-Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.

Project description:BackgroundUnderstanding the genetic basis of adaptive changes has been a major goal of evolutionary biology. In complex organisms without sequenced genomes, de novo transcriptome assembly using a longer read sequencing technology followed by expression profiling using short reads is likely to provide comprehensive identification of adaptive variation at the expression level and sequence polymorphisms in coding regions. We performed sequencing and de novo assembly of the bank vole heart transcriptome in lines selected for high metabolism and unselected controls.ResultsA single 454 Titanium run produced over million reads, which were assembled into 63,581 contigs. Searches against the SwissProt protein database and the ENSEMBL collection of mouse transcripts detected similarity to 11,181 and 14,051 genes, respectively. As judged by the representation of genes from the heart-related Gene Ontology categories and UniGenes detected in the mouse heart, our detection of the genes expressed in the heart was nearly complete (> 95% and almost 90% respectively). On average, 38.7% of the transcript length was covered by our sequences, with notably higher (45.0%) coverage of coding regions than of untranslated regions (24.5% of 5' and 32.7% of 3'UTRs). Lower sequence conservation between mouse and bank vole in untranslated regions was found to be partially responsible for poorer UTR representation. Our data might suggest a widespread transcription from noncoding genomic regions, a finding not reported in previous studies regarding transcriptomes in non-model organisms. We also identified over 19 thousand putative single nucleotide polymorphisms (SNPs). A much higher fraction of the SNPs than expected by chance exhibited variant frequency differences between selection regimes.ConclusionLonger reads and higher sequence yield per run provided by the 454 Titanium technology in comparison to earlier generations of pyrosequencing proved beneficial for the quality of assembly. An almost full representation of genes known to be expressed in the mouse heart was identified. Usage of the extensive genomic resources available for the house mouse, a moderately (20-40 mln years) divergent relative of the voles, enabled a comprehensive assessment of the transcript completeness. Transcript sequences generated in the present study allowed the identification of candidate SNPs associated with divergence of selection lines and constitute a valuable permanent resource forming a foundation for RNAseq experiments aiming at detection of adaptive changes both at the level of gene expression and sequence variants, that would facilitate studies of the genetic basis of evolutionary divergence.

Dataset Information

Mapping 3' transcript ends in the bank vole (Clethrionomys glareolus) mitochondrial genome with RNA-Seq.

Publications

Mapping 3' transcript ends in the bank vole (Clethrionomys glareolus) mitochondrial genome with RNA-Seq.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets