Identification of mosquito bloodmeals using mitochondrial cytochrome oxidase subunit I and cytochrome b gene sequences

Abstract.  Primer pairs were designed and protocols developed to selectively amplify segments of vertebrate mitochondrial cytochrome oxidase subunit 1 (COI) and cytochrome b (Cyt b ) mtDNA from the bloodmeals of mosquitoes (Diptera: Culicidae). The protocols use two pairs of nested COI primers and one pair of Cyt b primers to amplify short segments of DNA. Resultant sequences are then compared with sequences in GenBank, using the BLAST function, for putative host identification. Vertebrate DNA was amplified from 88% of our sample of 162 wild-caught, blood-fed mosquitoes from Oregon, U.S.A. and GenBank BLAST searches putatively identified 98% of the amplified sequences, including one amphibian, seven mammalian and 14 avian species. Criteria and caveats for putative identification of bloodmeals are discussed.     

Calibrating the molecular clock beyond cytochrome b: assessing the evolutionary rate of COI in birds

Estimating the age of species or their component lineages based on sequence data is crucial for many studies in avian evolutionary biology. Although calibrations of the molecular clock in birds have been performed almost exclusively using cytochrome b (cyt b), they are commonly extrapolated to other mitochondrial genes. The existence of a large, standardized cytochrome c oxidase subunit I (COI) library generated as a result of the DNA barcoding initiative provides the opportunity to obtain a calibration for this mitochondrial gene in birds. In this study we compare the evolutionary rate of COI relative to cyt b across ten different avian orders. We obtained divergence estimates for both genes from nearly 300 phylogenetically independent pairs of species through the analysis of almost 5000 public sequences. For each pair of species we calculated the difference in divergence between COI and cyt b. Our results indicate that COI evolves on average 14% slower than cyt b, but also reveal considerable variation both among and within avian orders, precluding the use of this value as a standard adjustment for the COI molecular clock for birds. Our findings suggest that this variation is partially explained by a clear negative relationship between the difference in divergence in these genes and the age of species. Distances for cyt b are higher than those for COI for closely related species, but the values become similar as the divergence between the species increases. This appears to be the result of a stronger pattern of negative time‐dependency in the rate of cyt b than in that of COI, a difference that could be related to lower functional constraints on a small number of sites in cyt b that allow it to initially accumulate mutations more rapidly than COI.     

Cryptic species and morphological plasticity in long-lived bivalves (Unionoida: Hyriidae) from inland Australia

Molecular (mitochondrial DNA, isozyme) and morphological diversity of freshwater mussels (Family Hyriidae) was examined at 21 sites encompassing four large river systems, across southwest Queensland, Australia. Evidence was found for two major morphological groups. One group, which occurred in every river system, closely matched a recognized species (Velesunio ambiguus) both morphologically and in a well-supported lineage within a mitochondrial phylogeny generated from partial cytochrome c oxidase subunit I (COI) sequences. The second group most closely matched Velesunio wilsonii in shell morphology but formed three deeply divergent mitochondrial DNA lineages. All four lineages occurred sympatrically in some areas and displayed corresponding fixed differences at nuclear allozyme loci, which suggests an absence of recent hybridization and the presence of separate species.     

Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

Type I STS markers are more informative than cytochrome B in phylogenetic reconstruction of the Mustelidae (Mammalia: Carnivora)

We compared the utility of five nuclear gene segments amplified with type I sequence-tagged site (STS) primers versus the complete mitochondrial cytochrome b (cyt b) gene in resolving phylogenetic relationships within the Mustelidae, a large and ecomorphologically diverse family of mammalian carnivores. Maximum parsimony and likelihood analyses of separate and combined data sets were used to address questions regarding the levels of homoplasy, incongruence, and information content within and among loci. All loci showed limited resolution in the separate analyses because of either a low amount of informative variation (nuclear genes) or high levels of homoplasy (cyt b). Individually or combined, the nuclear gene sequences had less homoplasy, retained more signal, and were more decisive, even though cyt b contained more potentially informative variation than all the nuclear sequences combined. We obtained a well-resolved and supported phylogeny when the nuclear sequences were combined. Maximum likelihood and Bayesian phylogenetic analyses of the total combined data (nuclear and mitochondrial DNA sequences) were able to better accommodate the high levels of homoplasy in the cyt b data than was an equally weighted maximum parsimony analysis. Furthermore, partition Bremer support analyses of the total combined tree showed that the relative support of the nuclear and mitochondrial genes differed according to whether or not the homoplasy in the cyt b gene was downweighted. Although the cyt b gene contributed phylogenetic signal for most major groupings, the nuclear gene sequences were more effective in reconstructing the deeper nodes of the combined tree in the equally weighted parsimony analysis, as judged by the variable-length bootstrap method. The total combined data supported the monophyly of the Lutrinae (otters), whereas the Melinae (badgers) and Mustelinae (weasels, martens) were both paraphyletic. The American badger, Taxidea taxus (Taxidiinae), was the most basal taxon. Because hundreds of type I STS primer sets spanning the complete genomes of the human and mouse have been published and thus represent many independently segregating loci, the potential utility of these markers for molecular systematics of mammals and other groups is enormous.     

Novel primers for complete mitochondrial cytochrome b gene sequencing in mammals

Sequence-based species identification relies on the extent and integrity of sequence data available in online databases such as GenBank. When identifying species from a sample of unknown origin, partial DNA sequences obtained from the sample are aligned against existing sequences in databases. When the sequence from the matching species is not present in the database, high-scoring alignments with closely related sequences might produce unreliable results on species identity. For species identification in mammals, the cytochrome b (cyt b) gene has been identified to be highly informative; thus, large amounts of reference sequence data from the cyt b gene are much needed. To enhance availability of cyt b gene sequence data on a large number of mammalian species in GenBank and other such publicly accessible online databases, we identified a primer pair for complete cyt b gene sequencing in mammals. Using this primer pair, we successfully PCR amplified and sequenced the complete cyt b gene from 40 of 44 mammalian species representing 10 orders of mammals. We submitted 40 complete, correctly annotated, cyt b protein coding sequences to GenBank. To our knowledge, this is the first single primer pair to amplify the complete cyt b gene in a broad range of mammalian species. This primer pair can be used for the addition of new cyt b gene sequences and to enhance data available on species represented in GenBank. The availability of novel and complete gene sequences as high-quality reference data can improve the reliability of sequence-based species identification.     

Comparison of evolutionary rates in the mitochondrial DNA cytochrome b gene and control region and their implications for phylogeny of the Cobitoidea (Teleostei: Cypriniformes)

It is widely accepted that mitochondrial DNA (mtDNA) control region evolves faster than protein encoding genes with few exceptions. In the present study, we sequenced the mitochondrial cytochrome b gene (cyt b) and control region (CR) and compared their rates in 93 specimens representing 67 species of loaches and some related taxa in the Cobitoidea (Order Cypriniformes). The results showed that sequence divergences of the CR were broadly higher than those of the cyt b (about 1.83 times). However, in considering only closely related species, CR sequence evolution was slower than that of cyt b gene (ratio of CR/cyt b is 0.78), a pattern that is found to be very common in Cypriniformes. Combined data of the cyt b and CR were used to estimate the phylogenetic relationship of the Cobitoidea by maximum parsimony, neighbor-joining, and Bayesian methods. With Cyprinus carpio and Danio rerio as outgroups, three analyses identified the same four lineages representing four subfamilies of loaches, with Botiinae on the basal-most clade. The phylogenetic relationship of the Cobitoidea was ((Catostomidae+Gyrinocheilidae)+(Botiinae+(Balitorinae+(Cobitinae+Nemacheilinae)))), which indicated that Sawada's Cobitidae (including Cobitinae and Botiinae) was not monophyletic. Our molecular phylogenetic analyses are in very close agreement with the phylogenetic results based on the morphological data proposed by Nalbant and Bianco, wherein these four subfamilies were elevated to the family level as Botiidae, Balitoridae, Cobitidae, and Nemacheilidae.     

Multigene Barcoding and Phylogeny of Geographically Widespread Muricids (Gastropoda: Neogastropoda) Along the Coast of China

The identification and phylogeny of muricids have been in a state of confusion for a long time due to the morphological convergence and plasticity. DNA-based identification and phylogeny methods often offer an analytically powerful addition or even an alternative. In this study, we employ a DNA barcoding method to identify 17 known and easily confused muricid species (120 individuals) from the whole China coast based on mitochondrial cytochrome c oxidase subunit I (COI) and 16S rRNA sequences, and nuclear ITS-1 and 28S rRNA sequences. The phylogeny of muricid subfamilies is also analysed based on all mitochondrial and nuclear sequences. The universal COI and 16S rRNA primers did not work broadly across the study group, necessitating the redesign of muricid specific COI and 16S rRNA primers in this paper. Our study demonstrates that COI gene is a suitable marker for barcoding muricids, which can distinguish all muricid species studied. Phylogenetic analysis of 16S rRNA, ITS-1 and 28S rRNA data also provide good support for the species resolution observed in COI data. The relationships of muricid subfamilies are resolved based on the separate and combined gene data that showed the monophyly of each the subfamilies Ergalataxinae, Rapaninae, Ocenebrinae and Muricinae, especially that Ergalataxinae did not fall within Rapaninae.     

Molecular markers to discriminate among four pest species of Helicoverpa (Lepidoptera: Noctuidae)

The four significant pest species in the Helicoverpa genus (H. armigera, H. assulta, H. punctigera and H. zea) are morphologically similar and can only be reliably distinguished through dissection of adult genitalia. Two partial regions of the mitochondrial DNA (mtDNA), the cytochrome oxidase subunit I (COI) and the cytochrome b (Cyt b) genes were amplified by PCR and digested with restriction endonucleases. The restriction patterns, generated by the endonucleases BstZ17I and HphI, demonstrated reliable differentiation of the four Helicoverpa pest species. This technique is fast, reliable and effective at distinguishing specimens irrespective of their life stages and offers support to conventional taxonomic differentiation based on morphological characters.     

Diversity and phylogeny of mitochondrial cytochrome B lineages from six morphospecies of avian Haemoproteus (Haemosporida: Haemoproteidae)

Species of Haemoproteus (Haemosporida: Haemoproteidae), avian haemosporidians, have traditionally been described based on morphology of their gametocytes and on limited experimental information on their vertebrate host specificity. We investigated to what extent the morphological species are represented by monophyletic groups based on DNA sequence data using 2 different fragment lengths of the cytochrome b (cyt. b) gene. Phylogenetic reconstructions of obtained cyt. b lineages from 6 morphospecies of Haemoproteus showed that all lineages formed monophyletic clusters matching the morphospecies. Comparing our data with a recently published study showed that this is not always the case; the morphospecies H. belopolskyi consists of 2 distinct clusters of lineages that apparently have converged in morphology. However, the overall broad congruence between the molecular and morphological clustering of lineages will facilitate the integration of the knowledge obtained by traditional and molecular parasitology. Mean between morphospecies variation was 10-fold higher than the within species variation (5.5% vs. 0.54%), suggesting that Haemoproteus lineages with a genetic differentiation >5% are expected to be morphologically differentiated in most cases. When investigate the utility of 2 different fragment sizes of the cyt. b gene, the partial, 479-bp, cyt. b protocol picked up all mitochondrial (mt)DNA lineages that are found when using the full cyt. b gene, 1073 bp, suggesting that this protocol is sufficient for identification of most mtDNA lineages. All of the mtDNA lineages were associated with unique alleles when amplification was possible at a nuclear locus, strengthening the hypothesis that the designation of lineages based on mtDNA is largely genome-wide representative. We, therefore, propose the use of a cyt. b fragment of this length as a standard gene fragment for a DNA bar-coding system for avian Haemoproteus species.     

Molecular phylogeny of avian genus Syrmaticus based on the mitochondrial cytochrome B gene and control region

Mitochondrial DNA cytochrome b (cyt b) and control region (CR) nucleotide sequences were used to study the molecular phylogeny of the genus Syrmaticus. We found that the substitution rates among the three codon positions of cyt b were heterogeneous and the transition-transversion ratio was highly biased. As to CR sequences of the genus, most variable sites were in the peripheral domains. All molecular phylogenetic trees based on the two genes showed that: 1) the Syrmaticus was monophyletic and included five species with the following cladistic relationship: (S. reevesii, (S. soemmerringii, (S. mikado, (S. humiae and S. ellioti)))). Using the TN genetic distance of cyt b, we inferred the divergence time of the five species according to putative molecular clock and found that values were largely in agreement with the geological scenarios. The origin and speciation processes of the studied group were inferred by combining molecular and biogeographical evidences.     

Molecular phylogeny of the chipmunks inferred from Mitochondrial cytochrome b and cytochrome oxidase II gene sequences

There are currently 25 recognized species of the chipmunk genus Tamias. In this study we sequenced the complete mitochondrial cytochrome b (cyt b) gene of 23 Tamias species. We analyzed the cyt b sequence and then analyzed a combined data set of cyt b along with a previous data set of cytochrome oxidase subunit II (COII) sequence. Maximum-likelihood was used to further test the fit of models of evolution to the cyt b data. Other sciurid cyt b sequence was added to examine the evolution of Tamias in the context of other sciurids. Relationships among Tamias species are discussed, particularly the possibility of a current sorting event among taxa of the southwestern United States and the extreme divergences among the three subgenera (Neotamias, Eutamias, and Tamias).     

Molecular phylogeny of glass sponges (Porifera, Hexactinellida): increased taxon sampling and inclusion of the mitochondrial protein-coding gene, cytochrome oxidase subunit I

Marine sponges of the class Hexactinellida (glass sponges) are among the most understudied groups of Porifera, and molecular approaches to investigating their evolution have only recently emerged. Although these first results appeared reliable as they largely corroborated morphology-based hypotheses, they were almost exclusively based on ribosomal RNA genes (rDNA) and should, therefore, be further tested with independent types of genetic data, such as protein-coding genes. To this end, we established the mitochondrial-encoded cytochrome oxidase subunit I gene (COI) as an additional marker, and conducted phylogenetic analyses on DNA- and amino-acid level, as well as a supermatrix analysis based on combined COI DNA and rDNA alignments. Furthermore, we increased taxon sampling compared to previous studies by adding seven additional species. The COI-based phylogenies were largely congruent with the rDNA-based phylogeny but suffered from poor bootstrap support for many nodes. However, addition of the COI sequences to the rDNA data set increased resolution of the overall molecular phylogeny. Thus, although obtaining COI sequences from glass sponges turned out to be quite challenging, this gene appears to be a valuable supplement to rDNA data for molecular evolutionary studies of this group. Some implications of our extended phylogeny for the evolution and systematics of Hexactinellida are discussed.     

Genetic diversity and phylogenetic analysis of blackbuck (Antilope cervicapra) in southern India

Blackbuck (Antilope cervicapra) is a threatened species endemic to the Indian subcontinent. Many populations of blackbuck are found in southern India. Populations of blackbuck are negatively affected in many places for various reasons, such as habitat destruction and poaching. Their range decreased sharply during the 20th century. There is very limited information available on the population dynamics of blackbuck in southern India. For the phylogenetic and genetic diversity analyses of blackbuck populations among different distribution ranges in southern India, we sequenced mt DNA of cytochrome b (Cyt b) for 120, cytochrome c oxidase subunit-1 (COI) for 137 and the control region (CR) for 137 fecal pellets from eleven different locations in southern India. We analyzed the genetic structure of three mitochondrial markers, the CR, Cyt b and the COI region, separately and in a combined dataset. The haplotype diversity and nucleotide diversity of CR were 0.969 and 0.047, respectively, and were higher than those of Cyt b and COI. A Bayesian phylogeny and an MJ network based on the CR and combined dataset (105 sequences) signified several distinct haplotype clusters within blackbuck, whereas no clusters were identified with the Cyt b and COI phylogenetic analyses. The analysis of molecular variance of the combined data set revealed 52.46% genetic variation within the population. Mismatch distribution analysis revealed that blackbuck populations underwent complex changes with analysis of the combined dataset in each population and analysis of each marker separately in the overall population. The results provide evidence that blackbuck in different geographic locations has a distinct population structure due to habitat fragmentation after the formation of the Western and Eastern Ghats.

     

AFLPs resolve phylogeny and reveal mitochondrial introgression within a species flock of African electric fish (Mormyroidea: Teleostei)

Estimating species phylogeny from a single gene tree can be especially problematic for studies of species flocks in which diversification has been rapid. Here we compare a phylogenetic hypothesis derived from cytochrome b (cyt b) sequences with another based on amplified fragment length polymorphisms (AFLP) for 60 specimens of a monophyletic riverine species flock of mormyrid electric fishes collected in Gabon, west-central Africa. We analyze the aligned cyt b sequences by Wagner parsimony and AFLP data generated from 10 primer combinations using neighbor-joining from a Nei-Li distance matrix, Wagner parsimony, and Dollo parsimony. The different analysis methods yield AFLP tree topologies with few conflicting nodes. Recovered basal relationships in the group are similar between cyt b and AFLP analyses, but differ substantially at many of the more derived nodes. More of the clades recovered with the AFLP characters are consistent with the morphological characters used to designate operational taxonomic units in this group. These results support our hypothesis that the mitochondrial gene tree differs from the overall species phylogeny due at least in part to mitochondrial introgession among lineages. Mapping the two forms of electric organ found in this group onto the AFLP tree suggests that posteriorly innervated electrocytes with nonpenetrating stalks have independently evolved from anteriorly innervated, penetrating-stalk electrocytes at least three times.     

Species identification in the taxonomically neglected, highly diverse, neotropical parasitoid wasp genus Notiospathius (Braconidae: Doryctinae) based on an integrative molecular and morphological approach

Various DNA sequence-based methods for species delineation have recently been developed to assess the species-richness of highly diverse, neglected invertebrate taxa. These methods, however, need to be tested under a variety of conditions, including the use of different markers and parameters. Here, we explored the species diversity of a species-rich group of braconid parasitoid wasps, the Neotropical genus Notiospathius, including 233 specimens from 10 different countries. We examined sequences of two mitochondrial (mt) (COI, cyt b) and one nuclear (wg) gene fragments. We analysed them separately as well as concatenating the mt data with the general mixed Yule-coalescent (GMYC) model for species delineation using different tree-building methods and parameters for reconstructing ultrametric trees. We evaluated the performance of GMYC analyses by comparing their species delineations with our morphospecies identifications. Reconstructing ultrametric trees with a relaxed lognormal clock rate using the program BEAST gave the most congruent results with morphology for the two mt markers. A tree obtained with wg using the programs MrBayes+Pathd8 had the fewest cases of incongruence with morphology, though the performance of this nuclear marker was considerably lower than that of COI and cyt b. Species delimitation using the coalescent prior to obtain ultrametric trees was morphologically more congruent with COI, whereas the Yule prior was more congruent with cyt b. The analyses concatenating the mt datasets failed to recover some species supported both by morphology and the separate analyses of the mt markers. The highest morphological congruence was obtained with the GMYC analysis on an ultrametric tree reconstructed with cyt b using the relaxed lognormal clock rate and the Yule prior, thus supporting the importance of using alternative markers when the information of the barcoding locus (COI) is not concordant with morphological evidence. Seventy-one species were delimited based on the congruence found among COI, cyt b and morphology. Both mt markers also revealed the existence of seven potential cryptic species. This high species richness from a scattered geographical sampling indicates that there is a remarkable number of Notiospathius species that remains undiscovered.     

Characterization of a new membrane-bound cytochrome c of Rhodopseudomonas capsulata

A cytochrome c (cyt. c) was solubilized with Triton-X-100 and co-purified with cytochrome c oxidase from membranes of chemotrophically grown cells of Rhodopseudomonas capsulata. Cyt. c and cytochrome oxidase were separated on Sephadex G-50 columns. Antibodies against cytochrome c2 from the same bacterium did not cross react with the membrane-bound cyt. c. The IEP of the membrane-bound cyt. c was found to be pH 8.2, the midpoint potential was 234 +/- 11 mV at pH 7.0. This cyt. c binds CO. The native cyt. c is a dimer with an apparent Mr of 25000 containing 2 mol heme per mol dimer, which is believed to function as an electron donor for the high-potential cytochrome c oxidase.     

A Molecular Phylogeny of the Octopoda (Mollusca: Cephalopoda) Evaluated in Light of Morphological Evidence

In this paper we examine the phylogenetic relationships of the Octopoda utilizing molecular sequence data from the cytochrome c oxidase subunit I (COI) gene and compare results from analyses of molecular data with classifications and phylogenies based on previous morphological studies. Partial COI sequences (657 bp, excluding primers) were obtained from 28 species representing most of the diversity in the Order Octopoda, along with a sequence from the established sister taxon to the Octopoda, Vampyroteuthis infernalis. Our results exhibit a number of basic differences from inferences based on standard morphological data. We attempt to resolve these differences based on our confidence in various morphological features. An important finding is the failure of the molecular data to support the monophyly of the Octopodidae. This family contains over 90% of the species in the Suborder Incirrata and has always been difficult to define. Statistical tests constraining Octopodidae monophyly by use of parsimony and maximum-likelihood techniques suggest that all incirrates may be derived from octopodids.     

Phylogeny of Glaucosomatidae inferred from molecular evidence

Most species of glaucosomatids (Teleostei: Glaucosomatidae) are endemic to Australia, except Glaucosoma buergeri that is widely distributed from Australia to Japan. This study elucidated phylogenetic relationships among glaucosomatids based on the morphological characters of the saccular‐otolith sagitta, in addition to molecular evidence of mitochondrial 16S rDNA, cytochrome oxidase I (COI) and cytochrome b (cyt b) sequences, and nuclear rhodopsin sequences. The topologies of individuals' phylogenetic trees, based on 16S rDNA, COI and cyt b sequences, were statistically indistinguishable from one another, and were only slightly different from a tree based on rhodopsin sequences. These molecular tree topologies, however, differed from species relationships in morphology‐based phylogenetic hypothesis proposed in previous studies. Specimens of G. buergeri from Australia and Taiwan showed differences in the sagitta and molecular differentiation at the four genes, suggesting a possible speciation event. Both molecular and morphological evidences indicate that Glaucosoma magnificum is the plesiomorphic sister species of other glaucosomatid species. Glaucosoma hebraicum is the sister species of a clade composed of G. buergeri and Glaucosoma scapulare. Molecular and morphological evidences also support the species status of G. hebraicum.     

High-resolution melting (HRM) of the cytochrome B gene: a powerful approach to identify blood-meal sources in Chagas disease Vectors

Methods to determine blood-meal sources of hematophagous Triatominae bugs (Chagas disease vectors) are serological or based on PCR employing species-specific primers or heteroduplex analysis, but these are expensive, inaccurate, or problematic when the insect has fed on more than one species. To solve those problems, we developed a technique based on HRM analysis of the mitochondrial gene cytochrome B (Cyt b). This technique recognized 14 species involved in several ecoepidemiological cycles of the transmission of Trypanosoma cruzi and it was suitable with DNA extracted from intestinal content and feces 30 days after feeding, revealing a resolution power that can display mixed feedings. Field samples were analyzed showing blood meal sources corresponding to domestic, peridomiciliary and sylvatic cycles. The technique only requires a single pair of primers that amplify the Cyt b gene in vertebrates and no other standardization, making it quick, easy, relatively inexpensive, and highly accurate.