Stalking the fourth domain in metagenomic data: searching for, discovering, and interpreting novel, deep branches in marker gene phylogenetic trees

Background

Most of our knowledge about the ancient evolutionary history of organisms has been derived from data associated with specific known organisms (i.e., organisms that we can study directly such as plants, metazoans, and culturable microbes). Recently, however, a new source of data for such studies has arrived: DNA sequence data generated directly from environmental samples. Such metagenomic data has enormous potential in a variety of areas including, as we argue here, in studies of very early events in the evolution of gene families and of species.

Methodology/Principal Findings

We designed and implemented new methods for analyzing metagenomic data and used them to search the Global Ocean Sampling (GOS) Expedition data set for novel lineages in three gene families commonly used in phylogenetic studies of known and unknown organisms: small subunit rRNA and the recA and rpoB superfamilies. Though the methods available could not accurately identify very deeply branched ss-rRNAs (largely due to difficulties in making robust sequence alignments for novel rRNA fragments), our analysis revealed the existence of multiple novel branches in the recA and rpoB gene families. Analysis of available sequence data likely from the same genomes as these novel recA and rpoB homologs was then used to further characterize the possible organismal source of the novel sequences.

Conclusions/Significance

Of the novel recA and rpoB homologs identified in the metagenomic data, some likely come from uncharacterized viruses while others may represent ancient paralogs not yet seen in any cultured organism. A third possibility is that some come from novel cellular lineages that are only distantly related to any organisms for which sequence data is currently available. If there exist any major, but so-far-undiscovered, deeply branching lineages in the tree of life, we suggest that methods such as those described herein currently offer the best way to search for them.     

Microbiome Analysis of Stool Samples from African Americans with Colon Polyps

Background

Colonic polyps are common tumors occurring in ~50% of Western populations with ~10% risk of malignant progression. Dietary agents have been considered the primary environmental exposure to promote colorectal cancer (CRC) development. However, the colonic mucosa is permanently in contact with the microbiota and its metabolic products including toxins that also have the potential to trigger oncogenic transformation.

Aim

To analyze fecal DNA for microbiota composition and functional potential in African Americans with pre-neoplastic lesions.

Materials & Methods

We analyzed the bacterial composition of stool samples from 6 healthy individuals and 6 patients with colon polyps using 16S ribosomal RNA-based phylogenetic microarray; the Human intestinal Tract Chip (HITChip) and 16S rRNA gene barcoded 454 pyrosequencing. The functional potential was determined by sequence-based metagenomics using 454 pyrosequencing.

Results

Fecal microbiota profiling of samples from the healthy and polyp patients using both a phylogenetic microarraying (HITChip) and barcoded 454 pyrosequencing generated similar results. A distinction between both sets of samples was only obtained when the analysis was performed at the sub-genus level. Most of the species leading to the dissociation were from the Bacteroides group. The metagenomic analysis did not reveal major differences in bacterial gene prevalence/abundances between the two groups even when the analysis and comparisons were restricted to available Bacteroides genomes.

Conclusion

This study reveals that at the pre-neoplastic stages, there is a trend showing microbiota changes between healthy and colon polyp patients at the sub-genus level. These differences were not reflected at the genome/functions levels. Bacteria and associated functions within the Bacteroides group need to be further analyzed and dissected to pinpoint potential actors in the early colon oncogenic transformation in a large sample size.     

Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage

Bacteria in the 16S rRNA clade SAR86 are among the most abundant uncultivated constituents of microbial assemblages in the surface ocean for which little genomic information is currently available. Bioinformatic techniques were used to assemble two nearly complete genomes from marine metagenomes and single-cell sequencing provided two more partial genomes. Recruitment of metagenomic data shows that these SAR86 genomes substantially increase our knowledge of non-photosynthetic bacteria in the surface ocean. Phylogenomic analyses establish SAR86 as a basal and divergent lineage of γ-proteobacteria, and the individual genomes display a temperature-dependent distribution. Modestly sized at 1.25–1.7 Mbp, the SAR86 genomes lack several pathways for amino-acid and vitamin synthesis as well as sulfate reduction, trends commonly observed in other abundant marine microbes. SAR86 appears to be an aerobic chemoheterotroph with the potential for proteorhodopsin-based ATP generation, though the apparent lack of a retinal biosynthesis pathway may require it to scavenge exogenously-derived pigments to utilize proteorhodopsin. The genomes contain an expanded capacity for the degradation of lipids and carbohydrates acquired using a wealth of tonB-dependent outer membrane receptors. Like the abundant planktonic marine bacterial clade SAR11, SAR86 exhibits metabolic streamlining, but also a distinct carbon compound specialization, possibly avoiding competition.     

Longleaf Pine (Pinus palustris P. Mill.) Restoration Using Herbicides: Overstory and Understory Vegetation Responses on a Coastal Plain Flatwoods Site in Florida,U.S.A.

A study was conducted on a Coastal Plain flatwoods site in Florida to determine the effects of common forestry herbicides on Longleaf pine seedling survival and growth and on the understory vegetation. Following removal of the overstory slash pine, five low‐rate herbicide treatments were applied over the top of planted Longleaf pine seedlings to provide short‐term understory vegetation control and accelerate seedling growth. The objective was to increase Longleaf pine growth by reducing the shrub competition while increasing the herbaceous ground cover. Despite causing reduction in seedling survival over the control treatment, imazapyr (0.21 ae kg/ha) resulted in the highest seedling growth (height and volume). The significant reduction of shrub cover, density, and height by imazapyr was believed to be responsible for the improved seedling growth in this treatment. Both hexazinone (0.56 ai kg/ha) and sulfometuron methyl (0.26 ai kg/ha) + hexazinone (0.56 ai kg/ha) treatments also reduced cover of Runner oak, a major shrub species, but the response was evident only 8 months after treatment. Although sulfometuron methyl (0.26 ai kg/ha) and sulfometuron methyl + hexazinone treatments did not result in any significant change in overall grass, forb, and shrub cover, both treatments resulted in greater Longleaf pine growth compared to the control. None of the herbicides significantly affected the major understory grasses and forbs. Overall, imazapyr provided the best desired results with significant increase in seedling growth and better control of shrub species with no significant effects on grass and other herbaceous species cover.     

Development of a mechanism-based high-throughput screen assay for leucine-rich repeat kinase 2—Discovery of LRRK2 inhibitors

LRRK2 is a large and complex protein that possesses kinase and GTPase activities and has emerged as the most relevant player in PD pathogenesis possibly through a toxic gain-of-function mechanism. Kinase activity is a critical component of LRRK2 function and represents a viable target for drug discovery. We now report the development of a mechanism-based TR-FRET assay for the LRRK2 kinase activity using full-length LRRK2. In this assay, PLK-peptide was chosen as the phosphoryl acceptor. A combination of steady-state kinetic studies and computer simulations was used to calculate the initial concentrations of ATP and PLK-peptide to generate a steady-state situation that favors the identification of ATP noncompetitive inhibitors. The assay was also run in the absence of GTP. Under these conditions, the assay was sensitive to inhibitors that directly interact with the kinase domain and those that modulate the kinase activity by directly interacting with other domains including the GTPase domain. The assay was optimized and used to robustly evaluate our compound library in a 384-well format. An inhibitor identified through the screen was further characterized as a noncompetitive inhibitor with both ATP and PLK-peptide and showed similar inhibition against LRRK2 WT and the mutant G2019S.     

Competition for 15N-labeled fertilizer in a pecan (Carya illinoensis K. Koch)-cotton (Gossypium hirsutum L.) alley cropping system in the southern United States

The degree of tree-crop competition for nitrogen (N) and its effect on fertilizer-use efficiency and N movement were examined in a pecan (Carya illinoensis K. Koch)-cotton (Gossypium hirsutum L.) alley cropping system. Assessment of competition was accomplished via the installation of a belowground polyethylene root barrier in half the number of plots in order to provide two treatments–barrier and non-barrier. The percentage of N derived from fertilizer (NDF) and fertilizer-use efficiency (UFN) were determined using ¹⁵N-enriched ammonium sulfate (5% atom enrichment) applied at 89.6 kg N ha^–1. In cotton, the barrier treatment resulted in higher leaf (38%), stem (66%), seed cotton (55%) and total (58%) biomass compared to the non-barrier treatment. Total N content in leaf, stem and seed cotton was 67% higher in barrier compared to non-barrier treatment. Percentage of NDF in cotton leaf and stem was significantly lower in barrier (15.8% and 17.3%, respectively) compared to non-barrier treatment (20.4% and 21.2%, respectively). For UFN, this trend was reversed, with plants in barrier treatment having a higher percentage of UFN. Root trenching did not affect pecan foliar N concentration, canopy N content, NDF or UFN. In soil, N recovery at 90–120 cm depth was lower in non-barrier treatment, indicating tree root uptake of fertilizer N. Although tree roots in non-barrier treatment had access to fertilizer N, competition was mainly for N already in the soil, since fertilizer was applied after major seasonal nutrient demands of the trees had been met. Overall, the alley cropping system in this study exhibits potential for efficient N cycling, given the apparent ability of pecan trees to intercept and uptake N fertilizer from deeper soil layers and return to surface soil via litterfall.     

A note on efficient computation of haplotypes via perfect phylogeny.

The problem of inferring haplotype phase from a population of genotypes has received a lot of attention recently. This is partly due to the observation that there are many regions on human genomic DNA where genetic recombination is rare (Helmuth, 2001; Daly et al., 2001; Stephens et al., 2001; Friss et al., 2001). A Haplotype Map project has been announced by NIH to identify and characterize populations in terms of these haplotypes. Recently, Gusfield introduced the perfect phylogeny haplotyping problem, as an algorithmic implication of the no-recombination in long blocks observation, together with the standard population-genetic assumption of infinite sites. Gusfield's solution based on matroid theory was followed by direct theta(nm2) solutions that use simpler techniques (Bafna et al., 2003; Eskin et al., 2003), and also bound the number of solutions to the PPH problem. In this short note, we address two questions that were left open. First, can the algorithms of Bafna et al. (2003) and Eskin et al. (2003) be sped-up to O(nm + m2) time, which would imply an O(nm) time-bound for the PPH problem? Second, if there are multiple solutions, can we find one that is most parsimonious in terms of the number of distinct haplotypes. We give reductions that suggests that the answer to both questions is "no." For the first problem, we show that computing the output of the first step (in either method) is equivalent to Boolean matrix multiplication. Therefore, the best bound we can presently achieve is O(nm(omega-1)), where omega < or = 2.52 is the exponent of matrix multiplication. Thus, any linear time solution to the PPH problem likely requires a different approach. For the second problem of computing a PPH solution that minimizes the number of distinct haplotypes, we show that the problem is NP-hard using a reduction from Vertex Cover (Garey and Johnson, 1979).     

Soil nitrogen dynamics as an indicator for longleaf pine restoration

Assessing the status of soil nutrients with their corresponding microbial communities provides important information about degraded soils during the restoration of coastal wet pine forests. Net nitrogen mineralization, nitrogen‐oxidizing bacteria (NOB), and soil microbial biomass were compared with patch‐derived volume along a 110‐year longleaf pine (Pinus palustris Mill.) chronosequence for identifying a trajectory and ecological benchmark during forest restoration. Net nitrogen mineralization rates decreased significantly in the maturing‐aged, pine patches, driven by a larger drop in net nitrification. Net nitrification and abundance of NOB were higher in young pine patches compared to soils from the maturing (86–110 years) pine patches. Gross nitrate fluxes followed the nonfungal portion of the soil microbial biomass along the chronosequence, declining in 64‐year‐old pine patches. Microbial biomass peaked in patches 17–34 years of age, but significantly declined in the older patches. Fungal biomass leveled off without decline. Ammonium was the major source of nitrogen within the maturing pine patches as well as the wetland patches, indicating that ammonium maintains longleaf pine during growth‐limiting conditions. Nitrate dominated during rapid tree growth, optimally in mesic conditions. The relative amounts of available ammonium to nitrate can be used to model nitrogen cycling in facultative‐wetland pine forests of the coastal United States as soils alternate between wet and mesic conditions. A key restoration benchmark occurred after 86 years of pine development when pine patch growth rates slowed, with lower numbers of NOB, when the nonfungal biomass leveled off, and net nitrification rates are at a minimum, during pine maturation.     

The Sorcerer II Global Ocean Sampling Expedition: Expanding the Universe of Protein Families

Metagenomics projects based on shotgun sequencing of populations of micro-organisms yield insight into protein families. We used sequence similarity clustering to explore proteins with a comprehensive dataset consisting of sequences from available databases together with 6.12 million proteins predicted from an assembly of 7.7 million Global Ocean Sampling (GOS) sequences. The GOS dataset covers nearly all known prokaryotic protein families. A total of 3,995 medium- and large-sized clusters consisting of only GOS sequences are identified, out of which 1,700 have no detectable homology to known families. The GOS-only clusters contain a higher than expected proportion of sequences of viral origin, thus reflecting a poor sampling of viral diversity until now. Protein domain distributions in the GOS dataset and current protein databases show distinct biases. Several protein domains that were previously categorized as kingdom specific are shown to have GOS examples in other kingdoms. About 6,000 sequences (ORFans) from the literature that heretofore lacked similarity to known proteins have matches in the GOS data. The GOS dataset is also used to improve remote homology detection. Overall, besides nearly doubling the number of current proteins, the predicted GOS proteins also add a great deal of diversity to known protein families and shed light on their evolution. These observations are illustrated using several protein families, including phosphatases, proteases, ultraviolet-irradiation DNA damage repair enzymes, glutamine synthetase, and RuBisCO. The diversity added by GOS data has implications for choosing targets for experimental structure characterization as part of structural genomics efforts. Our analysis indicates that new families are being discovered at a rate that is linear or almost linear with the addition of new sequences, implying that we are still far from discovering all protein families in nature.