Metabarcoding of marine zooplankton: prospects, progress and pitfalls

Metabarcoding (large-scale taxonomic identification of complex samples via analysis of one or few orthologous DNA regions, called barcodes) is revolutionizing analysis of biodiversity of marine zooplankton assemblages. Metabarcoding relies on high-throughput DNA sequencing (HTS) technologies, which yield millions of DNA sequences in parallel and allow large-scale analysis of environmental samples. Metabarcoding studies of marine zooplankton have used various regions of nuclear small- (18S) and large-subunit (28S) rRNA, which allow accurate classification of novel sequences and reliable amplification with consensus primers, but- due to their relatively conserved nature- may underestimate species diversity in a community. To discriminate species, more variable genes are needed. A limited number of metabarcoding studies have used mitochondrial cytochrome oxidase I (COI), which ensures detection of species-level diversity, but may require group-specific primers and thus result in inconsistent amplification success rates. Reference databases with sequences for accurately-identified species are critically needed to allow taxonomic designation of molecular operational taxonomic units (MOTU) and comparison with previous studies of zooplankton diversity. Potential and promising applications of metabarcoding include rapid detection of impacts of climate change, monitoring and assessment of ecosystem health, calculation of biotic indices, characterization of food webs and detection of introduced, non-indigenous species

miRanalyzer: an update on the detection and analysis of microRNAs in high-throughput sequencing experiments

We present a new version of miRanalyzer, a web server and stand-alone tool for the detection of known and prediction of new microRNAs in high-throughput sequencing experiments. The new version has been notably improved regarding speed, scope and available features. Alignments are now based on the ultrafast short-read aligner Bowtie (granting also colour space support, allowing mismatches and improving speed) and 31 genomes, including 6 plant genomes, can now be analysed (previous version contained only 7). Differences between plant and animal microRNAs have been taken into account for the prediction models and differential expression of both, known and predicted microRNAs, between two conditions can be calculated. Additionally, consensus sequences of predicted mature and precursor microRNAs can be obtained from multiple samples, which increases the reliability of the predicted microRNAs. Finally, a stand-alone version of the miRanalyzer that is based on a local and easily customized database is also available; this allows the user to have more control on certain parameters as well as to use specific data such as unpublished assemblies or other libraries that are not available in the web server. miRanalyzer is available at http://bioinfo2.ugr.es/miRanalyzer/miRanalyzer.php.The Ministry of Innovation and Science of the Spanish Government (BIO2010-20219 to M.H.); the Junta de Andalucia (P07FQM3163 to M.H.); the ‘Juan de la Cierva’ fellowship (to M.H.); the Department of Industry, Tourism and Trade of the Government of the Autonomous Community of the Basque Country (Etortek Research Programs 2009/20011 to A.M.A.); from the Innovation Technology Department of the Bizkaia County (to A.M.A.). Funding for open access charge: Department of Industry, Tourism and Trade of the Government of the Autonomous Community of the Basque Country (Etortek Research Programs 2009/2011 to A.M.A.)

The SAR11 Group of Alpha-Proteobacteria Is Not Related to the Origin of Mitochondria

Although free living, members of the successful SAR11 group of marine alpha-proteobacteria contain a very small and A+T rich genome, two features that are typical of mitochondria and related obligate intracellular parasites such as the Rickettsiales. Previous phylogenetic analyses have suggested that Candidatus Pelagibacter ubique, the first cultured member of this group, is related to the Rickettsiales+mitochondria clade whereas others disagree with this conclusion. In order to determine the evolutionary position of the SAR11 group and its relationship to the origin of mitochondria, we have performed phylogenetic analyses on the concatenation of 24 proteins from 5 mitochondria and 71 proteobacteria. Our results support that SAR11 group is not the sistergroup of the Rickettsiales+mitochondria clade and confirm that the position of this group in the alpha-proteobacterial tree is strongly affected by tree reconstruction artefacts due to compositional bias. As a consequence, genome reduction and bias toward a high A+T content may have evolved independently in the SAR11 species, which points to a different direction in the quest for the closest relatives to mitochondria and Rickettsiales. In addition, our analyses raise doubts about the monophyly of the newly proposed Pelagibacteraceae family

Site-specific time heterogeneity of the substitution process and its impact on phylogenetic inference

<p>Abstract</p> <p>Background</p> <p>Model violations constitute the major limitation in inferring accurate phylogenies. Characterizing properties of the data that are not being correctly handled by current models is therefore of prime importance. One of the properties of protein evolution is the variation of the relative rate of substitutions across sites and over time, the latter is the phenomenon called heterotachy. Its effect on phylogenetic inference has recently obtained considerable attention, which led to the development of new models of sequence evolution. However, thus far focus has been on the quantitative heterogeneity of the evolutionary process, thereby overlooking more qualitative variations.</p> <p>Results</p> <p>We studied the importance of variation of the site-specific amino-acid substitution process over time and its possible impact on phylogenetic inference. We used the CAT model to define an infinite mixture of substitution processes characterized by equilibrium frequencies over the twenty amino acids, a useful proxy for qualitatively estimating the evolutionary process. Using two large datasets, we show that qualitative changes in site-specific substitution properties over time occurred significantly. To test whether this unaccounted qualitative variation can lead to an erroneous phylogenetic tree, we analyzed a concatenation of mitochondrial proteins in which Cnidaria and Porifera were erroneously grouped. The progressive removal of the sites with the most heterogeneous CAT profiles across clades led to the recovery of the monophyly of Eumetazoa (Cnidaria+Bilateria), suggesting that this heterogeneity can negatively influence phylogenetic inference.</p> <p>Conclusion</p> <p>The time-heterogeneity of the amino-acid replacement process is therefore an important evolutionary aspect that should be incorporated in future models of sequence change.</p

Instability of Plastid DNA in the Nuclear Genome

Functional gene transfer from the plastid (chloroplast) and mitochondrial genomes to the nucleus has been an important driving force in eukaryotic evolution. Non-functional DNA transfer is far more frequent, and the frequency of such transfers from the plastid to the nucleus has been determined experimentally in tobacco using transplastomic lines containing, in their plastid genome, a kanamycin resistance gene (neo) readymade for nuclear expression. Contrary to expectations, non-Mendelian segregation of the kanamycin resistance phenotype is seen in progeny of some lines in which neo has been transferred to the nuclear genome. Here, we provide a detailed analysis of the instability of kanamycin resistance in nine of these lines, and we show that it is due to deletion of neo. Four lines showed instability with variation between progeny derived from different areas of the same plant, suggesting a loss of neo during somatic cell division. One line showed a consistent reduction in the proportion of kanamycin-resistant progeny, suggesting a loss of neo during meiosis, and the remaining four lines were relatively stable. To avoid genomic enlargement, the high frequency of plastid DNA integration into the nuclear genome necessitates a counterbalancing removal process. This is the first demonstration of such loss involving a high proportion of recent nuclear integrants. We propose that insertion, deletion, and rearrangement of plastid sequences in the nuclear genome are important evolutionary processes in the generation of novel nuclear genes. This work is also relevant in the context of transgenic plant research and crop production, because similar processes to those described here may be involved in the loss of plant transgenes

Phylogenomic Analysis of Odyssella thessalonicensis Fortifies the Common Origin of Rickettsiales, Pelagibacter ubique and Reclimonas americana Mitochondrion

Background: The evolution of the Alphaproteobacteria and origin of the mitochondria are topics of considerable debate. Most studies have placed the mitochondria ancestor within the Rickettsiales order. Ten years ago, the bacterium Odyssella thessalonicensis was isolated from Acanthamoeba spp., and the 16S rDNA phylogeny placed it within the Rickettsiales. Recently, the whole genome of O. thessalonicensis has been sequenced, and 16S rDNA phylogeny and more robust and accurate phylogenomic analyses have been performed with 65 highly conserved proteins. Methodology/Principal Findings: The results suggested that the O. thessalonicensis emerged between the Rickettsiales and other Alphaproteobacteria. The mitochondrial proteins of the Reclinomonas americana have been used to locate the phylogenetic position of the mitochondrion ancestor within the Alphaproteobacteria tree. Using the K tree score method, nine mitochondrion-encoded proteins, whose phylogenies were congruent with the Alphaproteobacteria phylogenomic tree, have been selected and concatenated for Bayesian and Maximum Likelihood phylogenies. The Reclinomonas americana mitochondrion is a sister taxon to the free-living bacteria Candidatus Pelagibacter ubique, and together, they form a clade that is deeply rooted in the Rickettsiales clade. Conclusions/Significance: The Reclinomonas americana mitochondrion phylogenomic study confirmed that mitochondri

Multigene Phylogeny of Choanozoa and the Origin of Animals

Animals are evolutionarily related to fungi and to the predominantly unicellular protozoan phylum Choanozoa, together known as opisthokonts. To establish the sequence of events when animals evolved from unicellular ancestors, and understand those key evolutionary transitions, we need to establish which choanozoans are most closely related to animals and also the evolutionary position of each choanozoan group within the opisthokont phylogenetic tree. Here we focus on Ministeria vibrans, a minute bacteria-eating cell with slender radiating tentacles. Single-gene trees suggested that it is either the closest unicellular relative of animals or else sister to choanoflagellates, traditionally considered likely animal ancestors. Sequencing thousands of Ministeria protein genes now reveals about 14 with domains of key significance for animal cell biology, including several previously unknown from deeply diverging Choanozoa, e.g. domains involved in hedgehog, Notch and tyrosine kinase signaling or cell adhesion (cadherin). Phylogenetic trees using 78 proteins show that Ministeria is not sister to animals or choanoflagellates (themselves sisters to animals), but to Capsaspora, another protozoan with thread-like (filose) tentacles. The Ministeria/Capsaspora clade (new class Filasterea) is sister to animals and choanoflagellates, these three groups forming a novel clade (filozoa) whose ancestor presumably evolved filose tentacles well before they aggregated as a periciliary collar in the choanoflagellate/sponge common ancestor. Our trees show ichthyosporean choanozoans as sisters to filozoa; a fusion between ubiquitin and ribosomal small subunit S30 protein genes unifies all holozoa (filozoa plus Ichthyosporea), being absent in earlier branching eukaryotes. Thus, several successive evolutionary innovations occurred among their unicellular closest relatives prior to the origin of the multicellular body-plan of animals

North African Influences and Potential Bias in Case-Control Association Studies in the Spanish Population

BACKGROUND: Despite the limited genetic heterogeneity of Spanish populations, substantial evidences support that historical African influences have not affected them uniformly. Accounting for such population differences might be essential to reduce spurious results in association studies of genetic factors with disease. Using ancestry informative markers (AIMs), we aimed to measure the African influences in Spanish populations and to explore whether these might introduce statistical bias in population-based association studies. METHODOLOGY/PRINCIPAL FINDINGS: We genotyped 93 AIMs in Spanish (from the Canary Islands and the Iberian Peninsula) and Northwest Africans, and conducted population and individual-based clustering analyses along with reference data from the HapMap, HGDP-CEPH, and other sources. We found significant differences for the Northwest African influence among Spanish populations from as low as ≈ 5% in Spanish from the Iberian Peninsula to as much as ≈ 17% in Canary Islanders, whereas the sub-Saharan African influence was negligible. Strikingly, the Northwest African ancestry showed a wide inter-individual variation in Canary Islanders ranging from 0% to 96%, reflecting the violent way the Islands were conquered and colonized by the Spanish in the XV century. As a consequence, a comparison of allele frequencies between Spanish samples from the Iberian Peninsula and the Canary Islands evidenced an excess of markers with significant differences. However, the inflation of p-values for the differences was adequately controlled by correcting for genetic ancestry estimates derived from a reduced number of AIMs. CONCLUSIONS/SIGNIFICANCE: Although the African influences estimated might be biased due to marker ascertainment, these results confirm that Northwest African genetic footprints are recognizable nowadays in the Spanish populations, particularly in Canary Islanders, and that the uneven African influences existing in these populations might increase the risk for false positives in association studies. Adjusting for population stratification assessed with a few dozen AIMs would be sufficient to control this effect