Complete Mitochondrial Genome of Concholepas Concholepas Inferred by 454 Pyrosequencing and Mtdna Expression in Two Mollusc Populations
Total Page:16
File Type:pdf, Size:1020Kb
Comparative Biochemistry and Physiology, Part D 8 (2013) 17–23 Contents lists available at SciVerse ScienceDirect Comparative Biochemistry and Physiology, Part D journal homepage: www.elsevier.com/locate/cbpd Complete mitochondrial genome of Concholepas concholepas inferred by 454 pyrosequencing and mtDNA expression in two mollusc populations Gustavo Núñez-Acuña, Andrea Aguilar-Espinoza, Cristian Gallardo-Escárate ⁎ Laboratory of Biotechnology and Aquatic Genomics, Interdisciplinary Center for Aquaculture Research, University of Concepción, PO Box 160-C, Concepción, Chile article info abstract Article history: Despite the great relevance of mitochondrial genome analysis in evolutionary studies, there is scarce information Received 2 August 2012 on how the transcripts associated with the mitogenome are expressed and their role in the genetic structuring of Received in revised form 20 October 2012 populations. This work reports the complete mitochondrial genome of the marine gastropod Concholepas Accepted 21 October 2012 concholepas, obtained by 454 pryosequencing, and an analysis of mitochondrial transcripts of two populations Available online 27 October 2012 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 anal- Keywords: fi Mitochondrial genome ysis of mitochondrial gene variability showed signi cant differences among populations. In terms of levels of rel- Concholepas concholepas ative abundance of transcripts associated with mitochondrion in the two populations (assessed by qPCR), the RNA-Seq genes associated with complexes III and IV of the mitochondrial genome had the highest levels of expression Gene 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 struc- ture and mRNA expression levels of C. concholepas. © 2012 Elsevier Inc. All rights reserved. 1. Introduction reported that there is a noncoding region in several species that can play a regulatory role in the duplication and transcription of the mito- In general, the genomic organization of the mtDNA in marine organ- chondrial genome, although it is difficult to study its homology in higher isms have been studied through long-PCR and random genome sequenc- taxa because it presents similarities only among closely related taxa ing with DNA nebulization and the creation of clone libraries with pUC or (Passamonti et al., 2011). Complete mitochondrial genome analysis in- Bluescript vectors (Burger et al., 2007). However, recent advances in creases the resolution of works on the molecular phylogeny of species, next-generation sequencing technology (NGS) have increased the veloc- allowing the identification of genetic reordering among the higher taxa ity of mitochondrial genome characterization of several vertebrate and of species (Rawlings et al., 2010). Mitochondrial genome analysis signifi- invertebrate marine species (Jex et al., 2008; Neira-Oviedo et al., 2011). cantly increases the resolution of philogenetic studies of lower taxa, com- The main advantage of using NGS platforms is the possibility of using bining the data of all the genes in a single set (Cunha et al., 2009). Finally, the data generated to identify specific genes and to infer its relative intraspecific analysis of the mitochondrial genome evidences possible abundance under different conditions of study (Jex et al., 2008; Neira- mtDNA polymorphism (single nucleotide polymorphisms, SNPs) associ- Oviedo et al., 2011), which is of particular interest in comparative geno- ated with different geographic locations (Xin et al., 2011). Nevertheless, mic ecology studies. to our knowledge there have been no studies relating mtDNA polymor- The mitochondrial genome is a macromolecule with a size that in the phism to patterns of mitochondrial gene expression at the biogeographic majority of metazoans ranges between 15 and 20 kb (Boyce et al., 1989; level. In this context, important questions remain unanswered, such as: Is Boore, 1999). In practically all metazoans the mitochondrion is circular, it possible to find local patterns of mitochondrial gene expression? Does except in some classes of cnydaria that have lineal mitochondrial DNA variation of mitochondrial gene expression imply local adaptations? chromosomes (Bridge et al., 1992). The mitochondrial genome generally Our model of study, Concholepas concholepas, is a marine gastropod has 13 enzymatic protein subunits involved in the process of oxidative found along the Chilean coast. It has become a commercially important phosphorylation, 22 tRNA and 2 ribosomal subunits. As well, it has been benthic resource for Chile (Leiva and Castilla, 2001). One of the main characteristics of these species at the ecological level is their long larval development stage in the water column, which can last from almost ⁎ Corresponding author at: Interdisciplinary Center for Aquaculture Research, Univer- three months to one year. This extensive stage of planktonic develop- sity of Concepción, PO Box 160-C, Concepción, Chile. Tel.: +56 41 2203422; fax: +56 41 2207310. ment results in a broad geographic distribution of the species from 6°S E-mail address: [email protected] (C. Gallardo-Escárate). to 56°S. Genetic studies based on the coding region of the cytochrome 1744-117X/$ – see front matter © 2012 Elsevier Inc. All rights reserved. http://dx.doi.org/10.1016/j.cbd.2012.10.004 18 G. Núñez-Acuña et al. / Comparative Biochemistry and Physiology, Part D 8 (2013) 17–23 oxidase I gene (COX1) do not show genetic-population structuring in accordance with the supplier's instructions. The RNA was purified among different C. concholepas groups over the range of their biogeo- and the concentration and purity were measured with a spectropho- graphic distribution (Cárdenas et al., 2009). Nevertheless, an unan- tometer (ND-1000, Nanodrop Technologies) and the integrity was swered question is whether there are genes in C. concholepas that visualized with electrophoresis in MOPS/formaldehyde agarose gels at show highly differentiated biogeographic patterns due to environmen- 1.2% staining with ethidium bromide at 0.001%, and then quantified tal factors (e.g. northern and southern Chile). New DNA molecular with a 2100 bioanalyzer (Agilent Technologies, USA). Selected extrac- markers have recently been described for C. concholepas as micro- tions presented 260/280 and 260/230 purity indices equal to or greater satellites (Cárdenas et al., 2006) and ESTs obtained by 454 pyro- than 2.0 and showed integral RNA in electrophoresis and bioanalyzer sequencing (Cárdenas et al., 2011). As well, there are now candidate measurements (RIN>8). Subsequently, sample pools of selected RNAs genes related to thermal, oxidative and immune responses that allow were prepared (one pool for each population group) and precipitated for testing hypotheses about local adaptation at the genomic level. overnight with 2× volumes of absolute ethanol and 0.1× volume of so- The objective of the present study was to describe for the first dium acetate 0.3 M at −80 °C. Double-stranded cDNA was synthesized time the mitochondrial genome for C. concholepas using 454 cDNA from the extracted RNA and pyrosequenced (full plate) in a 454 GS FLX pyrosequencing (RNA-seq). As well, we evaluated the levels of titanium platform by Macrogen Inc. (S. Korea). gene expression (in silico analysis) of all the mitochondrial genes in The raw data for both samples were trimmed and assembled sepa- two populations of C. concholepas (1000 km apart) based on data rately, and combined using CLC Genomics Workbench software (CLC from RNA-seq. Differences in gene expression were confirmed with bio, Denmark). The overlap settings used for these assemblies were real time qPCR analysis. mismatch cost=2, insert cost=3, minimum contig length=200 bp, similarity=0.8, and trimming quality score=0.05. The redundancy of 2. Material and methods singleton sequences was calculated using the Java application Duplicate Finder (http://bioinfotutlets.blogspot.com/). All contigs and singletons 2.1. Samples were annotated according to Gene Ontology terms using Blast2Go soft- ware (Conesa et al., 2005). From the annotated contigs, we selected Fifty adult individuals of C. concholepas (Gastropoda, Muricidae) were those that presented homology with genes present in the mitochondrial collected as samples from two locations: Caleta Chungungo, Punta de genome, and we assembled them using the mitochondrial genome of Choros (29°16′S–71°31′W) and Los Molinos, Valdivia (39°40′S–73°12′ the muricid Thais clavigera (GenBank accession no. DQ159954.1) as a ref- W) in February 2010. Shell length of samples were between 10.2 and erence (Ki et al., 2010). As well, for regions of the mitochondrial genome 12.7 cm. Approximately 100 mg of gill and muscle tissue was extracted that were not covered by EST, specific primers were designed and ampli- from the samples and conserved in 1 mL of RNAlater stabilization solu- fied by PCR and Sanger sequencing (Macrogen Inc., S. Korea) of the tion (Ambion®, USA) and stored at −80 °C until total RNA extraction. amplicons (Table