Gene capture coupled to high-throughput sequencing as a strategy for targeted metagenome exploration

Fungal dye-decolorizing peroxidase diversity: roles in either intra- or extracellular processes

Fungal dye-decolorizing peroxidases (DyPs) have found applications in the treatment of dye-contaminated industrial wastes or to improve biomass digestibility. Their roles in fungal biology are uncertain, although it has been repeatedly suggested that they could participate in lignin degradation and/or modification. Using a comprehensive set of 162 fully sequenced fungal species, we defined seven distinct fungal DyP clades on basis of a sequence similarity network. Sequences from one of these clades clearly diverged from all others, having on average the lower isoelectric points and hydropathy indices, the highest number of N-glycosylation sites, and N-terminal sequence peptides for secretion. Putative proteins from this clade are absent from brown-rot and ectomycorrhizal species that have lost the capability of degrading lignin enzymatically. They are almost exclusively present in white-rot and other saprotrophic Basidiomycota that digest lignin enzymatically, thus lending support for a specific role of DyPs from this clade in biochemical lignin modification. Additional nearly full-length fungal DyP genes were isolated from the environment by sequence capture by hybridization; they all belonged to the clade of the presumably secreted DyPs and to another related clade. We suggest focusing our attention on the presumably intracellular DyPs from the other clades, which have not been characterized thus far and could represent enzyme proteins with novel catalytic properties. KEY POINTS: • A fungal DyP phylogeny delineates seven main sequence clades. • Putative extracellular DyPs form a single clade of Basidiomycota sequences. • Extracellular DyPs are associated to white-rot fungi.

Revealing microbial species diversity using sequence capture by hybridization

Targeting small parts of the 16S rDNA phylogenetic marker by metabarcoding reveals microorganisms of interest but cannot achieve a taxonomic resolution at the species level, precluding further precise characterizations. To identify species behind operational taxonomic units (OTUs) of interest, even in the rare biosphere, we developed an innovative strategy using gene capture by hybridization. From three OTU sequences detected upon polyphenol supplementation and belonging to the rare biosphere of the human gut microbiota, we revealed 59 nearly full-length 16S rRNA genes, highlighting high bacterial diversity hidden behind OTUs while evidencing novel taxa. Inside each OTU, revealed 16S rDNA sequences could be highly distant from each other with similarities down to 85 %. We identified one new family belonging to the order Clostridiales, 39 new genera and 52 novel species. Related bacteria potentially involved in polyphenol degradation have also been identified through genome mining and our results suggest that the human gut microbiota could be much more diverse than previously thought.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

New Frontiers of Anaerobic Hydrocarbon Biodegradation in the Multi-Omics Era

The accumulation of petroleum hydrocarbons in the environment substantially endangers terrestrial and aquatic ecosystems. Many microbial strains have been recognized to utilize aliphatic and aromatic hydrocarbons under aerobic conditions. Nevertheless, most of these pollutants are transferred by natural processes, including rain, into the underground anaerobic zones where their degradation is much more problematic. In oxic zones, anaerobic microenvironments can be formed as a consequence of the intensive respiratory activities of (facultative) aerobic microbes. Even though aerobic bioremediation has been well-characterized over the past few decades, ample research is yet to be done in the field of anaerobic hydrocarbon biodegradation. With the emergence of high-throughput techniques, known as omics (e.g., genomics and metagenomics), the individual biodegraders, hydrocarbon-degrading microbial communities and metabolic pathways, interactions can be described at a contaminated site. Omics approaches provide the opportunity to examine single microorganisms or microbial communities at the system level and elucidate the metabolic networks, interspecies interactions during hydrocarbon mineralization. Metatranscriptomics and metaproteomics, for example, can shed light on the active genes and proteins and functional importance of the less abundant species. Moreover, novel unculturable hydrocarbon-degrading strains and enzymes can be discovered and fit into the metabolic networks of the community. Our objective is to review the anaerobic hydrocarbon biodegradation processes, the most important hydrocarbon degraders and their diverse metabolic pathways, including the use of various terminal electron acceptors and various electron transfer processes. The review primarily focuses on the achievements obtained by the current high-throughput (multi-omics) techniques which opened new perspectives in understanding the processes at the system level including the metabolic routes of individual strains, metabolic/electric interaction of the members of microbial communities. Based on the multi-omics techniques, novel metabolic blocks can be designed and used for the construction of microbial strains/consortia for efficient removal of hydrocarbons in anaerobic zones.

New insights into the pelagic microorganisms involved in the methane cycle in the meromictic Lake Pavin through metagenomics

Advances in metagenomics have given rise to the possibility of obtaining genome sequences from uncultured microorganisms, even for those poorly represented in the microbial community, thereby providing an important means to study their ecology and evolution. In this study, metagenomic sequencing was carried out at four sampling depths having different oxygen concentrations or environmental conditions in the water column of Lake Pavin. By analyzing the sequenced reads and matching the contigs to the proxy genomes of the closest cultivated relatives, we evaluated the metabolic potential of the dominant planktonic species involved in the methane cycle. We demonstrated that methane-producing communities were dominated by the genus Methanoregula while methane-consuming communities were dominated by the genus Methylobacter, thus confirming prior observations. Our work allowed the reconstruction of a draft of their core metabolic pathways. Hydrogenotrophs, the genes required for acetate activation in the methanogen genome, were also detected. Regarding methanotrophy, Methylobacter was present in the same areas as the non-methanotrophic, methylotrophic Methylotenera, which could suggest a relationship between these two groups. Furthermore, the presence of a large gene inventory for nitrogen metabolism (nitrate transport, denitrification, nitrite assimilation and nitrogen fixation, for instance) was detected in the Methylobacter genome.

Hybridization-based capture of pathogen mRNA enables paired host-pathogen transcriptional analysis

Dual transcriptional profiling of host and bacteria during infection is challenging due to the low abundance of bacterial mRNA. We report Pathogen Hybrid Capture (PatH-Cap), a method to enrich for bacterial mRNA and deplete bacterial rRNA simultaneously from dual RNA-seq libraries using transcriptome-specific probes. By addressing both the differential RNA content of the host relative to the infecting bacterium and the overwhelming abundance of uninformative structural RNAs (rRNA, tRNA) of both species in a single step, this approach enables analysis of very low-input RNA samples. By sequencing libraries before (pre-PatH-Cap) and after (post-PatH-Cap) enrichment, we achieve dual transcriptional profiling of host and bacteria, respectively, from the same sample. Importantly, enrichment preserves relative transcript abundance and increases the number of unique bacterial transcripts per gene in post-PatH-Cap libraries compared to pre-PatH-Cap libraries at the same sequencing depth, thereby decreasing the sequencing depth required to fully capture the transcriptional profile of the infecting bacteria. We demonstrate that PatH-Cap enables the study of low-input samples including single eukaryotic cells infected by 1-3 Pseudomonas aeruginosa bacteria and paired host-pathogen temporal gene expression analysis of Mycobacterium tuberculosis infecting macrophages. PatH-Cap can be applied to the study of a range of pathogens and microbial species, and more generally, to lowly-abundant species in mixed populations.

High-Level Abundances of Methanobacteriales and Syntrophobacterales May Help To Prevent Corrosion of Metal Sheet Piles

Iron sheet piles are widely used in flood protection, dike construction, and river bank reinforcement. Their corrosion leads to gradual deterioration and often makes replacement necessary. Natural deposit layers on these sheet piles can prevent degradation and significantly increase their life span. However, little is known about the mechanisms of natural protective layer formation. Here, we studied the microbially diverse populations of corrosion-protective deposit layers on iron sheet piles at the Gouderak pumping station in Zuid-Holland, the Netherlands. Deposit layers, surrounding sediment and top sediment samples were analyzed for soil physicochemical parameters, microbially diverse populations, and metabolic potential. Methanogens appeared to be enriched 18-fold in the deposit layers. After sequencing, metagenome assembly and binning, we obtained four nearly complete draft genomes of microorganisms (Methanobacteriales, two Coriobacteriales, and Syntrophobacterales) that were highly enriched in the deposit layers, strongly indicating a potential role in corrosion protection. Coriobacteriales and Syntrophobacterales could be part of a microbial food web degrading organic matter to supply methanogenic substrates. Methane-producing Methanobacteriales could metabolize iron, which may initially lead to mild corrosion but potentially stimulates the formation of a carbonate-rich protective deposit layer in the long term. In addition, Methanobacteriales and Coriobacteriales have the potential to interact with metal surfaces via direct interspecies or extracellular electron transfer. In conclusion, our study provides valuable insights into microbial populations involved in iron corrosion protection and potentially enables the development of novel strategies for in situ screening of iron sheet piles in order to reduce risks and develop more sustainable replacement practices.IMPORTANCE Iron sheet piles are widely used to reinforce dikes and river banks. Damage due to iron corrosion poses a significant safety risk and has significant economic impact. Different groups of microorganisms are known to either stimulate or inhibit the corrosion process. Recently, natural corrosion-protective deposit layers were found on sheet piles. Analyses of the microbial composition indicated a potential role for methane-producing archaea. However, the full metabolic potential of the microbial communities within these protective layers has not been determined. The significance of this work lies in the reconstruction of the microbial food web of natural corrosion-protective layers isolated from noncorroding metal sheet piles. With this work, we provide insights into the microbiological mechanisms that potentially promote corrosion protection in freshwater ecosystems. Our findings could support the development of screening protocols to assess the integrity of iron sheet piles to decide whether replacement is required.

Metabolic marker gene mining provides insight in global mcrA diversity and, coupled with targeted genome reconstruction, sheds further light on metabolic potential of the Methanomassiliicoccales

Over the past years, metagenomics has revolutionized our view of microbial diversity. Moreover, extracting near-complete genomes from metagenomes has led to the discovery of known metabolic traits in unsuspected lineages. Genome-resolved metagenomics relies on assembly of the sequencing reads and subsequent binning of assembled contigs, which might be hampered by strain heterogeneity or low abundance of a target organism. Here we present a complementary approach, metagenome marker gene mining, and use it to assess the global diversity of archaeal methane metabolism through the mcrA gene. To this end, we have screened 18,465 metagenomes for the presence of reads matching a database representative of all known mcrA proteins and reconstructed gene sequences from the matching reads. We use our mcrA dataset to assess the environmental distribution of the Methanomassiliicoccales and reconstruct and analyze a draft genome belonging to the ‘Lake Pavin cluster’, an uncultivated environmental clade of the Methanomassiliicoccales. Analysis of the ‘Lake Pavin cluster’ draft genome suggests that this organism has a more restricted capacity for hydrogenotrophic methylotrophic methanogenesis than previously studied Methanomassiliicoccales, with only genes for growth on methanol present. However, the presence of the soluble subunits of methyltetrahydromethanopterin:CoM methyltransferase (mtrAH) provide hypothetical pathways for methanol fermentation, and aceticlastic methanogenesis that await experimental verification. Thus, we show that marker gene mining can enhance the discovery power of metagenomics, by identifying novel lineages and aiding selection of targets for in-depth analyses. Marker gene mining is less sensitive to strain heterogeneity and has a lower abundance threshold than genome-resolved metagenomics, as it only requires short contigs and there is no binning step. Additionally, it is computationally cheaper than genome resolved metagenomics, since only a small subset of reads needs to be assembled. It is therefore a suitable approach to extract knowledge from the many publicly available sequencing projects.

Chloromethane formation and degradation in the fern phyllosphere

Chloromethane (CH₃Cl) is the most abundant halogenated trace gas in the atmosphere. It plays an important role in natural stratospheric ozone destruction. Current estimates of the global CH₃Cl budget are approximate. The strength of the CH₃Cl global sink by microbial degradation in soils and plants is under discussion. Some plants, particularly ferns, have been identified as substantial emitters of CH₃Cl. Their ability to degrade CH₃Cl remains uncertain. In this study, we investigated the potential of leaves from 3 abundant ferns (Osmunda regalis, Cyathea cooperi, Dryopteris filix-mas) to produce and degrade CH₃Cl by measuring their production and consumption rates and their stable carbon and hydrogen isotope signatures. Investigated ferns are able to degrade CH₃Cl at rates from 2.1 to 17 and 0.3 to 0.9μgg_dw^-1day^-1 for C. cooperi and D. filix-mas respectively, depending on CH₃Cl supplementation and temperature. The stable carbon isotope enrichment factor of remaining CH₃Cl was -39±13‰, whereas negligible isotope fractionation was observed for hydrogen (-8±19‰). In contrast, O. regalis did not consume CH₃Cl, but produced it at rates ranging from 0.6 to 128μgg_dw^-1day^-1, with stable isotope values of -97±8‰ for carbon and -202±10‰ for hydrogen, respectively. Even though the 3 ferns showed clearly different formation and consumption patterns, their leaf-associated bacterial diversity was not notably different. Moreover, we did not detect genes associated with the only known chloromethane utilization pathway "cmu" in the microbial phyllosphere of the investigated ferns. Our study suggests that still unknown CH₃Cl biodegradation processes on plants play an important role in global cycling of atmospheric CH₃Cl.

Comparing 16S rDNA amplicon sequencing and hybridization capture for pea aphid microbiota diversity analysis

OBJECTIVE

Targeted sequencing of 16S rDNA amplicons is routinely used for microbial community profiling but this method suffers several limitations such as bias affinity of universal primers and short read size. Gene capture by hybridization represents a promising alternative. Here we used a metagenomic extract from the pea aphid Acyrthosiphon pisum to compare the performances of two widely used PCR primer pairs with DNA capture, based on solution hybrid selection.

RESULTS

All methods produced an exhaustive description of the 8 bacterial taxa known to be present in this sample. In addition, the methods yielded similar quantitative results, with the number of reads strongly correlating with quantitative PCR controls. Both methods can thus be considered as qualitatively and quantitatively robust on such a sample with low microbial complexity.

Hybridization capture reveals microbial diversity missed using current profiling methods

BACKGROUND

Microorganisms comprise the majority of living organisms on our planet. For many years, exploration of the composition of microbial communities has been performed through the PCR-based study of the small subunit rRNA gene due to its high conservation across the domains of life. The application of this method has resulted in the discovery of many unexpected evolutionary lineages. However, amplicon sequencing is subject to numerous biases, with some taxa being missed, and is limited by the read length of second-generation sequencing platforms, which drastically reduces the phylogenetic resolution.

RESULTS

Here, we describe a hybridization capture strategy that allows the enrichment of 16S rRNA genes from metagenomic samples and enables an exhaustive identification and a complete reconstruction of the biomarker. Applying this approach to a microbial mock community and a soil sample, we demonstrated that hybridization capture is able to reveal greater microbial diversity than 16S rDNA amplicon sequencing and shotgun sequencing. The reconstruction of full-length 16S rRNA genes facilitated the improvement of phylogenetic resolution and the discovery of novel prokaryotic taxa.

CONCLUSIONS

Our results demonstrate that hybridization capture can lead to major breakthroughs in our understanding of microbial diversity, overcoming the limitations of conventional 16S rRNA gene studies. If applied to a broad range of environmental samples, this innovative approach could reveal the undescribed diversity of the still underexplored microbial communities and could provide a better understanding of ecosystem function.

Capture Hybridization of Long-Range DNA Fragments for High-Throughput Sequencing

Capture hybridization coupled with high-throughput sequencing (HTS) has become one of the most popular approaches to address some scientific problems not only for fundamental evolution but also for ecology and human disease in recent years. However, the technical problem of limited probe capture ability affects its widespread application. Here, we propose to capture hybridize long-range DNA fragments for HTS (termed LR-LCH). We provide a case of three amphibian samples to examine LR-LCH with 2 kb libraries and comparison of standard capture hybridization with 480 bp libraries. Capture sensitivity increased from an average 13.57% of standard capture hybridization to an average 19.80% of LR-LCH; capture efficiency also increased from an average 72.56% of standard capture hybridization to an average 97.71% of LR-LCH. These indicate that longer fragments in the library generally contain both relatively variable regions and relatively conservative regions. The divergent parts of target DNA are enriched along with conservative parts of DNA sequence that effectively captured during hybridization. We present a protocol that allows users to overcome the low capture sensitivity problem for high divergent regions.

Nutrient and acetate amendment leads to acetoclastic methane production and microbial community change in a non-producing Australian coal well

Coal mining is responsible for 11% of total anthropogenic methane emission thereby contributing considerably to climate change. Attempts to harvest coalbed methane for energy production are challenged by relatively low methane concentrations. In this study, we investigated whether nutrient and acetate amendment of a non-producing sub-bituminous coal well could transform the system to a methane source. We tracked cell counts, methane production, acetate concentration and geochemical parameters for 25 months in one amended and one unamended coal well in Australia. Additionally, the microbial community was analysed with 16S rRNA gene amplicon sequencing at 17 and 25 months after amendment and complemented by metagenome sequencing at 25 months. We found that cell numbers increased rapidly from 3.0 × 10⁴ cells ml^-1 to 9.9 × 10⁷ in the first 7 months after amendment. However, acetate depletion with concomitant methane production started only after 12-19 months. The microbial community was dominated by complex organic compound degraders (Anaerolineaceae, Rhodocyclaceae and Geobacter spp.), acetoclastic methanogens (Methanothrix spp.) and fungi (Agaricomycetes). Even though the microbial community had the functional potential to convert coal to methane, we observed no indication that coal was actually converted within the time frame of the study. Our results suggest that even though nutrient and acetate amendment stimulated relevant microbial species, it is not a sustainable way to transform non-producing coal wells into bioenergy factories.

Revealing large metagenomic regions through long DNA fragment hybridization capture

BACKGROUND

High-throughput DNA sequencing technologies have revolutionized genomic analysis, including the de novo assembly of whole genomes from single organisms or metagenomic samples. However, due to the limited capacity of short-read sequence data to assemble complex or low coverage regions, genomes are typically fragmented, leading to draft genomes with numerous underexplored large genomic regions. Revealing these missing sequences is a major goal to resolve concerns in numerous biological studies.

METHODS

To overcome these limitations, we developed an innovative target enrichment method for the reconstruction of large unknown genomic regions. Based on a hybridization capture strategy, this approach enables the enrichment of large genomic regions allowing the reconstruction of tens of kilobase pairs flanking a short, targeted DNA sequence.

RESULTS

Applied to a metagenomic soil sample targeting the linA gene, the biomarker of hexachlorocyclohexane (HCH) degradation, our method permitted the enrichment of the gene and its flanking regions leading to the reconstruction of several contigs and complete plasmids exceeding tens of kilobase pairs surrounding linA. Thus, through gene association and genome reconstruction, we identified microbial species involved in HCH degradation which constitute targets to improve biostimulation treatments.

CONCLUSIONS

This new hybridization capture strategy makes surveying and deconvoluting complex genomic regions possible through large genomic regions enrichment and allows the efficient exploration of metagenomic diversity. Indeed, this approach enables to assign identity and function to microorganisms in natural environments, one of the ultimate goals of microbial ecology.

Discovering Protein-Coding Genes from the Environment: Time for the Eukaryotes?

Eukaryotic microorganisms from diverse environments encompass a large number of taxa, many of them still unknown to science. One strategy to mine these organisms for genes of biotechnological relevance is to use a pool of eukaryotic mRNA directly extracted from environmental samples. Recent reports demonstrate that the resulting metatranscriptomic cDNA libraries can be screened by expression in yeast for a wide range of genes and functions from many of the different eukaryotic taxa. In combination with novel emerging high-throughput technologies, we anticipate that this approach should contribute to exploring the functional diversity of the eukaryotic microbiota.

The sequence capture by hybridization: a new approach for revealing the potential of mono-aromatic hydrocarbons bioattenuation in a deep oligotrophic aquifer

The formation water of a deep aquifer (853 m of depth) used for geological storage of natural gas was sampled to assess the mono-aromatic hydrocarbons attenuation potential of the indigenous microbiota. The study of bacterial diversity suggests that Firmicutes and, in particular, sulphate-reducing bacteria (Peptococcaceae) predominate in this microbial community. The capacity of the microbial community to biodegrade toluene and m- and p-xylenes was demonstrated using a culture-based approach after several hundred days of incubation. In order to reveal the potential for biodegradation of these compounds within a shorter time frame, an innovative approach named the solution hybrid selection method, which combines sequence capture by hybridization and next-generation sequencing, was applied to the same original water sample. The bssA and bssA-like genes were investigated as they are considered good biomarkers for the potential of toluene and xylene biodegradation. Unlike a PCR approach which failed to detect these genes directly from formation water, this innovative strategy demonstrated the presence of the bssA and bssA-like genes in this oligotrophic ecosystem, probably harboured by Peptococcaceae. The sequence capture by hybridization shows significant potential to reveal the presence of genes of functional interest which have low-level representation in the biosphere.

Sequence capture by hybridization reveals elusive hydrocarbon degradation potential

Highlight description of article (Ranchou‐Peyrouse et al.) to be published in MBT.

OCaPPI-Db: an oligonucleotide probe database for pathogen identification through hybridization capture

The detection and identification of bacterial pathogens involved in acts of bio- and agroterrorism are essential to avoid pathogen dispersal in the environment and propagation within the population. Conventional molecular methods, such as PCR amplification, DNA microarrays or shotgun sequencing, are subject to various limitations when assessing environmental samples, which can lead to inaccurate findings. We developed a hybridization capture strategy that uses a set of oligonucleotide probes to target and enrich biomarkers of interest in environmental samples. Here, we present Oligonucleotide Capture Probes for Pathogen Identification Database (OCaPPI-Db), an online capture probe database containing a set of 1,685 oligonucleotide probes allowing for the detection and identification of 30 biothreat agents up to the species level. This probe set can be used in its entirety as a comprehensive diagnostic tool or can be restricted to a set of probes targeting a specific pathogen or virulence factor according to the user’s needs.

Database URL: http://ocappidb.uca.works

Targeted Gene Capture by Hybridization to Illuminate Ecosystem Functioning

Microbial communities are extremely abundant and diverse on earth surface and play key role in the ecosystem functioning. Thus, although next-generation sequencing (NGS) technologies have greatly improved knowledge on microbial diversity, it is necessary to reduce the biological complexity to better understand the microorganism functions. To achieve this goal, we describe a promising approach, based on the solution hybrid selection (SHS) method for the selective enrichment in a target-specific biomarker from metagenomic and metatranscriptomic samples. The success of this method strongly depends on the determination of sensitive, specific, and explorative probes to assess the complete targeted gene repertoire. Indeed, in this method, RNA probes were used to capture large DNA or RNA fragments harboring biomarkers of interest that potentially allow to link structure and function of communities of interest.

Metagenomics of Thermophiles with a Focus on Discovery of Novel Thermozymes

Microbial populations living in environments with temperatures above 50°C (thermophiles) have been widely studied, increasing our knowledge in the composition and function of these ecological communities. Since these populations express a broad number of heat-resistant enzymes (thermozymes), they also represent an important source for novel biocatalysts that can be potentially used in industrial processes. The integrated study of the whole-community DNA from an environment, known as metagenomics, coupled with the development of next generation sequencing (NGS) technologies, has allowed the generation of large amounts of data from thermophiles. In this review, we summarize the main approaches commonly utilized for assessing the taxonomic and functional diversity of thermophiles through metagenomics, including several bioinformatics tools and some metagenome-derived methods to isolate their thermozymes.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

A metabarcoding framework for facilitated survey of endolithic phototrophs with tufA

Background

In spite of their ecological importance as primary producers and microbioeroders of marine calcium carbonate (CaCO₃) substrata, endolithic phototrophs spanning both prokaryotic (the cyanobacteria) and eukaryotic algae lack established molecular resources for their facilitated survey with high throughput sequencing. Here, the development of a metabarcoding framework for the elongation factor EF-Ttu (tufA) was tested on four Illumina-sequenced marine CaCO₃ microfloras for the characterization of their endolithic phototrophs, especially the abundant bioeroding Ostreobium spp. (Ulvophyceae). The framework consists of novel tufA degenerate primers and a comprehensive database enabling Operational Taxonomic Unit (OTU) identification at multiple taxonomic ranks with percent identity thresholds determined herein.

Results

The newly established tufA database comprises 4057 non-redundant sequences (from 1339 eukaryotic and prokaryotic phototrophs, and 2718 prokaryotic heterotrophs) including 27 classes in 10 phyla of phototrophic diversity summarized from data mining on GenBank^®, our barcoding of >150 clones produced from coral reef microfloras, and >300 eukaryotic phototrophs (>230 Ulvophyceae including >100 ‘Ostreobium’ spp., and >70 Florideophyceae, Phaeophyceae and miscellaneous taxa). Illumina metabarcoding with the newly designed primers resulted in 802 robust OTUs including 618 phototrophs and 184 heterotrophs (77 and 23 % of OTUs, respectively). Phototrophic OTUs belonged to 14 classes of phototrophs found in seven phyla, and represented ~98 % of all reads. The phylogenetic profiles of coral reef microfloras showed few OTUs in large abundance (proportion of reads) for the Chlorophyta (Ulvophyceae, i.e. Ostreobium and Phaeophila), the Rhodophyta (Florideophyceae) and Haptophyta (Coccolithophyceae), and a large diversity (richness) of OTUs in lower abundance for the Cyanophyta (Cyanophyceae) and the Ochrophyta (the diatoms, ‘Bacillariophyta’). The bioerosive ‘Ostreobium’ spp. represented four families in a large clade of subordinal divergence, i.e. the Ostreobidineae, and a fifth, phylogenetically remote family in the suborder Halimedineae (provisionally assigned as the ‘Pseudostreobiaceae’). Together they harbor 85–95 delimited cryptic species of endolithic microsiphons.

Conclusions

The novel degenerate primers allowed for amplification of endolithic phototrophs across a wide phylogenetic breadth as well as their recovery in very large proportions of reads (overall 98 %) and diversity (overall 77 % of OTUs). The established companion tufA database and determined identity thresholds allow for OTU identification at multiple taxonomic ranks to facilitate the monitoring of phototrophic assemblages via metabarcoding, especially endolithic communities rich in bioeroding Ulvophyceae, such as those harboring ‘Ostreobium’ spp., Phaeophila spp. and associated algal diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-016-0068-x) contains supplementary material, which is available to authorized users.

Distribution of Dehalococcoidia in the Anaerobic Deep Water of a Remote Meromictic Crater Lake and Detection of Dehalococcoidia-Derived Reductive Dehalogenase Homologous Genes

Here we describe the natural occurrence of bacteria of the class Dehalococcoidia (DEH) and their diversity at different depths in anoxic waters of a remote meromictic lake (Lake Pavin) using 16S rRNA gene amplicon sequencing and quantitative PCR. Detected DEH are phylogenetically diverse and the majority of 16S rRNA sequences have less than 91% similarity to previously isolated DEH 16S rRNA sequences. To predict the metabolic potential of detected DEH subgroups and to assess if they encode genes to transform halogenated compounds, we enriched DEH-affiliated genomic DNA by using a specific-gene capture method and probes against DEH-derived 16S rRNA genes, reductive dehalogenase genes and known insertion sequences. Two reductive dehalogenase homologous sequences were identified from DEH-enriched genomic DNA, and marker genes in the direct vicinity confirm that gene fragments were derived from DEH. The low sequence similarity with known reductive dehalogenase genes suggests yet-unknown catabolic potential in the anoxic zone of Lake Pavin.

Targeted capture in evolutionary and ecological genomics

The rapid expansion of next-generation sequencing has yielded a powerful array of tools to address fundamental biological questions at a scale that was inconceivable just a few years ago. Various genome-partitioning strategies to sequence select subsets of the genome have emerged as powerful alternatives to whole-genome sequencing in ecological and evolutionary genomic studies. High-throughput targeted capture is one such strategy that involves the parallel enrichment of preselected genomic regions of interest. The growing use of targeted capture demonstrates its potential power to address a range of research questions, yet these approaches have yet to expand broadly across laboratories focused on evolutionary and ecological genomics. In part, the use of targeted capture has been hindered by the logistics of capture design and implementation in species without established reference genomes. Here we aim to (i) increase the accessibility of targeted capture to researchers working in nonmodel taxa by discussing capture methods that circumvent the need of a reference genome, (ii) highlight the evolutionary and ecological applications where this approach is emerging as a powerful sequencing strategy and (iii) discuss the future of targeted capture and other genome-partitioning approaches in the light of the increasing accessibility of whole-genome sequencing. Given the practical advantages and increasing feasibility of high-throughput targeted capture, we anticipate an ongoing expansion of capture-based approaches in evolutionary and ecological research, synergistic with an expansion of whole-genome sequencing.

Metatranscriptomics of Soil Eukaryotic Communities

Functions expressed by eukaryotic organisms in soil can be specifically studied by analyzing the pool of eukaryotic-specific polyadenylated mRNA directly extracted from environmental samples. In this chapter, we describe two alternative protocols for the extraction of high-quality RNA from soil samples. Total soil RNA or mRNA can be converted to cDNA for direct high-throughput sequencing. Polyadenylated mRNA-derived full-length cDNAs can also be cloned in expression plasmid vectors to constitute soil cDNA libraries, which can be subsequently screened for functional gene categories. Alternatively, the diversity of specific gene families can also be explored following cDNA sequence capture using exploratory oligonucleotide probes.

Hunting Down Frame Shifts: Ecological Analysis of Diverse Functional Gene Sequences

Functional gene ecological analyses using amplicon sequencing can be challenging as translated sequences are often burdened with shifted reading frames. The aim of this work was to evaluate several bioinformatics tools designed to correct errors which arise during sequencing in an effort to reduce the number of frameshifts (FS). Genes encoding for alpha subunits of biphenyl (bphA) and benzoate (benA) dioxygenases were used as model sequences. FrameBot, a FS correction tool, was able to reduce the number of detected FS to zero. However, up to 44% of sequences were discarded by FrameBot as non-specific targets. Therefore, we proposed a de novo mode of FrameBot for FS correction, which works on a similar basis as common chimera identifying platforms and is not dependent on reference sequences. By nature of FrameBot de novo design, it is crucial to provide it with data as error free as possible. We tested the ability of several publicly available correction tools to decrease the number of errors in the data sets. The combination of maximum expected error filtering and single linkage pre-clustering proved to be the most efficient read processing approach. Applying FrameBot de novo on the processed data enabled analysis of BphA sequences with minimal losses of potentially functional sequences not homologous to those previously known. This experiment also demonstrated the extensive diversity of dioxygenases in soil. A script which performs FrameBot de novo is presented in the supplementary material to the study or available at https://github.com/strejcem/FBdenovo. The tool was also implemented into FunGene Pipeline available at http://fungene.cme.msu.edu/FunGenePipeline/.

Captured metagenomics: large-scale targeting of genes based on 'sequence capture' reveals functional diversity in soils

Microbial enzyme diversity is a key to understand many ecosystem processes. Whole metagenome sequencing (WMG) obtains information on functional genes, but it is costly and inefficient due to large amount of sequencing that is required. In this study, we have applied a captured metagenomics technique for functional genes in soil microorganisms, as an alternative to WMG. Large-scale targeting of functional genes, coding for enzymes related to organic matter degradation, was applied to two agricultural soil communities through captured metagenomics. Captured metagenomics uses custom-designed, hybridization-based oligonucleotide probes that enrich functional genes of interest in metagenomic libraries where only probe-bound DNA fragments are sequenced. The captured metagenomes were highly enriched with targeted genes while maintaining their target diversity and their taxonomic distribution correlated well with the traditional ribosomal sequencing. The captured metagenomes were highly enriched with genes related to organic matter degradation; at least five times more than similar, publicly available soil WMG projects. This target enrichment technique also preserves the functional representation of the soils, thereby facilitating comparative metagenomics projects. Here, we present the first study that applies the captured metagenomics approach in large scale, and this novel method allows deep investigations of central ecosystem processes by studying functional gene abundances.

Capturing prokaryotic dark matter genomes

Prokaryotes are the most diverse and abundant cellular life forms on Earth. Most of them, identified by indirect molecular approaches, belong to microbial dark matter. The advent of metagenomic and single-cell genomic approaches has highlighted the metabolic capabilities of numerous members of this dark matter through genome reconstruction. Thus, linking functions back to the species has revolutionized our understanding of how ecosystem function is sustained by the microbial world. This review will present discoveries acquired through the illumination of prokaryotic dark matter genomes by these innovative approaches.

MetCap: a bioinformatics probe design pipeline for large-scale targeted metagenomics

Background

Massive sequencing of genes from different environments has evolved metagenomics as central to enhancing the understanding of the wide diversity of micro-organisms and their roles in driving ecological processes. Reduced cost and high throughput sequencing has made large-scale projects achievable to a wider group of researchers, though complete metagenome sequencing is still a daunting task in terms of sequencing as well as the downstream bioinformatics analyses. Alternative approaches such as targeted amplicon sequencing requires custom PCR primer generation, and is not scalable to thousands of genes or gene families.

Results

In this study, we are presenting a web-based tool called MetCap that circumvents the limitations of amplicon sequencing of multiple genes by designing probes that are suitable for large-scale targeted metagenomics sequencing studies. MetCap provides a novel approach to target thousands of genes and genomic regions that could be used in targeted metagenomics studies. Automatic analysis of user-defined sequences is performed, and probes specifically designed for metagenome studies are generated. To illustrate the advantage of a targeted metagenome approach, we have generated more than 300,000 probes that match more than 400,000 publicly available sequences related to carbon degradation, and used these probes for target sequencing in a soil metagenome study. The results show high enrichment of target genes and a successful capturing of the majority of gene families. MetCap is freely available to users from: http://soilecology.biol.lu.se/metcap/.

Conclusion

MetCap is facilitating probe-based target enrichment as an easy and efficient alternative tool compared to complex primer-based enrichment for large-scale investigations of metagenomes. Our results have shown efficient large-scale target enrichment through MetCap-designed probes for a soil metagenome. The web service is suitable for any targeted metagenomics project that aims to study several genes simultaneously. The novel bioinformatics approach taken by the web service will enable researchers in microbial ecology to tap into the vast diversity of microbial communities using targeted metagenomics as a cost-effective alternative to whole metagenome sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0501-8) contains supplementary material, which is available to authorized users.

Solution hybrid selection capture for the recovery of functional full-length eukaryotic cDNAs from complex environmental samples

Eukaryotic microbial communities play key functional roles in soil biology and potentially represent a rich source of natural products including biocatalysts. Culture-independent molecular methods are powerful tools to isolate functional genes from uncultured microorganisms. However, none of the methods used in environmental genomics allow for a rapid isolation of numerous functional genes from eukaryotic microbial communities. We developed an original adaptation of the solution hybrid selection (SHS) for an efficient recovery of functional complementary DNAs (cDNAs) synthesized from soil-extracted polyadenylated mRNAs. This protocol was tested on the Glycoside Hydrolase 11 gene family encoding endo-xylanases for which we designed 35 explorative 31-mers capture probes. SHS was implemented on four soil eukaryotic cDNA pools. After two successive rounds of capture, >90% of the resulting cDNAs were GH11 sequences, of which 70% (38 among 53 sequenced genes) were full length. Between 1.5 and 25% of the cloned captured sequences were expressed in Saccharomyces cerevisiae. Sequencing of polymerase chain reaction-amplified GH11 gene fragments from the captured sequences highlighted hundreds of phylogenetically diverse sequences that were not yet described, in public databases. This protocol offers the possibility of performing exhaustive exploration of eukaryotic gene families within microbial communities thriving in any type of environment.

Exploiting the architecture and the features of the microsporidian genomes to investigate diversity and impact of these parasites on ecosystems

Fungal species play extremely important roles in ecosystems. Clustered at the base of the fungal kingdom are Microsporidia, a group of obligate intracellular eukaryotes infecting multiple animal lineages. Because of their large host spectrum and their implications in host population regulation, they influence food webs, and accordingly, ecosystem structure and function. Unfortunately, their ecological role is not well understood. Present also as highly resistant spores in the environment, their characterisation requires special attention. Different techniques based on direct isolation and/or molecular approaches can be considered to elucidate their role in the ecosystems, but integrating environmental and genomic data (for example, genome architecture, core genome, transcriptional and translational signals) is crucial to better understand the diversity and adaptive capacities of Microsporidia. Here, we review the current status of Microsporidia in trophic networks; the various genomics tools that could be used to ensure identification and evaluate diversity and abundance of these organisms; and how these tools could be used to explore the microsporidian life cycle in different environments. Our understanding of the evolution of these widespread parasites is currently impaired by limited sampling, and we have no doubt witnessed but a small subset of their diversity.

Target capture and massive sequencing of genes transcribed in Mytilus galloprovincialis

Next generation sequencing (NGS) allows fast and massive production of both genome and transcriptome sequence datasets. As the genome of the Mediterranean mussel Mytilus galloprovincialis is not available at present, we have explored the possibility of reducing the whole genome sequencing efforts by using capture probes coupled with PCR amplification and high-throughput 454-sequencing to enrich selected genomic regions. The enrichment of DNA target sequences was validated by real-time PCR, whereas the efficacy of the applied strategy was evaluated by mapping the 454-output reads against reference transcript data already available for M. galloprovincialis and by measuring coverage, SNPs, number of de novo sequenced introns, and complete gene sequences. Focusing on a target size of nearly 1.5 Mbp, we obtained a target coverage which allowed the identification of more than 250 complete introns, 10,741 SNPs, and also complete gene sequences. This study confirms the transcriptome-based enrichment of gDNA regions as a good strategy to expand knowledge on specific subsets of genes also in nonmodel organisms.

Phylogenomic data support a seventh order of Methylotrophic methanogens and provide insights into the evolution of Methanogenesis

Increasing evidence from sequence data from various environments, including the human gut, suggests the existence of a previously unknown putative seventh order of methanogens. The first genomic data from members of this lineage, Methanomassiliicoccus luminyensis and “Candidatus Methanomethylophilus alvus,” provide insights into its evolutionary history and metabolic features. Phylogenetic analysis of ribosomal proteins robustly indicates a monophyletic group independent of any previously known methanogenic order, which shares ancestry with the Marine Benthic Group D, the Marine Group II, the DHVE2 group, and the Thermoplasmatales. This phylogenetic position, along with the analysis of enzymes involved in core methanogenesis, strengthens a single ancient origin of methanogenesis in the Euryarchaeota and indicates further multiple independent losses of this metabolism in nonmethanogenic lineages than previously suggested. Genomic analysis revealed an unprecedented loss of the genes coding for the first six steps of methanogenesis from H₂/CO₂ and the oxidative part of methylotrophic methanogenesis, consistent with the fact that M. luminyensis and “Ca. M. alvus” are obligate H₂-dependent methylotrophic methanogens. Genomic data also suggest that these methanogens may use a large panel of methylated compounds. Phylogenetic analysis including homologs retrieved from environmental samples indicates that methylotrophic methanogenesis (regardless of dependency on H₂) is not restricted to gut representatives but may be an ancestral characteristic of the whole order, and possibly also of ancient origin in the Euryarchaeota. 16S rRNA and McrA trees show that this new order of methanogens is very diverse and occupies environments highly relevant for methane production, therefore representing a key lineage to fully understand the diversity and evolution of methanogenesis.

PhylOPDb: a 16S rRNA oligonucleotide probe database for prokaryotic identification

In recent years, high-throughput molecular tools have led to an exponential growth of available 16S rRNA gene sequences. Incorporating such data, molecular tools based on target-probe hybridization were developed to monitor microbial communities within complex environments. Unfortunately, only a few 16S rRNA gene-targeted probe collections were described. Here, we present PhylOPDb, an online resource for a comprehensive phylogenetic oligonucleotide probe database. PhylOPDb provides a convivial and easy-to-use web interface to browse both regular and explorative 16S rRNA-targeted probes. Such probes set or subset could be used to globally monitor known and unknown prokaryotic communities through various techniques including DNA microarrays, polymerase chain reaction (PCR), fluorescent in situ hybridization (FISH), targeted gene capture or in silico rapid sequence identification. PhylOPDb contains 74 003 25-mer probes targeting 2178 genera including Bacteria and Archaea. Database URL: http://g2im.u-clermont1.fr/phylopdb/