Rational Development of Liquid Biopsy Analysis in Renal Cell Carcinoma.

Renal cell carcinoma (RCC) is known for its variable clinical behavior and outcome, including heterogeneity in developing relapse or metastasis. Recent data highlighted the potential of somatic mutations as promising biomarkers for risk stratification in RCC. Likewise, the analysis of circulating tumor DNA (ctDNA) for such informative somatic mutations (liquid biopsy) is considered an important advance for precision oncology in RCC, allowing to monitor molecular disease evolution in real time. However, our knowledge about the utility of ctDNA analysis in RCC is limited, in part due to the lack of RCC-appropriate assays for ctDNA analysis. Here, by interrogating different blood compartments in xenograft models, we identified plasma cell-free (cf) DNA and extracellular vesicles (ev) DNA enriched for RCC-associated ctDNA. Additionally, we developed sensitive targeted sequencing and bioinformatics workflows capable of detecting somatic mutations in RCC-relevant genes with allele frequencies ≥ 0.5%. Applying this assay to patient-matched tumor and liquid biopsies, we captured tumor mutations in cf- and ev-DNA fractions isolated from the blood, highlighting the potentials of both fractions for ctDNA analysis. Overall, our study presents an RCC-appropriate sequencing assay and workflow for ctDNA analysis and provides a proof of principle as to the feasibility of detecting tumor-specific mutations in liquid biopsy in RCC patients.

A Distributed Whole Genome Sequencing Benchmark Study.

Population sequencing often requires collaboration across a distributed network of sequencing centers for the timely processing of thousands of samples. In such massive efforts, it is important that participating scientists can be confident that the accuracy of the sequence data produced is not affected by which center generates the data. A study was conducted across three established sequencing centers, located in Montreal, Toronto, and Vancouver, constituting Canada's Genomics Enterprise (www.cgen.ca). Whole genome sequencing was performed at each center, on three genomic DNA replicates from three well-characterized cell lines. Secondary analysis pipelines employed by each site were applied to sequence data from each of the sites, resulting in three datasets for each of four variables (cell line, replicate, sequencing center, and analysis pipeline), for a total of 81 datasets. These datasets were each assessed according to multiple quality metrics including concordance with benchmark variant truth sets to assess consistent quality across all three conditions for each variable. Three-way concordance analysis of variants across conditions for each variable was performed. Our results showed that the variant concordance between datasets differing only by sequencing center was similar to the concordance for datasets differing only by replicate, using the same analysis pipeline. We also showed that the statistically significant differences between datasets result from the analysis pipeline used, which can be unified and updated as new approaches become available. We conclude that genome sequencing projects can rely on the quality and reproducibility of aggregate data generated across a network of distributed sites.

Single-cell analysis of human adipose tissue identifies depot and disease specific cell types.

The complex relationship between metabolic disease risk and body fat distribution in humans involves cellular characteristics which are specific to body fat compartments. Here we show depot-specific differences in the stromal vascual fraction of visceral and subcutaneous adipose tissue by performing single-cell RNA sequencing of tissue specimen from obese individuals. We characterize multiple immune cells, endothelial cells, fibroblasts, adipose and hematopoietic stem cell progenitors. Subpopulations of adipose-resident immune cells are metabolically active and associated with metabolic disease status and those include a population of potential dysfunctional CD8+ T cells expressing metallothioneins. We identify multiple types of adipocyte progenitors that are common across depots, including a subtype enriched in individuals with type 2 diabetes. Depot-specific analysis reveals a class of adipocyte progenitors unique to visceral adipose tissue, which shares common features with beige preadipocytes. Our human single-cell transcriptome atlas across fat depots provides a resource to dissect functional genomics of metabolic disease.

Iron Acquisition in Mycobacterium avium subsp. paratuberculosis.

UNLABELLED: Mycobacterium avium subsp. paratuberculosis is a host-adapted pathogen that evolved from the environmental bacterium M. avium subsp. hominissuis through gene loss and gene acquisition. Growth of M. avium subsp. paratuberculosis in the laboratory is enhanced by supplementation of the media with the iron-binding siderophore mycobactin J. Here we examined the production of mycobactins by related organisms and searched for an alternative iron uptake system in M. avium subsp. paratuberculosis. Through thin-layer chromatography and radiolabeled iron-uptake studies, we showed that M. avium subsp. paratuberculosis is impaired for both mycobactin synthesis and iron acquisition. Consistent with these observations, we identified several mutations, including deletions, in M. avium subsp. paratuberculosis genes coding for mycobactin synthesis. Using a transposon-mediated mutagenesis screen conditional on growth without myobactin, we identified a potential mycobactin-independent iron uptake system on a M. avium subsp. paratuberculosis-specific genomic island, LSP(P)15. We obtained a transposon (Tn) mutant with a disruption in the LSP(P)15 gene MAP3776c for targeted study. The mutant manifests increased iron uptake as well as intracellular iron content, with genes downstream of the transposon insertion (MAP3775c to MAP3772c [MAP3775-2c]) upregulated as the result of a polar effect. As an independent confirmation, we observed the same iron uptake phenotypes by overexpressing MAP3775-2c in wild-type M. avium subsp. paratuberculosis. These data indicate that the horizontally acquired LSP(P)15 genes contribute to iron acquisition by M. avium subsp. paratuberculosis, potentially allowing the subsequent loss of siderophore production by this pathogen. IMPORTANCE: Many microbes are able to scavenge iron from their surroundings by producing iron-chelating siderophores. One exception is Mycobacterium avium subsp. paratuberculosis, a fastidious, slow-growing animal pathogen whose growth needs to be supported by exogenous mycobacterial siderophore (mycobactin) in the laboratory. Data presented here demonstrate that, compared to other closely related M. avium subspecies, mycobactin production and iron uptake are different in M. avium subsp. paratuberculosis, and these phenotypes may be caused by numerous deletions in its mycobactin biosynthesis pathway. Using a genomic approach, supplemented by targeted genetic and biochemical studies, we identified that LSP(P)15, a horizontally acquired genomic island, may encode an alternative iron uptake system. These findings shed light on the potential physiological consequence of horizontal gene transfer in M. avium subsp. paratuberculosis evolution.

BisQC: an operational pipeline for multiplexed bisulfite sequencing.

BACKGROUND: Bisulfite sequencing is the most efficient single nucleotide resolution method for analysis of methylation status at whole genome scale, but improved quality control metrics are needed to better standardize experiments. RESULTS: We describe BisQC, a step-by-step method for multiplexed bisulfite-converted DNA library construction, pooling, spike-in content, and bioinformatics. We demonstrate technical improvements for library preparation and bioinformatic analyses that can be done in standard laboratories. We find that decoupling amplification of bisulfite converted (bis) DNA from the indexing reaction is an advantage, specifically in reducing total PCR cycle number and pre-selecting high quality bis-libraries. We also introduce a progressive PCR method for optimal library amplification and size-selection. At the sequencing stage, we thoroughly test the benefits of pooling non-bis DNA library with bis-libraries and find that BisSeq libraries can be pooled with a high proportion of non-bis DNA libraries with minimal impact on BisSeq output. For informatics analysis, we propose a series of optimization steps including the utilization of the mitochondrial genome as a QC standard, and we assess the validity of using duplicate reads for coverage statistics. CONCLUSION: We demonstrate several quality control checkpoints at the library preparation, pre-sequencing, post-sequencing, and post-alignment stages, which should prove useful in determining sample and processing quality. We also determine that including a significant portion of non-bisulfite converted DNA with bisulfite converted DNA has a minimal impact on usable bisulfite read output.

Fusion of TTYH1 with the C19MC microRNA cluster drives expression of a brain-specific DNMT3B isoform in the embryonal brain tumor ETMR.

Embryonal tumors with multilayered rosettes (ETMRs) are rare, deadly pediatric brain tumors characterized by high-level amplification of the microRNA cluster C19MC. We performed integrated genetic and epigenetic analyses of 12 ETMR samples and identified, in all cases, C19MC fusions to TTYH1 driving expression of the microRNAs. ETMR tumors, cell lines and xenografts showed a specific DNA methylation pattern distinct from those of other tumors and normal tissues. We detected extreme overexpression of a previously uncharacterized isoform of DNMT3B originating at an alternative promoter that is active only in the first weeks of neural tube development. Transcriptional and immunohistochemical analyses suggest that C19MC-dependent DNMT3B deregulation is mediated by RBL2, a known repressor of DNMT3B. Transfection with individual C19MC microRNAs resulted in DNMT3B upregulation and RBL2 downregulation in cultured cells. Our data suggest a potential oncogenic re-engagement of an early developmental program in ETMR via epigenetic alteration mediated by an embryonic, brain-specific DNMT3B isoform.

Mutations in SETD2 and genes affecting histone H3K36 methylation target hemispheric high-grade gliomas.

Recurrent mutations affecting the histone H3.3 residues Lys27 or indirectly Lys36 are frequent drivers of pediatric high-grade gliomas (over 30% of HGGs). To identify additional driver mutations in HGGs, we investigated a cohort of 60 pediatric HGGs using whole-exome sequencing (WES) and compared them to 543 exomes from non-cancer control samples. We identified mutations in SETD2, a H3K36 trimethyltransferase, in 15% of pediatric HGGs, a result that was genome-wide significant (FDR = 0.029). Most SETD2 alterations were truncating mutations. Sequencing the gene in this cohort and another validation cohort (123 gliomas from all ages and grades) showed SETD2 mutations to be specific to high-grade tumors affecting 15% of pediatric HGGs (11/73) and 8% of adult HGGs (5/65) while no SETD2 mutations were identified in low-grade diffuse gliomas (0/45). Furthermore, SETD2 mutations were mutually exclusive with H3F3A mutations in HGGs (P = 0.0492) while they partly overlapped with IDH1 mutations (4/14), and SETD2-mutant tumors were found exclusively in the cerebral hemispheres (P = 0.0055). SETD2 is the only H3K36 trimethyltransferase in humans, and SETD2-mutant tumors showed a substantial decrease in H3K36me3 levels (P < 0.001), indicating that the mutations are loss-of-function. These data suggest that loss-of-function SETD2 mutations occur in older children and young adults and are specific to HGG of the cerebral cortex, similar to the H3.3 G34R/V and IDH mutations. Taken together, our results suggest that mutations disrupting the histone code at H3K36, including H3.3 G34R/V, IDH1 and/or SETD2 mutations, are central to the genesis of hemispheric HGGs in older children and young adults.

Mapping of determinants required for the function of the HIV-1 env nuclear retention sequence.

Control of HIV-1 RNA processing and transport are critical to the successful replication of the virus. In previous work, we identified a region within the HIV-1 env that is involved in mediating nuclear retention of unspliced viral RNA. To define this sequence further and identify elements required for function, deletion mutagenesis was carried out. Progressive 5' and 3' deletions map the nuclear retention sequence (NRS) within the intron between nts 8281 and 8381. While deletion of sequences comprising the 3'ss had no effect, removal of the 5'ss resulted in cytoplasmic accumulation of unspliced RNA. Sequence analysis determined that the region corresponding to the NRS is highly conserved among HIV-1 strains. To evaluate whether this NRS interacts with cellular factors, RNA electrophoretic mobility shift assays (REMSA) were performed. We show that the NRS specifically interacts with cellular factors present in HeLa nuclear extracts, and, by UV crosslinking, correlates with the binding of a 49-kDa protein. Immunoprecipitation of the UV crosslinked products determined that this 49-kDa protein corresponds to hnRNP C.

A genomics-guided approach for discovering and expressing cryptic metabolic pathways.

Genome analysis of actinomycetes has revealed the presence of numerous cryptic gene clusters encoding putative natural products. These loci remain dormant until appropriate chemical or physical signals induce their expression. Here we demonstrate the use of a high-throughput genome scanning method to detect and analyze gene clusters involved in natural-product biosynthesis. This method was applied to uncover biosynthetic pathways encoding enediyne antitumor antibiotics in a variety of actinomycetes. Comparative analysis of five biosynthetic loci representative of the major structural classes of enediynes reveals the presence of a conserved cassette of five genes that includes a novel family of polyketide synthase (PKS). The enediyne PKS (PKSE) is proposed to be involved in the formation of the highly reactive chromophore ring structure (or "warhead") found in all enediynes. Genome scanning analysis indicates that the enediyne warhead cassette is widely dispersed among actinomycetes. We show that selective growth conditions can induce the expression of these loci, suggesting that the range of enediyne natural products may be much greater than previously thought. This technology can be used to increase the scope and diversity of natural-product discovery.

The calicheamicin gene cluster and its iterative type I enediyne PKS.

The enediynes exemplify nature's ingenuity. We have cloned and characterized the biosynthetic locus coding for perhaps the most notorious member of the nonchromoprotein enediyne family, calicheamicin. This gene cluster contains an unusual polyketide synthase (PKS) that is demonstrated to be essential for enediyne biosynthesis. Comparison of the calicheamicin locus with the locus encoding the chromoprotein enediyne C-1027 reveals that the enediyne PKS is highly conserved among these distinct enediyne families. Contrary to previous hypotheses, this suggests that the chromoprotein and nonchromoprotein enediynes are generated by similar biosynthetic pathways.