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.

Dysregulation of the PRUNE2/PCA3 genetic axis in human prostate cancer: from experimental discovery to validation in two independent patient cohorts.

Background: We have previously shown that the long non-coding (lnc)RNA prostate cancer associated 3 (PCA3; formerly prostate cancer antigen 3) functions as a trans-dominant negative oncogene by targeting the previously unrecognized prostate cancer suppressor gene PRUNE2 (a homolog of the Drosophila prune gene), thereby forming a functional unit within a unique allelic locus in human cells. Here, we investigated the PCA3/PRUNE2 regulatory axis from early (tumorigenic) to late (biochemical recurrence) genetic events during human prostate cancer progression. Methods: The reciprocal PCA3 and PRUNE2 gene expression relationship in paired prostate cancer and adjacent normal prostate was analyzed in two independent retrospective cohorts of clinically annotated cases post-radical prostatectomy: a single-institutional discovery cohort (n=107) and a multi-institutional validation cohort (n=497). We compared the tumor gene expression of PCA3 and PRUNE2 to their corresponding expression in the normal prostate. We also serially examined clinical/pathological variables including time to disease recurrence. Results: We consistently observed increased expression of PCA3 and decreased expression of PRUNE2 in prostate cancer compared with the adjacent normal prostate across all tumor grades and stages. However, there was no association between the relative gene expression levels of PCA3 or PRUNE2 and time to disease recurrence, independent of tumor grades and stages. Conclusions: We concluded that upregulation of the lncRNA PCA3 and targeted downregulation of the protein-coding PRUNE2 gene in prostate cancer could be early (rather than late) molecular events in the progression of human prostate tumorigenesis but are not associated with biochemical recurrence. Further studies of PCA3/PRUNE2 dysregulation are warranted. Funding: We received support from the Human Tissue Repository and Tissue Analysis Shared Resource from the Department of Pathology of the University of New Mexico School of Medicine and a pilot award from the University of New Mexico Comprehensive Cancer Center. RP and WA were supported by awards from the Levy-Longenbaugh Donor-Advised Fund and the Prostate Cancer Foundation. EDN reports research fellowship support from the Brazilian National Council for Scientific and Technological Development (CNPq), Brazil, and the Associação Beneficente Alzira Denise Hertzog Silva (ABADHS), Brazil. This work has been funded in part by the NCI Cancer Center Support Grants (CCSG; P30) to the University of New Mexico Comprehensive Cancer Center (CA118100) and the Rutgers Cancer Institute of New Jersey (CA072720).

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.

Targeting a cell surface vitamin D receptor on tumor-associated macrophages in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive tumor with limited treatment options and poor prognosis. We applied the in vivo phage display technology to isolate peptides homing to the immunosuppressive cellular microenvironment of TNBC as a strategy for non-malignant target discovery. We identified a cyclic peptide (CSSTRESAC) that specifically binds to a vitamin D receptor, protein disulfide-isomerase A3 (PDIA3) expressed on the cell surface of tumor-associated macrophages (TAM), and targets breast cancer in syngeneic TNBC, non-TNBC xenograft, and transgenic mouse models. Systemic administration of CSSTRESAC to TNBC-bearing mice shifted the cytokine profile toward an antitumor immune response and delayed tumor growth. Moreover, CSSTRESAC enabled ligand-directed theranostic delivery to tumors and a mathematical model confirmed our experimental findings. Finally, in silico analysis showed PDIA3-expressing TAM in TNBC patients. This work uncovers a functional interplay between a cell surface vitamin D receptor in TAM and antitumor immune response that could be therapeutically exploited.

Concordance Analysis of ALK Gene Fusion Detection Methods in Patients with Non-Small-Cell Lung Cancer from Chile, Brazil, and Peru.

About 4% to 7% of the non-small-cell lung cancer patients have anaplastic lymphoma kinase (ALK) rearrangements, and specific targeted therapies improve patients' outcomes significantly. ALK gene fusions are detected by immunohistochemistry or fluorescent in situ hybridization as gold standards in South America. Next-generation sequencing-based assays are a reliable alternative, able to perform simultaneous detection of multiple events from a single sample. We analyzed 4240 non-small-cell lung cancer samples collected in 37 hospitals from Chile, Brazil, and Peru, where ALK rearrangements were determined as part of their standard of care (SofC) using either immunohistochemistry or fluorescent in situ hybridization. A subset of 1450 samples was sequenced with the Oncomine Focus Assay (OFA), and the concordance with the SofC tests was measured. An orthogonal analysis was performed using a real-time quantitative PCR echinoderm microtubule-associated protein-like 4-ALK fusion detection kit. ALK fusion prevalence is similar for Chile (3.67%; N = 2142), Brazil (4.05%; N = 1013), and Peru (4.59%; N = 675). Although a comparison between OFA and SofC assays showed similar sensitivity, OFA had significantly higher specificity and higher positive predictive value, which opens new opportunities for a more specific determination of ALK gene rearrangements.

Deciphering the total RNA content of extracellular vesicles.

Extracellular vesicles (EVs) have been recognized as relevant players in cell-cell communication. To fully explore their potential as carriers of biological information in clinical settings, protocols capable of dealing with minute amounts of proteins, lipids, and nucleic acids present in their cargo are a requirement. Here we delve into a protocol to decipher the total transcriptome of EVs, from undetectable amounts of EVs-derived RNA from clinical samples.

APOBEC-mediated DNA alterations: A possible new mechanism of carcinogenesis in EBV-positive gastric cancer.

Mechanisms of viral oncogenesis are diverse and include the off-target activity of enzymes expressed by the infected cells, which evolved to target viral genomes for controlling their infection. Among these enzymes, the single-strand DNA editing capability of APOBECs represent a well-conserved viral infection response that can also cause untoward mutations in the host DNA. Here we show, after evaluating somatic single-nucleotide variations and transcriptome data in 240 gastric cancer samples, a positive correlation between APOBEC3s mRNA-expression and the APOBEC-mutation signature, both increased in EBV+ tumors. The correlation was reinforced by the observation of APOBEC mutations preferentially occurring in the genomic loci of the most active transcripts. This EBV infection and APOBEC3 mutation-signature axis were confirmed in a validation cohort of 112 gastric cancer patients. Our findings suggest that APOBEC3 upregulation in EBV+ cancer may boost the mutation load, providing further clues to the mechanisms of EBV-induced gastric carcinogenesis. After further validation, this EBV-APOBEC axis may prove to be a secondary driving force in the mutational evolution of EBV+ gastric tumors, whose consequences in terms of prognosis and treatment implications should be vetted.

A transcriptome-based signature of pathological angiogenesis predicts breast cancer patient survival.

The specific genes and molecules that drive physiological angiogenesis differ from those involved in pathological angiogenesis, suggesting distinct mechanisms for these seemingly related processes. Unveiling genes and pathways preferentially associated with pathologic angiogenesis is key to understanding its mechanisms, thereby facilitating development of novel approaches to managing angiogenesis-dependent diseases. To better understand these different processes, we elucidated the transcriptome of the mouse retina in the well-accepted oxygen-induced retinopathy (OIR) model of pathological angiogenesis. We identified 153 genes changed between normal and OIR retinas, which represent a molecular signature relevant to other angiogenesis-dependent processes such as cancer. These genes robustly predict the survival of breast cancer patients, which was validated in an independent 1,000-patient test cohort (40% difference in 15-year survival; p = 2.56 x 10-21). These results suggest that the OIR model reveals key genes involved in pathological angiogenesis, and these may find important applications in stratifying tumors for treatment intensification or for angiogenesis-targeted therapies.

A non-functional galanin receptor-2 in a multiple sclerosis patient.

Multiple Sclerosis (MS) is an inflammatory neurodegenerative disease that affects approximately 2.5 million people globally. Even though the etiology of MS remains unknown, it is accepted that it involves a combination of genetic alterations and environmental factors. Here, after performing whole exome sequencing, we found a MS patient harboring a rare and homozygous single nucleotide variant (SNV; rs61745847) of the G-protein coupled receptor (GPCR) galanin-receptor 2 (GALR2) that alters an important amino acid in the TM6 molecular toggle switch region (W249L). Nuclear magnetic resonance imaging showed that the hypothalamus (an area rich in GALR2) of this patient exhibited an important volumetric reduction leading to an enlarged third ventricle. Ex vivo experiments with patient-derived blood cells (AKT phosphorylation), as well as studies in recombinant cell lines expressing the human GALR2 (calcium mobilization and NFAT mediated gene transcription), showed that galanin (GAL) was unable to stimulate cell signaling in cells expressing the variant GALR2 allele. Live cell confocal microscopy showed that the GALR2 mutant receptor was primarily localized to intracellular endosomes. We conclude that the W249L SNV is likely to abrogate GAL-mediated signaling through GALR2 due to the spontaneous internalization of this receptor in this patient. Although this homozygous SNV was rare in our MS cohort (1:262 cases), our findings raise the potential importance of impaired neuroregenerative pathways in the pathogenesis of MS, warrant future studies into the relevance of the GAL/GALR2 axis in MS and further suggest the activation of GALR2 as a potential therapeutic route for this disease.

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

Technical challenges of working with extracellular vesicles.

Extracellular Vesicles (EVs) are gaining interest as central players in liquid biopsies, with potential applications in diagnosis, prognosis and therapeutic guidance in most pathological conditions. These nanosized particles transmit signals determined by their protein, lipid, nucleic acid and sugar content, and the unique molecular pattern of EVs dictates the type of signal to be transmitted to recipient cells. However, their small sizes and the limited quantities that can usually be obtained from patient-derived samples pose a number of challenges to their isolation, study and characterization. These challenges and some possible options to overcome them are discussed in this review.

A total transcriptome profiling method for plasma-derived extracellular vesicles: applications for liquid biopsies.

Extracellular vesicles (EVs) are key mediators of intercellular communication. Part of their biological effects can be attributed to the transfer of cargos of diverse types of RNAs, which are promising diagnostic and prognostic biomarkers. EVs found in human biofluids are a valuable source for the development of minimally invasive assays. However, the total transcriptional landscape of EVs is still largely unknown. Here we develop a new method for total transcriptome profiling of plasma-derived EVs by next generation sequencing (NGS) from limited quantities of patient-derived clinical samples, which enables the unbiased characterization of the complete RNA cargo, including both small- and long-RNAs, in a single library preparation step. This approach was applied to RNA extracted from EVs isolated by ultracentrifugation from the plasma of five healthy volunteers. Among the most abundant RNAs identified we found small RNAs such as tRNAs, miRNAs and miscellaneous RNAs, which have largely unknown functions. We also identified protein-coding and long noncoding transcripts, as well as circular RNA species that were also experimentally validated. This method enables, for the first time, the full spectrum of transcriptome data to be obtained from minute patient-derived samples, and will therefore potentially allow the identification of cell-to-cell communication mechanisms and biomarkers.

Tissue-Associated Bacterial Alterations in Rectal Carcinoma Patients Revealed by 16S rRNA Community Profiling.

Sporadic and inflammatory forms of colorectal cancer (CRC) account for more than 80% of cases. Recent publications have shown mechanistic evidence for the involvement of gut bacteria in the development of both CRC-forms. Whereas, colon and rectal cancer have been routinely studied together as CRC, increasing evidence show these to be distinct diseases. Also, the common use of fecal samples to study microbial communities may reflect disease state but possibly not the tumor microenvironment. We performed this study to evaluate differences in bacterial communities found in tissue samples of 18 rectal-cancer subjects when compared to 18 non-cancer controls. Samples were collected during exploratory colonoscopy (non-cancer group) or during surgery for tumor excision (rectal-cancer group). High throughput 16S rRNA amplicon sequencing of the V4-V5 region was conducted on the Ion PGM platform, reads were filtered using Qiime and clustered using UPARSE. We observed significant increases in species richness and diversity in rectal cancer samples, evidenced by the total number of OTUs and the Shannon and Simpson indexes. Enterotyping analysis divided our cohort into two groups, with the majority of rectal cancer samples clustering into one enterotype, characterized by a greater abundance of Bacteroides and Dorea. At the phylum level, rectal-cancer samples had increased abundance of candidate phylum OD1 (also known as Parcubacteria) whilst non-cancer samples had increased abundance of Planctomycetes. At the genera level, rectal-cancer samples had higher abundances of Bacteroides, Phascolarctobacterium, Parabacteroides, Desulfovibrio, and Odoribacter whereas non-cancer samples had higher abundances of Pseudomonas, Escherichia, Acinetobacter, Lactobacillus, and Bacillus. Two Bacteroides fragilis OTUs were more abundant among rectal-cancer patients seen through 16S rRNA amplicon sequencing, whose presence was confirmed by immunohistochemistry and enrichment verified by digital droplet PCR. Our findings point to increased bacterial richness and diversity in rectal cancer, along with several differences in microbial community composition. Our work is the first to present evidence for a possible role of bacteria such as B. fragilis and the phylum Parcubacteria in rectal cancer, emphasizing the need to study tissue-associated bacteria and specific regions of the gastrointestinal tract in order to better understand the possible links between the microbiota and rectal cancer.

Towards a transcriptome-based theranostic platform for unfavorable breast cancer phenotypes.

Inflammatory breast carcinoma (IBC) is one of the most lethal forms of human breast cancer, and effective treatment for IBC is an unmet clinical need in contemporary oncology. Tumor-targeted theranostic approaches are emerging in precision medicine, but only a few specific biomarkers are available. Here we report up-regulation of the 78-kDa glucose-regulated protein (GRP78) in two independent discovery and validation sets of specimens derived from IBC patients, suggesting translational promise for clinical applications. We show that a GRP78-binding motif displayed on either bacteriophage or adeno-associated virus/phage (AAVP) particles or loop-grafted onto a human antibody fragment specifically targets orthotopic IBC and other aggressive breast cancer models in vivo. To evaluate the theranostic value, we used GRP78-targeting AAVP particles to deliver the human Herpes simplex virus thymidine kinase type-1 (HSVtk) transgene, obtaining simultaneous in vivo diagnosis through PET imaging and tumor treatment by selective activation of the prodrug ganciclovir at tumor sites. Translation of this AAVP system is expected simultaneously to image, monitor, and treat the IBC phenotype and possibly other aggressive (e.g., invasive and/or metastatic) subtypes of breast cancer, based on the inducible cell-surface expression of the stress-response chaperone GRP78, and possibily other cell-surface receptors in human tumors.

PRUNE2 is a human prostate cancer suppressor regulated by the intronic long noncoding RNA PCA3.

Prostate cancer antigen 3 (PCA3) is the most specific prostate cancer biomarker but its function remains unknown. Here we identify PRUNE2, a target protein-coding gene variant, which harbors the PCA3 locus, thereby classifying PCA3 as an antisense intronic long noncoding (lnc)RNA. We show that PCA3 controls PRUNE2 levels via a unique regulatory mechanism involving formation of a PRUNE2/PCA3 double-stranded RNA that undergoes adenosine deaminase acting on RNA (ADAR)-dependent adenosine-to-inosine RNA editing. PRUNE2 expression or silencing in prostate cancer cells decreased and increased cell proliferation, respectively. Moreover, PRUNE2 and PCA3 elicited opposite effects on tumor growth in immunodeficient tumor-bearing mice. Coregulation and RNA editing of PRUNE2 and PCA3 were confirmed in human prostate cancer specimens, supporting the medical relevance of our findings. These results establish PCA3 as a dominant-negative oncogene and PRUNE2 as an unrecognized tumor suppressor gene in human prostate cancer, and their regulatory axis represents a unique molecular target for diagnostic and therapeutic intervention.

Synchronous down-modulation of miR-17 family members is an early causative event in the retinal angiogenic switch.

Six members of the microRNA-17 (miR-17) family were mapped to three different chromosomes, although they share the same seed sequence and are predicted to target common genes, among which are those encoding hypoxia-inducible factor-1α (HIF1A) and VEGFA. Here, we evaluated the in vivo expression profile of the miR-17 family in the murine retinopathy of prematurity (ROP) model, whereby Vegfa expression is highly enhanced at the early stage of retinal neovascularization, and we found simultaneous reduction of all miR-17 family members at this stage. Using gene reporter assays, we observed binding of these miRs to specific sites in the 3' UTRs of Hif1a and Vegfa. Furthermore, overexpression of these miRs decreased HIF1A and VEGFA expression in vitro. Our data indicate that this miR-17 family elicits a regulatory synergistic down-regulation of Hif1a and Vegfa expression in this biological model. We propose the existence of a coordinated regulatory network, in which diverse miRs are synchronously regulated to target the Hif1a transcription factor, which in turn, potentiates and reinforces the regulatory effects of the miRs on Vegfa to trigger and sustain a significant physiological response.

Vascular ligand-receptor mapping by direct combinatorial selection in cancer patients.

Molecules differentially expressed in blood vessels among organs or between damaged and normal tissues, are attractive therapy targets; however, their identification within the human vasculature is challenging. Here we screened a peptide library in cancer patients to uncover ligand-receptors common or specific to certain vascular beds. Surveying ~2.35 x 10(6) motifs recovered from biopsies yielded a nonrandom distribution, indicating that systemic tissue targeting is feasible. High-throughput analysis by similarity search, protein arrays, and affinity chromatography revealed four native ligand-receptors, three of which were previously unrecognized. Two are shared among multiple tissues (integrin α4/annexin A4 and cathepsin B/apolipoprotein E3) and the other two have a restricted and specific distribution in normal tissue (prohibitin/annexin A2 in white adipose tissue) or cancer (RAGE/leukocyte proteinase-3 in bone metastases). These findings provide vascular molecular markers for biotechnology and medical applications.

Next-generation phage display: integrating and comparing available molecular tools to enable cost-effective high-throughput analysis.

BACKGROUND: Combinatorial phage display has been used in the last 20 years in the identification of protein-ligands and protein-protein interactions, uncovering relevant molecular recognition events. Rate-limiting steps of combinatorial phage display library selection are (i) the counting of transducing units and (ii) the sequencing of the encoded displayed ligands. Here, we adapted emerging genomic technologies to minimize such challenges. METHODOLOGY/PRINCIPAL FINDINGS: We gained efficiency by applying in tandem real-time PCR for rapid quantification to enable bacteria-free phage display library screening, and added phage DNA next-generation sequencing for large-scale ligand analysis, reporting a fully integrated set of high-throughput quantitative and analytical tools. The approach is far less labor-intensive and allows rigorous quantification; for medical applications, including selections in patients, it also represents an advance for quantitative distribution analysis and ligand identification of hundreds of thousands of targeted particles from patient-derived biopsy or autopsy in a longer timeframe post library administration. Additional advantages over current methods include increased sensitivity, less variability, enhanced linearity, scalability, and accuracy at much lower cost. Sequences obtained by qPhage plus pyrosequencing were similar to a dataset produced from conventional Sanger-sequenced transducing-units (TU), with no biases due to GC content, codon usage, and amino acid or peptide frequency. These tools allow phage display selection and ligand analysis at >1,000-fold faster rate, and reduce costs approximately 250-fold for generating 10(6) ligand sequences. CONCLUSIONS/SIGNIFICANCE: Our analyses demonstrates that whereas this approach correlates with the traditional colony-counting, it is also capable of a much larger sampling, allowing a faster, less expensive, more accurate and consistent analysis of phage enrichment. Overall, qPhage plus pyrosequencing is superior to TU-counting plus Sanger sequencing and is proposed as the method of choice over a broad range of phage display applications in vitro, in cells, and in vivo.