Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

Black Band Disease Microbial Community Variation on Corals in Three Regions of the Wider Caribbean

Black band disease (BBD) is a pathogenic consortium of microorganisms that primarily affects massive framework-building scleractinian corals on reefs worldwide. There has been considerable debate concerning the microbial community composition of BBD. The aim of this study was to utilize microbial profiling to assess overall patterns of variation in the BBD bacterial community with respect to geographic location, host coral species, time, and nutrient regime. Length heterogeneity polymerase chain reaction (LH-PCR) was employed to differentiate BBD communities based on the natural variation in the sequence lengths within hypervariable domains of the 16S rRNA gene. Analysis of LH-PCR profiles of 97 BBD samples using multivariate ordination methods and analysis of similarity revealed significant clustering with respect to geographic region when comparing BBD sampled from reefs near Lee Stocking Island in the Bahamas’ Exuma Chain, the Northern Florida Keys (NFK), and St. John in the US Virgin Islands. There was much variability in BBD community composition on a regional basis, between sites in the NFK, and in terms of coral host species. The observed differences among BBD microbial community profiles were driven primarily by variation in relative abundance of 313–316-bp amplicons, which correspond to cyanobacteria and α-proteobacteria. The results obtained in this study support previous reports of intrinsic variability and complexity of the BBD microbial community but also suggest that this variability has biogeographic patterns.     

Reduced variation in the yellow-achaete-scute region in natural populations of Drosophila melanogaster

Restriction map variation in 64 X chromosome lines extracted from three different populations of Drosophila melanogaster was investigated with seven six-nucleotide-recognizing restriction enzymes for a 106-kb region encompassing the yellow gene and the achaete-scute complex that is located in the region of reduced crossing over close to the telomere. Nine restriction site polymorphisms (out of 176 sites scored) and 19 length polymorphisms (15 insertions and 4 deletions) were detected. The estimated level of heterozygosity per nucleotide, H = 0.0003, is much lower than that reported for autosomal and sex-linked loci located in regions with normal levels of crossing over. The overall frequency of polymorphic restriction sites is reduced. Six out of nine restriction site polymorphisms are unique and the other three have frequencies less than 0.17. Some large insertions have reached relatively high frequencies, 0.08 to 0.17. Consistent with the theoretically predicted negative relationship between crossing over and the magnitude of linkage disequilibrium, an increase in the relative number of nonrandom associations was observed in the y-ac-sc region.     

Sequencing of pooled DNA samples (Pool-Seq) uncovers complex dynamics of transposable element insertions in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that parasitize genomes by semi-autonomously increasing their own copy number within the host genome. While TEs are important for genome evolution, appropriate methods for performing unbiased genome-wide surveys of TE variation in natural populations have been lacking. Here, we describe a novel and cost-effective approach for estimating population frequencies of TE insertions using paired-end Illumina reads from a pooled population sample. Importantly, the method treats insertions present in and absent from the reference genome identically, allowing unbiased TE population frequency estimates. We apply this method to data from a natural Drosophila melanogaster population from Portugal. Consistent with previous reports, we show that low recombining genomic regions harbor more TE insertions and maintain insertions at higher frequencies than do high recombining regions. We conservatively estimate that there are almost twice as many "novel" TE insertion sites as sites known from the reference sequence in our population sample (6,824 novel versus 3,639 reference sites, with on average a 31-fold coverage per insertion site). Different families of transposable elements show large differences in their insertion densities and population frequencies. Our analyses suggest that the history of TE activity significantly contributes to this pattern, with recently active families segregating at lower frequencies than those active in the more distant past. Finally, using our high-resolution TE abundance measurements, we identified 13 candidate positively selected TE insertions based on their high population frequencies and on low Tajima's D values in their neighborhoods.     

Comparison of morphological and molecular genetic quantification of relative abundance of arbuscular mycorrhizal fungi within roots

Nested PCR amplicons of ribosomal RNA genes have been used to identify individuals within assemblages of arbuscular mycorrhizal (AM) fungi in roots and to estimate their relative abundance. Microscopy has also been used to identify their relative abundance in roots, but only at low resolution, usually the genus level. We evaluated the robustness of using nested PCR amplicons of ribosomal RNA genes to estimate the relative abundance of undefined AM fungi in uniformly aged roots in comparison to visual estimates. The relative abundance of AM fungi was assessed as per cent root length colonised by morphotypes and relative sequence type abundance in clone libraries. Plants were grown in coastal soil to obtain assemblages of unknown AM fungi at two times (spring and autumn). Relative abundance of dominant genera of AM fungi in roots (Archaeospora and Glomus) based on an analysis of ribosomal RNA genes did not consistently correspond with relative abundance of morphotypes. This microscopic vs. molecular genetic comparison supports previous conclusions that there can be limitations in using nested PCR amplicons for quantifying the relative abundance of AM fungi in roots, with a sampling bias likely to be of significance. Both molecular genetic and morphological methods are used to estimate relative abundance of AM fungi as a precursor to understanding mycorrhizal function in field soils, but they are rarely verified using alternative approaches although this may be necessary.     

Genomic suppression subtractive hybridization as a tool to identify differences in mycorrhizal fungal genomes

Understanding the ecology of methanogens in natural and engineered environments is a prerequisite to predicting or managing methane emissions. In this study, a novel high-throughput fingerprint method was developed for determining methanogen diversity and relative abundance within environmental samples. The method described here, designated amplicon length heterogeneity PCR of the mcrA gene (LH-mcrA), is based on the natural length variation in the mcrA gene. The mcrA gene encodes the alpha-subunit of the methyl-coenzyme M reductase, which is involved in the terminal step of methane production by methanogens. The methanogenic communities from stored swine and dairy manures were distinct from each other. To validate the method, methanogenic communities in a plug flow-type bioreactor (PFBR) treating swine manure were characterized using LH-mcrA method and correlated to mcrA gene clone libraries. The diversity and relative abundance of the methanogenic groups were assessed. Methanobrevibacter, Methanosarcinaceae, Methanoculleus, Methanogenium, Methanocorpusculum and one unidentified group were assigned to particular LH-mcrA amplicons. Particular phylotypes related to Methanoculleus were predominant in the last compartment of the PFBR where the bulk of methane was produced. LH-mcrA method was found to be a reliable, fast and cost-effective alternative for diversity assessment of methanogenic communities in microbial systems.     

Using RNase sequence specificity to refine the identification of RNA-protein binding regions

Massively parallel pyrosequencing is a high-throughput technology that can sequence hundreds of thousands of DNA/RNA fragments in a single experiment. Combining it with immunoprecipitation-based biochemical assays, such as cross-linking immunoprecipitation (CLIP), provides a genome-wide method to detect the sites at which proteins bind DNA or RNA. In a CLIP-pyrosequencing experiment, the resolutions of the detected protein binding regions are partially determined by the length of the detected RNA fragments (CLIP amplicons) after trimming by RNase digestion. The lengths of these fragments usually range from 50-70 nucleotides. Many genomic regions are marked by multiple RNA fragments. In this paper, we report an empirical approach to refine the localization of protein binding regions by using the distribution pattern of the detected RNA fragments and the sequence specificity of RNase digestion. We present two regions to which multiple amplicons map as examples to demonstrate this approach.     

Application of denaturing high-performance liquid chromatography for rice variety identification and seed purity assessment

Recent DNA sequencing projects and the establishment of high-throughput assays have provided an abundance of sequence information and data on nucleotide polymorphisms in rice. Based on previously identified single-nucleotide polymorphisms (SNPs) and insertions/deletions, we employed denaturing high-performance liquid chromatography (DHPLC) to genotype rice varieties using a chromatographic pattern-based strategy. In this study, 12 amplicons harboring multiple and informative SNPs were screened. PCR products of the 12 amplicons from 47 rice varieties were analyzed by DHPLC and DNA sequencing. Each homozygous sample with a single peak pattern in the initial DHPLC analysis was mixed with zhenshan97 for a second DHPLC analysis. The 12 amplicons were found to be polymorphic across the hybrids, and mixed homozygous samples with 43 distinct DHPLC elution profiles detected. Sequence analysis confirmed that the distinct DHPLC patterns corresponded to different DNA sequences. A set of distinct characteristic profiles in six amplicons differentiated between all of the hybrids, inbred lines, and restorer lines and produced unique a fingerprint for these lines. In addition, we found that the DHPLC pattern of the hybrid was in accordance with the results obtained by DHPLC analysis of a mixed sample of the two parents. These results demonstrate that DHPLC can be efficiently applied for the rapid and automated identification of diverse rice varieties and could possibly be utilized for seed genetic purity testing on a high-throughput scale.     

Environmental factors influencing the distribution of rRNA from Verrucomicrobia in soil

The Verrucomicrobia constitute a newly discovered division of the Bacteria identified as a numerically abundant component of soil microbial communities in numerous sites around the world. The relative abundance of rRNA from Verrucomicrobia was investigated in the soil to examine the influence of specific environmental factors on the distribution of Verrucomicrobia and to better understand the distribution of this group in terrestrial ecosystems. The abundance of the verrucomicrobial rRNA was determined by using a novel oligonucleotide probe that is specific for verrucomicrobial 16S rRNA. The abundance of verrucomicrobial 16S rRNA in soil microbial communities was determined in relation to plant community composition and soil management history over a period of 2 years. Additional samples were analyzed to determine if verrucomicrobial rRNA relative abundance changes in relation to either soil depth or soil moisture content. The Verrucomicrobia composed 1.9+/-0.2% of the microbial community rRNA present in the 85 soil samples examined. The distribution of verrucomicrobial rRNA in the soil reveals that Verrucomicrobia are significantly affected by environmental characteristics that change in relation to time, soil history, and soil depth, and reveals that a statistically significant amount of the variation in verrucomicrobial rRNA abundance can be explained by changes in soil moisture content.     

Genome-Tagged Amplification (GTA): a PCR-based method to prepare sample-tagged amplicons from hundreds of individuals for next generation sequencing

Sampling the sequence of a relatively small fraction of the genome in large numbers of individuals is an important objective for population genetics and association genetics approaches. However, currently available ‘sequence capture’ methods either require expensive instrumentation or have problems dealing with high sample numbers and relatively small target sizes. We have developed Genome-Tagged Amplification (GTA) as a flexible PCR-based method for preparing pools of hundreds of amplicons from hundreds of samples for next generation sequencing. The method involves tagging of genomic DNA with barcode adapters at restriction sites, followed by PCR amplification from flanking DNA. It is freely scalable for both sample number and amplicon number and has no specialized equipment requirement. An optimized protocol is presented which provides a matrix of 96 × 192 combinations of samples x amplicons, corresponding to a complete 454 Titanium run. Initially, we used 454 sequencing; however, GTA could easily be adapted to Illumina sequencing platforms as read lengths have significantly increased in this system.     

New methods for finding common insertion sites and co-occurring common insertion sites in transposon- and virus-based genetic screens

Insertional mutagenesis screens in mice are used to identify individual genes that drive tumor formation. In these screens, candidate cancer genes are identified if their genomic location is proximal to a common insertion site (CIS) defined by high rates of transposon or retroviral insertions in a given genomic window. In this article, we describe a new method for defining CISs based on a Poisson distribution, the Poisson Regression Insertion Model, and show that this new method is an improvement over previously described methods. We also describe a modification of the method that can identify pairs and higher orders of co-occurring common insertion sites. We apply these methods to two data sets, one generated in a transposon-based screen for gastrointestinal tract cancer genes and another based on the set of retroviral insertions in the Retroviral Tagged Cancer Gene Database. We show that the new methods identify more relevant candidate genes and candidate gene pairs than found using previous methods. Identification of the biologically relevant set of mutations that occur in a single cell and cause tumor progression will aid in the rational design of single and combinatorial therapies in the upcoming age of personalized cancer therapy.     

Biases for detecting arbuscular mycorrhizal fungal mixture by terminal restriction fragment length polymorphism (T-RFLP)

Terminal restriction fragment length polymorphism (T-RFLP) analysis of amplified ribosomal RNA genes is used for profiling microbial communities and sometimes for species richness and relative abundance estimation in environmental samples. However, the T-RFLP fingerprint may be subject to biases during the procedure, influencing the detection of real community structures in the environment. To investigate possible sources of T-RFLP bias, 18S rRNA gene clones derived from two arbuscular mycorrhizal fungal sequences were combined in simple pairwise mixes to assess the effects of polymerase chain reaction cycle number, plant genomic DNA purification method and varying template ratio on the template-to-product ratio as measured by relative T-RF peak area. Varying cycle numbers indicated that amplification was still in the exponential phase at the cycle numbers lower than 18, so these small cycle numbers were used for the comparison of template-to-product quantities. Relative abundance estimated from T-RF peak ratios varied with different purification procedures, but the best results, closest to input ratios, were obtained by using phenol–chloroform purification. The presence of an excess of unpurified non-target plant genomic DNA generated a bias towards lower or overestimation of relative abundance. We conclude that a low number of amplification cycles and stringent DNA purification are necessary for accurate mixed sample analysis by T-RFLP.     

Sleeping Beauty transposition in the mouse genome is associated with changes in DNA methylation at the site of insertion

The Sleeping Beauty (SB) transposon (Tn) system is a nonviral gene delivery tool that has widespread application for transfer of therapeutic genes into the mammalian genome. To determine its utility as a gene delivery system, it was important to assess the epigenetic modifications associated with SB insertion into the genome and potential inactivation of the transgene. This study investigated the DNA methylation pattern of an SB Tn as well as the flanking genomic region at insertion sites in the mouse genome. The ubiquitous ROSA26 promoter and an initial part of the eGFP coding sequence in the SB Tn exhibited high levels of CpG methylation in transgenic mouse lines, irrespective of the chromosomal loci of the insertion sites. In contrast, no detectable CpG methylation in the endogenous mouse ROSA26 counterpart was observed in the same animals. Furthermore, significant hypomethylation was detected in neighboring chromosomal sequences of two unique SB Tn insertions compared to wild-type patterns. Taken together, these results suggest that SB Tn insertions into the mouse genome can be discriminated by DNA methylation machinery from an identical endogenous DNA sequence and can profoundly alter the DNA methylation status of the transgene cargo as well as flanking host genomic regions.     

Evidence that HIV budding in primary macrophages occurs through the exosome release pathway

Lipid rafts are specialized regions of cell membranes enriched in cholesterol and sphingolipids that are involved in immune activation and signaling. Studies in T-cells indicate that these membrane domains serve as sites for release of human immunodeficiency virus (HIV). By budding through lipid rafts in T-cells, HIV selectively incorporates raft markers and excludes non-raft proteins. This process has been well studied in T-cells, but it is unknown whether lipid rafts serve as budding sites for HIV in macrophages. Recently, we proposed a new model of retroviral biogenesis called the Trojan exosome hypothesis (Gould, S. J., Booth, A., and Hildreth, J. E. K. (2003) Proc. Natl. Acad. Sci. U. S. A. 100, 10592-10597). This model proposes that retroviruses coopt the existing cellular machinery for exosomal release. Here, we performed the first test designed to differentiate between the lipid raft hypothesis of retroviral biogenesis and the Trojan exosome hypothesis. Using macrophages, we examined the relative abundance of several host proteins on the cell surface, in lipid rafts, and on both HIV particles and exosomes derived from these cells. Our results show significant differences in the abundance of host proteins on the cell surface and in HIV. Moreover, our data demonstrate discordance in the abundance of some proteins in lipid rafts and in HIV. Finally, our data reveal a strong concordance between the host cell protein profile of exosomes and that of HIV. These results strongly support the Trojan exosome hypothesis and its prediction that retroviral budding represents exploitation of a pre-existing cellular pathway of intercellular vesicle trafficking.     

Population dynamics of an Ac-like transposable element in self- and cross-pollinating arabidopsis.

Theoretical models predict that the mating system should be an important factor driving the dynamics of transposable elements in natural populations due to differences in selective pressure on both element and host. We used a PCR-based approach to examine the abundance and levels of insertion polymorphism of Ac-III, a recently identified Ac-like transposon family, in natural populations of the selfing plant Arabidopsis thaliana and its close outcrossing relative, Arabidopsis lyrata. Although several insertions appeared to be ancient and shared between species, there is strong evidence for recent activity of this element family in both species. Sequences of the regions flanking insertions indicate that all Ac-III transposons segregating in natural populations are in noncoding regions and provide no evidence for local transposition events. Transposon display analysis suggests the presence of slightly higher numbers of insertion sites per individual but fewer total polymorphic insertions in the self-pollinating A. thaliana than A. lyrata. Element insertions appear to be segregating at significantly lower frequencies in A. lyrata than A. thaliana, which is consistent with a reduction in transposition rate, reduction in effective population size, or reduced efficacy of natural selection against element insertions in selfing populations.     

Global gene disruption in human cells to assign genes to phenotypes by deep sequencing

Insertional mutagenesis in a haploid background can disrupt gene function. We extend our earlier work by using a retroviral gene-trap vector to generate insertions in >98% of the genes expressed in a human cancer cell line that is haploid for all but one of its chromosomes. We apply phenotypic interrogation via tag sequencing (PhITSeq) to examine millions of mutant alleles through selection and parallel sequencing. Analysis of pools of cells, rather than individual clones enables rapid assessment of the spectrum of genes involved in the phenotypes under study. This facilitates comparative screens as illustrated here for the family of cytolethal distending toxins (CDTs). CDTs are virulence factors secreted by a variety of pathogenic Gram-negative bacteria responsible for tissue damage at distinct anatomical sites. We identify 743 mutations distributed over 12 human genes important for intoxication by four different CDTs. Although related CDTs may share host factors, they also exploit unique host factors to yield a profile characteristic for each CDT.     

Isolation and analysis of retroviral integration targets by solo long terminal repeat inverse PCR

Upon retroviral infection, the genomic RNA is reverse transcribed to make proviral DNA, which is then integrated into the host chromosome. Although the viral elements required for successful integration have been extensively characterized, little is known about the host DNA structure constituting preferred targets for proviral integration. In order to elucidate the mechanism for the target selection, comparison of host DNA sequences at proviral integration sites may be useful. To achieve simultaneous analysis of the upstream and downstream host DNA sequences flanking each proviral integration site, a Moloney murine leukemia virus-based retroviral vector was designed so that its integrated provirus could be removed by Cre-loxP homologous recombination, leaving a solo long terminal repeat (LTR). Taking advantage of the solo LTR, inverse PCR was carried out to amplify both the upstream and downstream cellular flanking DNA. The method called solo LTR inverse PCR, or SLIP, proved useful for simultaneously cloning the upstream and downstream flanking sequences of individual proviral integration sites from the polyclonal population of cells harboring provirus at different chromosomal sites. By the SLIP method, nucleotide sequences corresponding to 38 independent proviral integration targets were determined and, interestingly, atypical virus-host DNA junction structures were found in more than 20% of the cases. Characterization of retroviral integration sites using the SLIP method may provide useful insights into the mechanism for proviral integration and its target selection.     

Restriction-map variation with the yellow-achaete-scute region in five populations of Drosophila melanogaster

It has been proposed that the degree of recombination for a genomic region will affect the level of both nucleotide heterozygosity and the density of transposable elements. Both features of genomic diversity have been examined in a number of recent reports for regions undergoing relatively normal levels of recombination in Drosophila melanogaster. In this study the genomic variation associated with yellow-achaete- scute loci located at the tip of the X chromosome is examined by six- cutter restriction mapping. In this region, as usual for regions adjacent to telomeres, crossing-over is dramatically reduced, and published studies of visible mutants indicate extremely little restriction-map variation. Eight six-cutter restriction endonucleases were used to locate sequence variation in 14- and 16.5-kb regions in 109 lines sampled from North America, Africa, and Europe. The overall level of heterozygosity is estimated as 0.29%. Nine large insertions, all presumed to be transposable elements, were observed. Base-pair heterozygosity appears to be reduced compared with regions having normal levels of recombination. The estimated heterozygosity is much higher than reported in earlier studies of restriction-map variation among visible mutations in the complex. The incidence of large insertions is not elevated compared with that in other regions of the genome. This suggests that asymmetric synapsis and exchange is not an important mechanism for the elimination of transposable elements.      

Rapid amplification and cloning of Tn5 flanking fragments by inverse PCR

A simple approach is described to efficiently amplify DNA sequences flanking transposon Tn5 insertions. The method involves: (i) digestion with a restriction enzyme that cuts within Tn5; (ii) self-ligation under conditions favouring the production of monomeric circles; (iii) four parallel PCR reactions using primers designed to amplify left or right flanking sequences, and to distinguish target amplicons from non-specific products. This reveals the number of Tn5 insertions and the size of flanking genomic restriction fragments, without Southern blot analysis. The amplified product contains restriction sites that facilitate cohesive-end cloning. This rapid method is demonstrated using Tn5 and Tn5-Mob tagged DNA sequences involved in albicidin biosynthesis in Xanthomonas albilineans. It is generally applicable for efficient recovery of DNA sequences flanking transposon Tn5 derivatives in insertional mutagenesis studies.