Does nonneutral evolution shape observed patterns of DNA variation in animal mitochondrial genomes?

The development of genus-specific and species-specific real-time PCR assays for the authentication of Patagonian toothfish and Antarctic toothfish in commercial seafood products

BACKGROUND

The substitution or mislabeling of toothfish is an issue of significant concern for seafood authorities; it also reduces the effectiveness of marine conservation and management programs for its over-exploitation and illegal trafficking, boosting the need for identification methods.

RESULTS

Two species-specific real-time polymerase chain reaction (PCR) assays for the identification of Patagonian toothfish (Dissostichus eleginoides) and Antarctic toothfish (Dissostichus mawsoni) and a genus-specific real-time PCR assay for Dissostichus spp. identification were developed based on fragments of the 16S rRNA and COI (cytochrome c oxidase subunit I) genes. These methods were confirmed to be rapid, simple, and sensitive (absolute sensitivity of 0.0002 ng μL^-1 and relative sensitivity of 0.1 g kg^-1 with good specificity). These methods can be applied to processed and commercial fish products.

CONCLUSIONS

These approaches can be beneficial for protecting both consumers and producers from economic fraud and might also help protect toothfish from over-exploitation as well as combat illegal, unreported, and unregulated (IUU) fisheries. © 2021 Society of Chemical Industry.

Dispersal limitations and long-term persistence drive differentiation from haplotypes to communities within a tropical sky-island: Evidence from community metabarcoding

Neutral theory proposes that dispersal stochasticity is one of the main drivers of local diversity. Haplotypes-level genetic variation can now be efficiently sampled from across whole communities, thus making it possible to test neutral predictions from the genetic to species-level diversity, and higher. However, empirical data is still limited, with the few studies to date coming from temperate latitudes. Here, we focus on a tropical mountain within the Transmexican Volcanic Belt to evaluate spatially fine-scale patterns of arthropod community assembly to understand the role of dispersal limitation and landscape features as drivers of diversity. We sampled whole-communities of arthropods for eight orders at a spatial scale ranging from 50 m to 19 km, using whole community metabarcoding. We explored multiple hierarchical levels, from individual haplotypes to lineages at 0.5, 1.5, 3, 5, and 7.5% similarity thresholds, to evaluate patterns of richness, turnover, and distance decay of similarity with isolation-by-distance and isolation-by-resistance (costs to dispersal given by landscape features) approaches. Our results showed that distance and altitude influence distance decay of similarity at all hierarchical levels. This holds for arthropod groups of contrasting dispersal abilities, but with different strength depending on the spatial scale. Our results support a model where local-scale differentiation mediated by dispersal constraints, combined with long-term persistence of lineages, is an important driver of diversity within tropical sky islands.

Epigenetic Regulation of the Nuclear and Mitochondrial Genomes: Involvement in Metabolism, Development, and Disease

Our understanding of the interactions between the nuclear and mitochondrial genomes is becoming increasingly important as they are extensively involved in establishing early development and developmental progression. Evidence from various biological systems indicates the interdependency between the genomes, which requires a high degree of compatibility and synchrony to ensure effective cellular function throughout development and in the resultant offspring. During development, waves of DNA demethylation, de novo methylation, and maintenance methylation act on the nuclear genome and typify oogenesis and pre- and postimplantation development. At the same time, significant changes in mitochondrial DNA copy number influence the metabolic status of the developing organism in a typically cell-type-specific manner. Collectively, at any given stage in development, these actions establish genomic balance that ensures each developmental milestone is met and that the organism's program for life is established.

Elucidation of fecal inputs into the River Tagus catchment (Portugal) using source-specific mitochondrial DNA, HAdV, and phage markers

Determining the source of fecal contamination in a water body is important for the application of appropriate remediation measures. However, it has been suggested in the extant literature that this can best be achieved using a 'toolbox' of molecular- and culture-based methods. In response, this study deployed three indicators (Escherichia coli (EC), intestinal enterococci (IE) and somatic coliphages (SC)), one culture-dependent human marker (Bacteroides (GB-124) bacteriophage) and five culture-independent markers (human adenovirus (HAdV), human (HMMit), cattle (CWMit), pig (PGMit) and poultry (PLMit) mitochondrial DNA markers (mtDNA)) within the River Tagus catchment (n = 105). Water samples were collected monthly over a 13-month sampling campaign at four sites (impacted by significant specific human and non-human inputs and influenced by differing degrees of marine and freshwater mixing) to determine the dominant fecal inputs and assess geographical, temporal, and meteorological (precipitation, UV, temperature) fluctuations. Our results revealed that all sampling sites were not only highly impacted by fecal contamination but that this contamination originated from human and from a range of agricultural animal sources. HMMit was present in a higher percentage (83%) and concentration (4.20 log GC/100 mL) than HAdV (32%, 2.23 log GC/100 mL) and GB-124 bacteriophage with the latter being detected once. Animal mtDNA markers were detected, with CWMit found in 73% of samples with mean concentration of 3.74 log GC/100 mL. Correlation was found between concentrations of fecal indicators (EC, IE and SC), CWMit and season. Levels of CWMit were found to be related to physico-chemical parameters, such as temperature and UV radiation, possibly as a result of the increasing presence of livestock outside in warmer months. This study provides the first evaluation of such a source-associated 'toolbox' for monitoring surface water in Portugal, and the conclusions may inform future implementation of surveillance and remediation strategies for improving water quality.

The phylogenetic position of Arhaphe deviatica within Hemipteran insects: A potential model species for eco-devo studies of symbiosis

Insecta is known to be the most diverse group of species, exhibiting numerous forms of endosymbiotic associations. Molecular techniques have provided significant indicators for insect-microbe interactions. The present study aimed to register one of the true bugs of pentatomomorpha and clarify its taxonomic position through phylogenetic analysis of the partial 16S rRNA gene region. A maximum likelihood analysis retrieved a generally well-supported phylogeny based on Tamura 3-parameter model. Based on the partial mitochondrial 16S rRNA gene sequences, a phylogenetic study of suborder Heteroptera relationships within Hemipteras' order was constructed. Sequences of 221 bases of the 3' end of the gene from 28 species within 16 families were analyzed. This analysis and bootstrap confidence revealed two major clades comprising four suborders within Hemiptera, with a close relationship between Heteroptera + (Sternorrhyncha + (Auchenorrhycha + Coleorrhyncha)). Infraorder Pentatomomorpha is forming a sister group with a substantial bootstrap value to Cimicomorpha. Pyrrhocoroidea forms a sister relationship with Lygaeoidea + Coreoidea. There is a close relationship between Largidae and Pyrrhocoridae within Pyrrhocoroidea. The results show that the present species is firmly embedded in the genus Arhaphe with 94.35% sequence resemblance to its congeners. Besides, the recovered hemipteran species considered a potential model group for studying different symbionts. We propose both phylogenetic and ecological evolutionary developmental biology viewpoints for a more synthetic understanding of insect populations' molecular evolution.

Single-strand conformation polymorphism (SSCP) of mitochondrial genes helps to estimate genetic differentiation, demographic parameters and phylogeny of Glossina palpalis palpalis populations from West and Central Africa

A good understanding of tsetse fly population structure and migration is essential to optimize the control of sleeping sickness. This can be done by studying the genetics of tsetse fly populations. In this work, we estimated the genetic differentiation within and among geographically separated Glossina palpalis palpalis populations from Cameroon, the Democratic Republic of the Congo and Ivory Coast. We determined the demographic history of these populations and assessed phylogenetic relationships among individuals of this sub-species. A total of 418 tsetse flies were analysed: 258 were collected in four locations in Cameroon (Bipindi, Campo, Fontem and Bafia), 100 from Azaguié and Nagadoua in Ivory Coast and 60 from Malanga in the Democratic Republic of the Congo. We examined genetic variation at three mitochondrial loci: COI, COII-TLII, and 16S2. 34 haplotypes were found, of which 30 were rare, since each was present in <5% of the total number of individuals. No haplotype was shared among Cameroon, Ivory Coast and the Democratic Republic of the Congo populations. The fixation index F_ST of 0.88 showed a high genetic distance between Glossina palpalis palpalis populations from the three countries. That genetic distance was correlated to the geographic distance between populations. We also found that there is substantial gene flow between flies from locations separated by over 100 km in Cameroon and between flies from locations separated by over 200 km in Ivory Coast. Demographic parameters suggest that the tsetse flies from Fontem (Cameroon) had reduced in population size in the recent past. Phylogenetic analysis confirms that Glossina palpalis palpalis originating from the Democratic Republic of the Congo are genetically divergent from the two other countries as already published in previous studies.

Sequencing and characterization of the complete mitochondrial genome of Mishmi takin (Budorcas taxicolor taxicolor) and comparison with the other Caprinae species

Mishmi takin (Budorcas taxicolor taxicolor) is an endangered and least studies species. We sequenced and annotated a 16,584 bp long complete mitogenome of Mishmi takin (B. t. taxicolor). It was encoded by 37 genes: 13 Protein-Coding Genes (PCGs), two ribosomal RNA genes, 22 transfer RNA, and non-coding control region. All tRNA genes show a typical secondary cloverleaf structure, except tRNA^ser(AGY) in which dihydrouridine arm did not form a stable structure. We observed 82 bp nucleotide deletions in the control region of Mishmi takin. The overall analysis of 13 PCGs, two rRNA and 22 tRNA of Mishmi takin indicates its close relationship with the B. t. tibetana than the B. t. bedfordi. The overall genetic similarity indicated that all takin species is closely related to the P. nayaur than the C. ibex and Ovibos moschatus. This study provides a useful resource for further studies on conservation and population genetics of Budorcas.

Mitochondrial DNA Fitness Depends on Nuclear Genetic Background in Drosophila

Mitochondrial DNA (mtDNA) has been one of the most extensively studied molecules in ecological, evolutionary and clinical genetics. In its early application in evolutionary genetics, mtDNA was assumed to be a selectively neutral marker conferring negligible fitness consequences for its host. However, this dogma has been overturned in recent years due to now extensive evidence for non-neutral evolutionary dynamics. Since mtDNA proteins physically interact with nuclear proteins to provide the mitochondrial machinery for aerobic ATP production, among other cell functions, co-variation of the respective genes is predicted to affect organismal fitness. To test this hypothesis we used an mtDNA-nuclear DNA introgression model in Drosophila melanogaster to test the fitness of genotypes in perturbation-reperturbation population cages and in a non-competitive assay for female fecundity. Genotypes consisted of both conspecific and heterospecific mtDNA-nDNA constructs, with either D. melanogaster or D. simulans mtDNAs on two alternative D. melanogaster nuclear backgrounds, to investigate mitonuclear genetic interactions (G x G effects). We found considerable variation between nuclear genetic backgrounds on the selection of mtDNA haplotypes. In addition, there was variation in the selection on mtDNAs pre- and post- reperturbation, demonstrating overall poor repeatability of selection. There was a strong influence of nuclear background on non-competitive fecundity across all the mtDNA species types. In only one of the four cage types did we see a significant fecundity effect between genotypes that could help explain the respective change in genotype frequency over generational time. We discuss these results in the context of G x G interactions and the possible influence of stochastic environments on mtDNA-nDNA selection.

A case study characterizing animal fecal sources in surface water using a mitochondrial DNA marker

Water quality impairment by fecal waste in coastal watersheds is a public health issue. The present study provided evidence for the use of a mitochondrial (mtDNA) marker to detect animal fecal sources in surface water. The accurate identification of fecal pollution is based on the notion that fecal microorganisms preferentially inhabit a host animal's gut environment. In contrast, mtDNA host-specific markers are inherent to eukaryotic host cells, which offers the advantage by detecting DNA from the host rather than its fecal bacteria. The present study focused on sampling water presumably from non-point sources (NPS), which can increase bacterial and nitrogen concentrations to receiving water bodies. Stream sampling sites located within the Piscataqua River Watershed (PRW), New Hampshire, USA, were sampled from a range of sites that experienced nitrogen inputs such as sewer and septic systems and suburban runoff. Three mitochondrial (mtDNA) gene marker assays (human, bovine, and canine) were tested from surface water. Nineteen sites were sampled during an 18-month period. Analyses of the combined single and multiplex assay results showed that the proportion of occurrence was highest for bovine (15.6%; n = 77) compared to canine (5.6%; n = 70) and human (5.7%; n = 107) mtDNA gene markers. For the human mtDNA marker, there was a statistically significant relationship between presence vs. absence and land use (Fisher's test p = 0.0031). This result was evident particularly for rural suburban septic, which showed the highest proportion of presence (19.2%) compared to the urban sewered (3.3%), suburban sewered (0%), and agricultural (0%) as well as forested septic (0%) sites. Although further testing across varied land use is needed, our study provides evidence for using the mtDNA marker in large watersheds.

Codon usage bias and phylogenetic analysis of mitochondrial ND1 gene in pisces, aves, and mammals

The mitochondrially encoded NADH:ubiquinone oxidoreductase core subunit 1 (MT-ND1) gene is a subunit of the respiratory chain complex I and involved in the first step of the electron transport chain of oxidative phosphorylation (OXPHOS). To understand the pattern of compositional properties, codon usage and expression level of mitochondrial ND1 genes in pisces, aves, and mammals, we used bioinformatic approaches as no work was reported earlier. In this study, a perl script was used for calculating nucleotide contents and different codon usage bias parameters. The codon usage bias of MT-ND1 was low but the expression level was high as revealed from high ENC and CAI value. Correspondence analysis (COA) suggests that the pattern of codon usage for MT-ND1 gene is not same across species and that compositional constraint played an important role in codon usage pattern of this gene among pisces, aves, and mammals. From the regression equation of GC12 on GC3, it can be inferred that the natural selection might have played a dominant role while mutation pressure played a minor role in influencing the codon usage patterns. Further, ND1 gene has a discrepancy with cytochrome B (CYB) gene in preference of codons as evident from COA. The codon usage bias was low. It is influenced by nucleotide composition, natural selection, mutation pressure, length (number) of amino acids, and relative dinucleotide composition. This study helps in understanding the molecular biology, genetics, evolution of MT-ND1 gene, and also for designing a synthetic gene.

The relationship between mitochondrial DNA haplotype and the reproductive capacity of domestic pigs (Sus scrofa domesticus)

Background

The maternally inherited mitochondrial genome encodes key proteins of the electron transfer chain, which produces the vast majority of cellular ATP. Mitochondrial DNA (mtDNA) present in the mature oocyte acts as a template for all mtDNA that is replicated during development to meet the specific energy requirements of each tissue. Individuals that share a maternal lineage cluster into groupings known as mtDNA haplotypes. MtDNA haplotypes confer advantages and disadvantages to an organism and this affects its phenotype. In livestock, certain mtDNA haplotypes are associated with improved milk and meat quality, whilst, other species, mtDNA haplotypes have shown increased longevity, growth and susceptibility to diseases. In this work, we have set out to determine whether mtDNA haplotypes influence reproductive capacity. This has been undertaken using a pig model.

Results

To determine the genetic diversity of domestic pigs in Australia, we have sequenced the D-loop region of 368 pigs, and identified five mtDNA haplotypes (A to E). To assess reproductive capacity, we compared oocyte maturation, fertilization and development to blastocyst, and found that there were significant differences for maturation and fertilization amongst the haplotypes. We then determined that haplotypes C, D and E produced significantly larger litters. When we assessed the conversion of developmentally competent oocytes and their subsequent developmental stages to offspring, we found that haplotypes A and B had the lowest reproductive efficiencies. Amongst the mtDNA haplotypes, the number of mtDNA variants harbored at >25 % correlated with oocyte quality. MtDNA copy number for developmentally competent oocytes positively correlated with the level of the 16383delC variant. This variant is located in the conserved sequence box II, which is a regulatory region for mtDNA transcription and replication.

Conclusions

We have identified five mtDNA haplotypes in Australian domestic pigs indicating that genetic diversity is restricted. We have also shown that there are differences in reproductive capacity amongst the mtDNA haplotypes. We conclude that mtDNA haplotypes affect pig reproductive capacity and can be used as a marker to complement current selection methods to identify productive pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0375-4) contains supplementary material, which is available to authorized users.

Molecular evolution of mitochondrial coding genes in the oxidative phosphorylation pathway in malacostraca: purifying selection or accelerated evolution?

The mitochondrion is the energy-producing factory of eukaryotic cells, in which oxidative phosphorylation (OXPHOS) is the main pathway for the production of adenosine triphosphate (ATP) by cellular respiration. Because of their vital role in metabolism, mitochondrial proteins are predicted to evolve primarily under constant purifying selection. However, all mitochondrial coding genes of malacostraca had a significantly higher synonymous nt divergence (Ks) in this study. Complex I (NADH dehydrogenase) and complex V (ATP synthase) had a much higher ratio of non-synonymous to synonymous nt divergence (Ka/Ks) and non-synonymous diversity (πNS), whereas complex III (cytochrome bc₁ complex) and complex IV (cytochrome c oxidase) had a significantly lower Ka/Ks and non-synonymous diversity (πNS). The Ka/Ks, πNS, πS, and Ka results revealed that two types of mitochondrial genes, NADH dehydrogenase and ATP synthase, in malacostraca were consistent with accelerated evolution. Furthermore, two other types of mitochondrial genes, cytochrome bc₁ complex and cytochrome c oxidase, were consistent with purifying selection. Generally, the evolutionary pattern of all mitochondrial proteins of the OXPHOS pathway in malacostraca was not entirely consistent with purifying selection.

A Conventional Multiplex PCR Assay for the Detection of Toxic Gemfish Species (Ruvettus pretiosus and Lepidocybium flavobrunneum): A Simple Method To Combat Health Frauds

The meat of Ruvettus pretiosus and Lepidocybium flavobrunneum (gemfishes) contains high amounts of indigestible wax esters that provoke gastrointestinal disorders. Although some countries have banned the sale of these species, mislabeling cases have been reported in sushi catering. This work developed a simple conventional multiplex PCR, which discriminates the two toxic gemfishes from other potentially replaced species, such as tunas, cod, and sablefish. A common degenerate forward primer and three species-specific reverse primers were designed to amplify cytochrome oxidase subunit I (COI) gene regions of different lengths (479, 403, and 291 bp) of gemfishes, tunas, and sablefish, respectively. A primer pair was designed to amplify a fragment (193 bp) of the cytb gene of cod species. Furthermore, a primer pair targeting the 16S rRNA gene was intended as common positive control (115 bp). The method developed in this study, by producing the expected amplicon for all of the DNA samples tested (reference and commercial), provides a rapid and reliable response in identifying the two toxic species to combat health frauds.

FMiR: A Curated Resource of Mitochondrial DNA Information for Fish

Mitochondrial genome sequences have been widely used for evolutionary and phylogenetic studies. Among vertebrates, fish are an important, diverse group, and their mitogenome sequences are growing rapidly in public repositories. To facilitate mitochondrial genome analysis and to explore the valuable genetic information, we developed the Fish Mitogenome Resource (FMiR) database to provide a workbench for mitogenome annotation, species identification and microsatellite marker mining. The microsatellites are also known as simple sequence repeats (SSRs) and used as molecular markers in studies on population genetics, gene duplication and marker assisted selection. Here, easy-to-use tools have been implemented for mining SSRs and for designing primers to identify species/habitat specific markers. In addition, FMiR can analyze complete or partial mitochondrial genome sequence to identify species and to deduce relational distances among sequences across species. The database presently contains curated mitochondrial genomes from 1302 fish species belonging to 297 families and 47 orders reported from saltwater and freshwater ecosystems. In addition, the database covers information on fish species such as conservation status, ecosystem, family, distribution and occurrence downloaded from the FishBase and IUCN Red List databases. Those fish information have been used to browse mitogenome information for the species belonging to a particular category. The database is scalable in terms of content and inclusion of other analytical modules. The FMiR is running under Linux operating platform on high performance server accessible at URL http://mail.nbfgr.res.in/fmir.

Mitochondrial DNA Fragmentation as a Molecular Tool to Monitor Thermal Processing of Plant-Derived, Low-Acid Foods, and Biomaterials

Cycle threshold (Ct) increase, quantifying plant-derived DNA fragmentation, was evaluated for its utility as a time-temperature integrator. This novel approach to monitoring thermal processing of fresh, plant-based foods represents a paradigm shift. Instead of using quantitative polymerase chain reaction (qPCR) to detect pathogens, identify adulterants, or authenticate ingredients, this rapid technique was used to quantify the fragmentation of an intrinsic plant mitochondrial DNA (mtDNA) gene over time-temperature treatments. Universal primers were developed which amplified a mitochondrial gene common to plants (atp1). These consensus primers produced a robust qPCR signal in 10 vegetables, 6 fruits, 3 types of nuts, and a biofuel precursor. Using sweet potato (Ipomoea batatas) puree as a model low-acid product and simple linear regression, Ct value was highly correlated to time-temperature treatment (R(2) = 0.87); the logarithmic reduction (log CFU/mL) of the spore-forming Clostridium botulinum surrogate, Geobacillus stearothermophilus (R(2) = 0.87); and cumulative F-value (min) in a canned retort process (R(2) = 0.88), all comparisons conducted at 121 °C. D121 and z-values were determined for G. stearothermophilus ATCC 7953 and were 2.71 min and 11.0 °C, respectively. D121 and z-values for a 174-bp universal plant amplicon were 11.3 min and 9.17 °C, respectively, for mtDNA from sweet potato puree. We present these data as proof-of-concept for a molecular tool that can be used as a rapid, presumptive method for monitoring thermal processing in low-acid plant products.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

Ecotype evolution in Glossina palpalis subspecies, major vectors of sleeping sickness

Background

The role of environmental factors in driving adaptive trajectories of living organisms is still being debated. This is even more important to understand when dealing with important neglected diseases and their vectors.

Methodology/Principal Findings

In this paper, we analysed genetic divergence, computed from seven microsatellite loci, of 614 tsetse flies (Glossina palpalis gambiensis and Glossina palpalis palpalis, major vectors of animal and human trypanosomes) from 28 sites of West and Central Africa. We found that the two subspecies are so divergent that they deserve the species status. Controlling for geographic and time distances that separate these samples, which have a significant effect, we found that G. p. gambiensis from different landscapes (Niayes of Senegal, savannah and coastal environments) were significantly genetically different and thus represent different ecotypes or subspecies. We also confirm that G. p. palpalis from Ivory Coast, Cameroon and DRC are strongly divergent.

Conclusions/Significance

These results provide an opportunity to examine whether new tsetse fly ecotypes might display different behaviour, dispersal patterns, host preferences and vectorial capacities. This work also urges a revision of taxonomic status of Glossina palpalis subspecies and highlights again how fast ecological divergence can be, especially in host-parasite-vector systems.

The role of environmental factors in driving adaptive trajectories of living organisms is still being debated. This is even more important to understand when dealing with important and /or neglected diseases and their vectors. In this paper, we analysed genetic divergence, computed from several genetic markers, of 614 tsetse flies (Glossina palpalis gambiensis and Glossina palpalis palpalis, major vectors of animal and human trypanosomes) from 28 sites of West and Central Africa. We found that the two subspecies are so divergent that they deserve the species status. We found that G. p. gambiensis from different landscapes (Niayes of Senegal, savannah and coastal environments) were significantly genetically different, and thus represent different adaptive entities or even subspecies. We also confirm that G. p. palpalis from Ivory Coast, Cameroon and DRC are strongly divergent. These results provide an opportunity to examine whether these different types of tsetse fly might display different behaviour, dispersal patterns, host preferences and vectorial capacities. This work also urges a revision of taxonomic status of Glossina palpalis subspecies and highlights again how fast ecological divergence can be, especially in host-parasite-vector systems.

Patterns of positive selection in seven ant genomes

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide data sets in bees and flies with the same methodology to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

Assortative mating between two sympatric closely-related specialists: inferred from molecular phylogenetic analysis and behavioral data

Host plant shifting of phytophagous insects can lead to the formation of host associated differentiation and ultimately speciation. In some cases, host plant specificity alone acts as a nearly complete pre-mating isolating barrier among insect populations. We here test whether effective pre-mating isolation and host-independent behavioral isolation have evolved under the condition of extreme host specilization using two sympatric flea beetles with incomplete post-mating isolation under laboratory conditions. Phylogenetic analysis and coalescent simulation results showed that there is a limited interspecific gene flow, indicating effctive isolation between these species. Three types of mating tests in the absence of host plant cues showed that strong host-independent behavioral isolation has evolved between them. We conclude that almost perfect assortative mating between these two extreme host specialists results from a combination of reduced encounter rates due to differential host preference and strong sexual isolation.

Applications of DNA barcoding to fish landings: authentication and diversity assessment

DNA barcoding methodologies are being increasingly applied not only for scientific purposes but also for diverse real-life uses. Fisheries assessment is a potential niche for DNA barcoding, which serves for species authentication and may also be used for estimating within-population genetic diversity of exploited fish. Analysis of single-sequence barcodes has been proposed as a shortcut for measuring diversity in addition to the original purpose of species identification. Here we explore the relative utility of different mitochondrial sequences (12S rDNA, COI, cyt b, and D-Loop) for application as barcodes in fisheries sciences, using as case studies two marine and two freshwater catches of contrasting diversity levels. Ambiguous catch identification from COI and cyt b was observed. In some cases this could be attributed to duplicated names in databases, but in others it could be due to mitochondrial introgression between closely related species that may obscure species assignation from mtDNA. This last problem could be solved using a combination of mitochondrial and nuclear genes. We suggest to simultaneously analyze one conserved and one more polymorphic gene to identify species and assess diversity in fish catches.

The maintenance of mitochondrial genetic variation by negative frequency-dependent selection

Mitochondrial genes generally show high levels of standing genetic variation, which is puzzling given the accumulating evidence for phenotypic effects of mitochondrial genetic variation. Negative frequency-dependent selection, where the relative fitness of a genotype is inversely related to its frequency in a population, provides a potent and potentially general process that can maintain mitochondrial polymorphism. We assessed the change in mitochondrial haplotype frequencies over 10 generations of experimental evolution in 180 seed beetle populations in the laboratory, where haplotypes competed for propagation to subsequent generations. We found that haplotypes consistently increased in frequency when they were initially rare and decreased in frequency when initially common. Our results have important implications for the use of mtDNA haplotype frequency data to infer population level processes and they revive the general hypothesis that negative frequency-dependent selection, presumably caused by habitat heterogeneity, may commonly promote polymorphism in ecologically relevant life history genes.

Animal mitochondria, positive selection and cyto-nuclear coevolution: insights from pulmonates

Pulmonate snails have remarkably high levels of mtDNA polymorphism within species and divergence between species, making them an interesting group for the study of mutation and selection on mitochondrial genomes. The availability of sequence data from most major lineages – collected largely for studies of phylogeography - provides an opportunity to perform several tests of selection that may provide general insights into the evolutionary forces that have produced this unusual pattern. Several protein coding mtDNA datasets of pulmonates were analyzed towards this direction. Two different methods for the detection of positive selection were used, one based on phylogeny, and the other on the McDonald-Kreitman test. The cyto-nuclear coevolution hypothesis, often implicated to account for the high levels of mtDNA divergence of some organisms, was also addressed by assessing the divergence pattern exhibited by a nuclear gene. The McDonald-Kreitman test indicated multiple signs of positive selection in the mtDNA genes, but was significantly biased when sequence divergence was high. The phylogenetic method identified five mtDNA datasets as affected by positive selection. In the nuclear gene, the McDonald-Kreitman test provided no significant results, whereas the phylogenetic method identified positive selection as likely present. Overall, our findings indicate that: 1) slim support for the cyto-nuclear coevolution hypothesis is present, 2) the elevated rates of mtDNA polymorphims and divergence in pulmonates do not appear to be due to pervasive positive selection, 3) more stringent tests show that spurious positive selection is uncovered when distant taxa are compared and 4) there are significant examples of positive selection acting in some cases, so it appears that mtDNA evolution in pulmonates can escape from strict deleterious evolution suggested by the Muller’s ratchet effect.

Reconstructing the phylogenetic history of long-term effective population size and life-history traits using patterns of amino acid replacement in mitochondrial genomes of mammals and birds

The nearly neutral theory, which proposes that most mutations are deleterious or close to neutral, predicts that the ratio of nonsynonymous over synonymous substitution rates (dN/dS), and potentially also the ratio of radical over conservative amino acid replacement rates (Kr/Kc), are negatively correlated with effective population size. Previous empirical tests, using life-history traits (LHT) such as body-size or generation-time as proxies for population size, have been consistent with these predictions. This suggests that large-scale phylogenetic reconstructions of dN/dS or Kr/Kc might reveal interesting macroevolutionary patterns in the variation in effective population size among lineages. In this work, we further develop an integrative probabilistic framework for phylogenetic covariance analysis introduced previously, so as to estimate the correlation patterns between dN/dS, Kr/Kc, and three LHT, in mitochondrial genomes of birds and mammals. Kr/Kc displays stronger and more stable correlations with LHT than does dN/dS, which we interpret as a greater robustness of Kr/Kc, compared with dN/dS, the latter being confounded by the high saturation of the synonymous substitution rate in mitochondrial genomes. The correlation of Kr/Kc with LHT was robust when controlling for the potentially confounding effects of nucleotide compositional variation between taxa. The positive correlation of the mitochondrial Kr/Kc with LHT is compatible with previous reports, and with a nearly neutral interpretation, although alternative explanations are also possible. The Kr/Kc model was finally used for reconstructing life-history evolution in birds and mammals. This analysis suggests a fairly large-bodied ancestor in both groups. In birds, life-history evolution seems to have occurred mainly through size reduction in Neoavian birds, whereas in placental mammals, body mass evolution shows disparate trends across subclades. Altogether, our work represents a further step toward a more comprehensive phylogenetic reconstruction of the evolution of life-history and of the population-genetics environment.

Tough adults, frail babies: an analysis of stress sensitivity across early life-history stages of widely introduced marine invertebrates

All ontogenetic stages of a life cycle are exposed to environmental conditions so that population persistence depends on the performance of both adults and offspring. Most studies analysing the influence of abiotic conditions on species performance have focussed on adults, while studies covering early life-history stages remain rare. We investigated the responses of early stages of two widely introduced ascidians, Styela plicata and Microcosmus squamiger, to different abiotic conditions. Stressors mimicked conditions in the habitats where both species can be found in their distributional ranges and responses were related to the selection potential of their populations by analysing their genetic diversity. Four developmental stages (egg fertilisation, larval development, settlement, metamorphosis) were studied after exposure to high temperature (30°C), low salinities (26 and 22‰) and high copper concentrations (25, 50 and 100 µg/L). Although most stressors effectively led to failure of complete development (fertilisation through metamorphosis), fertilisation and larval development were the most sensitive stages. All the studied stressors affected the development of both species, though responses differed with stage and stressor. S. plicata was overall more resistant to copper, and some stages of M. squamiger to low salinities. No relationship was found between parental genetic composition and responses to stressors. We conclude that successful development can be prevented at several life-history stages, and therefore, it is essential to consider multiple stages when assessing species' abilities to tolerate stress. Moreover, we found that early development of these species cannot be completed under conditions prevailing where adults live. These populations must therefore recruit from elsewhere or reproduce during temporal windows of more benign conditions. Alternatively, novel strategies or behaviours that increase overall reproductive success might be responsible for ensuring population survival.

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Evolutionary history of frogfishes (Teleostei: Lophiiformes: Antennariidae): a molecular approach

Fishes of the family Antennariidae (order Lophiiformes) are primarily shallow-water benthic forms found in nearly all tropical and subtropical oceans and seas of the world, with some taxa extending into temperate waters. Despite an earlier attempt based on morphology, no previous hypothesis of intergeneric relationships of the Antennariidae exists. To resolve phylogenetic relationships within the Antennariidae, and to test the validity of species groups within Antennarius, DNA sequences from the mitochondrial 16S and cytochrome oxidase c subunit 1 (COI) genes, and nuclear recombination activating gene 2 (RAG2), for 25 described and four undescribed antennariid species, representing 10 of 12 known genera and one undescribed genus, were unambiguously aligned and analyzed using Bayesian and maximum likelihood methods. The markers were partitioned and analyzed for substitution saturation and only the third codon position of COI (COI-3) was found to have reached saturation. However, analysis of both datasets, one with the saturated data and one without, differed only slightly. All molecular analyses recovered two major clades, one comprised of Fowlerichthys, Antennarius, Histrio, and Antennatus; and another containing Rhycherus, Antennariidae gen. et sp. nov., Kuiterichthys, Phyllophryne, Echinophryne, Tathicarpus, Lophiocharon, and Histiophryne. Evidence is presented to illustrate a correlation between phylogeny, geographic distribution, and reproductive life history. The results of these analyses provide the first hypothesis of evolutionary relationships within the Antennariidae.

The whereabouts of an ancient wanderer: global phylogeography of the solitary ascidian Styela plicata

Genetic tools have greatly aided in tracing the sources and colonization history of introduced species. However, recurrent introductions and repeated shuffling of populations may have blurred some of the genetic signals left by ancient introductions. Styela plicata is a solitary ascidian distributed worldwide. Although its origin remains unclear, this species is believed to have spread worldwide by travelling on ship's hulls. The goals of this study were to infer the genetic structure and global phylogeography of S. plicata and to look for present-day and historical genetic patterns. Two genetic markers were used: a fragment of the mitochondrial gene Cytochrome Oxidase subunit I (COI) and a fragment of the nuclear gene Adenine Nucleotide Transporter/ADP-ATP Translocase (ANT). A total of 368 individuals for COI and 315 for ANT were sequenced from 17 locations worldwide. The levels of gene diversity were moderate for COI to high for ANT. The Mediterranean populations showed the least diversity and allelic richness for both markers, while the Indian, Atlantic and Pacific Oceans had the highest gene and nucleotide diversities. Network and phylogenetic analyses with COI and ANT revealed two groups of alleles separated by 15 and 4 mutational steps, respectively. The existence of different lineages suggested an ancient population split. However, the geographic distributions of these groups did not show any consistent pattern, indicating different phylogeographic histories for each gene. Genetic divergence was significant for many population-pairs irrespective of the geographic distance among them. Stochastic introduction events are reflected in the uneven distribution of COI and ANT allele frequencies and groups among many populations. Our results confirmed that S. plicata has been present in all studied oceans for a long time, and that recurrent colonization events and occasional shuffling among populations have determined the actual genetic structure of this species.

Searching for evidence of selection in avian DNA barcodes

The barcode of life project has assembled a tremendous number of mitochondrial cytochrome c oxidase I (COI) sequences. Although these sequences were gathered to develop a DNA-based system for species identification, it has been suggested that further biological inferences may also be derived from this wealth of data. Recurrent selective sweeps have been invoked as an evolutionary mechanism to explain limited intraspecific COI diversity, particularly in birds, but this hypothesis has not been formally tested. In this study, I collated COI sequences from previous barcoding studies on birds and tested them for evidence of selection. Using this expanded data set, I re-examined the relationships between intraspecific diversity and interspecific divergence and sampling effort, respectively. I employed the McDonald-Kreitman test to test for neutrality in sequence evolution between closely related pairs of species. Because amino acid sequences were generally constrained between closely related pairs, I also included broader intra-order comparisons to quantify patterns of protein variation in avian COI sequences. Lastly, using 22 published whole mitochondrial genomes, I compared the evolutionary rate of COI against the other 12 protein-coding mitochondrial genes to assess intragenomic variability. I found no conclusive evidence of selective sweeps. Most evidence pointed to an overall trend of strong purifying selection and functional constraint. The COI protein did vary across the class Aves, but to a very limited extent. COI was the least variable gene in the mitochondrial genome, suggesting that other genes might be more informative for probing factors constraining mitochondrial variation within species.

Conserved PCR primer set designing for closely-related species to complete mitochondrial genome sequencing using a sliding window-based PSO algorithm

Background

Complete mitochondrial (mt) genome sequencing is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. For long template sequencing, i.e., like the entire mtDNA, it is essential to design primers for Polymerase Chain Reaction (PCR) amplicons which are partly overlapping each other. The presented chromosome walking strategy provides the overlapping design to solve the problem for unreliable sequencing data at the 5′ end and provides the effective sequencing. However, current algorithms and tools are mostly focused on the primer design for a local region in the genomic sequence. Accordingly, it is still challenging to provide the primer sets for the entire mtDNA.

Methodology/Principal Findings

The purpose of this study is to develop an integrated primer design algorithm for entire mt genome in general, and for the common primer sets for closely-related species in particular. We introduce ClustalW to generate the multiple sequence alignment needed to find the conserved sequences in closely-related species. These conserved sequences are suitable for designing the common primers for the entire mtDNA. Using a heuristic algorithm particle swarm optimization (PSO), all the designed primers were computationally validated to fit the common primer design constraints, such as the melting temperature, primer length and GC content, PCR product length, secondary structure, specificity, and terminal limitation. The overlap requirement for PCR amplicons in the entire mtDNA is satisfied by defining the overlapping region with the sliding window technology. Finally, primer sets were designed within the overlapping region. The primer sets for the entire mtDNA sequences were successfully demonstrated in the example of two closely-related fish species. The pseudo code for the primer design algorithm is provided.

Conclusions/Significance

In conclusion, it can be said that our proposed sliding window-based PSO algorithm provides the necessary primer sets for the entire mt genome amplification and sequencing.

Application of mitochondrial genes sequences for measuring the genetic diversity of Arabian oryx

Arabian oryx (Oryx leucoryx) had faced extinction in the wild more than three decades ago and was saved by the prudent efforts of captive breeding programs. A clear understanding of the molecular diversity of contemporary Arabian oryx population is important for the long term success of captive breeding and reintroduction of this potentially endangered species. We have sequenced the segments of mitochondrial DNA including12S rRNA, 16S rRNA, cytochrome b (Cyt-b) and control region (CR) genes of 24 captive-bred and reintroduced animals. Although the sequences of 12S rRNA, 16S rRNA and Cyt-b were found to be identical for all the samples, typical sequence variations in the CR gene were observed in the form of 7 haplotypes. One of these haplotypes has been reported earlier while the remaining 6 haplotypes are novel and represent different lineages from the founders. The haplotype and nucleotide diversities were found to be 0.789 and 0.009 respectively. The genetic distances among the 7 mtDNA haplotypes varied from 0.001 to 0.017. These findings are of potential relevance to the management of captive breeding programs for the conservation of Arabian oryx.

High proportions of deleterious polymorphisms in constrained human genes

Previous studies on human mitochondrial genomes showed that the ratio of intra-specific diversities at nonsynonymous-to-synonymous positions was two to ten times higher than the ratio of interspecific divergences at these positions, suggesting an excess of slightly deleterious nonsynonymous polymorphisms. However, such an overabundance of nonsynonymous single nucleotide polymorphisms (SNPs) was not found in human nuclear genomes. Here, genome-wide estimates using >14,000 human-chimp nuclear genes and 1 million SNPs from four human genomes showed a significant proportion of deleterious nonsynonymous SNPs (∼ 15%). Importantly, this study reveals a negative correlation between the magnitude of selection pressure and the proportion of deleterious SNPs on human genes. The proportion of deleterious amino acid replacement polymorphisms is 3.5 times higher in genes under high purifying selection compared with that in less constrained genes (28% vs. 8%). These results are explained by differences in the extent of contribution of mildly deleterious mutations to diversity and substitution.

Nucleotide diversity of a ND5 fragment confirms that population expansion is the most suitable explanation for the mtDNA haplotype polymorphism of Drosophila subobscura

Results from mitochondria (mt) DNA restriction site analyses (RSAs) have revealed that wild populations of Drosophila subobscura are formed by two common (I and II) and some rare, often endemic, low-frequency haplotypes. In the study reported here, we analysed nucleotide diversity in a 942-bp fragment of the mtDNA ND5 gene in 48 D. subobscura individuals captured from three populations that showed haplotypes I, II or the less common ones, as well as in one additional individual belonging to D. guanche that was taken as an outgroup. RSAs and sequencing results were compared. The two approaches yielded similar nucleotide variability parameters, suggesting a consistency in the results obtained from mtDNA dynamics in natural populations of D. subobscura. Patterns of polymorphism at ND5 are most consistent with the hypothesis of population expansion after a bottleneck that may have occurred since the last glaciation or which may occur seasonally after the summer and winter. However, we cannot rule out that selection has a role in maintaining the two major haplotypes at intermediate frequencies in worldwide populations of D. subobscura.

Analysis of Canis mitochondrial DNA demonstrates high concordance between the control region and ATPase genes

Background

Phylogenetic studies of wild Canis species have relied heavily on the mitochondrial DNA control region (mtDNA CR) to infer species relationships and evolutionary lineages. Previous analyses of the CR provided evidence for a North American evolved eastern wolf (C. lycaon), that is more closely related to red wolves (C. rufus) and coyotes (C. latrans) than grey wolves (C. lupus). Eastern wolf origins, however, continue to be questioned. Therefore, we analyzed mtDNA from 89 wolves and coyotes across North America and Eurasia at 347 base pairs (bp) of the CR and 1067 bp that included the ATPase6 and ATPase8 genes. Phylogenies and divergence estimates were used to clarify the evolutionary history of eastern wolves, and regional comparisons of nonsynonomous to synonomous substitutions (dN/dS) at the ATPase6 and ATPase8 genes were used to elucidate the potential role of selection in shaping mtDNA geographic distribution.

Results

We found high concordance across analyses between the mtDNA regions studied. Both had a high percentage of variable sites (CR = 14.6%; ATP = 9.7%) and both phylogenies clustered eastern wolf haplotypes monophyletically within a North American evolved lineage apart from coyotes. Divergence estimates suggest the putative red wolf sequence is more closely related to coyotes (D_xyCR = 0.01982 ± 0.00494 SD; D_xyATP = 0.00332 ± 0.00097 SD) than the eastern wolf sequences (D_xyCR = 0.03047 ± 0.00664 SD; D_xyATP = 0.00931 ± 0.00205 SD). Neutrality tests on both genes were indicative of the population expansion of coyotes across eastern North America, and dN/dS ratios suggest a possible role for purifying selection in the evolution of North American lineages. dN/dS ratios were higher in European evolved lineages from northern climates compared to North American evolved lineages from temperate regions, but these differences were not statistically significant.

Conclusions

These results demonstrate high concordance between coding and non-coding regions of mtDNA, and provide further evidence that the eastern wolf possessed distinct mtDNA lineages prior to recent coyote introgression. Purifying selection may have influenced North American evolved Canis lineages, but detection of adaptive selection in response to climate is limited by the power of current statistical tests. Increased sampling and development of alternative analytical tools will be necessary to disentangle demographic history from processes of natural selection.

Correcting for purifying selection: an improved human mitochondrial molecular clock

There is currently no calibration available for the whole human mtDNA genome, incorporating both coding and control regions. Furthermore, as several authors have pointed out recently, linear molecular clocks that incorporate selectable characters are in any case problematic. We here confirm a modest effect of purifying selection on the mtDNA coding region and propose an improved molecular clock for dating human mtDNA, based on a worldwide phylogeny of > 2000 complete mtDNA genomes and calibrating against recent evidence for the divergence time of humans and chimpanzees. We focus on a time-dependent mutation rate based on the entire mtDNA genome and supported by a neutral clock based on synonymous mutations alone. We show that the corrected rate is further corroborated by archaeological dating for the settlement of the Canary Islands and Remote Oceania and also, given certain phylogeographic assumptions, by the timing of the first modern human settlement of Europe and resettlement after the Last Glacial Maximum. The corrected rate yields an age of modern human expansion in the Americas at approximately 15 kya that-unlike the uncorrected clock-matches the archaeological evidence, but continues to indicate an out-of-Africa dispersal at around 55-70 kya, 5-20 ky before any clear archaeological record, suggesting the need for archaeological research efforts focusing on this time window. We also present improved rates for the mtDNA control region, and the first comprehensive estimates of positional mutation rates for human mtDNA, which are essential for defining mutation models in phylogenetic analyses.

Conflicting mitochondrial and nuclear phylogenies for the widely disjunct Emys (Testudines: Emydidae) species complex, and what they tell us about biogeography and hybridization

Understanding the mechanisms by which widely disjunct members of a clade came to occupy their current distribution is one of the fundamental challenges of biogeography. Here, we used data from 7 nuclear and 1 mitochondrial gene to examine the phylogenetic and biogeographic history of Emys, a clade of turtles that is broadly disjunct in western and eastern North America and Europe. We found strong disagreement between mitochondrial and nuclear gene trees, with mitochondrial DNA supporting the monophyly of the North American taxa (marmorata + blandingii) to the exclusion of the European orbicularis, and nuclear genes supporting the monophyly of (blandingii + orbicularis) to the exclusion of marmorata. We used fossil-calibrated molecular chronograms, in combination with supporting evidence from the fossil record and paleoclimatology, to identify a potential example of ancient hybridization and mitochondrial gene capture 12 million years ago, which explains this discrepancy. Based on the weight of evidence, we argue that the invasion of Eurasia by Emys orbicularis occurred about 16 Ma via a trans-Beringian land bridge. The case of Emys emphasizes how single-gene trees can be strongly affected by population processes, including hybridization, and that the effects of these processes can persist through long periods of evolutionary history. Given the chaotic state of the current taxonomy of these turtles, our work also emphasizes the care that should be used in implementing taxonomic changes based on 1 or a few gene trees and the importance of taking a conservative approach in renaming or splitting higher taxa based on apparent nonmonophyly.

Inverse relationship between longevity and evolutionary rate of mitochondrial proteins in mammals and birds

Recently, an unexpected, positive correlation between the rate of evolution of mitochondrial proteins and longevity was reported. Here we re-analyze this relationship in various mammalian lineages using a bayesian phylogenetic analysis of amino-acid sequences, allowing for variable evolutionary rates across sites and species. A negative relationship between protein evolutionary rate and species longevity is reported for all oxidative phosphorylation complexes. A detailed analysis of the cytochrome b in 528 mammals reinforced this result, which contradicts previous publications. Reconducting the analysis in birds yielded similar results. We explain the discrepancy between this and previous reports by our improved taxon sampling and more appropriate methodology: unlike distance-based methods, the tree-based bayesian approach can take into account the high variation of substitution rate across amino-acid sites, and the resulting multiple substitution events. We discuss how our analysis contradicts Rottenberg's rationale, but does not dismiss his proposal of a longevity-dependent selective pressure on mitochondrial mutation rate in mammals and birds. This is because his interpretation invokes adaptation as the single evolutionary force at work, disregarding the effects of mutation, genetic drift, and purifying selection.

Domestic wastewater influent profiling using mitochondrial real-time PCR for source tracking animal contamination

Real-time PCR amplifying mammalian and avian mitochondrial DNA (mtDNA) was developed to characterize domestic and light industrial wastewater influent from two municipal wastewater treatment facilities (WWTF) over a 24-week period. Influent samples were assayed with species-specific primers and dual-labeled probes for human, bovine, swine, dog, cat, Canada goose and white-tailed deer to detect and quantify eukaryotic mtDNA contributors to local municipal wastewaters. Human (mean=9.6 x 10(4) copies/ml) and dog (mean=5.3 x 10(2) copies/ml) mtDNA were detected in all 24 samples (12 samples/site). Bovine and swine mtDNA were detected sporadically and at lower levels than human mtDNA, means=3.0 x 10(4) and 9.5 x 10(2) copies/ml, respectively. Domestic cat, Canada goose and white-tailed deer were detected only once in 24 samples. Mitochondrial DNA concentrations were compared to other bacterial, chemical and spectrophotometric parameters. Human mtDNA was positively correlated with ammonia concentration (P=0.01) and initial OD(600) reading (P=0.02) at one WWTF. Bovine mtDNA was positively correlated with biological oxygen demand (BOD) (P=0.02), final DNA concentration (P=0.03), initial and final humic acid concentrations (P=0.01, P=0.01), and final OD(600) (P=0.03) at one WWTF and total suspended solids (TSS) (P=0.04, P=0.09) at both facilities. Fecal coliforms were not positively or negatively correlated with mtDNA concentrations of any species assayed. For source tracking purposes, a combination of human (10(5) copies/ml) and dog mtDNA signal (10(2) copies/ml) could be indicative of municipal domestic wastewater contamination of environmental waters.

Bayesian, maximum parsimony and UPGMA models for inferring the phylogenies of antelopes using mitochondrial markers

This investigation was aimed to compare the inference of antelope phylogenies resulting from the 16S rRNA, cytochrome-b (cyt-b) and d-loop segments of mitochondrial DNA using three different computational models including Bayesian (BA), maximum parsimony (MP) and unweighted pair group method with arithmetic mean (UPGMA). The respective nucleotide sequences of three Oryx species (Oryx leucoryx, Oryx dammah and Oryx gazella) and an out-group (Addax nasomaculatus) were aligned and subjected to BA, MP and UPGMA models for comparing the topologies of respective phylogenetic trees. The 16S rRNA region possessed the highest frequency of conserved sequences (97.65%) followed by cyt-b (94.22%) and d-loop (87.29%). There were few transitions (2.35%) and none transversions in 16S rRNA as compared to cyt-b (5.61% transitions and 0.17% transversions) and d-loop (11.57% transitions and 1.14% transversions) while comparing the four taxa. All the three mitochondrial segments clearly differentiated the genus Addax from Oryx using the BA or UPGMA models. The topologies of all the gamma-corrected Bayesian trees were identical irrespective of the marker type. The UPGMA trees resulting from 16S rRNA and d-loop sequences were also identical (Oryx dammah grouped with Oryx leucoryx) to Bayesian trees except that the UPGMA tree based on cyt-b showed a slightly different phylogeny (Oryx dammah grouped with Oryx gazella) with a low bootstrap support. However, the MP model failed to differentiate the genus Addax from Oryx. These findings demonstrate the efficiency and robustness of BA and UPGMA methods for phylogenetic analysis of antelopes using mitochondrial markers.

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.

Extreme primary and secondary protein structure variability in the chimeric male-transmitted cytochrome c oxidase subunit II protein in freshwater mussels: evidence for an elevated amino acid substitution rate in the face of domain-specific purifying selection

BACKGROUND

Freshwater unionoidean bivalves, and species representing two marine bivalve orders (Mytiloida and Veneroida), exhibit a mode of mtDNA inheritance involving distinct maternal (F) and paternal (M) transmission routes concomitant with highly divergent gender-associated mtDNA genomes. Additionally, male unionoidean bivalves have a approximately 550 bp 3' coding extension to the cox2 gene (Mcox2e), that is apparently absent from all other metazoan taxa.

RESULTS

Our molecular sequence analyses of MCOX2e indicate that both the primary and secondary structures of the MCOX2e region are evolving much faster than other regions of the F and M COX2-COX1 gene junction. The near N-terminus approximately 2/3 of the MCOX2e region contains an interspecifically variable number of predicted transmembrane helices (TMH) and interhelical loops (IHL) whereas the C-terminus approximately 1/3 is relatively conserved and hydrophilic while containing conserved functional motifs. MCOX2e displays an overall pattern of purifying selection that leads to the preservation of TMH/IHL and C-terminus tail sub-regions. However, 14 amino acid positions in the MCOX2e TMH/IHL sub-region might be targeted by diversifying selection, each representing a site where there exists interspecific variation for the constituent amino acids residing in a TMH or IHL.

CONCLUSION

Our results indicate that Mcox2e is unique to unionoidean bivalves, likely the result of a single insertion event that took place over 65 MYA and that MCOX2e is functional. The predicted TMH number, length and position variability likely stems from substitution-based processes rather than the typically implicated insertion/deletion events. MCOX2e has relatively high rates of primary and secondary structure evolution, with some amino acid residues potentially subjected to site-specific positive selection, yet an overall pattern of purifying selection leading to the preservation of the TMH/IHL and hydrophilic C-terminus tail subregions. The more conserved C-terminus tail (relative to the TMH/IHL sub-region of MCOX2e) is likely biologically active because it contains functional motifs. The rapid evolution of primary and secondary structure in MCOX2e, combined with the action of both positive and purifying selection, provide supporting evidence for the hypothesis that MCOX2e has a novel reproductive function within unionoidean bivalves. All tolled, our data indicate that unionoidean bivalve MCOX2 is the first reported chimeric animal mtDNA-encoded protein.