High-throughput monitoring of wild bee diversity and abundance via mitogenomics

1. Bee populations and other pollinators face multiple, synergistically acting threats, which have led to population declines, loss of local species richness and pollination services, and extinctions. However, our understanding of the degree, distribution and causes of declines is patchy, in part due to inadequate monitoring systems, with the challenge of taxonomic identification posing a major logistical barrier. Pollinator conservation would benefit from a high-throughput identification pipeline. 2. We show that the metagenomic mining and resequencing of mitochondrial genomes (mitogenomics) can be applied successfully to bulk samples of wild bees. We assembled the mitogenomes of 48 UK bee species and then shotgun-sequenced total DNA extracted from 204 whole bees that had been collected in 10 pan-trap samples from farms in England and been identified morphologically to 33 species. Each sample data set was mapped against the 48 reference mitogenomes. 3. The morphological and mitogenomic data sets were highly congruent. Out of 63 total species detections in the morphological data set, the mitogenomic data set made 59 correct detections (93�7% detection rate) and detected six more species (putative false positives). Direct inspection and an analysis with species-specific primers suggested that these putative false positives were most likely due to incorrect morphological IDs. Read frequency significantly predicted species biomass frequency (R2 = 24�9%). Species lists, biomass frequencies, extrapolated species richness and community structure were recovered with less error than in a metabarcoding pipeline. 4. Mitogenomics automates the onerous task of taxonomic identification, even for cryptic species, allowing the tracking of changes in species richness and istributions. A mitogenomic pipeline should thus be able to contain costs, maintain consistently high-quality data over long time series, incorporate retrospective taxonomic revisions and provide an auditable evidence trail. Mitogenomic data sets also provide estimates of species counts within samples and thus have potential for tracking population trajectories.

Phylogeny of Pelargonium (Geraniaceae) based on DNA sequences from three genomes

Phylogenetic hypotheses for the largely South African genus Pelargonium L'Hér. (Geraniaceae) were derived based on DNA sequence data from nuclear, chloroplast and mitochondrial encoded regions. The datasets were unequally represented and comprised cpDNA trnL-F sequences for 152 taxa, nrDNA ITS sequences for 55 taxa, and mtDNA nad1 b/c exons for 51 taxa. Phylogenetic hypotheses derived from the separate three datasets were overall congruent. A single hypothesis synthesising the information in the three datasets was constructed following a total evidence approach and implementing dataset specific stepmatrices in order to correct for substitution biases. Pelargonium was found to consist of five main clades, some with contrasting evolutionary patterns with respect to biogeographic distributions, dispersal capacity, pollination biology and karyological diversification. The five main clades are structured in two (subgeneric) clades that correlate with chromosome size. One of these clades includes a "winter rainfall clade" containing more than 70% of all currently described Pelargonium species, and all restricted to the South African Cape winter rainfall region. Apart from (woody) shrubs and small herbaceous rosette subshrubs, this clade comprises a large "xerophytic" clade including geophytes, stem and leaf succulents, harbouring in total almost half of the genus. This clade is considered to be the result of in situ proliferation, possibly in response to late-Miocene and Pliocene aridification events. Nested within it is a radiation comprising c. 80 species from the geophytic Pelargonium section Hoarea, all characterised by the possession of (a series of) tunicate tubers.

Wolbachia infection and dramatic intraspecific mitochondrial dna divergence in a fig wasp

Mitochondria and Wolbachia are maternally inherited genomes that exhibit strong linkage disequilibrium in many organisms. We surveyed Wolbachia infections in 187 specimens of the fig wasp species, Ceratosolen solmsi, and found an infection prevalence of 89.3%. DNA sequencing of 20 individuals each from Wolbachia-infected and -uninfected subpopulations revealed extreme mtDNA divergence (up to 9.2% and 15.3% in CO1 and cytochrome b, respectively) between infected and uninfected wasps. Further, mtDNA diversity was significantly reduced within the infected group. Our sequencing of a large part of the mitochondrial genome from both Wolbachia-infected and -uninfected individuals revealed that high sequence divergence is common throughout the mitochondrial genome. These patterns suggest a partial selective sweep of mitochondria subsequent to the introduction of Wolbachia into C. solsmi, by hybrid introgression from a related species.

Metal bioaccumulation and cellular fractionation in an epigeic earthworm (Lumbricus rubellus): the interactive influences of population exposure histories, site-specific geochemistry and mitochondrial genotype

Subcellular fractionation techniques were used to describe temporal changes (at intervals from T0 to T70 days) in the Pb, Zn and P partitioning profiles of Lumbricus rubellus populations from one calcareous (MDH) and one acidic (MCS) geographically isolated Pb/Zn-mine sites and one reference site (CPF). MDH and MCS individuals were laboratory maintained on their native field soils; CPF worms were exposed to both MDH and MCS soils. Site-specific differences in metal partitioning were found: notably, the putatively metal-adapted populations, MDH and MCS, preferentially partitioned higher proportions of their accumulated tissue metal burdens into insoluble CaPO4-rich organelles compared with naive counterparts, CPF. Thus, it is plausible that efficient metal immobilization is a phenotypic trait characterising metal tolerant ecotypes. Mitochondrial cytochrome oxidase II (COII) genotyping revealed that the populations indigenous to mine and reference soils belong to distinct genetic lineages, differentiated by 13%, with 7 haplotypes within the reference site lineage but fewer (3 and 4, respectively) in the lineage common to the two mine sites. Collectively, these observations raise the possibility that site-related genotype differences could influence the toxico-availability of metals and, thus, represent a potential confounding variable in field-based eco-toxicological assessments.

Land plants and DNA barcodes: short-term and long-term goals

Land plants have had the reputation of being problematic for DNA barcoding for two general reasons: (i) the standard DNA regions used in algae, animals and fungi have exceedingly low levels of variability and (ii) the typically used land plant plastid phylogenetic markers (e.g. rbcL, trnL-F, etc.) appear to have too little variation. However, no one has assessed how well current phylogenetic resources might work in the context of identification (versus phylogeny reconstruction). In this paper, we make such an assessment, particularly with two of the markers commonly sequenced in land plant phylogenetic studies, plastid rbcL and internal transcribed spacers of the large subunits of nuclear ribosomal DNA (ITS), and find that both of these DNA regions perform well even though the data currently available in GenBank/EBI were not produced to be used as barcodes and BLAST searches are not an ideal tool for this purpose. These results bode well for the use of even more variable regions of plastid DNA (such as, for example, psbA-trnH) as barcodes, once they have been widely sequenced. In the short term, efforts to bring land plant barcoding up to the standards being used now in other organisms should make swift progress. There are two categories of DNA barcode users, scientists in fields other than taxonomy and taxonomists. For the former, the use of mitochondrial and plastid DNA, the two most easily assessed genomes, is at least in the short term a useful tool that permits them to get on with their studies, which depend on knowing roughly which species or species groups they are dealing with, but these same DNA regions have important drawbacks for use in taxonomic studies (i.e. studies designed to elucidate species limits). For these purposes, DNA markers from uniparentally (usually maternally) inherited genomes can only provide half of the story required to improve taxonomic standards being used in DNA barcoding. In the long term, we will need to develop more sophisticated barcoding tools, which would be multiple, low-copy nuclear markers with sufficient genetic variability and PCR-reliability; these would permit the detection of hybrids and permit researchers to identify the 'genetic gaps' that are useful in assessing species limits.

Determination of mitochondrial genetic diversity in mammals

Mitochondrial DNA (mtDNA) is one of the most Popular population genetic markers. Its relevance as an indicator Of Population size and history has recently been questioned by several large-scale studies in animals reporting evidence for recurrent adaptive evolution, at least in invertebrates. Here we focus on mammals, a more restricted taxonomic group for which the issue of mtDNA near neutrality is crucial. By analyzing the distribution of mtDNA diversity across species and relating 4 to allozyme diversity, life-history traits, and taxonomy, we show that (i) mtDNA in mammals (toes not reject the nearly neutral model; (ii) mtDNA diversity, however, is unrelated to any of the 14 life-history and ecological variables that we analyzed, including body mass, geographic range, and The World Conservation Union (IUCN) categorization; (iii) mtDNA diversity is highly variable between mammalian orders and families; (iv) this taxonomic effect is most likely explained by variations of mutation rate between lineages. These results are indicative of a strong stochasticity of effective population size in mammalian species. They Suggest that, even in the absence of selection, mtDNA genetic diversity is essentially unpredictable, knowing species biology, and probably uncorrelated to species abundance.

Simple sequence repeats reveal uneven distribution of genetic diversity in chloroplast genomes of Brassica oleracea L. and (n=9) wild relatives

Diversity in the chloroplast genome of 171 accessions representing the Brassica 'C' (n = 9) genome, including domesticated and wild B. oleracea and nine inter-fertile related wild species, was investigated using six chloroplast SSR (microsatellite) markers. The lack of diversity detected among 105 cultivated and wild accessions of B. oleracea contrasted starkly with that found within its wild relatives. The vast majority of B. oleracea accessions shared a single haplotype, whereas as many as six haplotypes were detected in two wild species, B. villosa Biv. and B. cretica Lam.. The SSRs proved to be highly polymorphic across haplotypes, with calculated genetic diversity values (H) of 0.23-0.87. In total, 23 different haplotypes were detected in C genome species, with an additional five haplotypes detected in B. rapa L. (A genome n = 10) and another in B. nigra L. (B genome, n = 8). The low chloroplast diversity of B. oleracea is not suggestive of multiple domestication events. The predominant B. oleracea haplotype was also common in B. incana Ten. and present in low frequencies in B. villosa, B. macrocarpa Guss, B. rupestris Raf. and B. cretica. The chloroplast SSRs reveal a wealth of diversity within wild Brassica species that will facilitate further evolutionary and phylogeographic studies of this important crop genus.

Stocking may increase mitochondrial DNA diversity but fails to halt the decline of endangered Atlantic salmon populations

Over the last 50 years, Spanish Atlantic salmon (Salmo salar) populations have been in decline. In order to bolster these populations, rivers were stocked with fish of northern European origin during the period 1974-1996, probably also introducing the furunculosis-inducing pathogen, Aeromonas salmonicida. Here we assess the relative importance of processes influencing mitochondrial (mt)DNA variability in these populations from 1948 to 2002. Genetic material collected over this period from four rivers in northern Spain (Cantabria) was used to detect variability at the mtDNA ND1 gene. Before stocking, a single haplotype was found at high frequency (0.980). Following stocking, haplotype diversity (h) increased in all rivers (mean h before stocking was 0.041, and 0.245 afterwards). These increases were due principally to the dramatic increase in frequency of a previously very low frequency haplotype, reported at higher frequencies in northern European populations proximate to those used to stock Cantabrian rivers. Genetic structuring increased after stocking: among-river differentiation was low before stocking (1950s/1960s Phi(ST) = -0.00296-0.00284), increasing considerably at the height of stocking (1980s Phi(ST) = 0.18932) and decreasing post-stocking (1990s/2002 Phi(ST) = 0.04934-0.03852). Gene flow from stocked fish therefore seems to have had a substantial role in increasing mtDNA variability. Additionally, we found significant differentiation between individuals that had probably died from infectious disease and apparently healthy, angled fish, suggesting a possible role for pathogen-driven selection of mtDNA variation. Our results suggest that stocking with non-native fish may increase genetic diversity in the short term, but may not reverse population declines.

Rare and fleeting: an example of interspecific recombination in animal mitochondrial DNA

Recombination is thought to occur only rarely in animal mitochondrial DNA ( mtDNA). However, detection of mtDNA recombination requires that cells become heteroplasmic through mutation, intramolecular recombination or ' leakage' of paternal mtDNA. Interspecific hybridization increases the probability of detecting mtDNA recombinants due to higher levels of sequence divergence and potentially higher levels of paternal leakage. During a study of historical variation in Atlantic salmon ( Salmo salar) mtDNA, an individual with a recombinant haplotype containing sequence from both Atlantic salmon and brown trout ( Salmo trutta) was detected. The individual was not an F1 hybrid but it did have an unusual nuclear genotype which suggested that it was a later-generation backcross. No other similar recombinant haplotype was found from the same population or three neighbouring Atlantic salmon populations in 717 individuals collected during 1948 - 2002. Interspecific recombination may increase mtDNA variability within species and can have implications for phylogenetic studies.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes

Homeobox genes encode DNA-binding proteins, many of which are implicated in the control of embryonic development. Evolutionarily, most homeobox genes fall into two related clades: the ANTP and the PRD classes. Some genes in ANTP class, notably Hox, ParaHox, and NK genes, have an intriguing arrangement into physical clusters. To investigate the evolutionary history of these gene clusters, we examined homeobox gene chromosomal locations in the cephalochordate amphioxus, Branchiostoma floridae. We deduce that 22 amphioxus ANTP class homeobox genes localize in just three chromosomes. One contains the Hox cluster plus AmphiEn, AmphiMnx, and AmphiDll. The ParaHox cluster resides in another chromosome, whereas a third chromosome contains the NK type homeobox genes, including AmphiMsx and ArnphiTlx. By comparative analysis we infer that clustering of ANTP class homeobox genes evolved just once, during a series of extensive cis-duplication events of genes early in animal evolution. A trans-duplication event occurred later to yield the Hox and ParaHox gene clusters on different chromosomes. The results obtained have implications for understanding the origin of homeobox gene clustering, the diversification of the ANTP class of homeobox genes, and the evolution of animal genomes.

The poliovirus 2C cis-acting replication element-mediated uridylylation of VPg is not required for synthesis of negative-sense genomes

Nucleotides in the terminal loop of the poliovirus 2C cis-acting replication element (2C(CRE)), a 61 nt structured RNA, function as the template for the addition of two uridylate (U) residues to the viral protein VPg. This uridylylation reaction leads to the formation of VPgpUpU, which is used by the viral RNA polymerase as a nucleotide-peptide primer for genome replication. Although VPg primes both positive- and negative-strand replication, the specific requirement for 2C(CRE)-mediated uridylylation for one or both events has not been demonstrated. We have used a cell-free in vitro translation and replication reaction to demonstrate that 2C(CRE) is not required for the initiation of the negative-sense strand, which is synthesized in the absence of 2C(CRE)-mediated VPgpUpU formation. We propose that the 3' poly(A) tail could serve as the template for the formation of a VPg-poly(U) primer that functions in the initiation of negative-sense strands.

Do mitochondriotropic antioxidants prevent chlorinative stress-induced mitochondrial and cellular injury?

Reactive chlorine species such as hypochlorous acid ( HOCl) are cytotoxic oxidants generated by activated neutrophils at the sites of chronic inflammation. Since mitochondria are key mediators of apoptosis and necrosis, we hypothesized that mitochondriotropic antioxidants could limit HOCl-mediated intracellular oxidative injury to human fetal liver cells, preserve mitochondrial function, and prevent cell death. In this current study, we show that recently developed mitochondria-targeted antioxidants ( MitoQ and SS31) significantly protected against HOCl-induced mitochondrial damage and cell death at concentrations >= 25 nM. Our study highlights the potential application of mitochondria-specific targeted antioxidants for the prevention of cellular dysfunction and cell death under conditions of chlorinative stress, as occurs during chronic inflammation.

The pro-inflammatory oxidant hypochlorous acid induces Bax-dependent mitochondrial permeabilisation and cell death through AIF-/EndoG-dependent pathways

At sites of chronic inflammation, such as in the inflamed rheumatoid joint, activated neutrophils release hydrogen peroxide (H2O2) and the enzyme myeloperoxidase to catalyse the formation of hypochlorous acid (HOCl). 3-chlorotyrosine, a marker of HOCl in vivo, has been observed in synovial fluid proteins from rheumatoid arthritis patients. However the mechanisms of HOCl-induced cytotxicity are unknown. We determined the molecular mechanisms by which HOCl induced cell death in human mesenchymal progenitor cells (MPCs) differentiated into a chondrocytic phenotype as a model of human cartilage cells and show that HOCl induced rapid Bax conformational change, mitochondrial permeability and release of intra-mitochondrial pro-apoptotic proteins which resulted in nuclear translocation of AIF and EndoG. siRNA-mediated knockdown of Bax substantially prevented mitochondrial permeability, release of intra-mitochondrial pro-apoptotic proteins. Cell death was inhibited by siRNA-mediated knockdown of Bax, AIF or EndoG. Although we observed several biochemical markers of apoptosis, caspase activation was not detected either by western blotting, fluorescence activity assays or by using caspase inhibitors to inhibit cell death. This was further supported by findings that (1) in vitro exposure of recombinant human caspases to HOCl caused significant inhibition of caspase activity and (2) the addition of HOCl to staurosporine-treated MPCs inhibited the activity of cellular caspases. Our results show for the first time that HOCl induced Bax-dependent mitochondrial permeability which led to cell death without caspase activity by processes involving AIF/EndoG-dependent pathways. Our study provides a novel insight into the potential mechanisms of cell death in the inflamed human joint. (c) 2006 Elsevier Inc. All rights reserved.