Search

A long-standing challenge in human microbiome research is achieving the taxonomic and functional resolution needed to generate testable hypotheses about the gut microbiota’s impact on health and disease. With a growing number of live microbial interventions in clinical development, this challenge is renewed by a need to understand the pharmacokinetics and pharmacodynamics of therapeutic candidates. While short-read sequencing of the bacterial 16S rRNA gene has been the standard for microbiota profiling, recent improvements in the fidelity of long-read sequencing underscores the need for a re-evaluation of the value of distinct microbiome-sequencing approaches. We leveraged samples from participants enrolled in a phase 1b clinical trial of a novel live biotherapeutic product to perform a comparative analysis of short-read and long-read amplicon and metagenomic sequencing approaches to assess their utility for generating clinical microbiome data. Across all methods, overall community taxonomic profiles were comparable and relationships between samples were conserved. Comparison of ubiquitous short-read 16S rRNA amplicon profiling to long-read profiling of the 16S-ITS-23S rRNA amplicon showed that only the latter provided strain-level community resolution and insight into novel taxa. All methods identified an active ingredient strain in treated study participants, though detection confidence was higher for long-read methods. Read coverage from both metagenomic methods provided evidence of active-ingredient strain replication in some treated participants. Compared to short-read metagenomics, approximately twice the proportion of long reads were assigned functional annotations. Finally, compositionally similar bacterial metagenome-assembled genomes (MAGs) were recovered from short-read and long-read metagenomic methods, although a greater number and more complete MAGs were recovered from long reads. Despite higher costs, both amplicon and metagenomic long-read approaches yielded added microbiome data value in the form of higher confidence taxonomic and functional resolution and improved recovery of microbial genomes compared to traditional short-read methodologies.

A longstanding challenge in human microbiome research is achieving the taxonomic and functional resolution needed to generate testable hypotheses about the gut microbiome’s impact on health and disease. More recently, this challenge has extended to a need for in-depth understanding of the pharmacokinetics and pharmacodynamics of clinical microbiome-based interventions. Whole genome metagenomic sequencing provides high taxonomic resolution and information on metagenome functional capacity, but the required deep sequencing is costly. For this reason, short-read sequencing of the bacterial 16S ribosomal RNA (rRNA) gene is the standard for microbiota profiling, despite its poor taxonomic resolution. The recent falling costs and improved fidelity of long-read sequencing warrant an evaluation of this approach for clinical microbiome analysis. We used samples from participants enrolled in a Phase 1b clinical trial of a novel live biotherapeutic product to perform a comparative analysis of short-read and long-read amplicon and metagenomic sequencing approaches to assess their value for generating informative and actionable clinical microbiome data. Comparison of ubiquitous short-read 16S rRNA amplicon profiling to long-read profiling of the 16S-ITS-23S rRNA amplicon showed that only the latter provided strain-level community resolution and insight into novel taxa. Across all methods, overall community taxonomic profiles were comparable and relationships between samples were conserved, highlighting the accuracy of modern microbiome analysis pipelines. All methods identified an active ingredient strain in treated study participants, though detection confidence was higher for long-read methods. Read coverage from both metagenomic methods provided evidence of active ingredient strain replication in some treated participants. Compared to short-read metagenomics, approximately twice the proportion of long reads were assigned functional annotations (63% vs. 34%). Finally, similar bacterial metagenome-assembled genomes (MAGs) were recovered across short-read and long-read metagenomic methods, although MAGs recovered from long reads were more complete. Overall, despite higher costs, long-read microbiome characterization provides added scientific value for clinical microbiome research in the form of higher taxonomic and functional resolution and improved recovery of microbial genomes compared to traditional short-read methodologies.Data SummaryAll supporting data, code and protocols have been provided within the article or as supplementary data files. Two supplementary figures and four supplementary tables are available with the online version of this article. Sequencing data are accessible in the National Center for Biotechnology Information (NCBI) database under BioProject accession number PRJNA754443. The R code and additional data files used for analysis and figure generation are accessible in a GitHub repository (https://github.com/jeanette-gehrig/Gehrig_et_al_sequencing_comparison).Impact StatementAccurate sequencing and analysis are essential for informative microbiome profiling, which is critical for the development of novel microbiome-targeted therapeutics. Recent improvements in long-read sequencing technology provide a promising, but more costly, alternative to ubiquitous short-read sequencing. To our knowledge, a direct comparison of the informational value of short-read and HiFi long-read sequencing approaches has not been reported for clinical microbiome samples. Using samples from participants in a Phase 1b trial of a live biotherapeutic product, we compare microbiome profiles generated from short-read and long-read sequencing for both amplicon-based 16S ribosomal RNA profiling and metagenomic sequencing. Though overall taxonomic profiles were similar across methods, only long-read amplicon sequencing provided strain-level resolution, and long-read metagenomic sequencing resulted in a significantly greater proportion of functionally annotated genes. Detection of a live biotherapeutic active ingredient strain in treated participants was achieved with all methods, and both metagenomic methods provided evidence of active replication of this strain in some participants. Similar taxonomies were recovered through metagenomic assemblies of short and long reads, although assemblies were more complete with long reads. Overall, we show the utility of long-read microbiome sequencing in direct comparison to commonly used short-read methods for clinically relevant microbiome profiling.

Genomic DNA was amplified about 5 billion-fold from single, flow-sorted bacterial cells by the multiple displacement amplification (MDA) reaction, using φ 29 DNA polymerase. A 662-bp segment of the 16S rRNA gene could be accurately sequenced from the amplified DNA. MDA methods enable new strategies for studying nonculturable microorganisms.

Fundamental to most genetic analysis is availability of genomic DNA of adequate quality and quantity. Because DNA yield from human samples is frequently limiting, much effort has been invested in developing methods for whole genome amplification (WGA) by random or degenerate oligonucleotide-primed PCR. However, existing WGA methods like degenerate oligonucleotide-primed PCR suffer from incomplete coverage and inadequate average DNA size. We describe a method, termed multiple displacement amplification (MDA), which provides a highly uniform representation across the genome. Amplification bias among eight chromosomal loci was less than 3-fold in contrast to 4–6 orders of magnitude for PCR-based WGA methods. Average product length was >10 kb. MDA is an isothermal, strand-displacing amplification yielding about 20–30 μg product from as few as 1–10 copies of human genomic DNA. Amplification can be carried out directly from biological samples including crude whole blood and tissue culture cells. MDA-amplified human DNA is useful for several common methods of genetic analysis, including genotyping of single nucleotide polymorphisms, chromosome painting, Southern blotting and restriction fragment length polymorphism analysis, subcloning, and DNA sequencing. MDA-based WGA is a simple and reliable method that could have significant implications for genetic studies, forensics, diagnostics, and long-term sample storage.

Estrogen signaling is mediated by ERα and -β. ERs are converted from an inactive form to a transcriptionally active state through conformational changes induced by ligand and estrogen-responsive element (ERE) sequences. We show here that ERα and ERβ bind to an ERE independently from ER ligands. We found that although the binding affinity of ERβ for an ERE is 2-fold lower than that of ERα, both ERs use the same nucleotides for DNA contacts. We show that both EREs and ligands are independent modulators of ER conformation. Specifically, the ERE primarily determines the receptor-DNA affinity, whereas the structure of the ER ligand dictates the affinity of ER for particular cofactors. We found that the ligand-dependent cofactor transcriptional intermediary factor-2, through a distinct surface, also interacts with ERα preferentially and independently of ligand. The extent of interaction, however, is dependent upon the ER-ERE affinity. In transfected cells, ERα is more transcriptionally active than ERβ. The ERE sequence, however, determines the potency of gene induction when either ER subtype binds to an agonist. Antagonists prevent ERs from inducing transcription independently from ERE sequences. Thus, ERE- and ligand-induced structural changes are independent determinants for the recruitment of cofactors and transcriptional responses. The ability of ERα to differentially recruit a cofactor could contribute to ER subtype-specific gene responses.

Estrogen-induced signaling mediated by estrogen receptors (ERs) is also affected by aberrant ERs that act as constitutively active or dominant negative modulators. Variant ERs can contribute to carcinogenesis and to the loss of estrogen responsiveness, rendering antiestrogen therapy ineffective. Determining target gene response during co-synthesis of different ER species is difficult, because dimers formed in the presence of more than one ER species are a heterogenous population of homo- or heterodimers. We engineered a homofusion ERalpha as a prototype single-chain receptor by genetically conjugating two ER monomers into a covalently fused single-chain protein to obtain a homogeneous population. This permits analysis of symmetrical or asymmetrical mutations that simulate variant homo- and heterodimers. Although a monomer, the homofusion receptor exhibited similar biochemical and functional properties to the dimeric ERalpha. We used activation function-2 (AF2) defective mutants as a model in either one or both receptor domains for a dominant-negative phenotype by suppressing the reporter activity induced by the WT receptor. When co-expressed with ERalpha, the fusion variant deficient in both AF2 functions suppressed the reporter activity effectively induced by ERalpha. These results show the utility of fusion receptors as models for generation of receptor-based agonists and antagonists.

Human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT), nucleocapsid protein (NC), genomic RNA, and the growing DNA strand all influence the copying of the HIV-1 RNA genome into DNA. A detailed understanding of these activities is required to understand the process of reverse transcription. HIV-1 viral DNA is initiated from a tRNA 3 Lys primer bound to the viral genome at the primer binding site. The U3 and R regions of the RNA genome are the first sequences to be copied. The TAR hairpin, a structure found within the R region of the viral genome, is the site of increased RT pausing, RNase H activity, and RT dissociation. Template RNA was digested approximately 17 bases behind the site where polymerase paused at the base of TAR. In most template RNAs, this was the only cleavage made by the RT responsible for initiating polymerization. If the RT that initiated DNA synthesis dissociated from the base of the TAR hairpin and an RT rebound at the end of the primer, there was competition between the polymerase and RNase H activities. After the complete heteroduplex was formed, there were additional RNase H cleavages that did not involve polymerization. Levels of NC that prevented TAR DNA self-priming did not protect genomic RNA from RNase H digestion. RNase H digestion of the 100-bp heteroduplex produced a 14-base RNA from the 5′ end of the RNA that remained annealed to the 3′ end of the minus-strand strong-stop DNA only if NC was present in the reaction.

The estrogen receptor (ER) is a transcription factor that binds to a specific DNA sequence found in the regulatory regions of estrogen-responsive genes, called the estrogen response element (ERE). Many genes that contain EREs have been identified, and most of these EREs contain one or more changes from the core consensus sequence, a 13-nucleotide segment with 10 nucleotides forming an inverted repeat. A number of genes have multiple copies of these imperfect EREs. In order to understand why natural EREs have developed in this manner, we have attempted to define the basic sequence requirements for ER binding. To this end, we measured the binding of homodimeric ER to a variety of nonconsensus EREs. We discovered that an ERE containing even a single change from the consensus may be unable to bind ER. However, an ERE with two changes from the consensus may be capable of binding avidly to ER in the context of certain flanking sequences. We found that changes in the sequences flanking a nonconsensus ERE can greatly alter ER-ERE affinity, either positively or negatively. Careful study of sequences flanking a series of EREs made it possible to develop rules that predict whether ER binds to a given natural ERE and also to predict the relative amounts of binding when comparing two EREs.

Chicken ovalbumin upstream promoter-transcription factor (COUP-TF) was identified as a low abundance protein in bovine uterus that co-purified with estrogen receptor (ER) in a ligand-independent manner and was separated from the ER by its lower retention on estrogen response element (ERE)-Sepharose. In gel mobility shift assays, COUP-TF bound as an apparent dimer to ERE and ERE half-sites. COUP-TF bound to an ERE half-site with high affinity,K d = 1.24 nm. In contrast, ER did not bind a single ERE half-site. None of the class II nuclear receptors analyzed, i.e. retinoic acid receptor, retinoid X receptor, thyroid receptor, peroxisome proliferator-activated receptor, or vitamin D receptor, were constituents of the COUP-TF·DNA binding complex detected in gel mobility shift assays. Direct interaction of COUP-TF with ER was indicated by GST “pull-down” and co-immunoprecipitation assays. The nature of the ER ligand influenced COUP-TF-ERE half-site binding. When ER was liganded by the antiestrogen 4-hydroxytamoxifen (4-OHT), COUP-TF-half-site interaction decreased. Conversely, COUP-TF transcribed and translated in vitroenhanced the ERE binding of purified estradiol (E2)-liganded ER but not 4-OHT-liganded ER. Co-transfection of ER-expressing MCF-7 human breast cancer cells with an expression vector for COUP-TFI resulted in a dose-dependent inhibition of E2-induced expression of a luciferase reporter gene under the control of three tandem copies of EREc38. The ability of COUP-TF to bind specifically to EREs and half-sites, to interact with ER, and to inhibit E2-induced gene expression suggests COUP-TF regulates ER action by both direct DNA binding competition and through protein-protein interactions.

Conformational differences between type I antiestrogen-liganded estrogen receptor and estradiol (E2)-liganded estrogen receptor (ER) are thought to be responsible for differentiating agonist versus antagonist ER activity at individual genes. To examine the impact of ER ligand on estrogen-response element (ERE) binding kinetics and receptor conformation, we quantitated the effect of site-directed, ER-specific antibodies raised against synthetic peptides corresponding to the DNA-binding domain of human ER on ER-ERE binding in vitro. Although 4-hydroxytamoxifen-liganded-ER (4-OHT-ER) and E2-ER bind a consensus ERE with equal high affinity, the stoichiometry of 4-OHT-ER-ERE binding at saturation is approximately 50% lower than that of E2-ER binding to all ERE sequences tested. In contrast, the ERE binding stoichiometry of tamoxifen aziridine-liganded ER (TAz-ER) is identical to that of E2-ER: one receptor dimer bound per ERE. The difference in binding stoichiometry is caused by dissociation of one molecule of 4-OHT from the ER as the dimeric receptor binds DNA. Addition of low concentrations of ER-specific polyclonal antibodies AT3A or AT3B prevented 4-OHT ligand dissociation, yielding an increase in specific 4-OHT-ER-ERE binding to a level equal to that of E2-ER- or TAz-ER-ERE binding. However, higher amounts of AT3A or AT3B inhibited specific ERE binding of both 4-OHT- and E2-ER. We conclude that differences in ER conformation when liganded with 4-OHT versus E2 are revealed by these antibodies and that such differences in receptor conformation may influence subsequent interaction of the receptor with other proteins necessary for transactivation.

Background Rolling circle amplification of ligated probes is a simple and sensitive means for genotyping directly from genomic DNA. SNPs and mutations are interrogated with open circle probes (OCP) that can be circularized by DNA ligase when the probe matches the genotype. An amplified detection signal is generated by exponential rolling circle amplification (ERCA) of the circularized probe. The low cost and scalability of ligation/ERCA genotyping makes it ideally suited for automated, high throughput methods. Results A retrospective study using human genomic DNA samples of known genotype was performed for four different clinically relevant mutations: Factor V Leiden, Factor II prothrombin, and two hemochromatosis mutations, C282Y and H63D. Greater than 99% accuracy was obtained genotyping genomic DNA samples from hundreds of different individuals. The combined process of ligation/ERCA was performed in a single tube and produced fluorescent signal directly from genomic DNA in less than an hour. In each assay, the probes for both normal and mutant alleles were combined in a single reaction. Multiple ERCA primers combined with a quenched-peptide nucleic acid (Q-PNA) fluorescent detection system greatly accellerated the appearance of signal. Probes designed with hairpin structures reduced misamplification. Genotyping accuracy was identical from either purified genomic DNA or genomic DNA generated using whole genome amplification (WGA). Fluorescent signal output was measured in real time and as an end point. Conclusions Combining the optimal elements for ligation/ERCA genotyping has resulted in a highly accurate single tube assay for genotyping directly from genomic DNA samples. Accuracy exceeded 99 % for four probe sets targeting clinically relevant mutations. No genotypes were called incorrectly using either genomic DNA or whole genome amplified sample.

Background Single nucleotide polymorphisms (SNPs) are the foundation of powerful complex trait and pharmacogenomic analyses. The availability of large SNP databases, however, has emphasized a need for inexpensive SNP genotyping methods of commensurate simplicity, robustness, and scalability. We describe a solution-based, microtiter plate method for SNP genotyping of human genomic DNA. The method is based upon allele discrimination by ligation of open circle probes followed by rolling circle amplification of the signal using fluorescent primers. Only the probe with a 3' base complementary to the SNP is circularized by ligation. Results SNP scoring by ligation was optimized to a 100,000 fold discrimination against probe mismatched to the SNP. The assay was used to genotype 10 SNPs from a set of 192 genomic DNA samples in a high-throughput format. Assay directly from genomic DNA eliminates the need to preamplify the target as done for many other genotyping methods. The sensitivity of the assay was demonstrated by genotyping from 1 ng of genomic DNA. We demonstrate that the assay can detect a single molecule of the circularized probe. Conclusions Compatibility with homogeneous formats and the ability to assay small amounts of genomic DNA meets the exacting requirements of automated, high-throughput SNP scoring.

Understanding how viral components collaborate to convert the human immunodeficiency virus type 1 genome from single-stranded RNA into double-stranded DNA is critical to the understanding of viral replication. Not only must the correct reactions be carried out, but unwanted side reactions must be avoided. After minus-strand strong stop DNA (−sssDNA) synthesis, degradation of the RNA template by the RNase H domain of reverse transcriptase (RT) produces single-stranded DNA that has the potential to self-prime at the imperfectly base-paired TAR hairpin, making continued DNA synthesis impossible. Although nucleocapsid protein (NC) interferes with −sssDNA self-priming in reverse transcription reactions in vitro, NC alone did not prevent self-priming of a synthetic −sssDNA oligomer. NC did not influence DNA bending and therefore cannot inhibit self-priming at hairpins by directly blocking hairpin formation. Using DNA oligomers as a model for genomic RNA fragments, we found that a 17-base DNA fragment annealed to the 3′ end of the −sssDNA prevented self-priming in the presence of NC. This implies that to avoid self-priming, an RNA-DNA hybrid that is more thermodynamically stable than the hairpin must remain within the hairpin region. This suggests that NC prevents self-priming by generating or stabilizing a thermodynamically favored RNA-DNA heteroduplex instead of the kinetically favored TAR hairpin. In support of this idea, sequence changes that increased base pairing in the DNA TAR hairpin resulted in an increase in self-priming in vitro. We present a model describing the role of NC-dependent inhibition of self-priming in first-strand transfer.

Estrogen-inducible genes contain an enhancer called the estrogen response element (ERE), a double-stranded inverted repeat. The estrogen receptor (ER) is generally thought to bind to the double-stranded ERE. However, some reports provide evidence that an ER homodimer can bind a single strand of the ERE and suggest that single-stranded ERE binding is the preferred binding mode for ER. Since these two models describe quite different mechanisms of receptor action, we have attempted to reconcile the observations. Analyzing DNA structure by nuclease sensitivity, we found that two identical molecules of a single strand of DNA containing the ERE sequence can partially anneal in an antiparallel manner. Bimolecular annealing produces double-stranded inverted repeats, with adjacent unannealed tails. The amount of annealing correlates exactly with the ability of ER to bind bimolecular EREs. Either strand of an ERE could anneal to itself in a way that would bind ER. We conclude that ER binds only the annealed double-stranded ERE both in vitro and in vivo.

Quantitative DNase I footprinting assays were employed to simultaneously measure the amount of estrogen receptor (ER) bound to each site in constructs containing multiple estrogen response elements (EREs). These assays revealed identical, high affinity ER-ERE binding, K d of approximately 0.25 nM, for estradiol-liganded ER (E 2 -ER), 4-hydroxytamoxifen liganded ER (4-OHT-ER), tamoxifen aziridine liganded ER (TAz-ER), and unliganded dimeric ER, for each ERE in constructs containing up to four tandem EREs. Increasing concentrations of ER resulted in the same pattern of occupancy for each ERE, whether or not the site was located near other EREs. Similarly, the presence or absence of E 2 , 4-OHT, or TAz ligand did not change ER-ERE interaction. Since activated ER-ERE binding affinity is identical, whether ER is liganded or unliganded, ligand cannot regulate ER-ERE binding affinity. These results support the hypothesis that ligand-dependent conformational changes primarily determine how ER interacts with components of the transcription initiation complex that mediate gene transactivation. In addition, footprint assays revealed that, following ER binding, an AT-rich site adjacent to the ERE becomes hypersensitive to DNase I digestion. This sequence may be easily or intrinsically bent, assisting in recruiting ER to ERE sites.