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.

High-throughput genetic manipulation of multicellular organisms using a machine-vision guided embryonic microinjection robot.

Microinjection is a technique used for transgenesis, mutagenesis, cell labeling, cryopreservation, and in vitro fertilization in multiple single and multicellular organisms. Microinjection requires specialized skills and involves rate-limiting and labor-intensive preparatory steps. Here, we constructed a machine-vision guided generalized robot that fully automates the process of microinjection in fruit fly (Drosophila melanogaster) and zebrafish (Danio rerio) embryos. The robot uses machine learning models trained to detect embryos in images of agar plates and identify specific anatomical locations within each embryo in 3D space using dual view microscopes. The robot then serially performs a microinjection in each detected embryo. We constructed and used three such robots to automatically microinject tens of thousands of Drosophila and zebrafish embryos. We systematically optimized robotic microinjection for each species and performed routine transgenesis with proficiency comparable to highly skilled human practitioners while achieving up to 4× increases in microinjection throughput in Drosophila. The robot was utilized to microinject pools of over 20,000 uniquely barcoded plasmids into 1,713 embryos in 2 days to rapidly generate more than 400 unique transgenic Drosophila lines. This experiment enabled a novel measurement of the number of independent germline integration events per successfully injected embryo. Finally, we showed that robotic microinjection of cryoprotective agents in zebrafish embryos significantly improves vitrification rates and survival of cryopreserved embryos post-thaw as compared to manual microinjection. We anticipate that the robot can be used to carry out microinjection for genome-wide manipulation and cryopreservation at scale in a wide range of organisms.

Comparative analysis of genome-wide DNA methylation identifies patterns that associate with conserved transcriptional programs in osteosarcoma.

Osteosarcoma is an aggressive tumor of the bone that primarily affects young adults and adolescents. Osteosarcoma is characterized by genomic chaos and heterogeneity. While inactivation of tumor protein p53 (TP53) is nearly universal other high frequency mutations or structural variations have not been identified. Despite this genomic heterogeneity, key conserved transcriptional programs associated with survival have been identified across human, canine and induced murine osteosarcoma. The epigenomic landscape, including DNA methylation, plays a key role in establishing transcriptional programs in all cell types. The role of epigenetic dysregulation has been studied in a variety of cancers but has yet to be explored at scale in osteosarcoma. Here we examined genome-wide DNA methylation patterns in 24 human and 44 canine osteosarcoma samples identifying groups of highly correlated DNA methylation marks in human and canine osteosarcoma samples. We also link specific DNA methylation patterns to key transcriptional programs in both human and canine osteosarcoma. Building on previous work, we built a DNA methylation-based measure for the presence and abundance of various immune cell types in osteosarcoma. Finally, we determined that the underlying state of the tumor, and not changes in cell composition, were the main driver of differences in DNA methylation across the human and canine samples. SIGNIFICANCE: Genome wide comparison of DNA methylation patterns in osteosarcoma across two species lays the ground work for the exploration of DNA methylation programs that help establish conserved transcriptional programs in the context of varied mutational landscapes.

The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol.

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

Multiscale, multi-perspective imaging assisted robotic microinjection of 3D biological structures.

Microinjection is a widely used technique employed by biologists with applications in transgenesis, cryopreservation, mutagenesis, labeling/dye injection and in-vitro fertilization. However, microinjection is an extremely laborious manual procedure, which makes it a critical bottleneck in the field and thus ripe for automation. Here, we present a computer-guided robot that automates the targeted microinjection of Drosophila melanogaster and zebrafish (Danio rerio) embryos, two important model organisms in biological research. The robot uses a series of cameras to image an agar plate containing embryos at multiple magnifications and perspectives. This imaging is combined with machine learning and computer vision algorithms to pinpoint a location on the embryo for targeted microinjection with microscale precision. We demonstrate the utility of this microinjection robot to successfully microinject Drosophila melanogaster and zebrafish embryos. Results obtained indicate that the robotic microinjection approach can significantly increase the throughput of microinjection as compared to manual microinjection while maintaining survival rates comparable to human operators. In the future, this robotic platform can be used to perform high throughput microinjection experiments and can be extended to automatically microinject a host of organisms such as roundworms (Caenorhabditis elegans), mosquito (Culicidae) embryos, sea urchins (Echinoidea) and frog (Xenopus) oocytes.

AR gene rearrangement analysis in liquid biopsies reveals heterogeneity in lethal prostate cancer.

Castration-resistant prostate cancer (CRPC) is driven by AR gene aberrations that arise during androgen receptor (AR)-targeted therapy. AR amplification and mutations have been profiled in circulating tumor cells (CTCs), but whether AR gene rearrangements can be assessed in CTCs is unknown. In this study, we leveraged CRPC cell lines with defined AR gene rearrangements to develop and validate a CTC DNA analysis approach that utilized whole genome amplification and targeted DNA-sequencing of AR and other genes important in CRPC. We tested the utility of this approach by analyzing matched CTC DNA and plasma cell-free DNA (cfDNA) from a case series of ten CRPC patients. One of ten CTC samples and two of ten cfDNA samples were positive for AR gene rearrangements. All AR gene rearrangements were discordant between matched liquid biopsy samples. One patient harbored separate AR gene rearrangements in CTC DNA and cfDNA, but concordant AR amplification and AR T878A mutation. This patient also displayed concordant loss of TP53 and PTEN, but the loss of RB1 in cfDNA only. The overall frequency of discordant alterations in these genes between matched CTC DNA and cfDNA was high. This study establishes the technical feasibility of analyzing structural rearrangements, mutations, and copy number variants in AR and other CRPC genes using two different sources of DNA from a single blood sample. Paired CTC DNA and cfDNA analysis may have utility for capturing the heterogeneity of genetic alterations in CRPC patients.

Dissecting and tuning primer editing by proofreading polymerases.

Proofreading polymerases have 3' to 5' exonuclease activity that allows the excision and correction of mis-incorporated bases during DNA replication. In a previous study, we demonstrated that in addition to correcting substitution errors and lowering the error rate of DNA amplification, proofreading polymerases can also edit PCR primers to match template sequences. Primer editing is a feature that can be advantageous in certain experimental contexts, such as amplicon-based microbiome profiling. Here we develop a set of synthetic DNA standards to report on primer editing activity and use these standards to dissect this phenomenon. The primer editing standards allow next-generation sequencing-based enzymological measurements, reveal the extent of editing, and allow the comparison of different polymerases and cycling conditions. We demonstrate that proofreading polymerases edit PCR primers in a concentration-dependent manner, and we examine whether primer editing exhibits any sequence specificity. In addition, we use these standards to show that primer editing is tunable through the incorporation of phosphorothioate linkages. Finally, we demonstrate the ability of primer editing to robustly rescue the drop-out of taxa with 16S rRNA gene-targeting primer mismatches using mock communities and human skin microbiome samples.

A rapid, cost-effective tailed amplicon method for sequencing SARS-CoV-2.

BACKGROUND: The global COVID-19 pandemic has led to an urgent need for scalable methods for clinical diagnostics and viral tracking. Next generation sequencing technologies have enabled large-scale genomic surveillance of SARS-CoV-2 as thousands of isolates are being sequenced around the world and deposited in public data repositories. A number of methods using both short- and long-read technologies are currently being applied for SARS-CoV-2 sequencing, including amplicon approaches, metagenomic methods, and sequence capture or enrichment methods. Given the small genome size, the ability to sequence SARS-CoV-2 at scale is limited by the cost and labor associated with making sequencing libraries. RESULTS: Here we describe a low-cost, streamlined, all amplicon-based method for sequencing SARS-CoV-2, which bypasses costly and time-consuming library preparation steps. We benchmark this tailed amplicon method against both the ARTIC amplicon protocol and sequence capture approaches and show that an optimized tailed amplicon approach achieves comparable amplicon balance, coverage metrics, and variant calls to the ARTIC v3 approach. CONCLUSIONS: The tailed amplicon method we describe represents a cost-effective and highly scalable method for SARS-CoV-2 sequencing.

Inhalation exposure to cigarette smoke and inflammatory agents induces epigenetic changes in the lung.

Smoking-related lung tumors are characterized by profound epigenetic changes including scrambled patterns of DNA methylation, deregulated histone acetylation, altered gene expression levels, distorted microRNA profiles, and a global loss of cytosine hydroxymethylation marks. Here, we employed an enhanced version of bisulfite sequencing (RRBS/oxRRBS) followed by next generation sequencing to separately map DNA epigenetic marks 5-methyl-dC and 5-hydroxymethyl-dC in genomic DNA isolated from lungs of A/J mice exposed whole-body to environmental cigarette smoke for 10 weeks. Exposure to cigarette smoke significantly affected the patterns of cytosine methylation and hydroxymethylation in the lungs. Differentially hydroxymethylated regions were associated with inflammatory response/disease, organismal injury, and respiratory diseases and were involved in regulation of cellular development, function, growth, and proliferation. To identify epigenetic changes in the lung associated with exposure to tobacco carcinogens and inflammation, A/J mice were intranasally treated with the tobacco carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), the inflammatory agent lipopolysaccharide (LPS), or both. NNK alone caused minimal epigenetic alterations, while exposure either to LPS or NNK/LPS in combination led to increased levels of global cytosine methylation and formylation, reduced cytosine hydroxymethylation, decreased histone acetylation, and altered expression levels of multiple genes. Our results suggest that inflammatory processes are responsible for epigenetic changes contributing to lung cancer development.

Diverse AR Gene Rearrangements Mediate Resistance to Androgen Receptor Inhibitors in Metastatic Prostate Cancer.

PURPOSE: Prostate cancer is the second leading cause of male cancer deaths. Castration-resistant prostate cancer (CRPC) is a lethal stage of the disease that emerges when endocrine therapies are no longer effective at suppressing activity of the androgen receptor (AR) transcription factor. The purpose of this study was to identify genomic mechanisms that contribute to the development and progression of CRPC. EXPERIMENTAL DESIGN: We used whole-genome and targeted DNA-sequencing approaches to identify mechanisms underlying CRPC in an aggregate cohort of 272 prostate cancer patients. We analyzed structural rearrangements at the genome-wide level and carried out a detailed structural rearrangement analysis of the AR locus. We used genome engineering to perform experimental modeling of AR gene rearrangements and long-read RNA sequencing to analyze effects on expression of AR and truncated AR variants (AR-V). RESULTS: AR was among the most frequently rearranged genes in CRPC tumors. AR gene rearrangements promoted expression of diverse AR-V species. AR gene rearrangements occurring in the context of AR amplification correlated with AR overexpression. Cell lines with experimentally derived AR gene rearrangements displayed high expression of tumor-specific AR-Vs and were resistant to endocrine therapies, including the AR antagonist enzalutamide. CONCLUSIONS: AR gene rearrangements are an important mechanism of resistance to endocrine therapies in CRPC.

Measuring sequencer size bias using REcount: a novel method for highly accurate Illumina sequencing-based quantification.

Quantification of DNA sequence tags from engineered constructs such as plasmids, transposons, or other transgenes underlies many functional genomics measurements. Typically, such measurements rely on PCR followed by next-generation sequencing. However, PCR amplification can introduce significant quantitative error. We describe REcount, a novel PCR-free direct counting method. Comparing measurements of defined plasmid pools to droplet digital PCR data demonstrates that REcount is highly accurate and reproducible. We use REcount to provide new insights into clustering biases due to molecule length across different Illumina sequencers and illustrate the impacts on interpretation of next-generation sequencing data and the economics of data generation.