Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

The role of mitochondria in sex- and age-specific gene expression in a species without sex chromosomes

Mitochondria perform an array of functions, many of which involve interactions with gene products encoded by the nucleus. These mitochondrial functions, particularly those involving energy production, can be expected to differ between sexes and across ages. Here, we measured mitochondrial effects on sex- and age-specific gene expression in parental and reciprocal F1 hybrids between allopatric populations of Tigriopus californicus with over 20% mitochondrial DNA divergence. Because the species lacks sex chromosomes, sex-biased mitochondrial effects are not confounded by the effects of sex chromosomes. Results revealed pervasive sex differences in mitochondrial effects, including effects on energetics and aging involving nuclear interactions throughout the genome. Using single-individual RNA sequencing, sex differences were found to explain more than 80% of the variance in gene expression. Males had higher expression of mitochondrial genes and mitochondrially targeted proteins (MTPs) involved in oxidative phosphorylation (OXPHOS), while females had elevated expression of non-OXPHOS MTPs, indicating strongly sex-dimorphic energy metabolism at the whole organism level. Comparison of reciprocal F1 hybrids allowed insights into the nature of mito-nuclear interactions, showing both mitochondrial effects on nuclear expression, and nuclear effects on mitochondrial expression. While based on a small set of crosses, sex-specific increases in mitochondrial expression with age were associated with longer life. Network analyses identified nuclear components of strong mito-nuclear interactions and found them to be sexually dimorphic. These results highlight the profound impact of mitochondria and mito-nuclear interactions on sex- and age-specific gene expression.

The role of mitochondria in sex- and age-specific gene expression in a species without sex chromosomes

Mitochondria perform an array of functions, many of which involve interactions with gene products encoded by the nucleus. These mitochondrial functions, particularly those involving energy production, can be expected to differ between sexes and across ages. Here we measured mitochondrial effects on sex- and age-specific gene expression in parental and reciprocal F1 hybrids between allopatric populations of Tigriopus californicus with over 20% mitochondrial DNA divergence. Because the species lacks sex chromosomes, sex-biased mitochondrial effects are not confounded by the effects of sex chromosomes. Using single-individual RNA sequencing, sex differences were found to explain more than 80% of the variance in gene expression. Males had higher expression of mitochondrial genes and mitochondrially targeted proteins (MTPs) involved in oxidative phosphorylation (OXPHOS), while females had elevated expression of non-OXPHOS MTPs, indicating strongly sex-dimorphic energy metabolism at the whole organism level. Comparison of reciprocal F1 hybrids allowed insights into the nature of mito-nuclear interactions, showing both mitochondrial effects on nuclear expression, as well as nuclear effects on mitochondrial expression. Across both sexes, increases in mitochondrial expression with age were associated with longer life. Network analyses identified nuclear components of strong mito-nuclear interactions, and found them to be sexually dimorphic. These results highlight the profound impact of mitochondria and mito-nuclear interactions on sex- and age-specific gene expression.

Disproportionate role of nuclear-encoded proteins in organismal and mitochondrial thermal performance in a copepod

Determining the mechanisms by which organisms evolve thermal tolerance is crucial to predicting how populations may respond to changes in local temperature regimes. Although evidence of relationships between mitochondrial background and thermal adaptation have been found, the presence of both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded proteins warrants experiments aimed at parsing out the relative role of each genome in thermal adaptation. We investigated the relative role of mtDNA-encoded products in thermal tolerance between two divergent populations of Tigriopus californicus using first-generation (F1) hybrids that vary in maternally inherited mtDNA but are heterozygous for population-specific alleles across nuclear loci. We tested two measures of thermal tolerance, (1) survivorship to acute thermal stress and (2) thermal stability of mitochondrial performance in Complex I-fueled ATP synthesis, both across a range of increasing temperatures. We found that the southern population (San Diego, CA, USA) outperformed the northern population (Strawberry Hill, OR, USA) in survivorship, and that both reciprocal F1 hybrid crosses had intermediate survival. Mitochondria from the San Diego population displayed greater stability in ATP synthesis with increasing temperatures compared with those from Strawberry Hill. Interestingly, hybrids from both cross directions had synthesis profiles that were very similar to that of Strawberry Hill. Taken together, these results suggest that the relative role of the mtDNA in these phenotypes is negligible compared with that of elements encoded by nuclear DNA in this system.

Quality of DNA extracted from freshwater fish scales and mucus and its application in genetic diversity studies of Perca fluviatilis and Rutilus rutilus

Studies on genetic diversity require biological material containing a reliable source of DNA that can be extracted and analyzed. Recently, non-invasive sampling has become a preferred sampling method of biological material. The suitability of a less invasive approach that involves obtaining samples by swabbing the fish skin (including live, non-anesthetized fish) should be considered. In this study, we compared the efficiency of DNA extraction, amplification, and sequencing of mtDNA fragments of two fish species Perca fluviatilis and Rutilus rutilus based on DNA collected from the scales and mucus using the modified Aljanabi and Martinez method. The results revealed a higher quality of DNA extracted from the mucus; however, the mean DNA concentration obtained from the scales of both fish species was higher. We verified the method suitable for amplification and sequencing of mtDNA fragments of both fish species using newly designed markers (D-loop, ATP6) and examined the potential risk of intraspecific cross-contamination. The DNA sequence alignment analysis revealed identical sequences attributed to the same individual when DNA, extracted from two different sources (scales and mucus), was used. We demonstrated that the quantity and quality of DNA extracted from the scales and mucus using the proposed method were high enough to carry out genetic diversity studies based on sampling of live fish with the possibility to release it after collecting samples.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Hypoxia tolerance, but not low pH tolerance, is associated with a latitudinal cline across populations of Tigriopus californicus

Intertidal organisms must tolerate daily fluctuations in environmental parameters, and repeated exposure to co-occurring conditions may result in tolerance to multiple stressors correlating. The intertidal copepod Tigriopus californicus experiences diurnal variation in dissolved oxygen levels and pH as the opposing processes of photosynthesis and cellular respiration lead to coordinated highs during the day and lows at night. While environmental parameters with overlapping spatial gradients frequently result in correlated traits, less attention has been given to exploring temporally correlated stressors. We investigated whether hypoxia tolerance correlates with low pH tolerance by separately testing the hypoxia and low pH stress tolerance separately of 6 genetically differentiated populations of T. californicus. We independently checked for similarities in tolerance for each of the two stressors by latitude, sex, size, and time since collection as predictors. We found that although hypoxia tolerance correlated with latitude, low pH tolerance did not, and no predictor was significant for both stressors. We concluded that temporally coordinated exposure to low pH and low oxygen did not result in populations developing equivalent tolerance for both. Although climate change alters several environmental variables simultaneously, organisms' abilities to tolerate these changes may not be similarly coupled.

Extreme mitochondrial DNA divergence underlies genetic conflict over sex determination

Cytoplasmic male sterility (CMS) is a form of genetic conflict over sex determination that results from differences in modes of inheritance between genomic compartments.^1-3 Indeed, maternally transmitted (usually mitochondrial) genes sometimes enhance their transmission by suppressing the male function in a hermaphroditic organism to the detriment of biparentally inherited nuclear genes. Therefore, these hermaphrodites become functionally female and may coexist with regular hermaphrodites in so-called gynodioecious populations.³ CMS has been known in plants since Darwin's times⁴ but is previously unknown in the animal kingdom.^5-8 We relate the first observation of CMS in animals. It occurs in a freshwater snail population, where some individuals appear unable to sire offspring in controlled crosses and show anatomical, physiological, and behavioral characters consistent with a suppression of the male function. Male sterility is associated with a mitochondrial lineage that underwent a spectacular acceleration of DNA substitution rates, affecting the entire mitochondrial genome-this acceleration concerns both synonymous and non-synonymous substitutions and therefore results from increased mitogenome mutation rates. Consequently, mitochondrial haplotype divergence within the population is exceptionally high, matching that observed between snail taxa that diverged 475 million years ago. This result is reminiscent of similar accelerations in mitogenome evolution observed in plant clades where gynodioecy is frequent,^9,10 both being consistent with arms-race evolution of genome regions implicated in CMS.^11,12 Our study shows that genomic conflicts can trigger independent evolution of similar sex-determination systems in plants and animals and dramatically accelerate molecular evolution.

Frequent Paternal Mitochondrial Inheritance and Rapid Haplotype Frequency Shifts in Copepod Hybrids

Mitochondria are assumed to be maternally inherited in most animal species, and this foundational concept has fostered advances in phylogenetics, conservation, and population genetics. Like other animals, mitochondria were thought to be solely maternally inherited in the marine copepod Tigriopus californicus, which has served as a useful model for studying mitonuclear interactions, hybrid breakdown, and environmental tolerance. However, we present PCR, Sanger sequencing, and Illumina Nextera sequencing evidence that extensive paternal mitochondrial DNA (mtDNA) transmission is occurring in inter-population hybrids of T. californicus. PCR on four types of crosses between three populations (total sample size of 376 F1 individuals) with 20% genome-wide mitochondrial divergence showed 2% to 59% of F1 hybrids with both paternal and maternal mtDNA, where low and high paternal leakage values were found in different cross directions of the same population pairs. Sequencing methods further verified nucleotide similarities between F1 mtDNA and paternal mtDNA sequences. Interestingly, the paternal mtDNA in F1s from some crosses inherited haplotypes that were uncommon in the paternal population. Compared to some previous research on paternal leakage, we employed more rigorous methods to rule out contamination and false detection of paternal mtDNA due to non-functional nuclear mitochondrial DNA fragments. Our results raise the potential that other animal systems thought to only inherit maternal mitochondria may also have paternal leakage, which would then affect the interpretation of past and future population genetics or phylogenetic studies that rely on mitochondria as uniparental markers.

Massive gene rearrangement in mitogenomes of phytoseiid mites

Mitochondrial (mt) gene sequences have been widely used to infer phylogeny in animals. The relative order of mt genes in the mitogenome can also be a useful marker for evolution, but the propensity of mt gene rearrangements vary tremendously among taxa. Ticks and mites in Acari exemplify this trend as some families retain the ancestral arthropod gene order, while others show highly divergent gene orders. Mites in Phytoseiidae, many of which are effective biological control agents, show some of the most divergent gene orders. However, the diversity of mitogenome order within this family is little known. We thus sequenced three mt genomes of phytoseiid mites from two of the most speciose genera: Amblyseius swirskii (Athias-Henriot), Amblyseius tsugawai (Ehara) and Neoseiulus womersleyi (Schicha). We find differences in mt GC skew and nucleotide composition, especially between N. womersleyi and the two Amblyseius species. Each species within Phytoseiidae (including three previously available sequences) present a unique gene order. Phytoseiid mitogenomes show some of the highest numbers of breakpoints when compared to the ancestral arthropod order (up to 33), as well as high numbers of breakpoints within the family (14-30). This suggests a history of massive, ongoing mitogenome rearrangements in the family. Phylogenetic analyses of mt sequences confirm that the degree of gene rearrangements follows phylogenetic relatedness. We discuss possible causes for the high degree of mt gene rearrangement within phytoseiid mites as well as selection in the mt and nuclear genome tied to the independent evolution of many diverse feeding strategies in the family. Finally, we suggest N. womersleyi should be used instead of the synonym Amblyseius pseudolongispinosus.

Pervasive Mitonuclear Coadaptation Underlies Fast Development in Interpopulation Hybrids of a Marine Crustacean

Cellular energy production requires coordinated interactions between genetic components from the nuclear and mitochondrial genomes. This coordination results in coadaptation of interacting elements within populations. Interbreeding between divergent gene pools can disrupt coadapted loci and result in hybrid fitness breakdown. While specific incompatible loci have been detected in multiple eukaryotic taxa, the extent of the nuclear genome that is influenced by mitonuclear coadaptation is not clear in any species. Here, we used F₂ hybrids between two divergent populations of the copepod Tigriopus californicus to examine mitonuclear coadaptation across the nuclear genome. Using developmental rate as a measure of fitness, we found that fast-developing copepods had higher ATP synthesis capacity than slow developers, suggesting variation in developmental rates is at least partly associated with mitochondrial dysfunction. Using Pool-seq, we detected strong biases for maternal alleles across 7 (of 12) chromosomes in both reciprocal crosses in high-fitness hybrids, whereas low-fitness hybrids showed shifts toward the paternal population. Comparison with previous results on a different hybrid cross revealed largely different patterns of strong mitonuclear coadaptation associated with developmental rate. Our findings suggest that functional coadaptation between interacting nuclear and mitochondrial components is reflected in strong polygenic effects on this life-history phenotype, and reveal that molecular coadaptation follows independent evolutionary trajectories among isolated populations.

Faraway, so close. The comparative method and the potential of non-model animals in mitochondrial research

Inference from model organisms has been the engine for many discoveries in life science, but indiscriminate generalization leads to oversimplifications and misconceptions. Model organisms and inductive reasoning are irreplaceable: there is no other way to tackle the complexity of living systems. At the same time, it is not advisable to infer general patterns from a restricted number of species, which are very far from being representative of the diversity of life. Not all models are equal. Some organisms are suitable to find similarities across species, other highly specialized organisms can be used to focus on differences. In this opinion piece, we discuss the dominance of the mechanistic/reductionist approach in life sciences and make a case for an enhanced application of the comparative approach to study processes in all their various forms across different organisms. We also enlist some rising animal models in mitochondrial research, to exemplify how non-model organisms can be chosen in a comparative framework. These taxa often do not possess implemented tools and dedicated methods/resources. However, because of specific features, they have the potential to address still unanswered biological questions. Finally, we discuss future perspectives and caveats of the comparative method in the age of 'big data'. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Variation in developmental temperature alters adulthood plasticity of thermal tolerance in Tigriopus californicus

In response to environmental change, organisms rely on both genetic adaptation and phenotypic plasticity to adjust key traits that are necessary for survival and reproduction. Given the accelerating rate of climate change, plasticity may be particularly important. For organisms in warming aquatic habitats, upper thermal tolerance is likely to be a key trait, and many organisms express plasticity in this trait in response to developmental or adulthood temperatures. Although plasticity at one life stage may influence plasticity at another life stage, relatively little is known about this possibility for thermal tolerance. Here, we used locally adapted populations of the copepod Tigriopus californicus to investigate these potential effects in an intertidal ectotherm. We found that low latitude populations had greater critical thermal maxima (CT_max) than high latitude populations, and variation in developmental temperature altered CT_max plasticity in adults. After development at 25°C, CT_max was plastic in adults, whereas no adulthood plasticity in this trait was observed after development at 20°C. This pattern was identical across four populations, suggesting that local thermal adaptation has not shaped this effect among these populations. Differences in the capacities to maintain ATP synthesis rates and to induce heat shock proteins at high temperatures, two likely mechanisms of local adaptation in this species, were consistent with changes in CT_max owing to phenotypic plasticity, which suggests that there is likely mechanistic overlap between the effects of plasticity and adaptation. Together, these results indicate that developmental effects may have substantial impacts on upper thermal tolerance plasticity in adult ectotherms.

Age Dependent Dysfunction of Mitochondrial and ROS Metabolism Induced by Mitonuclear Mismatch

Mitochondrial and nuclear genomes have to coevolve to ensure the proper functioning of the different mitochondrial complexes that are assembled from peptides encoded by both genomes. Mismatch between these genomes is believed to be strongly selected against due to the consequent impairments of mitochondrial functions and induction of oxidative stress. Here, we used a Drosophila model harboring an incompatibility between a mitochondrial tRNA^tyr and its nuclear-encoded mitochondrial tyrosine synthetase to assess the cellular mechanisms affected by this incompatibility and to test the relative contribution of mitonuclear interactions and aging on the expression of impaired phenotypes. Our results show that the mitochondrial tRNA mutation caused a decrease in mitochondrial oxygen consumption in the incompatible nuclear background but no effect with the compatible nuclear background. Mitochondrial DNA copy number increased in the incompatible genotype but that increase failed to rescue mitochondrial functions. The flies harboring mismatch between nuclear and mitochondrial genomes had almost three times the relative mtDNA copy number and fifty percent higher rate of hydrogen peroxide production compared to other genome combinations at 25 days of age. We also found that aging exacerbated the mitochondrial dysfunctions. Our results reveal the tight interactions linking mitonuclear mismatch to mitochondrial dysfunction, mitochondrial DNA regulation, ROS production and aging.

Novel microRNAs are associated with population divergence in transcriptional response to thermal stress in an intertidal copepod

The role of gene expression in adaptation to differing thermal environments has been assayed extensively. Yet, in most natural systems, analyses of gene expression reveal only one level of the complexity of regulatory machineries. MicroRNAs (miRNAs) are small noncoding RNAs which are key components of many gene regulatory networks, and they play important roles in a variety of cellular pathways by modulating post-transcriptional quantities of mRNA available for protein synthesis. The characterization of miRNA loci and their regulatory dynamics in nonmodel systems are still largely understudied. In this study, we examine the role of miRNAs in response to high thermal stress in the intertidal copepod Tigriopus californicus. Allopatric populations of this species show varying levels of local adaptation with respect to thermal regimes, and previous studies showed divergence in gene expression between populations from very different thermal environments. We examined the transcriptional response to temperature stress in two populations separated by only 8 km by utilizing RNA-seq to quantify both mRNA and miRNA levels. Using the currently available genome sequence, we first describe the repertoire of miRNAs in T. californicus and assess the degree to which transcriptional response to temperature stress is governed by miRNA activity. The two populations showed large differences in the number of genes involved, the magnitude of change in commonly used genes and in the number of miRNAs involved in transcriptional modulation during stress. Our results suggest that an increased level of regulatory network complexity may underlie improved survivorship under thermal stress in one of the populations.

Mitochondrial genomes of the key zooplankton copepods Arctic Calanus glacialis and North Atlantic Calanus finmarchicus with the longest crustacean non-coding regions

We determined the nearly complete mitochondrial genomes of the Arctic Calanus glacialis and its North Atlantic sibling Calanus finmarchicus, which are key zooplankton components in marine ecosystems. The sequenced part of C. glacialis mitogenome is 27,342 bp long and consists of two contigs, while for C. finmarchicus it is 29,462 bp and six contigs, what makes them the longest reported copepod mitogenomes. The typical set of metazoan mitochondrial genes is present in these mitogenomes, although the non-coding regions (NCRs) are unusually long and complex. The mitogenomes of the closest species C. glacialis and C. finmarchicus, followed by the North Pacific C. sinicus, are structurally similar and differ from the much more typical of deep-water, Arctic C. hyperboreus. This evolutionary trend for the expansion of NCRs within the Calanus mitogenomes increases mitochondrial DNA density, what resulted in its similar density to the nuclear genome. Given large differences in the length and structure of C. glacialis and C. finmarchicus mitogenomes, we conclude that the species are genetically distinct and thus cannot hybridize. The molecular resources presented here: the mitogenomic and rDNA sequences, and the database of repetitive elements should facilitate the development of genetic markers suitable in pursuing evolutionary research in copepods.

A simple and efficient total genomic DNA extraction method for individual zooplankton

Molecular approaches are widely applied in species identification and taxonomic studies of minute zooplankton. One of the most focused zooplankton nowadays is from Subclass Copepoda. Accurate species identification of all life stages of the generally small sized copepods through molecular analysis is important, especially in taxonomic and systematic assessment of harpacticoid copepod populations and to understand their dynamics within the marine community. However, total genomic DNA (TGDNA) extraction from individual harpacticoid copepods can be problematic due to their small size and epibenthic behavior. In this research, six TGDNA extraction methods done on individual harpacticoid copepods were compared. The first new simple, feasible, efficient and consistent TGDNA extraction method was designed and compared with the commercial kit and modified available TGDNA extraction methods. The newly described TGDNA extraction method, "Incubation in PCR buffer" method, yielded good and consistent results based on the high success rate of PCR amplification (82%) compared to other methods. Coupled with its relatively consistent and economical method the "Incubation in PCR buffer" method is highly recommended in the TGDNA extraction of other minute zooplankton species.

DNA Barcoding of Metazoan Zooplankton Copepods from South Korea

Copepods, small aquatic crustaceans, are the most abundant metazoan zooplankton and outnumber every other group of multicellular animals on earth. In spite of ecological and biological importance in aquatic environment, their morphological plasticity, originated from their various lifestyles and their incomparable capacity to adapt to a variety of environments, has made the identification of species challenging, even for expert taxonomists. Molecular approaches to species identification have allowed rapid detection, discrimination, and identification of cryptic or sibling species based on DNA sequence data. We examined sequence variation of a partial mitochondrial cytochrome C oxidase I gene (COI) from 133 copepod individuals collected from the Korean Peninsula, in order to identify and discriminate 94 copepod species covering six copepod orders of Calanoida, Cyclopoida, Harpacticoida, Monstrilloida, Poecilostomatoida and Siphonostomatoida. The results showed that there exists a clear gap with ca. 20 fold difference between the averages of within-specific sequence divergence (2.42%) and that of between-specific sequence divergence (42.79%) in COI, suggesting the plausible utility of this gene in delimitating copepod species. The results showed, with the COI barcoding data among 94 copepod species, that a copepod species could be distinguished from the others very clearly, only with four exceptions as followings: Mesocyclops dissimilis–Mesocyclops pehpeiensis (0.26% K2P distance in percent) and Oithona davisae–Oithona similis (1.1%) in Cyclopoida, Ostrincola japonica–Pseudomyicola spinosus (1.5%) in Poecilostomatoida, and Hatschekia japonica–Caligus quadratus (5.2%) in Siphonostomatoida. Thus, it strongly indicated that COI may be a useful tool in identifying various copepod species and make an initial progress toward the construction of a comprehensive DNA barcode database for copepods inhabiting the Korean Peninsula.

Mitochondrial genome of the intertidal acorn barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Sessilia): Gene order comparison and phylogenetic consideration within Sessilia

The complete mitochondrial genome of the intertidal barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Maxillopoda: Sessilia) is presented. The genome is a circular molecule of 15,200 bp, which encodes 13 PCGs, 2 ribosomal RNA genes, and 22 transfer RNA genes. All non-coding regions are 591 bp in length, with the longest one speculated as the control region (389 bp), which is located between srRNA and trnK. The overall A+T content of the mitochondrial genome of T. serrata is 65.4%, which is lowest among all the eight mitochondrial genomes reported from sessile barnacles. There are variations of initiation and stop codons in the reported sessile barnacle mitochondrial genomes. Large-scale gene rearrangements are found in these genomes as compared to the pancrustacean ground pattern. ML and Bayesian analyses of all 15 complete mitochondrial genomes available from Maxillopoda lead to identical phylogenies. The phylogenetic tree based on mitochondrial PCGs shows that Argulus americanus (Branchiura) cluster with Armillifer armillatus (Pentastomida), distinct from all ten species from Cirripedia. Within the order Sessilia, Amphibalanus amphitrite (Balanidae) clusters with Striatobalanus amaryllis (Archaeobalanidae), and Nobia grandis (Pyrgomatidae). However, the two Megabalanus (Balanidae) are separated from the above grouping, resulting in non-monophyly of the family Balanidae. Moreover, the two Megabalanus have large-scale rearrangements as compared to the gene order shared by former three species. Therefore, both phylogenetic analysis using PCG sequences and gene order comparison suggest that Balanidae is not a monophyletic group. Given the limited taxa and moderate support values of the internal branches, the non-monophyly of the family Balanidae requires further verification.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities

BACKGROUND

Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform.

RESULTS

The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations.

CONCLUSIONS

Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.

Phylogenetic information content of Copepoda ribosomal DNA repeat units: ITS1 and ITS2 impact

The utility of various regions of the ribosomal repeat unit for phylogenetic analysis was examined in 16 species representing four families, nine genera, and two orders of the subclass Copepoda (Crustacea). Fragments approximately 2000 bp in length containing the ribosomal DNA (rDNA) 18S and 28S gene fragments, the 5.8S gene, and the internal transcribed spacer regions I and II (ITS1 and ITS2) were amplified and analyzed. The DAMBE (Data Analysis in Molecular Biology and Evolution) software was used to analyze the saturation of nucleotide substitutions; this test revealed the suitability of both the 28S gene fragment and the ITS1/ITS2 rDNA regions for the reconstruction of phylogenetic trees. Distance (minimum evolution) and probabilistic (maximum likelihood, Bayesian) analyses of the data revealed that the 28S rDNA and the ITS1 and ITS2 regions are informative markers for inferring phylogenetic relationships among families of copepods and within the Cyclopidae family and associated genera. Split-graph analysis of concatenated ITS1/ITS2 rDNA regions of cyclopoid copepods suggested that the Mesocyclops, Thermocyclops, and Macrocyclops genera share complex evolutionary relationships. This study revealed that the ITS1 and ITS2 regions potentially represent different phylogenetic signals.

Prediction of the peptidomes of Tigriopus californicus and Lepeophtheirus salmonis (Copepoda, Crustacea)

Transcriptome mining is a powerful method for crustacean peptide discovery, especially when large sequence datasets are available and an appropriate reference is extant. Recently, a 206,041-sequence transcriptome for the copepod Calanus finmarchicus was mined for peptide-encoding transcripts, with ones for 17 families/subfamilies identified. Here, the deduced Calanus pre/preprohormones were used as templates for peptide discovery in the copepods Tigriopus californicus and Lepeophtheirus salmonis; large transcriptome shotgun assembly datasets are publicly accessible for both species. Sixty-five Tigriopus and 17 Lepeophtheirus transcripts, encompassing 22 and 13 distinct peptide families/subfamilies, respectively, were identified, with the structures of 161 and 70 unique mature peptides predicted from the deduced precursors. The identified peptides included members of the allatostatin A, allatostatin C, bursicon α, bursicon β, CAPA/periviscerokinin/pyrokinin, crustacean cardioactive peptide, crustacean hyperglycemic hormone/ion transport peptide, diuretic hormone 31, FLRFamide, leucokinin, myosuppressin, neuroparsin, neuropeptide F, orcokinin, and tachykinin-related peptide families, most of which possess novel structures, though isoforms from other copepods are known. Of particular note was the discovery of novel isoforms of adipokinetic hormone-corazonin-like peptide, allatotropin, corazonin, eclosion hormone and intocin, peptide families previously unidentified in copepods. In addition, Tigriopus precursors for two previously unknown peptide groups were discovered, one encoding GSEFLamides and the other DXXRLamides; precursors for the novel FXGGXamide family were identified from both Tigriopus and Lepeophtheirus. These data not only greatly expand the catalog of known copepod peptides, but also provide strong foundations for future functional studies of peptidergic signaling in members of this ecologically important crustacean subclass.

The mitochondrial genome impacts respiration but not fermentation in interspecific Saccharomyces hybrids

In eukaryotes, mitochondrial DNA (mtDNA) has high rate of nucleotide substitution leading to different mitochondrial haplotypes called mitotypes. However, the impact of mitochondrial genetic variant on phenotypic variation has been poorly considered in microorganisms because mtDNA encodes very few genes compared to nuclear DNA, and also because mitochondrial inheritance is not uniparental. Here we propose original material to unravel mitotype impact on phenotype: we produced interspecific hybrids between S. cerevisiae and S. uvarum species, using fully homozygous diploid parental strains. For two different interspecific crosses involving different parental strains, we recovered 10 independent hybrids per cross, and allowed mtDNA fixation after around 80 generations. We developed PCR-based markers for the rapid discrimination of S. cerevisiae and S. uvarum mitochondrial DNA. For both crosses, we were able to isolate fully isogenic hybrids at the nuclear level, yet possessing either S. cerevisiae mtDNA (Sc-mtDNA) or S. uvarum mtDNA (Su-mtDNA). Under fermentative conditions, the mitotype has no phenotypic impact on fermentation kinetics and products, which was expected since mtDNA are not necessary for fermentative metabolism. Alternatively, under respiratory conditions, hybrids with Sc-mtDNA have higher population growth performance, associated with higher respiratory rate. Indeed, far from the hypothesis that mtDNA variation is neutral, our work shows that mitochondrial polymorphism can have a strong impact on fitness components and hence on the evolutionary fate of the yeast populations. We hypothesize that under fermentative conditions, hybrids may fix stochastically one or the other mt-DNA, while respiratory environments may increase the probability to fix Sc-mtDNA.

Elevated oxidative damage is correlated with reduced fitness in interpopulation hybrids of a marine copepod

Aerobic energy production occurs via the oxidative phosphorylation pathway (OXPHOS), which is critically dependent on interactions between the 13 mitochondrial DNA (mtDNA)-encoded and approximately 70 nuclear-encoded protein subunits. Disruptive mutations in any component of OXPHOS can result in impaired ATP production and exacerbated oxidative stress; in mammalian systems, such mutations are associated with ageing as well as numerous diseases. Recent studies have suggested that oxidative stress plays a role in fitness trade-offs in life-history evolution and functional ecology. Here, we show that outcrossing between populations with divergent mtDNA can exacerbate cellular oxidative stress in hybrid offspring. In the copepod Tigriopus californicus, we found that hybrids that showed evidence of fitness breakdown (low fecundity) also exhibited elevated levels of oxidative damage to DNA, whereas those with no clear breakdown did not show significantly elevated damage. The extent of oxidative stress in hybrids appears to be dependent on the degree of genetic divergence between their respective parental populations, but this pattern requires further testing using multiple crosses at different levels of divergence. Given previous evidence in T. californicus that hybridization disrupts nuclear/mitochondrial interactions and reduces hybrid fitness, our results suggest that such negative intergenomic epistasis may also increase the production of damaging cellular oxidants; consequently, mtDNA evolution may play a significant role in generating postzygotic isolating barriers among diverging populations.

A disproportionate role for mtDNA in Dobzhansky-Muller incompatibilities?

Evolution in allopatric populations can lead to incompatibilities that result in reduced hybrid fitness and ultimately reproductive isolation upon secondary contact. The Dobzhansky-Muller (DM) model nicely accounts for the evolution of such incompatibilities. Although DM incompatibilities were originally conceived as resulting of interactions between nuclear genes, recent studies have documented cases where incompatibilities have arisen between nuclear and mitochondrial genomes (mtDNA). Although mtDNA comprises only a tiny component (typically <<0.01%) of an organism's genetic material, several features of mtDNA may lead to a disproportionate contribution to the evolution of hybrid incompatibilities: (i) essentially all functions of mtDNA require interaction with nuclear gene products. All mtDNA-encoded proteins are components of the oxidative phosphorylation (OXPHOS) system and all mtDNA-encoded RNAs are part of the mitochondrial protein synthetic machinery; both processes require interaction with nuclear-encoded proteins for function. (ii) Transcription and replication of mtDNA also involve mitonuclear interactions as nuclear-encoded proteins must bind to regulatory motifs in the mtDNA to initiate these processes. (iii) Although features of mtDNA vary amongst taxa, metazoan mtDNA is typically characterized by high nucleotide substitution rates, lack of recombination and reduced effective population sizes that collectively lead to increased chance fixation of mildly deleterious mutations. Combined, these features create an evolutionary dynamic where rapid mtDNA evolution favours compensatory nuclear gene evolution, ultimately leading to co-adaptation of mitochondrial and nuclear genomes. When previously isolated lineages hybridize in nature or in the lab, intergenomic co-adaptation is disrupted and hybrid breakdown is observed; the role of intergenomic co-adaptation in hybrid breakdown and speciation will generally be most pronounced when rates of mtDNA evolution are high or when restricted gene flow results in significant population differentiation.

In vivo quantification reveals extensive natural variation in mitochondrial form and function in Caenorhabditis briggsae

We have analyzed natural variation in mitochondrial form and function among a set of Caenorhabditis briggsae isolates known to harbor mitochondrial DNA structural variation in the form of a heteroplasmic nad5 gene deletion (nad5Δ) that correlates negatively with organismal fitness. We performed in vivo quantification of 24 mitochondrial phenotypes including reactive oxygen species level, membrane potential, and aspects of organelle morphology, and observed significant among-isolate variation in 18 traits. Although several mitochondrial phenotypes were non-linearly associated with nad5Δ levels, most of the among-isolate phenotypic variation could be accounted for by phylogeographic clade membership. In particular, isolate-specific mitochondrial membrane potential was an excellent predictor of clade membership. We interpret this result in light of recent evidence for local adaptation to temperature in C. briggsae. Analysis of mitochondrial-nuclear hybrid strains provided support for both mtDNA and nuclear genetic variation as drivers of natural mitochondrial phenotype variation. This study demonstrates that multicellular eukaryotic species are capable of extensive natural variation in organellar phenotypes and highlights the potential of integrating evolutionary and cell biology perspectives.

Naturally occurring mitochondrial DNA haplotypes exhibit metabolic differences: insight into functional properties of mitochondria

Linking the mitochondrial genotype and the organismal phenotype is of paramount importance in evolution of mitochondria. In this study, we determined the differences in catalytic properties of mitochondria dictated by divergences in the siII and siIII haplogroups of Drosophila simulans using introgressions of siII mtDNA type into the siIII nuclear background. We used a novel in situ method (permeabilized fibers) that allowed us to accurately measure the consumption of oxygen by mitochondria in constructed siII-introgressed flies and in siIII-control flies. Our results showed that the catalytic capacity of the electron transport system is not impaired by introgressions, suggesting that the functional properties of mitochondria are tightly related to the mtDNA haplogroup and not to the nuclear DNA or to the mito-nuclear interactions. This is the first study, to our knowledge, that demonstrates a naturally occurring haplogroup can confer specific functional differences in aspects of mitochondrial metabolism. This study illustrates the importance of mtDNA changes on organelle evolution and highlights the potential bioenergetic and metabolic impacts that divergent mitochondrial haplogroups may have upon a wide variety of species including humans.

Genomic resources for sea lice: analysis of ESTs and mitochondrial genomes

Sea lice are common parasites of both farmed and wild salmon. Salmon farming constitutes an important economic market in North America, South America, and Northern Europe. Infections with sea lice can result in significant production losses. A compilation of genomic information on different genera of sea lice is an important resource for understanding their biology as well as for the study of population genetics and control strategies. We report on over 150,000 expressed sequence tags (ESTs) from five different species (Pacific Lepeophtheirus salmonis (49,672 new ESTs in addition to 14,994 previously reported ESTs), Atlantic L. salmonis (57,349 ESTs), Caligus clemensi (14,821 ESTs), Caligus rogercresseyi (32,135 ESTs), and Lernaeocera branchialis (16,441 ESTs)). For each species, ESTs were assembled into complete or partial genes and annotated by comparisons to known proteins in public databases. In addition, whole mitochondrial (mt) genome sequences of C. clemensi (13,440 bp) and C. rogercresseyi (13,468 bp) were determined and compared to L. salmonis. Both nuclear and mtDNA genes show very high levels of sequence divergence between these ectoparastic copepods suggesting that the different species of sea lice have been in existence for 37-113 million years and that parasitic association with salmonids is also quite ancient. Our ESTs and mtDNA data provide a novel resource for the study of sea louse biology, population genetics, and control strategies. This genomic information provides the material basis for the development of a 38K sea louse microarray that can be used in conjunction with our existing 44K salmon microarray to study host-parasite interactions at the molecular level. This report represents the largest genomic resource for any copepod species to date.

Observing copepods through a genomic lens

Background

Copepods outnumber every other multicellular animal group. They are critical components of the world's freshwater and marine ecosystems, sensitive indicators of local and global climate change, key ecosystem service providers, parasites and predators of economically important aquatic animals and potential vectors of waterborne disease. Copepods sustain the world fisheries that nourish and support human populations. Although genomic tools have transformed many areas of biological and biomedical research, their power to elucidate aspects of the biology, behavior and ecology of copepods has only recently begun to be exploited.

Discussion

The extraordinary biological and ecological diversity of the subclass Copepoda provides both unique advantages for addressing key problems in aquatic systems and formidable challenges for developing a focused genomics strategy. This article provides an overview of genomic studies of copepods and discusses strategies for using genomics tools to address key questions at levels extending from individuals to ecosystems. Genomics can, for instance, help to decipher patterns of genome evolution such as those that occur during transitions from free living to symbiotic and parasitic lifestyles and can assist in the identification of genetic mechanisms and accompanying physiological changes associated with adaptation to new or physiologically challenging environments. The adaptive significance of the diversity in genome size and unique mechanisms of genome reorganization during development could similarly be explored. Genome-wide and EST studies of parasitic copepods of salmon and large EST studies of selected free-living copepods have demonstrated the potential utility of modern genomics approaches for the study of copepods and have generated resources such as EST libraries, shotgun genome sequences, BAC libraries, genome maps and inbred lines that will be invaluable in assisting further efforts to provide genomics tools for copepods.

Summary

Genomics research on copepods is needed to extend our exploration and characterization of their fundamental biological traits, so that we can better understand how copepods function and interact in diverse environments. Availability of large scale genomics resources will also open doors to a wide range of systems biology type studies that view the organism as the fundamental system in which to address key questions in ecology and evolution.

Repeated parallel evolution of minimal rRNAs revealed from detailed comparative analysis

The concept of a minimal ribosomal RNA-containing ribosome, a structure with a minimal set of elements capable of providing protein biosynthesis, is essential for understanding this fundamental cellular process. Nematodes and trypanosomes have minimal mitochondrial rRNAs and detailed reconstructions of their secondary structures indicate that certain conserved helices have been lost in these taxa. In contrast, several recent studies on acariform mites have argued that minimal rRNAs may evolve via shortening of secondary structure elements but not the loss of these elements as shown for trypanosomes and nematodes. Based on extensive structural analysis of chelicerate arthropods, we demonstrate that extremely short rRNAs of acariform mites share certain structural modifications with nematodes and trypanosomes: loss of helices of the GTPase region and divergence in the evolutionarily conserved connecting loop between helices H1648 and H1764 of the large subunit rRNA. These highly concerted parallel modifications indicate that minimal rRNAs were generated under the strong selection that favored or tolerated reductions of helices in particular locations while maintaining the functionality of the rRNA molecules throughout evolution. We also discuss potential evolution of minimal rRNAs and atypical transfer RNAs.

Distinctive mitochondrial genome of Calanoid copepod Calanus sinicus with multiple large non-coding regions and reshuffled gene order: useful molecular markers for phylogenetic and population studies

Background

Copepods are highly diverse and abundant, resulting in extensive ecological radiation in marine ecosystems. Calanus sinicus dominates continental shelf waters in the northwest Pacific Ocean and plays an important role in the local ecosystem by linking primary production to higher trophic levels. A lack of effective molecular markers has hindered phylogenetic and population genetic studies concerning copepods. As they are genome-level informative, mitochondrial DNA sequences can be used as markers for population genetic studies and phylogenetic studies.

Results

The mitochondrial genome of C. sinicus is distinct from other arthropods owing to the concurrence of multiple non-coding regions and a reshuffled gene arrangement. Further particularities in the mitogenome of C. sinicus include low A + T-content, symmetrical nucleotide composition between strands, abbreviated stop codons for several PCGs and extended lengths of the genes atp6 and atp8 relative to other copepods. The monophyletic Copepoda should be placed within the Vericrustacea. The close affinity between Cyclopoida and Poecilostomatoida suggests reassigning the latter as subordinate to the former. Monophyly of Maxillopoda is rejected. Within the alignment of 11 C. sinicus mitogenomes, there are 397 variable sites harbouring three 'hotspot' variable sites and three microsatellite loci.

Conclusion

The occurrence of the circular subgenomic fragment during laboratory assays suggests that special caution should be taken when sequencing mitogenomes using long PCR. Such a phenomenon may provide additional evidence of mitochondrial DNA recombination, which appears to have been a prerequisite for shaping the present mitochondrial profile of C. sinicus during its evolution. The lack of synapomorphic gene arrangements among copepods has cast doubt on the utility of gene order as a useful molecular marker for deep phylogenetic analysis. However, mitochondrial genomic sequences have been valuable markers for resolving phylogenetic issues concerning copepods. The variable site maps of C. sinicus mitogenomes provide a solid foundation for population genetic studies.

Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus

The accumulation of genetic incompatibilities between isolated populations is thought to lead to the evolution of intrinsic postzygotic isolation. The molecular basis for these mechanisms, however, remains poorly understood. The intertidal copepod Tigriopus californicus provides unique opportunities for addressing mechanistic questions regarding the early stages of speciation; hybrids between highly divergent populations are fertile and viable, but exhibit reduced fitness at the F(2) or later generations. Given the current scarcity of genomic information in taxa at incipient stages of reproductive isolation, we utilize high-throughout 454 pyrosequencing to characterize a substantial fraction of protein-coding regions (the transcriptome) of T. californicus. Our sequencing effort was divided equally between two divergent populations in order to estimate levels of divergence and to reveal patterns of selection across the transcriptome. Assembly of sequences generated over 40,000 putatively unique transcripts (unigenes) for each population, 19,622 of which were orthologous between populations. BLAST searches of public databases determined protein identity and functional features for 15,402 and 12,670 unigenes, respectively. Based on rates of nonsynonymous and synonymous substitutions in 5897 interpopulation orthologs (those >150 bp and with at least 2X coverage), we identified 229 potential targets of positive selection. Many of these genes are predicted to be involved in several metabolic processes, and to function in hydrolase, peptidase and binding activities. The library of T. californicus coding regions, annotated with their predicted functions and level of divergence, will serve as an invaluable resource for elucidating molecular mechanisms underlying the early stages of speciation.

Genetic architecture of metabolic rate: environment specific epistasis between mitochondrial and nuclear genes in an insect

The extent to which mitochondrial DNA (mtDNA) variation is involved in adaptive evolutionary change is currently being reevaluated. In particular, emerging evidence suggests that mtDNA genes coevolve with the nuclear genes with which they interact to form the energy producing enzyme complexes in the mitochondria. This suggests that intergenomic epistasis between mitochondrial and nuclear genes may affect whole-organism metabolic phenotypes. Here, we use crossed combinations of mitochondrial and nuclear lineages of the seed beetle Callosobruchus maculatus and assay metabolic rate under two different temperature regimes. Metabolic rate was affected by an interaction between the mitochondrial and nuclear lineages and the temperature regime. Sequence data suggests that mitochondrial genetic variation has a role in determining the outcome of this interaction. Our genetic dissection of metabolic rate reveals a high level of complexity, encompassing genetic interactions over two genomes, and genotype × genotype × environment interactions. The evolutionary implications of these results are twofold. First, because metabolic rate is at the root of life histories, our results provide insights into the complexity of life-history evolution in general, and thermal adaptation in particular. Second, our results suggest a mechanism that could contribute to the maintenance of nonneutral mtDNA polymorphism.

The complete mitochondrial genome of the cyclopoid copepod Paracyclopina nana: a highly divergent genome with novel gene order and atypical gene numbers

In this paper, we describe the complete mitogenome of the cyclopoid copepod Paracyclopina nana with emphasis on the highly rearranged gene order and high divergence against published copepod mitogenomes. The P. nana mtDNA is 15,981 bp in length (70.9% AT) and consists of 37 genes (12 protein-coding genes, 2 rRNAs, 23 tRNAs) that are atypical for metazoan mitogenomes. Unusually, it contains an extra tRNA (tRNA-Ala) but it does not contain the ATPase 8 gene. The P. nana mitogenome has a long putative control region with high AT content (1351 bp, 77.0% AT). The Cyt b was considerably short in length, compared to other crustaceans. Compared to typical mitogenomes of arthropods and copepods, the gene order of the P. nana mitogenome is highly rearranged with a novel gene structure. In addition, P. nana has highly divergent mt genes (mostly less than 50%), judged by amino acid substitution. We present the first complete mitogenome sequence from a cyclopoid copepod, thereby increasing our understanding of copepod and crustacean evolution from the mitochondrial point of view.