A comparison of random sequence reads versus 16S rDNA sequences for estimating the biodiversity of a metagenomic library

The construction of metagenomic libraries has permitted the study of microorganisms resistant to isolation and the analysis of 16S rDNA sequences has been used for over two decades to examine bacterial biodiversity. Here, we show that the analysis of random sequence reads (RSRs) instead of 16S is a suitable shortcut to estimate the biodiversity of a bacterial community from metagenomic libraries. We generated 10,010 RSRs from a metagenomic library of microorganisms found in human faecal samples. Then searched them using the program BLASTN against a prokaryotic sequence database to assign a taxon to each RSR. The results were compared with those obtained by screening and analysing the clones containing 16S rDNA sequences in the whole library. We found that the biodiversity observed by RSR analysis is consistent with that obtained by 16S rDNA. We also show that RSRs are suitable to compare the biodiversity between different metagenomic libraries. RSRs can thus provide a good estimate of the biodiversity of a metagenomic library and, as an alternative to 16S, this approach is both faster and cheaper.

Algebraic approach to optimal clone selection applied in metagenomic projects

Due to the wide diversity of unknown organisms in the environment, 99% of them cannot be grown in traditional culture medium in laboratories. Therefore, metagenomics projects are proposed to study microbial communities present in the environment, from molecular techniques, especially the sequencing. Thereby, for the coming years it is expected an accumulation of sequences produced by these projects. Thus, the sequences produced by genomics and metagenomics projects present several challenges for the treatment, storing and analysis such as: the search for clones containing genes of interest. This work presents the OCI Metagenomics, which allows defines and manages dynamically the rules of clone selection in metagenomic libraries, thought an algebraic approach based on process algebra. Furthermore, a web interface was developed to allow researchers to easily create and execute their own rules to select clones in genomic sequence database. This software has been tested in metagenomic cosmid library and it was able to select clones containing genes of interest. Copyright 2010 ACM.

Novel inorganic pyrophosphatase from soil metagenomic and family and subfamily prediction

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Construction and analysis of a metagenomic library from an enhanced biological phosphorus removal biomass

The enhanced biological phosphorus removal (EBPR) process is regularly used for the treatment of wastewater, but suffers from erratic performance. Successful EBPR relies on the growth of bacteria called polyphosphate-accumulating organisms (PAOs), which store phosphorus intracellularly as polyphosphate, thus removing it from wastewater. Metabolic models have been proposed which describe the measured chemical transformations, however genetic evidence is lacking to confirm these hypotheses. The aim of this research was to generate a metagenomic library from biomass enriched in PAOs as determined by phenotypic data and fluorescence in situ hybridisation (FISH) using probes specific for the only described PAO to date, Candidatus Accumulibacter phosphatis. DNA extraction methods were optimised and two fosmid libraries were constructed which contained 93 million base pairs of metagenomic data. Initial screening of the library for 16S rRNA genes revealed fosmids originating from a range of non-pure-cultured wastewater bacteria. The metagenomic libraries constructed will provide the ability to link phylogenetic and metabolic data for bacteria involved in nutrient removal from wastewater. Keywords DNA extraction; EBPR; metagenomic library; 16S rRNA gene.

Análise metagenômica da microbiota de ambientes aquáticos do estado do Rio Grande do Norte - Brasil

The screening for genes in metagenomic libraries from soil creates opportunities to explore the enormous genetic and metabolic diversity of microorganisms. Rivers are ecosystems with high biological diversity, but few were examined using the metagenomic approach. With this objective, a metagenomic library was constructed from DNA soil samples collected at three different points along the Jundiaí-river (Rio Grande do Norte-Brazil). The points sampled are from open area, rough terrain and with the direct incidence of sunlight. This library was analyzed functionally and based in sequence. For functional analysis Luria-Bertani solid medium (LB) with NaCl concentration varied from 0.17M to 0.85M was used for functional analysis. Positives clones resistant to hypersaline medium were obtained. The recombinant DNAs were extracted and transformed into Escherichia coli strain DH10B and survival curves were obtained for quantification of abiotic stress resistance. The sequences of clones were obtained and submitted to the BLASTX tool. Some clones were found to hypothetical proteins of microorganisms from both Archaea and Bacteria division. One of the clones showed a complete ORF with high similarity to glucose-6-phosphate isomerase which participates in the synthesis of glycerol pathway and serves as a compatible solute to balance the osmotic pressure inside and outside of cells. Subsequently, in order to identify genes encoding osmolytes or enzymes related halotolerance, environmental DNA samples from the river soil, from the water column of the estuary and ocean were collected and pyrosequenced. Sequences of osmolytes and enzymes of different microorganisms were obtained from the UniProt and used as RefSeqs for homology identification (TBLASTN) in metagenomic databases. The sequences were submitted to HMMER for the functional domains identification. Some enzymes were identified: alpha-trehalose-phosphate synthase, L-ectoina synthase (EctC), transaminase L-2 ,4-diaminobutyric acid (EctB), L-2 ,4-diaminobutyric acetyltransferase (EctA), L-threonine 3 dehydrogenase (sorbitol pathway), glycerol-3-phosphate dehydrogenase, inositol 3-phosphate dehydrogenase, chaperones, L-proline, glycine betaine binding ABC transporter, myo-inositol-1-phosphate synthase protein of proline simportadora / PutP sodium-and trehalose-6-phosphate phosphatase These proteins are commonly related to saline environments, however the identification of them in river environment is justified by the high salt concentration in the soil during prolonged dry seasons this river. Regarding the richness of the microbiota the river substrate has an abundance of halobacteria similar to the sea and more than the estuary. These data confirm the existence of a specialized response against salt stress by microorganisms in the environment of the Jundiaí river

Prospecção de genes de interesse biotecnológico : uma abordagem metagenômica

The total number of prokaryotic cells on Earth has been estimated at 4 to 6x1030 and only about 1% of microorganisms present in the environment can be cultivated by standard techniques of cultivation and plating. Therefore, it is a huge biological and genetic pool that can be exploited, for the identification and characterization of genes with biotechnological potential. Within this perspective, the metagenomics approach was applied in this work. Functional screening methods were performed aiming to identify new genes related to DNA repair and / or oxidative stress resistance, hydrocarbon degradation and hydrolytic activities (lipase, amylase and protease). Metagenomic libraries were built utilizing DNA extracted from soil samples collected in João Câmara RN. The libraries were analyzed functionally using specific substrate containing solid medium (hydrolytic activity), supplemented with H2O2 (DNA repair and / or resistance to oxidative stress) and liquid medium supplemented with light Arabian oil (activity, degradation of hydrocarbons). After confirmation of activity and exclusion of false-positive results, 49 clones were obtained, being 2 positive for amylase activity, 22 resistant to oxidative stress generated by H2O2 and 25 clones active for hydrocarbons degradation. Analysis of the sequences showed hypothetical proteins, dienelactona hydrolase, DNA polymerase, acetyltransferase, phosphotransferase, methyltransferase, endonucleases, among other proteins. The sequence data obtained matched with the functions tested, highlighting the success of metagenomics approaches combined with functional screening methods, leading to very promising results

Diversidade bacteriana em um latossolo sob cultivo intensivo e floresta através da análise metagenômica

Pós-graduação em Microbiologia Agropecuária - FCAV

Prospecção de genes com interesse biotecnológico

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

Expressão heteróloga de biossurfactantes identificados em bibliotecas metagenômicas

The microorganisms have a vast genetic diversity and they are present throughout the biosphere, however, only about 1% of the species can be cultivated by traditional cultivation techniques. Within this diversity there is a huge pool genetic and biological being explored. The metagenomics has enabled direct access to microbial genome derived from environmental samples using independent methods of cultivation. The methodology enables to obtain functional information about the proteins, as well as identify potential products with biotechnological interest and new industrially exploitable biological resources, such as new solutions to environmental impacts. Oil-contaminated areas are characterized by a large accumulation of hydrocarbons and surfactants may be used for bioremediation. Thus, the metagenomic approach was used in this study in order to select genes involved in the degradation and hydrocarbon emulsification. In a previous work, the environmental DNA (eDNA) was extracted from soil samples collected from two different areas (Caatinga and Saline River) of Rio Grande do Norte (Brazil), the metagenomic libraries were constructed and functionally analyzed. The clone able to degrade the oil was evaluated for the ability to synthesize biosurfactants. The sequence analysis revealed an ORF with 897 bp, 298 amino acids and a protein with around 34 kDa. The search for homology in GenBank revealed sequence similarity with a hypothetical protein of representatives Halobacteriaceae family, who were recently shown as strains producing biosurfactants. The presence of the inserted coding sequence and the acquired phenotype was confirmed. Primers were designed and the ORF amplified by PCR. The ORF was subcloned into pETDuet-1 expression vector for subsequent purification of the protein of interest containing a histidine tail. The tests performed to confirm the biosurfactant activity and the ability of hydrocarbon degradation showed positive results. The immunodetection test (western blot) using the monoclonal AntiHis® confirmed the presence of the environmental protein. This study was the first to report a possible protein with biosurfactant activity obtained from a metagenomic approach

Análise metagenômica da microbiota de ambientes aquáticos do estado do Rio Grande do Norte - Brasil

The screening for genes in metagenomic libraries from soil creates opportunities to explore the enormous genetic and metabolic diversity of microorganisms. Rivers are ecosystems with high biological diversity, but few were examined using the metagenomic approach. With this objective, a metagenomic library was constructed from DNA soil samples collected at three different points along the Jundiaí-river (Rio Grande do Norte-Brazil). The points sampled are from open area, rough terrain and with the direct incidence of sunlight. This library was analyzed functionally and based in sequence. For functional analysis Luria-Bertani solid medium (LB) with NaCl concentration varied from 0.17M to 0.85M was used for functional analysis. Positives clones resistant to hypersaline medium were obtained. The recombinant DNAs were extracted and transformed into Escherichia coli strain DH10B and survival curves were obtained for quantification of abiotic stress resistance. The sequences of clones were obtained and submitted to the BLASTX tool. Some clones were found to hypothetical proteins of microorganisms from both Archaea and Bacteria division. One of the clones showed a complete ORF with high similarity to glucose-6-phosphate isomerase which participates in the synthesis of glycerol pathway and serves as a compatible solute to balance the osmotic pressure inside and outside of cells. Subsequently, in order to identify genes encoding osmolytes or enzymes related halotolerance, environmental DNA samples from the river soil, from the water column of the estuary and ocean were collected and pyrosequenced. Sequences of osmolytes and enzymes of different microorganisms were obtained from the UniProt and used as RefSeqs for homology identification (TBLASTN) in metagenomic databases. The sequences were submitted to HMMER for the functional domains identification. Some enzymes were identified: alpha-trehalose-phosphate synthase, L-ectoina synthase (EctC), transaminase L-2 ,4-diaminobutyric acid (EctB), L-2 ,4-diaminobutyric acetyltransferase (EctA), L-threonine 3 dehydrogenase (sorbitol pathway), glycerol-3-phosphate dehydrogenase, inositol 3-phosphate dehydrogenase, chaperones, L-proline, glycine betaine binding ABC transporter, myo-inositol-1-phosphate synthase protein of proline simportadora / PutP sodium-and trehalose-6-phosphate phosphatase These proteins are commonly related to saline environments, however the identification of them in river environment is justified by the high salt concentration in the soil during prolonged dry seasons this river. Regarding the richness of the microbiota the river substrate has an abundance of halobacteria similar to the sea and more than the estuary. These data confirm the existence of a specialized response against salt stress by microorganisms in the environment of the Jundiaí river

Metagenomic Anlaysis of microsymbiont selection by the legume plant host

Rhizobium leguminosarum bv.viciae is able to establish nitrogen-fixing symbioses with legumes of the genera Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants (Laguerre et al., Young et al.) provided evidence that different plant hosts are able to select different rhizobial genotypes among those available in a given soil. However, these studies were necessarily limited by the paucity of relevant biodiversity markers. We have now reappraised this problem with the help of genomic tools. A well-characterized agricultural soil (INRA Bretennieres) was used as source of rhizobia. Plants of Pisum sativum, Lens culinaris, Vicia sativa and V. faba were used as traps. Isolates from 100 nodules were pooled, and DNA from each pool was sequenced (BGI-Hong Kong; Illumina Hiseq 2000, 500 bp PE libraries, 100 bp reads, 12 Mreads). Reads were quality filtered (FastQC, Trimmomatic), mapped against reference R. leguminosarum genomes (Bowtie2, Samtools), and visualized (IGV). An important fraction of the filtered reads were not recruited by reference genomes, suggesting that plant isolates contain genes that are not present in the reference genomes. For this study, we focused on three conserved genomic regions: 16S-23S rDNA, atpD and nodDABC, and a Single Nucleotide Polymorphism (SNP) analysis was carried out with meta / multigenomes from each plant. Although the level of polymorphism varied (lowest in the rRNA region), polymorphic sites could be identified that define the specific soil population vs. reference genomes. More importantly, a plant-specific SNP distribution was observed. This could be confirmed with many other regions extracted from the reference genomes (data not shown). Our results confirm at the genomic level previous observations regarding plant selection of specific genotypes. We expect that further, ongoing comparative studies on differential meta / multigenomic sequences will identify specific gene components of the plant-selected genotypes

Comparative Metagenomic Analysis Of Coral Microbial Communities Using A Reference-independent Approach.

By comparing the SEED and Pfam functional profiles of metagenomes of two Brazilian coral species with 29 datasets that are publicly available, we were able to identify some functions, such as protein secretion systems, that are overrepresented in the metagenomes of corals and may play a role in the establishment and maintenance of bacteria-coral associations. However, only a small percentage of the reads of these metagenomes could be annotated by these reference databases, which may lead to a strong bias in the comparative studies. For this reason, we have searched for identical sequences (99% of nucleotide identity) among these metagenomes in order to perform a reference-independent comparative analysis, and we were able to identify groups of microbial communities that may be under similar selective pressures. The identification of sequences shared among the metagenomes was found to be even better for the identification of groups of communities with similar niche requirements than the traditional analysis of functional profiles. This approach is not only helpful for the investigation of similarities between microbial communities with high proportion of unknown reads, but also enables an indirect overview of gene exchange between communities.

New Fefe-hydrogenase Genes Identified In A Metagenomic Fosmid Library From A Municipal Wastewater Treatment Plant As Revealed By High-throughput Sequencing.

A fosmid metagenomic library was constructed with total community DNA obtained from a municipal wastewater treatment plant (MWWTP), with the aim of identifying new FeFe-hydrogenase genes encoding the enzymes most important for hydrogen metabolism. The dataset generated by pyrosequencing of a fosmid library was mined to identify environmental gene tags (EGTs) assigned to FeFe-hydrogenase. The majority of EGTs representing FeFe-hydrogenase genes were affiliated with the class Clostridia, suggesting that this group is the main hydrogen producer in the MWWTP analyzed. Based on assembled sequences, three FeFe-hydrogenase genes were predicted based on detection of the L2 motif (MPCxxKxxE) in the encoded gene product, confirming true FeFe-hydrogenase sequences. These sequences were used to design specific primers to detect fosmids encoding FeFe-hydrogenase genes predicted from the dataset. Three identified fosmids were completely sequenced. The cloned genomic fragments within these fosmids are closely related to members of the Spirochaetaceae, Bacteroidales and Firmicutes, and their FeFe-hydrogenase sequences are characterized by the structure type M3, which is common to clostridial enzymes. FeFe-hydrogenase sequences found in this study represent hitherto undetected sequences, indicating the high genetic diversity regarding these enzymes in MWWTP. Results suggest that MWWTP have to be considered as reservoirs for new FeFe-hydrogenase genes.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes. Methods: Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries. Results: Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues. Conclusion: To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.