Variability and bias in microbiome metagenomic sequencing: an interlaboratory study comparing experimental protocols

Several studies have documented the significant impact of methodological choices in microbiome analyses. The myriad of methodological options available complicate the replication of results and generally limit the comparability of findings between independent studies that use differing techniques and measurement pipelines. Here we describe the Mosaic Standards Challenge (MSC), an international interlaboratory study designed to assess the impact of methodological variables on the results. The MSC did not prescribe methods but rather asked participating labs to analyze 7 shared reference samples (5 × human stool samples and 2 × mock communities) using their standard laboratory methods. To capture the array of methodological variables, each participating lab completed a metadata reporting sheet that included 100 different questions regarding the details of their protocol. The goal of this study was to survey the methodological landscape for microbiome metagenomic sequencing (MGS) analyses and the impact of methodological decisions on metagenomic sequencing results. A total of 44 labs participated in the MSC by submitting results (16S or WGS) along with accompanying metadata; thirty 16S rRNA gene amplicon datasets and 14 WGS datasets were collected. The inclusion of two types of reference materials (human stool and mock communities) enabled analysis of both MGS measurement variability between different protocols using the biologically-relevant stool samples, and MGS bias with respect to ground truth values using the DNA mixtures. Owing to the compositional nature of MGS measurements, analyses were conducted on the ratio of Firmicutes: Bacteroidetes allowing us to directly apply common statistical methods. The resulting analysis demonstrated that protocol choices have significant effects, including both bias of the MGS measurement associated with a particular methodological choices, as well as effects on measurement robustness as observed through the spread of results between labs making similar methodological choices. In the analysis of the DNA mock communities, MGS measurement bias was observed even when there was general consensus among the participating laboratories. This study was the result of a collaborative effort that included academic, commercial, and government labs. In addition to highlighting the impact of different methodological decisions on MGS result comparability, this work also provides insights for consideration in future microbiome measurement study design.

Evaluation of Bacteria and Fungi DNA Abundance in Human Tissues

Whereas targeted and shotgun sequencing approaches are both powerful in allowing the study of tissue-associated microbiota, the human: microorganism abundance ratios in tissues of interest will ultimately determine the most suitable sequencing approach. In addition, it is possible that the knowledge of the relative abundance of bacteria and fungi during a treatment course or in pathological conditions can be relevant in many medical conditions. Here, we present a qPCR-targeted approach to determine the absolute and relative amounts of bacteria and fungi and demonstrate their relative DNA abundance in nine different human tissue types for a total of 87 samples. In these tissues, fungi genomes are more abundant in stool and skin samples but have much lower levels in other tissues. Bacteria genomes prevail in stool, skin, oral swabs, saliva, and gastric fluids. These findings were confirmed by shotgun sequencing for stool and gastric fluids. This approach may contribute to a more comprehensive view of the human microbiota in targeted studies for assessing the abundance levels of microorganisms during disease treatment/progression and to indicate the most informative methods for studying microbial composition (shotgun versus targeted sequencing) for various samples types.

Abstract A26: Comparison of DNA extraction kits for metagenomic studies in feces

DNA extraction is a critical step in metagenomic studies since it requires the proper lysis of cell walls, allowing the recovery of nucleic acids that would represent the total diversity of communities. Recent studies indicated that the use of different DNA extraction methods in the same sample yielded distinct microbial profiles, highlighting the importance of standardization in DNA extraction methodologies to accurately measure the human microbiome. These issues are even more critical when limiting amounts of samples are available, or when studies target communities with reduced abundance, such as fungi, which represents only 0.1% of total gut microbiota, as has been poorly explored. Here we investigated the performance of three commercial kits for DNA extraction from human fecal samples: QIAamp DNA Stool Mini Kit (Qiagen), E.Z.N.A. Stool DNA Kit (Omega), and ZymoBIOMICS DNA Kit (Zymo), all with or without adding an enzymatic extraction with MetaPolyzyme (MAC4L - Sigma), in an attempt to improve the microbiome community profiling. After extraction, DNA-quality and yield were evaluated using Qubit and real-time PCR, V4-V5 16S rRNA and ITS2 followed by DNA sequencing, as well as further analysis using shotgun DNA sequencing. No significant differences were observed in DNA yield as evaluated by Qubit; an average of 40.3 ng/ul of DNA was obtained for all samples and kits. For both fungi and bacteria analysis, the addition of MAC4L in the extraction protocol apparently did not change the alpha diversity (observed and Shannon) for both individuals. Our analysis for samples extracted with different kits, with or without MAC4L, indicated the percentage of reads from bacteria to vary between 85-90%, whereas fungi were between 0.29-0.49%. Human DNA was also found, as expected, and relative amounts varied from 0.1 to 0.3%. qPCR showed that for all kits up to 50ng of DNA/15uL reaction can be safely used, with no significant DNA polymerase inhibition. These results showed similar efficacy for the three kits in representing bacteria and fungi from feces, and maybe because all these kits had a bead-beading step we had no clear benefits from adding extra enzymatic digestion steps with MAC4L. Finally, we decided to use the Omega kit in our lab, since it combines a simple extraction protocol and is reliable and cost-effective.

Citation Format: Isabella Kuniko T. Takenaka, Bruno S. Moda, Thaís F. Bartelli, Maria G. Amorim, Gabriela E. Albuquerque, Andrew M. Thomas, Emmanuel Dias-Neto, Diana N. Nunes. Comparison of DNA extraction kits for metagenomic studies in feces [abstract]. In: Proceedings of the AACR Special Conference on the Microbiome, Viruses, and Cancer; 2020 Feb 21-24; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2020;80(8 Suppl):Abstract nr A26.

Genomics and epidemiology for gastric adenocarcinomas (GE4GAC): a Brazilian initiative to study gastric cancer

Gastric cancer (GC) is the fifth most common type of cancer worldwide with high incidences in Asia, Central, and South American countries. This patchy distribution means that GC studies are neglected by large research centers from developed countries. The need for further understanding of this complex disease, including the local importance of epidemiological factors and the rich ancestral admixture found in Brazil, stimulated the implementation of the GE4GAC project. GE4GAC aims to embrace epidemiological, clinical, molecular and microbiological data from Brazilian controls and patients with malignant and pre-malignant gastric disease. In this letter, we summarize the main goals of the project, including subject and sample accrual and current findings.

Msc6p is required for mitochondrial translation initiation in the absence of formylated Met-tRNAfMet

Mitochondrial translation normally requires formylation of the initiator tRNA-met, a reaction catalyzed by the enzyme formyltransferase, Fmt1p and MTFMT in Saccharomyces cerevisiae and human mitochondria, respectively. Yeast fmt1 mutants devoid of Fmt1p, however, can synthesize all mitochondrial gene products by initiating translation with a non-formylated methionyl-tRNA. Yeast synthetic respiratory-deficient fmt1 mutants have uncovered several factors suggested to play a role in translation initiation with non-formylated methionyl-tRNA. Here, we present evidence that Msc6p, a member of the pentatricopeptide repeat (PPR) motif family, is another essential factor for mitochondrial translation in fmt1 mutants. The PPR motif is characteristic of RNA-binding proteins found in chloroplasts and plant and fungal mitochondria, and is generally involved in RNA stability and transport. Moreover, in the present study, we show that the respiratory deficiency of fmt1msc6 double mutants can be rescued by overexpression of the yeast mitochondrial initiation factor mIF-2, encoded by IFM1. The role of Msc6p in translational initiation is further supported by pull-down assays showing that it transiently interacts with mIF-2. Altogether, our data indicate that Msc6p is an important factor in mitochondrial translation with an auxiliary function related to the mIF-2-dependent formation of the initiation complex.

Structural and thermodynamic studies of the tobacco calmodulin-like rgs-CaM protein

The tobacco calmodulin-like protein rgs-CaM is involved in host defense against virus and is reported to possess an associated RNA silencing suppressor activity. Rgs-CaM is also believed to act as an antiviral factor by interacting and targeting viral silencing suppressors for autophagic degradation. Despite these functional data, calcium interplay in the modulation of rgs-CaM is still poorly understood. Here we show that rgs-CaM displays a prevalent alpha-helical conformation and possesses three functional Ca2+-binding sites. Using computational modeling and molecular dynamics simulation, we demonstrate that Ca2+ binding to rgs-CaM triggers expansion of its tertiary structure with reorientation of alpha-helices within the EF-hands. This conformational change leads to the exposure of a large negatively charged region that may be implicated in the electrostatic interactions between rgs-CaM and viral suppressors. Moreover, the kd values obtained for Ca2+ binding to the three functional sites are not within the affinity range of a typical Ca2+ sensor.