Diversity of zooplankton assemblages at Station L4, Western English Channel, measured using next generation DNA sequencing

Background: Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel Next Generation Sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify the richness and diversity of a mixed zooplankton assemblage from a productive monitoring site in the Western English Channel. Methodology/Principle Findings: Plankton WP2 replicate net hauls (200 µm) were taken at the Western Channel Observatory long-term monitoring station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,042 sequences were obtained for all samples. The sequences clustered in to 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 138 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 75 taxonomic groups. Conclusions: The percentage of OTUs assigned to major eukaryotic taxonomic groups broadly aligns between the metagenetic and morphological analysis and are dominated by Copepoda. However, the metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for estimating diversity and species richness of zooplankton communities.

Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

Evolutionary origin of Plasmodium and other Apicomplexa based on rRNA genes.

We have explored the evolutionary history of the Apicomplexa and two related protistan phyla, Dinozoa and Ciliophora, by comparing the nucleotide sequences of small subunit ribosomal RNA genes. We conclude that the Plasmodium lineage, to which the malarial parasites belong, diverged from other apicomplexan lineages (piroplasmids and coccidians) several hundred million years ago, perhaps even before the Cambrian. The Plasmodium radiation, which gave rise to several species parasitic to humans, occurred approximately 129 million years ago; Plasmodium parasitism of humans has independently arisen several times. The origin of apicomplexans (Plasmodium), dinoflagellates, and ciliates may be > 1 billion years old, perhaps older than the three multicellular kingdoms of animals, plants, and fungi. Digenetic parasitism independently evolved several times in the Apicomplexa.

Quantification of microbial communities in subsurface marine sediments of the Black Sea and off Namibia

Organic-rich subsurface marine sediments were taken by gravity coring up to a depth of 10 m below seafloor at six stations from the anoxic Black Sea and the Benguela upwelling system off Namibia during the research cruises Meteor 72-5 and 76-1, respectively. The quantitative microbial community composition at various sediment depths was analyzed using total cell counting, catalyzed reporter deposition fluorescence in situ hybridization (CARD FISH) and quantitative real-time PCR (Q-PCR). Total cell counts decreased with depths from 10(9) to 10(10) cells/mL at the sediment surface to 10(7)-10(9) cells/mL below one meter depth. Based on CARD FISH and Q-PCR analyses overall similar proportions of Bacteria and Archaea were found. The down-core distribution of prokaryotic and eukaryotic small subunit ribosomal RNA genes (16S and 18S rRNA) as well as functional genes involved in different biogeochemical processes was quantified using Q-PCR. Crenarchaeota and the bacterial candidate division JS-1 as well as the classes Anaerolineae and Caldilineae of the phylum Chloroflexi were highly abundant. Less abundant but detectable in most of the samples were Eukarya as well as the metal and sulfate-reducing Geobacteraceae (only in the Benguela upwelling influenced sediments). The functional genes cbbL, encoding for the large subunit of RuBisCO, the genes dsrA and aprA, indicative of sulfate-reducers as well as the mcrA gene of methanogens were detected in the Benguela upwelling and Black Sea sediments. Overall, the high organic carbon content of the sediments goes along with high cell counts and high gene copy numbers, as well as an equal abundance of Bacteria and Archaea.

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5'ETS molecule using three distinct methods and located the acceptor site between two known 5'ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5'ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5'ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5'ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5' ETS molecule using three distinct methods and located the acceptor site between two known 5' ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5' ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5' ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5' ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

Functions and intramolecular interactions of ribosomal RNA in decoding of termination codons

Two regions in the 3$\prime$ domain of 16S rRNA (the RNA of the small ribosomal subunit) have been implicated in decoding of termination codons. Using segment-directed PCR random mutagenesis, I isolated 33 translational suppressor mutations in the 3$\prime$ domain of 16S rRNA. Characterization of the mutations by both genetic and biochemical methods indicated that some of the mutations are defective in UGA-specific peptide chain termination and that others may be defective in peptide chain termination at all termination codons. The studies of the mutations at an internal loop in the non-conserved region of helix 44 also indicated that this structure, in a non-conserved region of 16S rRNA, is involved in both peptide chain termination and assembly of 16S rRNA.^ With a suppressible trpA UAG nonsense mutation, a spontaneously arising translational suppressor mutation was isolated in the rrnB operon cloned into a pBR322-derived plasmid. The mutation caused suppression of UAG at two codon positions in trpA but did not suppress UAA or UGA mutations at the same trpA positions. The specificity of the rRNA suppressor mutation suggests that it may cause a defect in UAG-specific peptide chain termination. The mutation is a single nucleotide deletion (G2484$\Delta$) in helix 89 of 23S rRNA (the large RNA of the large ribosomal subunit). The result indicates a functional interaction between two regions of 23S rRNA. Furthermore, it provides suggestive in vivo evidence for the involvement of the peptidyl-transferase center of 23S rRNA in peptide chain termination. The $\Delta$2484 and A1093/$\Delta$2484 (double) mutations were also observed to alter the decoding specificity of the suppressor tRNA lysT(U70), which has a mutation in its acceptor stem. That result suggests that there is an interaction between the stem-loop region of helix 89 of 23S rRNA and the acceptor stem of tRNA during decoding and that the interaction is important for the decoding specificity of tRNA.^ Using gene manipulation procedures, I have constructed a new expression vector to express and purify the cellular protein factors required for a recently developed, realistic in vitro termination assay. The gene for each protein was cloned into the newly constructed vector in such a way that expression yielded a protein with an N-terminal affinity tag, for specific, rapid purification. The amino terminus was engineered so that, after purification, the unwanted N-terminal tag can be completely removed from the protein by thrombin cleavage, yielding a natural amino acid sequence for each protein. I have cloned the genes for EF-G and all three release factors into this new expression vector and the genes for all the other protein factors into a pCAL-n expression vector. These constructs will allow our laboratory group to quickly and inexpensively purify all the protein factors needed for the new in vitro termination assay. (Abstract shortened by UMI.) ^

Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15–25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

Three small nucleolar RNAs that are involved in ribosomal RNA precursor processing

Three small nucleolar RNAs (snoRNAs), E1, E2 and E3, have been described that have unique sequences and interact directly with unique segments of pre-rRNA in vivo. In this report, injection of antisense oligodeoxynucleotides into Xenopus laevis oocytes was used to target the specific degradation of these snoRNAs. Specific disruptions of pre-rRNA processing were then observed, which were reversed by injection of the corresponding in vitro-synthesized snoRNA. Degradation of each of these three snoRNAs produced a unique rRNA maturation phenotype. E1 RNA depletion shut down 18 rRNA formation, without overaccumulation of 20S pre-rRNA. After E2 RNA degradation, production of 18S rRNA and 36S pre-rRNA stopped, and 38S pre-rRNA accumulated, without overaccumulation of 20S pre-rRNA. E3 RNA depletion induced the accumulation of 36S pre-rRNA. This suggests that each of these snoRNAs plays a different role in pre-rRNA processing and indicates that E1 and E2 RNAs are essential for 18S rRNA formation. The available data support the proposal that these snoRNAs are at least involved in pre-rRNA processing at the following pre-rRNA cleavage sites: E1 at the 5′ end and E2 at the 3′ end of 18S rRNA, and E3 at or near the 5′ end of 5.8S rRNA.

Phylogenetic relationships of Trypanosoma chelodina and Trypanosoma binneyi from Australian tortoises and platypuses inferred from small subunit rRNA analyses

Trypanosome infections are often difficult to detect by conventional microscopy and their pleomorphy often confounds differential diagnosis. Molecular techniques are now being used to diagnose infections and to determine phylogenetic relationships between species. Complete small subunit rRNA gene sequences were determined for isolates of Trypanosoma chelodina from the Brisbane River tortoise (Emydura signata), the saw-shelled tortoise (Elseya latisternum), and the eastern snake-necked tortoise (Chelodina longicollis) from southeast Queensland, Australia. Partial sequence data were also obtained for T. binneyi from a platypus (Ornithorhynchus anatinus) from Tasmania. Phylogenetic relationships between T. chelodina, T. binneyi and other species were examined by maximum parsimony and likelihood methods. The Australian tortoise and platypus trypanosomes did not exhibit any close phylogenetic relationships with those of mammals, reptiles or amphibians, but were closely related to each other, and to fish trypanosomes. This contra-indicates their co-evolution with their vertebrate hosts but does not exclude co-evolution with different groups of invertebrate vectors, notably insects and leeches.

Identification of Biomphalaria havanensis and Biomphalaria obstructa populations from Cuba using polymerase chain reaction and restriction fragment length polymorphism of the ribosomal RNA intergenic spacer

In Cuba, several Biomphalaria species have been reported such as B. orbignyi, B. schrammi, B. helophila, B. havanensis and B. peregrina; only the latter three are considered as potential hosts of Schistosoma mansoni. The specific identification of Biomphalaria species is based on anatomical and morphological characters of genital organs and shells. The correct identification of these snails is complicated by the high variation in these characters, similarity among species and in some cases by the small size of the snails. In this paper, we reported the classical morphological identification, the use of PCR and RFLP analysis of the internal transcribed spacer region of the ribosomal RNA genes for molecular identification of seven snail populations from different localities in Cuba. Using morphological and molecular analysis, we showed that among the studied Cuban Biomphalaria populations only B. havanensis and B. obstructa species were found.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.

Study of the nucleolar cycle and ribosomal RNA distribution during meiosis in triatomines (Triatominae, Heteroptera)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The characterization of mouse carboxypeptidase N small subunit gene structure and presence in developing embryos

Carboxypeptidase N (CPN) is a plasma zinc metalloprotease, which consists of two enzymatically active small subunits and two large subunits that protect the protein from degradation. CPN cleaves carboxy-terminal arginines and lysines from peptides found in the bloodstream such as complement anaphylatoxins, kinins, and creatine kinase MM. In this study, the mouse CPN small subunit (CPN1) coding region, gene structure, and chromosomal location were characterized and the expression of CPN1 was investigated in mouse embryos at different stages of development. The CPN1 gene, which was approximately 29 kb in length, contained nine exons and localized to mouse chromosome 19D2. The fifth and sixth exons of CPN1 encoded the amino acids necessary for substrate binding and catalytic activity. CPN1 RNA was expressed predominately in adult liver and contained a 1371 bp open reading frame encoding 457 amino acids. In the mouse embryo, CPN1 RNA was observed at 8.5 days post coitus (dpc), while its protein was detected at 10.5 dpc. In situ hybridization of the fetal liver detected CPN1 RNA in erythroid progenitor cells at 10.5, 13.5, and 16.5 dpc and in hepatocytes at 16.5 dpc. This was compared to the expression of the complement component C3, the parent molecule of complement anaphylatoxin C3a. Consistently throughout the experiments, CPN1 message and protein preceded the expression of C3. To obtain a better understanding of the biological significance of CPN1 in vivo, studies were initiated to produce a genetically engineered mouse in which the CPN1 gene was ablated. To facilitate this project a targeting vector was constructed by removing the functionally important fifth and sixth exons of the CPN1 gene. Collectively, these studies have: (1) provided important detailed information regarding the structure and organization of the murine CPN1 gene, (2) yielded insights into the developmental expression of mouse CPN1 in relationship to C3 expression, and (3) set the stage for the generation of a CPN1 “knock-out” mouse, which can be used to determine the biological significance of CPN1 in both normal and diseased conditions. ^

A hyperphosphorylated form of the large subunit of RNA polymerase II is associated with splicing complexes and the nuclear matrix.

A hyperphosphorylated form of the largest subunit of RNA polymerase II (pol IIo) is associated with the pre-mRNA splicing process. Pol IIo was detected in association with a subset of small nuclear ribonucleoprotein particle and Ser-Arg protein splicing factors and also with pre-mRNA splicing complexes assembled in vitro. A subpopulation of pol IIo was localized to nuclear "speckle" domains enriched in splicing factors, indicating that it may also be associated with RNA processing in vivo. Moreover, pol IIo was retained in a similar pattern following in situ extraction of cells and was quantitatively recovered in the nuclear matrix fraction. The results implicate nuclear matrix-associated hyperphosphorylated pol IIo as a possible link in the coordination of transcription and splicing processes.