Functional metagenomics reveals novel β-galactosidases not predictable from gene sequences

Abstract

The techniques of metagenomics have allowed researchers to access the genomic potential of uncultivated microbes, but there remain significant barriers to determination of gene function based on DNA sequence alone. Functional metagenomics, in which DNA is cloned and expressed in surrogate hosts, can overcome these barriers, and make important contributions to the discovery of novel enzymes. In this study, a soil metagenomic library carried in an IncP cosmid was used for functional complementation for β-galactosidase activity in both Sinorhizobium meliloti (α-Proteobacteria) and Escherichia coli (γ-Proteobacteria) backgrounds. One β-galactosidase, encoded by six overlapping clones that were selected in both hosts, was identified as a member of glycoside hydrolase family 2. We could not identify ORFs obviously encoding possible β-galactosidases in 19 other sequenced clones that were only able to complement S. meliloti. Based on low sequence identity to other known glycoside hydrolases, yet not β-galactosidases, three of these ORFs were examined further. Biochemical analysis confirmed that all three encoded β-galactosidase activity. Lac36W_ORF11 and Lac161_ORF7 had conserved domains, but lacked similarities to known glycoside hydrolases. Lac161_ORF10 had neither conserved domains nor similarity to known glycoside hydrolases. Bioinformatic and structural modeling implied that Lac161_ORF10 protein represented a novel enzyme family with a five-bladed propeller glycoside hydrolase domain. By discovering founding members of three novel β-galactosidase families, we have reinforced the value of functional metagenomics for isolating novel genes that could not have been predicted from DNA sequence analysis alone.

About the Authors:

Jiujun Cheng

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

Tatyana Romantsov

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

Katja Engel

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

Andrew C. Doxey

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

David R. Rose

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

Josh D. Neufeld

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

Trevor C. Charles

* E-mail: [email protected]

Affiliation: Department of Biology, University of Waterloo, Waterloo, ON, Canada

ORCID http://orcid.org/0000-0002-0344-5932

Introduction

Soils harbour the greatest genetic diversity of any habitats on Earth [1]. Our knowledge of microorganisms comprising soil communities is hampered by cultivation challenges for many microorganisms in these communities [2], although improvements in cultivation methods are addressing this bottleneck [3]. The genomes of metabolically versatile soil microbes are potential sources of biocatalysts for use in various industrial processes. Limited knowledge of links between sequence and function prevent rapid progress in bioinformatics-based systems biology. Metagenomics can be used to explore the collective genetic constituency of environmental microbes, including those that are difficult to culture through conventional microbiological techniques. Sequence-based and function-based strategies are used in metagenomics, depending on the main objectives of the particular study. Sequence-based metagenomics identifies genes by sequence similarity to known database sequences. However, it is difficult, if not impossible, to reliably predict the function of truly novel genes without experimental evidence. Functional screening strategies are based on phenotypic detection of the desired activity, heterologous complementation of host strains, and induced gene expression [2, 4, 5]. These experimental activities have identified novel genes showing little similarity to genes of known function [6–9]. In addition, heterologous complementation screening strategies facilitate simultaneous screening of millions of metagenomic clones. Most functional screens are performed in Escherichia coli of the γ-Proteobacteria. Because gene expression is often host-dependent [10], multi-host systems increase the likelihood of successful gene expression [5, 11–19].

Glycoside hydrolases (GH) hydrolyze the glycosidic linkages of glycosides and oligosaccharides, and are classified into 131 families based on the similarity of amino acid sequences [20]; http://www.cazy.org/Glycoside-Hydrolases.html). The β-galactosidase (EC 3.2.1.23) enzymes are grouped within GH1, GH2, GH35, GH42, and GH59 families. The hydrolytic activity of β-galactosidase can contribute to various applications, such as reducing the lactose content of dairy products [21], producing bioethanol from cheese whey [22], and as the basis of biosensors for detection of lactose [23]. The associated transgalactosylation activity of β-galactosidases can be used to synthesize galactosylated products [24]. Functional screening of metagenomic libraries resulted in discovery of a GH43 enzyme acting on multiple substrates including lactose [25], cold-adapted or thermostable β-galactosidases [26–28], GH1/GH2 [29–33], two novel β-galactosidases without any similarity to known GHs [7], glycosyltransferase family 4 and β-glycosidase with β-galactosidase activity respectively [34, 35].

In this study, we demonstrate the value of metagenomic cosmid libraries for enzyme discovery. Using lactose as the sole carbon source to support growth of Sinorhizobium meliloti, we identified three new β-galactosidases from a cornfield soil library and characterized the biochemical properties of these novel enzymes. In doing this, we revealed novel protein sequence space associations with β-galactosidase activity and substrate specificity.

Materials and methods

Bacterial strains, plasmids, cosmids, and growth conditions

Several bacterial strains, plasmids, and cosmids were used in this study (Table 1). All E. coli strains were grown at 37°C in LB medium (1% tryptone, 0.5% yeast extract, 0.5% NaCl, pH 7.0). All S. meliloti strains were grown at 30°C in LB supplemented with 2.5 mM CaCl2 and 2.5 mM MgSO4 (LBmc) [36]. Antibiotics were used at the following final concentrations: streptomycin (Sm, 100 μg/ml for E. coli, 200 μg/ml for S. meliloti), neomycin (Nm, 200 μg/ml), rifampicin (Rif, 100 μg/ml), kanamycin (Km, 50 μg/ml), tetracycline (Tc, 20 μg/ml for E. coli, 10 μg/ml for S. meliloti), gentamicin (Gm, 10 μg/ml for E. coli), chloramphenicol (Cm, 10 μg/ml for E. coli).

[Figure omitted. See PDF.]

Table 1. Bacterial strains, plasmids, and cosmids.

https://doi.org/10.1371/journal.pone.0172545.t001

Functional screening of β-galactosidases

The procedure for functional selection and screen of β-galactosidases is illustrated in Fig 1. A previously pooled metagenomic library of corn field soil DNA (12AC; [18]) was used for screening β-galactosidases. The 12AC library contains ~7.9 × 104 cosmid clones with average sizes of 33.4 kb. The cosmid DNA was isolated from the pooled library clones using GeneJET Plasmid Miniprep Kit (Thermo Scientific). E. coli DH5α (lacZYA) was grown in LB to an OD600 of 0.6. Cells were collected by centrifugation at 4°C and at 12,300 × g for 20 min and then washed three times with cold 10% glycerol. Cells were gently suspended in 2 ml of ice-cold 10% glycerol to about 3 × 1010 cells/ml. Electrocompetent cell volumes of 40 μl were mixed with 1 μl of cosmid DNA (45 ng) in a cold 1.5-ml microcentrifuge tube on ice, and then transferred to a cold electroporation cuvette (0.2 cm, Bio-Rad). Electroporation was performed using a Gene Pulser (Bio-Rad; C = 25 μF; PC = 200 ohm; V = 3.0 kV). Liquid SOC medium (1 ml) was added to the cuvette after one pulse. Electroporated cells were transferred to a 1.5-ml microcentrifuge tube and incubated at 37°C in a water bath for 30 min, inverting the tube every 5 min. The tube was then shaken at 200 rpm for 30 min at 37°C. Following concentration by centrifugation, cells were spread on LB Tc plates, and incubated overnight at 37°C. Multiple electroporations were performed to obtain the desired numbers of recombinant E. coli DH5α clones. Transformants were pooled, then stored at -75°C after addition of DMSO (7% final concentration).

[Figure omitted. See PDF.]

Fig 1. Functional metagenomics selection/screen and characterization of novel β-galactosidases.

https://doi.org/10.1371/journal.pone.0172545.g001

The pooled E. coli DH5α cosmid clones were washed three times with 0.85% NaCl and then spread on defined M9 medium [43] supplemented with L-arginine (50 μg/ml), thiamine (10 μg/ml), tetracycline (15 μg/ml), lactose (15 mM) as the sole carbon source, and the chromogenic substrate 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-gal; 36 μg/ml). Plates were incubated at 37°C for 1–3 days. Positive blue colonies were streak purified one time on M9-lactose plates. The Lac+ clones were inoculated in 3 ml of LB Tc medium and grown overnight at 37°C. Cosmid DNA was isolated using the GeneJET Plasmid Miniprep Kit (Thermo Scientific), digested simultaneously with EcoRI-BamHI-HindIII, then resolved on 1% agarose gels. Cosmids representative of distinct restriction patterns were retransformed into E. coli DH5α, and then spread on M9-lactose to confirm the Lac+ phenotype.

To screen for Lac+ clones in S. meliloti, 12AC cosmids were conjugated from E. coli DH5α into S. meliloti RmF728 (lacEFGZ1K1) [38]) with helper plasmid pRK600. The pooled 12AC library clones of 0.25 ml were mixed with 2 ml each of overnight-grown S. meliloti RmF728 and E. coli DH5α (pRK600). Cells were collected by centrifugation at 12,300 × g for 3 min, washed twice with 2 ml of 0.85% NaCl, and then suspended in 0.5 ml of 0.85% NaCl. Mixed cells were spotted on LB agar and incubated overnight at 30°C. Following collection of the mating spot in a 2-ml microtube and washing twice with 0.85% NaCl, the conjugation mixture was serially diluted and plated on defined M9 medium (Nm Tc) that was supplemented with biotin (0.3 μg/ml), thiamine (10 μg/ml), X-gal (36 μg/ml), and lactose (15 mM) as the sole carbon source. Lac+ colonies were streak purified once on M9 lactose plates. Cosmids were then transferred from S. meliloti to E. coli DH5α (RifR) via conjugation with the helper plasmid pRK600. E. coli DH5α carrying the empty cosmid pJC8 was used as a negative control. Lac+ cosmid DNA was prepared and analyzed by EcoRI-BamHI-HindIII digestion as described previously.

The Lac+ phenotypes of S. meliloti strains were also verified by assaying β-galactosidase activity. Thirty-nine random Lac+ strains were grown overnight in LBmc, washed twice with 0.85% NaCl, and then subcultured (1:200 dilution) in M9 lactose medium. Following growth for 48 h, β-galactosidase activity was measured using o-nitrophenyl β-galactoside (ONPG) as described previously [44]. Briefly, whole cells were permeabilized with SDS and incubated with ONPG for 1 h at room temperature. Following termination of the reaction with Na2CO3, absorption at 420 nm was recorded. Specific activity was represented as (1000 × A420)/(time × OD600 × culture volume).

Cloning, expression and purification

The KOD Xtreme DNA polymerase (Novagen) was used for all PCR amplifications with several different primers (S1 Table). PCR amplifications consisted of one cycle of 94°C for 5 min, 30 cycles of 94°C for 30 s, 50–57°C for 30 s, 68°C for 30 s to 3 min, and incubation at 68°C for 10 min. The Lac161_ORF10 was PCR amplified using primers lac161NdeI and lac161HindIII, and cloned into the NdeI-HindIII sites in pET-30a(+) to obtain pTR5. To clone putative GH genes with a C-terminal His tag in a broad-host-range plasmid, a 0.37-kb DNA fragment containing the NdeI site to the end of T7 terminator from pET-30b(+) was amplified using primers JC226 and JC227, and inserted into the NdeI-NheI sites in pSRKGm [42] to obtain plasmid pJC98. The Lac161_ORF7 was obtained by PCR amplification using primers JC220 and JC221, and then inserted into the NdeI-SalI sites in pJC98, yielding pJC102. Lac36W_ORF11 was PCR amplified using primer pair JC212 and JC213, and then cloned into the NdeI-XhoI sites in pJC98 to obtain plasmid pJC97. Plasmids were verified by restriction enzyme digestion analysis.

The expression plasmids pTR5, pJC97, and pJC102 were introduced into E. coli BL21(DE3)pLysS using the CaCl2 transformation method. Gene overexpression was induced by adding 0.1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 20°C for 16 h. Cell pellets were suspended in a lysis buffer (100 mM potassium phosphate; pH 7.4, 5 mM MgSO4, 30 μg/ml DNase, 1 mg/ml lysozyme, 2 mM β-mercaptoethanol, and 0.5 mM phenylmethylsulfonyl fluoride), incubated on ice for 30 min, and then disrupted by three passes through EmulsiFlex-C3 (Avestin Inc. Ottawa, Ontario) pressure cell at an internal cell pressure of 1.6 × 108 Pa. His-tagged proteins were purified from supernatants of cell extracts under native conditions using Co2+-NTA affinity chromatography (Clontech Laboratories). Purified proteins were dialyzed twice at 4°C against 50 mM potassium phosphate and 10 mM Tris-HCl (pH 7.4).

Characterization of β-galactosidases

Enzyme activities were measured using a Glucose Oxidase Activity Assay Kit (Sigma-Aldrich) for quantifying the amount of glucose produced upon adding different concentrations of enzyme and lactose. Assays in triplicate were carried out in 96-well microtiter plates The following buffers were used for profiling optimal pH: 100 mM 2-(N-morpholino)ethanesulfonic acid (Mes, pH 5.5–6.5, 100 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (pH 7.0–7.5), and 100 mM Tris-HCL (pH 8.0–9.0). Temperature ranged from 20-55ºC. Kinetic parameters were determined with lactose (0.5–15 mM), 100 mM MES buffer of pH 6.5 (Lac161_ORF10) or pH 6.0 (Lac161_ORF7 and Lac36W_ORF11) and enzyme (0–5 mM). Reactions were incubated at 37°C (Lac161_ORF10), 42°C (Lac36W_ ORF11), and 50°C (Lac161_ORF3) for 30 min and terminated with Tris-HCl (pH 7) to a final concentration of 1 M. Aliquots of glucose oxidase/peroxidase reagent (125 μl) were added to each well and left to develop at 37°C for 30 min. Absorbance was measured at 450 nm and compared with a standard glucose curve to determine the amount of glucose released. All reactions were performed in triplicate.

Bioinformatic analysis

The cosmids were individually indexed and sequenced on Illumina HiSeq 2000 platform. Raw sequence data were assembled as described previously [45]. Open reading frames were annotated using MetaGeneMark [46]. Functions of proteins were predicted by BLAST analysis against NCBI non-redundant protein sequences, Pfam [47], and CAZy analysis toolkit [48]. Transmembrane helices were predicted by the TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM). Signal peptides were predicted using SignalP 4.0 [49]. Conserved protein domains were searched against NCBI Conserved Domain Database and analyzed with CDTree [50]. Protein structure was predicted with Phyre 2.0 [51]. Taxonomic affiliations of cosmid inserts were assigned based on compositional classifier PhyloPythiaS [52]. Complete sequences of metagenomic Lac+ cosmids have been deposited in GenBank with accession numbers KF255992-KF255994 and KF796593-KF796611.

Protein homology search against metagenomic datasets

SSEARCH36 [53] was used to search 158 metagenomes (32 aquatic, 76 human gut, 50 soil) for homologs to Lac161_ORF7, Lac161_ORF10, and Lac36W_ORF11, with an E-value threshold of 0.01. The database of metagenomes was compiled based on the set of aquatic and human gut metagenomes [54, 55], as well as a variety of soil metagenomes obtained from the MG-RAST server (http://metagenomics.anl.gov/). Accession numbers for all datasets are available in S4 Table. For comparison, and to estimate a background level of protein abundance using a housekeeping gene, all metagenomes were also searched for metagenomic homologs of the rpoB protein using HMMer (http://hmmer.org/) as implemented in MetAnnotate [56].

Results

Functional screening of β-galactosidases

Cosmid clones expressing β-galactosidase genes were screened in previously pooled metagenomic library 12AC [18]. Functional β-galactosidase enzymes can hydrolyse lactose (galactose-β-1,4-glucose) to galactose and glucose, facilitating the growth of bacterial hosts (lac) on M9 minimal media when lactose is used as the sole carbon source [18]. Because both the library host E. coli HB101 (lacY1) and surrogate S. meliloti RmF728 (lac) are resistant to streptomycin, which would affect selection of transconjugants in S. meliloti, 12AC cosmids were transferred from E. coli HB101 to DH5α (lacZYA) via electroporation. We obtained ~8.2 × 105 recombinant clones (TcR) of E. coli DH5α, which was ~10 fold greater than the number of original cosmid clones in the library. A total of 161 blue colonies were recovered on the selection/screen medium following spreading the E. coli DH5α clones on M9-lactose plate with X-gal. EcoR1-HindIII-BamHI restriction enzyme digest demonstrated that these 161 clones represented 17 different banding patterns.

We used S. meliloti from the α-Proteobacteria as a soil-dwelling surrogate host for screening in an effort to expand the range of recovered β-galactosidase-encoding clones. S. meliloti strain RmF728 is a derivative of the well studied Rm1021, and has been modified to carry a genomic deletion that removes the lactose metabolism genes [38]. The 12AC cosmids were transferred from E. coli DH5α to S. meliloti RmF728 via en masse triparental conjugation, and 1052 Lac+ colonies that were recovered on M9-lactose medium demonstrated reliable growth after streak purification. The colony color of these clones on M9-lactose containing X-gal ranged from white to varying shades of blue, suggesting that β-galactosidases in those white clones could hydrolyze lactose but not X-Gal. The measurement of β-galactosidase activities of 39 randomly selected S. meliloti clones grown in M9 lactose medium confirmed that the ability to grow on lactose as sole carbon source was due to cosmid clone-encoded β-galactosidase activity (S2 Table).

Each of the Lac+ clones was transferred from S. meliloti by triparental conjugation to E. coli DH5α (RifR). Electrophoretic comparison of 291 randomly chosen cosmids digested with EcoRI-HindIII-BamHI demonstrated 208 distinct patterns (71.5%), which suggested that the use of S. meliloti as a surrogate host for this screen yielded a greater diversity of β-galactosidase genes than when E. coli was used. There was some overlap with the clones isolated by complementation of E. coli DH5α, with four restriction enzyme digestion patterns common to both screens. In general, the clones demonstrating a Lac+ phenotype in both E. coli and S. meliloti exhibited higher activity (S2 Table).

Sequencing and annotation of Lac+ cosmids

We chose 3 distinct E. coli and 22 distinct S. meliloti Lac+ cosmids randomly for high-throughput sequencing (Table 2). Complete insert sequences were successfully assembled from 22 of 25 clones and these sequences and annotated ORFs of the other 22 cosmids were deposited in GenBank (Table 2). The Lac100B, Lac112W, and Lac224 clones from S. meliloti assembled only partially. Cloned metagenomic inserts were predicted to originate from multiple bacterial phyla (Cytophaga, Thermomicrobia, Verrucomicrobia, α-, β-, γ- and δ-Proteobacteria) and were GC rich overall (53% to 71%, 64% average).

[Figure omitted. See PDF.]

Table 2. 12AC metagenomic clones complementing E. coli DH5α (lacZYA) and S. meliloti RmF728 (lacEFGZ1K1) grown in M9-lactose medium.

https://doi.org/10.1371/journal.pone.0172545.t002

The E. coli Lac+ clones LacEc1, LacEc104, and LacEc123, and S. meliloti clones Lac24B, Lac36B, and Lac35B were predicted to originate from Serratia of the γ-Proteobacteria (Table 2). These clones overlapped over a segment of 15,344 bases (Fig A in S1 File). The 5’ region (positions 3–3,458) exhibited 93% identity to a chromosomal region (positions 2,604,251–2,607,707) of Serratia marcescens subsp. marcescens Db11 chromosome (GenBank accession HG326223), but the 3’ region (positions 6,652–15,344) of cloned DNA matched best to another region (positions 2,6143,370–2,623,056, 93% identity) of strain Db11. Eleven of 13 ORFs predicted in the overlapping region were 89–98% identical to the clustered orthologs (Fig B in S1 File). The second ORF (Lac35B, GenBank AGW45499) encodes a β-galactosidase (EC 3.2.1.23) with conserved domains of GH2 (Fig 1; Fig B in S1 File). The enzyme matched to the predicted β-galactosidase (SMDB11_2462) of S. marcescens subsp. marcescens Db11 with 98% amino acid sequence identity. Additionally, the annotated β-galactosidase also shares 66% amino acid sequence identity to the well characterized β-galactosidase LacZ (GenBank, BAE76126) of E. coli K12 substr. W3110. The amino acid residues important for catalytic function in E. coli LacZ [57] are conserved in the annotated β-galactosidase at Glu415, His417, Glu460, Tyr502, and Glu536 (S2 File).

Expression of the gene encoding the GH2 β-galactosidase from the cosmid clones in both E. coli and S. meliloti suggested a functional promoter(s) upstream of the gene. There were two regions homologous to the conserved -35 and -10 sites of RpoD promoters of E. coli [58] and S. meliloti [59] (Fig C in S1 File) in the 102-base intergenic region between the 3’ end (position 23) of an ORF (Lac35B_ORF9, GenBank AGW45500) encoding a two-component response regulator CitA, probably involved in Mg-citrate transport, and the β-galactosidase gene (Lac35B_ORF10, GenBank AGW45499) in Lac35B. These two putative promoters could drive expression of the β-galactosidase gene in E. coli and S. meliloti. Unlike the E. coli lac operon, there was no LacI homolog predicted in the cloned metagenomic DNA of those six cosmids. In addition, expression of the gene encoding β-galactosidase was neither dramatically inhibited by 15 mM glucose nor stimulated by addition of 0.4 mM IPTG in M9 medium, as determined on X-gal plates (S3 File).

Because the lactose permease LacY in E. coli DH5α [60] and ABC-type transporter LacEFGK1 for lactose in S. meliloti RmF728 are deleted [38], complementation would require a lactose transporter be encoded within the overlapping region (Fig B in S1 File). We detected an ABC-type transporter system consisting of periplasmic solute-binding protein, permease, and ATP-binding protein (ORF19-ORF17; GenBank AGW45491-AGW45493), but the transporter is probably involved in metal ion uptake. However, ORF21 (GenBank AGW45496) is predicted to encode a major facilitator transporter (IPR020846) with 14 transmembrane helices. This protein belongs to the same major facilitator superfamily as E. coli lactose permease (LacY) and may be functional as a lactose transporter when expressed in E. coli DH5α and S. meliloti RmF728.

The Lac+ clones Lac20, Lac71, and Lac172 isolated in S. meliloti shared a region of 14,707 bp (Fig A and B in S4 File) that was 93% identical to a segment (positions 2,578,724–2,593,427) of the S. marcescens WW4 chromosome (GenBank accession CP003959). The 14 annotated ORFs within this region exhibited 85–100% amino acid sequence identities to the clustered orthologs (Fig B in S4 File). These data suggest that the cloned DNA in Lac20, Lac71, and Lac172 originated from γ-Proteobacteria. One of the two major facilitator transporters (Lac20_ORF31, GenBank accession AHN97675; Lac20_ORF33, GenBank accession AHN97677) might be involved in lactose uptake in S. meliloti. We were unable to identify an ORF encoding a predicted β-galactosidase based on protein sequence homology.

Examination of the annotated ORFs of the remaining 13 Lac+ cosmids from S. meliloti (Table 2) also did not suggest any candidate that resembled known β-galactosidases. Based on a protein sequence comparison to the CAZy database, which showed low level similarity to proteins carrying known GHs, but not β-galactosidases, we chose Lac36W_ORF11 (GenBank accession AGW45517), Lac161_ORF7 (GenBank accession AGW45552), and Lac161_ORF10 (GenBank accession AGW45555) for biochemical characterization of putative β-galactosidase activity. We amplified the selected ORFs with PCR and cloned the amplicons into expression vectors, generating C-terminal His tags for overexpression in E. coli and subsequent affinity purification. Following processing, the resulting affinity-purified proteins were assayed for β-galactosidase activity.

Biochemical characterization of Lac36W_ORF11

Cosmid Lac36W exhibited β-galactosidase activity in S. meliloti (S2 Table), but not in E. coli. The cosmid contained a metagenomic DNA fragment of 34,259 bp with 67.3% GC (GenBank accession KF255993). The cloned DNA was assigned taxonomically to Xanthomonas of the ɣ-Proteobacteria (Table 2).

Protein sequence searches of the predicted 33 ORFs against the CAZy database suggested that Lac36W_ORF11 (GenBank accession AGW45517) showed sequence similarity to the protein ERE_21070 of Eubacterium rectale M104/1 (GenBank accession CBK94002), which has a Glyco_hydro_53 domain (Genbank accession CBK94002). Lac36W_ORF11 had a signal peptide of 21 amino acids predicted by SignalP 4.0 [49] and two domains (Fig 2): a 7 transmembrane region of receptor with diverse intracellular signaling module (7TMR-DISM-7TM; PF07695), and diguanylate cyclase domain (DGC or GGDEF, PF00990). We predict that the 7TMR-DISM-7TM domain might function as a lactose binding domain like other 7TM-containing proteins [61] and the C-terminal region may act as a β-galactosidase, although it exhibited no similarity to the endo-β-1,4-galactanase domain in the protein ERE_21070 and other known GH family members. Therefore, we cloned the entire Lac36W_ORF11 and expressed it in E. coli. Purified Lac36W_ORF11 protein was able to hydrolyze lactose to galactose and glucose (Table 3). The enzyme maintained 60% activity in the pH range of 6.5–8.0 (Fig 3E) and still kept ~50% activity at 50°C (Fig 3F). Because there was no similarity to any known GH domain and carbohydrate binding module (CBM), we proposed that Lac36W_ORF11 (GenBank, AGW45517) is a novel β-galactosidase.

[Figure omitted. See PDF.]

Fig 2. Conserved domains (E < 0.01) in β-galactosidases isolated from 12AC metagenomic library clones.

https://doi.org/10.1371/journal.pone.0172545.g002

[Figure omitted. See PDF.]

Fig 3. Biochemical characterization of novel β-galactosidases.

pH profiles of Lac161_ORF10 (A), Lac161_ORF7 (C), Lac36W_ORF11 (E). Temperature profiles of Lac161_ORF10 (B), Lac161_ORF7 (D), Lac36W_ORF11 (F).

https://doi.org/10.1371/journal.pone.0172545.g003

[Figure omitted. See PDF.]

Table 3. Biochemical characterization of novel β-galactosidases from 12AC Lac+ metagenomic clones complementing S. meliloti RmF728 (lac).

The β-galactosidase activity of purified proteins was assayed using lactose as substrate. Released glucose was quantified with a glucose oxidase activity kit.

https://doi.org/10.1371/journal.pone.0172545.t003

The Lac36W_ORF11 was situated within a putative operon, flanked by Lac36W_ORF12, immediately downstream, and Lac36W_ORF10, immediately upstream. Lac36W_ORF12 encodes a putative methionine-S-sulfoxide reductase and is located 107-bp downstream of the Lac36W_ORF11 (Fig A in S5 File). Lac36W_ORF10 was located 5 bp upstream of Lac36W_ORF11 and encodes a hypothetical protein (DUF2007). The nature of the promoter for this predicted operon and its basis for function in S. meliloti, but not E. coli, is unclear. The Lac36W_ORF14 is predicted to encode a transcriptional regulator (LysR-like) that may play a role in regulation of the operon. There were no ORFs encoding homologs to known transporters in the cloned 34-kb DNA fragment and, as a result, the genetic basis for uptake of lactose is unknown.

Biochemical characterization of Lac161_ORF7

Cosmid Lac161 complemented the Lac- phenotype of S. meliloti RmF728 but could not complement E. coli DH5α. The cosmid carried an insert of 35,906 bp with 59.1% GC content (GenBank accession KF255994). The metagenomic DNA was assigned taxonomically to Chthoniobacter of the phylum Verrucomicrobia (Table 2).

Among the 29 annotated ORFs, Lac161_ORF7 (GenBank accession AGW45552) was predicted to be a membrane-bound dehydrogenase protein with three domains (Fig 2): a putative membrane-bound dehydrogenase domain Piru-Ver-Nterm,

TIGR02604), Heat repeat 2 (PF13646), and a cytochrom_C (PF00034) with a putative heme-binding motif CxxCH (TIGER02603). The Heat_2 and Cytochrom_C domains may be involved in intracellular transport and electron transfer. In addition, Lac161_ORF7 was homologous to several proteins that were annotated as probable GHs, such as HVO_B0215 (GenBank accession ADE01485.1; CBM16, CAZy) of Haloferax volcanii DS2. Further sequence alignment analysis did not show any similarity to known GHs or CBMs. To determine whether the gene product exhibited any GH activity, the Lac161_ORF7 was cloned and expressed. Purified ORF7 protein was able to hydrolyze lactose with a Km of 1.8 mM, which is the lowest of the three β-galactosidases studied in this work (Table 3). In addition, the Km value of Lac161_ORF7 is similar to the reported Km (2.0) of E. coli LacZ [62]. The ORF7 protein was most active at the same pH of 6.0 as Lac36W_ORF11 (Table 3; Fig 3C and 3E), but the highest activity of Lac161_ORF7 was observed at 50°C (Fig 3D). In addition, Lac161_ORF7 had the highest kcat/Km among the β-galactosidases identified in this study. These results implied that Lac161_ORF7 (GenBank accession AGW45552) represents a novel β-galactosidase family.

Biochemical characterization of Lac161_ORF10

Protein sequence comparison with the Pfam database suggested that Lac161_ORF10 (GenBank accession AGW45555) grouped with a family of proteins of unknown function DUF377 (PF04041; Fig 1), some of which have been predicted to be β-fructosidases (GH32 and GH68), or α-L-arabinase and β-xylosidase (GH43 and GH62) [63]. Because of this observation, the Lac161_ORF10 was overexpressed and purified. The resulting gene product was able to hydrolyze lactose with a Km of 3.2 mM, similar to the values of Lac36W_ORF11 and Lac161_ORF7 (Table 3). The optimal pH and temperature of β-galactosidase activity was 6.5 and 37°C, respectively (Fig 3A and 3B). In order to further investigate the range of substrate specificity, four other disaccharides were tested as substrates. When sucrose (glucose-β-1,2-fructose) was added, no glucose was released, suggesting that Lac161_ORF10 was not a β-fructofuranosidase (or invertase, GH32). Additionally, the ORF10 protein was unable to catalyze hydrolysis of xyloside (xylose-β-1,4-xylose, often associated with GH43), maltose (glucose-α-1,4-glucoside, often associated with GH64), and cellobiose (glucose-β-1,4 glucoside, often associated with GH1). Sequence analysis and activity assays therefore suggested that Lac161_ORF10 (Genbank accession AGW45555) is also a new β-galactosidase, like Lac36W_ORF11 (GenBank accession AGW45517) and Lac161_ORF7 (GenBank accession AGW45552) proteins.

Lac161_ORF7 and ORF_10 encoding the two novel β-galactosidases might form one operon along with Lac161_ORF8 and Lac161_ORF9 (Fig B in S4 File). The Lac161_ORF8 encodes a hypothetical protein (GenBank accession AGW45553) homologous to an enolase superfamily including o-succinylbenzoate synthase (cd03320). Lac161_ORF9 encodes a hypothetical protein (GenBank accession AGW45554) with a similar domain to methane oxygenase (PF14100). The reason for gene expression in S. meliloti but not E. coli is not yet known.

Bioinformatic and structural modeling of Lac161_ORF10

A search of Lac161_ORF10 (GenBank accession AGW45555) against the NCBI Conserved Domain Database (CDD; [50] revealed no significant hits to characterized protein domains (E < 0.01). However, a glycoside hydrolase superfamily domain (GH43_62_32_68 superfamily, cl14647) was detected over region 130–234 as the top-scoring CDD hit overall (E = 0.10). More specifically, the match corresponds to a GH_J clan domain, which includes GH32 and GH68 enzymes. The presence of a GH_J domain within Lac161_ORF10 is further supported by the domain architectures of related sequences. The top 10 homologs of Lac161_ORF10 detected by BLAST were mainly from Bacteroides (S3 Table), and all possessed this domain over the aligning region (E < 0.01). According to CDTree, Lac161_ORF10 represented a highly distinct branch within the GH_J sequence cluster (Fig 4A), which provided some explanation for the observed weak similarity to existing CDD domains.

[Figure omitted. See PDF.]

Fig 4. Bioinformatic characterization of a putative glycosyl hydrolase domain in Lac161_ORF10.

(A) The NCBI Conserved Domain Database (CDD) predicts Lac161_ORF10 as a divergent member of the GH_J clan of glycosyl hydrolases. (B) Structural model of Lac161_ORF10 generated by Phyre 2.0, with a predicted cluster of 8 ligand-binding residues highlighted in yellow. The putative binding site was predicted by 3dLigandSite based on the Phyre model with PDB ID 1vkd (chain A) as the template. A NAG ligand is shown in red, which approximates the location of a lactose molecule in Lac161_ORF10. (C) An alignment of Lac161_ORF10 with the most similar members of the CDD’s GH_J sequence cluster (Genbank accession numbers are included on the right of the tree). The most conserved columns are coloured light blue. A predicted active site feature (D197) is highlighted in yellow, and is consistent with 3dLigandSite’s predicted cluster of ligand-binding residues.

https://doi.org/10.1371/journal.pone.0172545.g004

Proteins within the GH_J superfamily, including GH32 and GH68, all posses a five-bladed propeller fold, and share a funnel-shaped active site typically composed of a catalytic nucleophile (e.g., Asp) and proton donor (e.g., Glu) acting as the general acid/base as well as a RDP motif [64] involved in stabilizing the transition state (Fig 4B). Our analysis suggests that Lac161_ORF10 also shared some of these characteristics.

Using Phyre 2.0 [51], a structural model of Lac161_ORF10 was generated. Phyre predicted a five-bladed propeller fold for Lac161_ORF10 (Fig 4B) with high confidence (99.9%) based on the template PDB 1vkd_A, a predicted glycoside hydrolase from Thermotoga maritima (Tmari_1232). Both Tmari_1232 and Lac161_ORF10 are members of the Pfam DUF377 family, further supporting the model. We then analyzed potential active sites using a sequence and structure-based approach. According to the CDD sequence alignment, Lac161_ORF10 possesses a KDP motif (residues 196–198) that aligns to the active site RDP motif in the reference 1y9m structure (Fig 4C). Ligand-binding sites were also predicted in the structural model using 3dLigandSite [65]. This revealed a predicted cluster of eight residues, including the previously identified D-197 residue, as forming the putative active site (Fig 4B). However, alternate alignments and putative active sites from those reported above are possible given the structural repetition of five-bladed propellers. Ultimately, Lac161_ORF10 (GenBank, AGW45555) appears to represent a novel family of β-galactosidase with a GH_J-like five-bladed propeller glycoside hydrolase domain, and an active site similar in composition to other members of this superfamily.

Metagenome abundance

In order to investigate the distribution of sequences similar to the newly described β-galactosidase sequences throughout different metagenomes, we performed protein homology searches with these sequences against collections of aquatic, human gut, and soil metagenomic databases (Fig 5). The results were normalized against rpoB gene abundance. Homologs to each of the three genes were represented in all three habitats. However, Lac36W_ORF11 in the human gut was by far of greatest relative abundance. Lac36W_ORF11 was also high in soil, but not as high as in human gut. Although of overall lower relative abundance, Lac161_ORF10 was also of greater abundance in human gut than in soil or aquatic metagenomes. Lac161_ORF7 exhibited a distinct profile, being extremely rare in the human gut, at low levels in aquatic samples, and at higher levels in soil metagenomes. It will be of interest to determine whether these homologs are also functional β-galactosidases.

[Figure omitted. See PDF.]

Fig 5. Protein homology searches of novel β-galactosidase sequences of Lac161_ORF10, Lac161_ORF7 and Lac36W_ORF11 against aquatic, human gut, and soil metagenomic databases, normalized to the rpoB gene.

https://doi.org/10.1371/journal.pone.0172545.g005

Discussion

Metagenomics provides unprecedented access to the genomic potential of uncultivated microbial communities. Despite enormous progress resulting from developments in high throughput sequencing, the potential for novel enzyme discovery remains highest using a functional metagenomics approach, in which genes are isolated based on their function rather than by DNA sequence similarity to already known genes [9]. Using such an approach, we have discovered genes encoding novel types of lactose hydrolyzing enzymes. The enzymes encoded by these genes were biochemically similar to known enzymes, although they would not have been easily predicted by their sequences without knowing that they were carried on a segment of DNA that encoded β-galactosidase activity. These results demonstrate the importance of sequence-agnostic functional screens for the discovery of enzymes of novel origin, and suggest that further implementation of this strategy will contribute to fundamental knowledge about the relationship between sequence and protein function, improve the resolution of sequence based metagenomics, and expand the repertoire of novel enzymes available for industrial applications.

This work follows on other metagenomic functional screening efforts that have discovered β-galactosidases of GH1 [29, 31], GH2 [31–33], GH3 [66] GH42 [26, 27, 31–33], GH43 [25, 30, 66], GT4 [35], glycosidase [34] and two new GH members [7]. Some of the Lac+ clones recovered from S. meliloti can hydrolyze lactose but not X-Gal, suggestiong that X-Gal was not a good fit with the active sites of those enzymes. The same property of other metagenomic β-galacctosidases has been reported [28]. Here we have highlighted the application of functional metagenomics for mining novel enzymes from soil microbial communities. Although the functional metagenomics strategy has potential for expanding the availability of enzymes that can be further developed for biotech applications, it is perhaps just as important to apply such strategies to the expansion of knowledge that will inform functional interpretation of DNA sequence. This in turn could impact on the ability to derive metabolic information from genome sequence, even from uncultivated microorganisms.

We suggest that the use of a diversity of surrogate hosts for functional metagenomic screening has the potential to substantially extend the breadth of gene discovery. By discovering founding members of three novel β-galactosidase families, we have reinforced the value of functional metagenomics for isolating novel genes that could not have been predicted from DNA sequence analysis alone. Activity-based screening of metagenomic library clones for biocatalysts is dependent on the expression of genes of interest and presence of accessory components required for the enzyme activity in the surrogate hosts [5, 11]. Multi-host-systems have been developed to improve functional screening [12–14, 16, 17, 19, 67]. In the present work, functional screening of the corn field soil library (12AC) for the ability to complement β-galactosidase mutants resulted in a greater number of distinct clones using S. meliloti than the more widely used E. coli. In addition, three novel β-galactosidase genes were identified only in S. meliloti. These data emphasize the importance of developing of multi-host systems for functional screening.

Supporting information

[Figure omitted. See PDF.]

S1 File. Lac+ clones isolated from E. coli (LacEc1, LacEc104 and LacEc123) and S. meliloti (Lac24B, Lac35B and Lac36B).

(A) An overlapping region of 15,344 bp was present in those cosmids. (B) A β-galactosidase of family GH2 (ORF10, solid box), and putative lactose transporter (ORF21, dash lined box) were predicted in Lac35B. The regions encoding orthologs in γ-Proteobacteria Serratia marcescens subsp. marcescens Db11 chromosome (GenBank HG326223; 2,623,056–2,604,251 nt) were highlighted. (C) Putative RpoD promoters (P) active in both E. coli and S. meliloti were located upstream of the β-galactosidase gene. The same enzyme was encoded by LacEc1_ORF31, LacEc104_ORF20, LacEc123_ORF13, Lac24B_ORF9, and Lac36B_ORF3 respectively.

https://doi.org/10.1371/journal.pone.0172545.s001

(PDF)

S2 File. Sequence alignment of β-glactosidases LacZ of E. coli K12 substr, W3110 (Genbank, BAE76126) and 12AC metagenomic clone LacEc1 (β-Gal; LacEc1_ORF31; GenBank, KF96609).

Conserved amino acids Glu415, His417, Glu460, Tyr502 and Glu536 at the active sites of LacEc1_ORF31 were highlighted.

https://doi.org/10.1371/journal.pone.0172545.s002

(PDF)

S3 File. β-Gaclactosidase activity of Lac+ clones isolated from E. coli DH5α.

Empty cosmid pJC8 was used as a negative control. X-Gal was used as chromogenic substrate. (A) M9-glucose (15 mM), (B) M9-lactose (15 mM), (C) M9-glycerol (30 mM) + 0.4 mM IPTG.

https://doi.org/10.1371/journal.pone.0172545.s003

(PDF)

S4 File. Lac+ clone Lac20, Lac71 and Lac172 isolated from S. meliloti.

(A) An overlapping region of 14,707 bp was present in those cosmids. (B) The major facilitator transporter(s) (solid box) in the region might be involved in lactose uptake. The hypothetical protein(s) (dash lined box) might be a β-galactosidase. Orthologs in γ-Proteobacteria Serratia marcescens WW4 chromosome (GenBank CP003959; 2,578,724–2,593,247 nt) were highlighted.

https://doi.org/10.1371/journal.pone.0172545.s004

(PDF)

S5 File. A DNA fragment carrying genes encoding β-galactosidases in Lac+ cosmids Lac36W and Lac161.

(A) A gene locus from cosmid Lac36W (GenBank, KF255993). Lac36W_07, cytosine/adenosine deaminase; Lac36W_08, hypothetical protein; Lac36W_09, glutaminyl-tRNA synthetase; Lac36W_10, hypothetical protein; Lac36W_11, β-galactosidase; Lac36W_12, methionine-S-sulfoxide reductase; Lac36W_13, hypothetical protein; Lac36W_14, LysR family transcriptional regulator. The locations of potential promoter regions (P) were showed. (B) A gene locus from cosmid Lac161 (GenBank, KF255994). Lac161_06, histidine kinase; Lac161_07, β-galactosidase; Lac161_08, hypothetical protein; Lac161_09, hypothetical protein; Lac161_10, β-galactosidase; Lac161_11, hypothetical protein; Lac161_12, host specificity protein. The positions of potential promoter regions (P) are shown.

https://doi.org/10.1371/journal.pone.0172545.s005

(PDF)

S1 Table. DNA oligonucleotides used in this study with restriction recognition sites underlined.

https://doi.org/10.1371/journal.pone.0172545.s006

(PDF)

S2 Table. β-Galactosidase activities of random 12AC lac+ clones in S. meliloti RmF728.

https://doi.org/10.1371/journal.pone.0172545.s007

(PDF)

S3 Table. Top twenty homologs of Lac161_ORF10 detected by a BlastP search of the NCBI nr database.

https://doi.org/10.1371/journal.pone.0172545.s008

(PDF)

S4 Table. Detected abundance of three novel beta-galactosidases in a variety of metagenomic datasets.

https://doi.org/10.1371/journal.pone.0172545.s009

(PDF)

Acknowledgments

We are grateful to Julia Hanchard and Shirley Wong for technical assistance. This work was financially supported by a Strategic Projects grant and Discovery Grants from the Natural Sciences and Engineering Research Council of Canada (NSERC).

Author Contributions

1. Conceptualization: JDN DRR TCC.

2. Formal analysis: JC ACD.

3. Funding acquisition: JDN TCC.

4. Investigation: JC TR KE.

5. Methodology: JC ACD JDN TCC.

6. Project administration: KE.

7. Supervision: DRR JDN TCC.

8. Visualization: JC ACD TCC.

9. Writing – original draft: JC.

10. Writing – review & editing: JC JDN TCC.

Citation: Cheng J, Romantsov T, Engel K, Doxey AC, Rose DR, Neufeld JD, et al. (2017) Functional metagenomics reveals novel β-galactosidases not predictable from gene sequences. PLoS ONE 12(3): e0172545. https://doi.org/10.1371/journal.pone.0172545

References

1. Curtis TP, Sloan WT, Scannell JW. Estimating prokaryotic diversity and its limits. Proc Natl Acad Sci USA. 2002;99(16):10494–9. pmid:12097644

2. Simon C, Daniel R. Metagenomic analyses: past and future trends. Appl Environ Microbiol. 2011;77(4):1153–61. PubMed Central PMCID: PMCPMC3067235. pmid:21169428

3. Shade A, Hogan CS, Klimowicz AK, Linske M, McManus PS, Handelsman J. Culturing captures members of the soil rare biosphere. Environ Microbiol. 2012;14(9):2247–52. pmid:22788977

4. André I, Potocki-Veronese G, Barbe S, Moulis C, Remaud-Simeon M. CAZyme discovery and design for sweet dreams. Curr Opin Chem Biol. 2014;19:17–24. pmid:24780275

5. Taupp M, Mewis K, Hallam SJ. The art and design of functional metagenomic screens. Curr Opin Biotechnol. 2011;22(3):465–72. pmid:21440432

6. Ferrer M, Beloqui A, Timmis KN, Golyshin PN. Metagenomics for mining new genetic resources of microbial communities. J Mol Microbiol Biotechnol. 2009;16(1–2):109–23. pmid:18957866

7. Beloqui A, Nechitaylo TY, Lopez-Cortes N, Ghazi A, Guazzaroni ME, Polaina J, et al. Diversity of glycosyl hydrolases from cellulose-depleting communities enriched from casts of two earthworm species. Appl Environ Microbiol. 2010;76(17):5934–46. pmid:20622123

8. Iqbal HA, Feng Z, Brady SF. Biocatalysts and small molecule products from metagenomic studies. Curr Opin Chem Biol. 2012;16(1–2):109–16. PubMed Central PMCID: PMCPMC3335775. pmid:22455793

9. Ufarté L, Potocki-Veronese G, Laville É. Discovery of new protein families and functions: new challenges in functional metagenomics for biotechnologies and microbial ecology. Frontiers in microbiology. 2015;6:563. PubMed Central PMCID: PMCPMC4456863. pmid:26097471

10. Gabor EM, Alkema WBL, Janssen DB. Quantifying the accessibility of the metagenome by random expression cloning techniques. Environ Microbiol. 2004;6(9):879–86. pmid:15305913

11. Martinez A, Kolvek SJ, Yip CLT, Hopke J, Brown KA, Macneil IA, et al. Genetically modified bacterial strains and novel bacterial artificial chromosome shuttle vectors for constructing environmental libraries and detecting heterologous natural products in multiple expression hosts. Appl Environ Microbiol. 2004;70(4):2452–63. PubMed Central PMCID: PMCPMC383137. pmid:15066844

12. Wang C, Meek DJ, Panchal P, Boruvka N, Archibald FS, Driscoll BT, et al. Isolation of poly-3-hydroxybutyrate metabolism genes from complex microbial communities by phenotypic complementation of bacterial mutants. Appl Environ Microbiol. 2006;72(1):384–91. PubMed Central PMCID: PMCPMC1352230. pmid:16391068

13. Aakvik T, Degnes K, Dahlsrud R, Schmidt F, Dam R, Yu L, et al. A plasmid RK2-based broad-host-range cloning vector useful for transfer of metagenomic libraries to a variety of bacterial species. FEMS Microbiol Lett. 2009;296:149–58. pmid:19459950

14. Craig JW, Chang F-Y, Kim JH, Obiajulu SC, Brady SF. Expanding small-molecule functional metagenomics through parallel screening of broad-host-range cosmid environmental DNA libraries in diverse proteobacteria. Appl Environ Microbiol. 2010;76(5):1633–41. PubMed Central PMCID: PMCPMC2832356. pmid:20081001

15. Hao Y, Winans SC, Glick BR, Charles TC. Identification and characterization of new LuxR/LuxI-type quorum sensing systems from metagenomic libraries. Environ Microbiol. 2010;12(1):105–17. PubMed Central PMCID: PMCPMC2905786. pmid:19735279

16. Li Y, Wexler M, Richardson DJ, Bond PL, Johnston AW. Screening a wide host-range, waste-water metagenomic library in tryptophan auxotrophs of Rhizobium leguminosarum and of Escherichia coli reveals different classes of cloned trp genes. Environ Microbiol. 2005;7(12):1927–36. pmid:16309391

17. Ly MA, Liew EF, Le NB, Coleman NV. Construction and evaluation of pMycoFos, a fosmid shuttle vector for Mycobacterium spp. with inducible gene expression and copy number control. J Microbiol Methods. 2011;86(3):320–6. pmid:21689690

18. Cheng J, Pinnell L, Engel K, Neufeld JD, Charles TC. Versatile broad-host-range cosmids for construction of high quality metagenomic libraries. J Microbiol Methods. 2014;99:27–34. pmid:24495694

19. Biver S, Steels S, Portetelle D, Vandenbol M. Bacillus subtilis as a tool for screening soil metagenomic libraries for antimicrobial activities. J Microbiol Biotechnol. 2013;23(6):850–5. Epub 2013/05/17. pmid:23676924

20. Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B. The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res. 2009;37(Database issue):D233–8. PubMed Central PMCID: PMCPMC2686590. pmid:18838391

21. Harju M, Kallioinen H, Tossavainen O. Lactose hydrolysis and other conversions in dairy products: Technological aspects. Int Dairy J. 2012;22(2):104–9.

22. Guimaraes PMR, Teixeira JA, Domingues L. Fermentation of lactose to bio-ethanol by yeasts as part of integrated solutions for the valorisation of cheese whey. Biotechnol Adv. 2010;28(3):375–84. pmid:20153415

23. Marrakchi M, Dzyadevych SV, Lagarde F, Martelet C, Jaffrezic-Renault N. Conductometric biosensor based on glucose oxidase and beta-galactosidase for specific lactose determination in milk. Mat Sci Eng C-Bio S. 2008;28(5–6):872–5.

24. Gosling A, Stevens GW, Barber AR, Kentish SE, Gras SL. Recent advances refining galactooligosaccharide production from lactose. Food Chem. 2010;121(2):307–18.

25. Ferrer M, Ghazi A, Beloqui A, Vieites JM, López-Cortés N, Marín-Navarro J, et al. Functional metagenomics unveils a multifunctional glycosyl hydrolase from the family 43 catalysing the breakdown of plant polymers in the calf rumen. PLoS ONE. 2012;7(6):e38134. PubMed Central PMCID: PMCPMC3382598. pmid:22761666

26. Wang K, Li G, Yu SQ, Zhang CT, Liu YH. A novel metagenome-derived beta-galactosidase: gene cloning, overexpression, purification and characterization. Appl Microbiol Biotechnol. 2010;88(1):155–65. pmid:20614217

27. Zhang X, Li H, Li CJ, Ma T, Li G, Liu YH. Metagenomic approach for the isolation of a thermostable beta-galactosidase with high tolerance of galactose and glucose from soil samples of Turpan Basin. BMC Microbiol. 2013;13:237. pmid:24156692

28. Erich S, Kuschel B, Schwarz T, Ewert J, Bohmer N, Niehaus F, et al. Novel high-performance metagenome beta-galactosidases for lactose hydrolysis in the dairy industry. J Biotechnol. 2015;210:27–37. pmid:26122513

29. Gupta R, Govil T, Capalash N, Sharma P. Characterization of a glycoside hydrolase family 1 βGupta R, Govil from hot spring metagenome with transglycosylation activity. Applied Biochemistry and Biotechnology. 2012;168(6):1681–93. pmid:23015191

30. Wierzbicka-Wos A, Bartasun P, Cieslinski H, Kur J. Cloning and characterization of a novel cold-active glycoside hydrolase family 1 enzyme with beta-glucosidase, beta-fucosidase and beta-galactosidase activities. BMC Biotechnol. 2013;13:22. PubMed Central PMCID: PMC3605331. pmid:23497058

31. Vester JK, Glaring MA, Stougaard P. Discovery of novel enzymes with industrial potential from a cold and alkaline environment by a combination of functional metagenomics and culturing. Microb Cell Fact. 2014;13:72. PubMed Central PMCID: PMCPMC4035831. pmid:24886068

32. Maruthamuthu M, Jimenez DJ, Stevens P, van Elsas JD. A multi-substrate approach for functional metagenomics-based screening for (hemi)cellulases in two wheat straw-degrading microbial consortia unveils novel thermoalkaliphilic enzymes. BMC Genomics. 2016;17:86. PubMed Central PMCID: PMCPMC4730625. pmid:26822785

33. Li Y, Wang Z, Xu X, Jin L. A Ca-alginate particle co-immobilized with Phanerochaete chrysosporium cells and the combined cross-linked enzyme aggregates from Trametes versicolor. Bioresour Technol. 2015;198:464–9. pmid:26413897

34. Matsuzawa T, Yaoi K. Screening, identification, and characterization of a novel saccharide-stimulated beta-glycosidase from a soil metagenomic library. Appl Microbiol Biotechnol. 2016.

35. Wang SD, Guo GS, Li L, Cao LC, Tong L, Ren GH, et al. Identification and characterization of an unusual glycosyltransferase-like enzyme with beta-galactosidase activity from a soil metagenomic library. Enzyme Microb Technol. 2014;57:26–35. pmid:24629264

36. Finan TM, Hartweig E, Lemieux K, Bergman K, Walker GC, Signer ER. General transduction in Rhizobium meliloti. J Bacteriol. 1984;159(1):120–4. pmid:6330024

37. Meade HM, Long SR, Ruvkun GB, Brown SE, Ausubel FM. Physical and genetic characterization of symbiotic and auxotrophic mutants of Rhizobium meliloti induced by transposon Tn5 mutagenesis. J Bacteriol. 1982;149(1):114–22. Epub 1982/01/01. PubMed Central PMCID: PMCPmc216598. pmid:6274841

38. Charles TC, Finan TM. Analysis of a 1600-kilobase Rhizobium meliloti megaplasmid using defined deletions generated in vivo. Genetics. 1991;127(1):5–20. pmid:1849856

39. Hanahan D. Studies on transformation of Escherichia coli with plasmids. J Mol Biol. 1983;166(4):557–80. Epub 1983/06/05. pmid:6345791

40. Boyer HW, Roulland-Dussoix D. A complementation analysis of the restriction and modification of DNA in Escherichia coli. J Mol Biol. 1969;41(3):459–72. Epub 1969/05/14. pmid:4896022

41. Schäfer A, Tauch A, Jäger W, Kalinowski J, Thierbach G, Pühler A. Small mobilizable multi-purpose cloning vectors derived from the Escherichia coli plasmids pK18 and pK19: selection of defined deletions in the chromosome of Corynebacterium glutamicum. Gene. 1994;145(1):69–73. pmid:8045426

42. Khan SR, Gaines J, Roop RM, Farrand SK. Broad-host-range expression vectors with tightly regulated promoters and their use to examine the influence of TraR and TraM expression on Ti plasmid quorum sensing. Appl Environ Microbiol. 2008;74(16):5053–62. PubMed Central PMCID: PMCPMC2519271. pmid:18606801

43. Cheng J, Sibley CD, Zaheer R, Finan TM. A Sinorhizobium meliloti minE mutant has an altered morphology and exhibits defects in legume symbiosis. Microbiology. 2007;153(2):375–87.

44. Cowie A, Cheng J, Sibley CD, Fong Y, Zaheer R, Patten CL, et al. An integrated approach to functional genomics: construction of a novel reporter gene fusion library for Sinorhizobium meliloti. Appl Environ Microbiol. 2006;72(11):7156–67. pmid:16963549

45. Lam KN, Hall MW, Engel K, Vey G, Cheng J, Neufeld JD, et al. Evaluation of a pooled strategy for high-throughput sequencing of cosmid clones from metagenomic libraries. PLoS One. 2014;9(6):e98968. PubMed Central PMCID: PMC4049660. pmid:24911009

46. Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010;38(12):e132. PubMed Central PMCID: PMCPMC2896542. pmid:20403810

47. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(Database issue):D222–30. PubMed Central PMCID: PMC3965110. pmid:24288371

48. Park BH, Karpinets TV, Syed MH, Leuze MR, Uberbacher EC. CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. Glycobiology. 2010;20(12):1574–84. pmid:20696711

49. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nature methods. 2011;8(10):785–6. pmid:21959131

50. Marchler-Bauer A, Zheng C, Chitsaz F, Derbyshire MK, Geer LY, Geer RC, et al. CDD: conserved domains and protein three-dimensional structure. Nucleic Acids Res. 2013;41(Database issue):D348–52. PubMed Central PMCID: PMC3531192. pmid:23197659

51. Kelley LA, Sternberg MJ. Protein structure prediction on the Web: a case study using the Phyre server. Nature protocols. 2009;4(3):363–71. pmid:19247286

52. Patil KR, Roune L, McHardy AC. The PhyloPythiaS web werver for taxonomic assignment of metagenome sequences. PLoS ONE. 2012;7(6):e38581. pmid:22745671

53. Pearson WR, Lipman DJ. Improved tools for biological sequence comparison. Proc Natl Acad Sci. 1988;85(8):2444–8. PubMed Central PMCID: PMCPMC280013. pmid:3162770

54. Doxey AC, Kurtz DA, Lynch MDJ, Sauder LA, Neufeld JD. Aquatic metagenomes implicate Thaumarchaeota in global cobalamin production. ISME J. 2015;9(2):461–71. PubMed Central PMCID: PMCPMC4303638. pmid:25126756

55. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464(7285):59–65. PubMed Central PMCID: PMCPMC3779803. pmid:20203603

56. Petrenko P, Lobb B, Kurtz DA, Neufeld JD, Doxey AC. MetAnnotate: function-specific taxonomic profiling and comparison of metagenomes. BMC Biology. 2015;13(1):92. PubMed Central PMCID: PMCPMC4636000.

57. Jacobson RH, Zhang XJ, DuBose RF, Matthews BW. Three-dimensional structure of beta-galactosidase from E. coli. Nature. 1994;369(6483):761–6. pmid:8008071

58. Lisser S, Margalit H. Compilation of E. coli mRNA promoter sequences. Nucleic Acids Res. 1993;21(7):1507–16. Epub 1993/04/11. PubMed Central PMCID: PMCPmc309355. pmid:8479900

59. MacLellan SR, MacLean AM, Finan TM. Promoter prediction in the rhizobia. Microbiology. 2006;152(Pt 6):1751–63. Epub 2006/06/01. pmid:16735738

60. Meselson M, Yuan R. DNA restriction enzyme from E. coli. Nature. 1968;217(5134):1110–4. Epub 1968/03/23. pmid:4868368

61. Anantharaman V, Aravind L. Application of comparative genomics in the identification and analysis of novel families of membrane-associated receptors in bacteria. BMC Genomics. 2003;4(1):34. Epub 2003/08/14. PubMed Central PMCID: PMCPmc212514. pmid:12914674

62. Wallenfels K, Malhotra OP. Galactosidases. Adv Carbohydr Chem. 1961;16:239–98. pmid:14004681

63. Naumoff DG. Beta-fructosidase superfamily: homology with some alpha-L-arabinases and beta-D-xylosidases. Proteins. 2001;42(1):66–76. pmid:11093261

64. Lammens W, Le Roy K, Schroeven L, Van Laere A, Rabijns A, Van den Ende W. Structural insights into glycoside hydrolase family 32 and 68 enzymes: functional implications. J Exp Bot. 2009;60(3):727–40. Epub 2009/01/09. pmid:19129163

65. Wass MN, Kelley LA, Sternberg MJ. 3DLigandSite: predicting ligand-binding sites using similar structures. Nucleic Acids Res. 2010;38(Web Server issue):W469–73. Epub 2010/06/02. PubMed Central PMCID: PMCPmc2896164. pmid:20513649

66. Ramirez-Escudero M, Del Pozo MV, Marin-Navarro J, Gonzalez B, Golyshin PN, Polaina J, et al. Structural and functional characterization of a ruminal beta-glycosidase defines a novel subfamily of glycoside hydrolase family 3 with permuted domain topology. J Biol Chem. 2016;291(46):24200–14. PubMed Central PMCID: PMCPMC5104943. pmid:27679487

67. Wexler M, Johnston AW. Wide host-range cloning for functional metagenomics. Methods Mol Biol (Clifton, NJ). 2010;668:77–96. Epub 2010/09/11.