Metabolic liver cancer: associations of rare and common germline variants in one-carbon metabolism and DNA methylation genes.

Animal studies implicate one-carbon metabolism and DNA methylation genes in hepatocellular carcinoma (HCC) development in the setting of metabolic perturbations. Using human samples, we investigated the associations between common and rare variants in these closely related biochemical pathways and risk for metabolic HCC development in a multicenter international study. We performed targeted exome sequencing of 64 genes among 556 metabolic HCC cases and 643 cancer-free controls with metabolic conditions. Multivariable logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs), adjusting for multiple comparisons. Gene-burden tests were used for rare variant associations. Analyses were performed in the overall sample and among non-Hispanic whites. The results show that among non-Hispanic whites, presence of rare functional variants in ABCC2 was associated with 7-fold higher risk of metabolic HCC (OR = 6.92, 95% CI: 2.38-20.15, P = 0.0004), and this association remained significant when analyses were restricted to functional rare variants observed in ≥2 participants (cases 3.2% versus controls 0.0%, P = 1.02 × 10-5). In the overall multiethnic sample, presence of rare functional variants in ABCC2 was nominally associated with metabolic HCC (OR = 3.60, 95% CI: 1.52-8.58, P = 0.004), with similar nominal association when analyses were restricted to functional rare variants observed in ≥2 participants (cases 2.9% versus controls 0.2%, P = 0.006). A common variant in PNPLA3 (rs738409[G]) was associated with higher HCC risk in the overall sample (P = 6.36 × 10-6) and in non-Hispanic whites (P = 0.0002). Our findings indicate that rare functional variants in ABCC2 are associated with susceptibility to metabolic HCC in non-Hispanic whites. PNPLA3-rs738409 is also associated with metabolic HCC risk.

ACE2 and TMPRSS2 SARS-CoV-2 infectivity genes: deep mutational scanning and characterization of missense variants.

The human angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) proteins play key roles in the cellular internalization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the coronavirus responsible for the coronavirus disease of 2019 (COVID-19) pandemic. We set out to functionally characterize the ACE2 and TMPRSS2 protein abundance for variant alleles encoding these proteins that contained non-synonymous single-nucleotide polymorphisms (nsSNPs) in their open reading frames (ORFs). Specifically, a high-throughput assay, deep mutational scanning (DMS), was employed to test the functional implications of nsSNPs, which are variants of uncertain significance in these two genes. Specifically, we used a 'landing pad' system designed to quantify the protein expression for 433 nsSNPs that have been observed in the ACE2 and TMPRSS2 ORFs and found that 8 of 127 ACE2, 19 of 157 TMPRSS2 isoform 1 and 13 of 149 TMPRSS2 isoform 2 variant proteins displayed less than ~25% of the wild-type protein expression, whereas 4 ACE2 variants displayed 25% or greater increases in protein expression. As a result, we concluded that nsSNPs in genes encoding ACE2 and TMPRSS2 might potentially influence SARS-CoV-2 infectivity. These results can now be applied to DNA sequence data for patients infected with SARS-CoV-2 to determine the possible impact of patient-based DNA sequence variation on the clinical course of SARS-CoV-2 infection.

Molecular classification and identification of an aggressive signature in low-grade B-cell lymphomas.

Non-follicular low-grade B-cell lymphomas (LGBCL) are biologically diverse entities that share clinical and histologic features that make definitive pathologic categorization challenging. While most patients with LGBCL have an indolent course, some experience aggressive disease, highlighting additional heterogeneity across these subtypes. To investigate the potential for shared biology across subtypes, we performed RNA sequencing and applied machine learning approaches that identified five clusters of patients that grouped independently of subtype. One cluster was characterized by inferior outcome, upregulation of cell cycle genes, and increased tumor immune cell content. Integration of whole exome sequencing identified novel LGBCL mutations and enrichment of TNFAIP3 and BCL2 alterations in the poor survival cluster. Building on this, we further refined a transcriptomic signature associated with early clinical failure in two independent cohorts. Taken together, this study identifies unique clusters of LGBCL defined by novel gene expression signatures and immune profiles associated with outcome across diagnostic subtypes.

All2: A tool for selecting mosaic mutations from comprehensive multi-cell comparisons.

Accurate discovery of somatic mutations in a cell is a challenge that partially lays in immaturity of dedicated analytical approaches. Approaches comparing a cell's genome to a control bulk sample miss common mutations, while approaches to find such mutations from bulk suffer from low sensitivity. We developed a tool, All2, which enables accurate filtering of mutations in a cell without the need for data from bulk(s). It is based on pair-wise comparisons of all cells to each other where every call for base pair substitution and indel is classified as either a germline variant, mosaic mutation, or false positive. As All2 allows for considering dropped-out regions, it is applicable to whole genome and exome analysis of cloned and amplified cells. By applying the approach to a variety of available data, we showed that its application reduces false positives, enables sensitive discovery of high frequency mutations, and is indispensable for conducting high resolution cell lineage tracing.

Correction: All2: A tool for selecting mosaic mutations from comprehensive multi-cell comparisons.

[This corrects the article DOI: 10.1371/journal.pcbi.1009487.].

Inherited variants at 3q13.33 and 3p24.1 are associated with risk of diffuse large B-cell lymphoma and implicate immune pathways.

We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10-13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10-12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.

SLCO1B1: Application and Limitations of Deep Mutational Scanning for Genomic Missense Variant Function.

SLCO1B1 (solute carrier organic anion transporter family member 1B1) is an important transmembrane hepatic uptake transporter. Genetic variants in the SLCO1B1 gene have been associated with altered protein folding, resulting in protein degradation and decreased transporter activity. Next-generation sequencing (NGS) of pharmacogenes is being applied increasingly to associate variation in drug response with genetic sequence variants. However, it is difficult to link variants of unknown significance with functional phenotypes using "one-at-a-time" functional systems. Deep mutational scanning (DMS) using a "landing pad cell-based system" is a high-throughput technique designed to analyze hundreds of gene open reading frame (ORF) missense variants in a parallel and scalable fashion. We have applied DMS to analyze 137 missense variants in the SLCO1B1 ORF obtained from the Exome Aggregation Consortium project. ORFs containing these variants were fused to green fluorescent protein and were integrated into "landing pad" cells. Florescence-activated cell sorting was performed to separate the cells into four groups based on fluorescence readout indicating protein expression at the single cell level. NGS was then performed and SLCO1B1 variant frequencies were used to determine protein abundance. We found that six variants not previously characterized functionally displayed less than 25% and another 12 displayed approximately 50% of wild-type protein expression. These results were then functionally validated by transporter studies. Severely damaging variants identified by DMS may have clinical relevance for SLCO1B1-dependent drug transport, but we need to exercise caution since the relatively small number of severely damaging variants identified raise questions with regard to the application of DMS to intrinsic membrane proteins such as organic anion transporter protein 1B1. SIGNIFICANCE STATEMENT: The functional implications of a large numbers of open reading frame (ORF) "variants of unknown significance" (VUS) in transporter genes have not been characterized. This study applied deep mutational scanning to determine the functional effects of VUS that have been observed in the ORF of SLCO1B1(s olute carrier organic anion transporter family member 1B1). Several severely damaging variants were identified, studied, and validated. These observations have implications for both the application of deep mutational scanning to intrinsic membrane proteins and for the clinical effect of drugs and endogenous compounds transported by SLCO1B1.

SCELLECTOR: ranking amplification bias in single cells using shallow sequencing.

BACKGROUND: The study of mosaic mutation is important since it has been linked to cancer and various disorders. Single cell sequencing has become a powerful tool to study the genome of individual cells for the detection of mosaic mutations. The amount of DNA in a single cell needs to be amplified before sequencing and multiple displacement amplification (MDA) is widely used owing to its low error rate and long fragment length of amplified DNA. However, the phi29 polymerase used in MDA is sensitive to template fragmentation and presence of sites with DNA damage that can lead to biases such as allelic imbalance, uneven coverage and over representation of C to T mutations. It is therefore important to select cells with uniform amplification to decrease false positives and increase sensitivity for mosaic mutation detection. RESULTS: We propose a method, Scellector (single cell selector), which uses haplotype information to detect amplification quality in shallow coverage sequencing data. We tested Scellector on single human neuronal cells, obtained in vitro and amplified by MDA. Qualities were estimated from shallow sequencing with coverage as low as 0.3× per cell and then confirmed using 30× deep coverage sequencing. The high concordance between shallow and high coverage data validated the method. CONCLUSION: Scellector can potentially be used to rank amplifications obtained from single cell platforms relying on a MDA-like amplification step, such as Chromium Single Cell profiling solution.

Identification of missing variants by combining multiple analytic pipelines.

BACKGROUND: After decades of identifying risk factors using array-based genome-wide association studies (GWAS), genetic research of complex diseases has shifted to sequencing-based rare variants discovery. This requires large sample sizes for statistical power and has brought up questions about whether the current variant calling practices are adequate for large cohorts. It is well-known that there are discrepancies between variants called by different pipelines, and that using a single pipeline always misses true variants exclusively identifiable by other pipelines. Nonetheless, it is common practice today to call variants by one pipeline due to computational cost and assume that false negative calls are a small percent of total. RESULTS: We analyzed 10,000 exomes from the Alzheimer's Disease Sequencing Project (ADSP) using multiple analytic pipelines consisting of different read aligners and variant calling strategies. We compared variants identified by using two aligners in 50,100, 200, 500, 1000, and 1952 samples; and compared variants identified by adding single-sample genotyping to the default multi-sample joint genotyping in 50,100, 500, 2000, 5000 and 10,000 samples. We found that using a single pipeline missed increasing numbers of high-quality variants correlated with sample sizes. By combining two read aligners and two variant calling strategies, we rescued 30% of pass-QC variants at sample size of 2000, and 56% at 10,000 samples. The rescued variants had higher proportions of low frequency (minor allele frequency [MAF] 1-5%) and rare (MAF < 1%) variants, which are the very type of variants of interest. In 660 Alzheimer's disease cases with earlier onset ages of ≤65, 4 out of 13 (31%) previously-published rare pathogenic and protective mutations in APP, PSEN1, and PSEN2 genes were undetected by the default one-pipeline approach but recovered by the multi-pipeline approach. CONCLUSIONS: Identification of the complete variant set from sequencing data is the prerequisite of genetic association analyses. The current analytic practice of calling genetic variants from sequencing data using a single bioinformatics pipeline is no longer adequate with the increasingly large projects. The number and percentage of quality variants that passed quality filters but are missed by the one-pipeline approach rapidly increased with sample size.

Proteomic detection of immunoglobulin light chain variable region peptides from amyloidosis patient biopsies.

Immunoglobulin light chain (LC) amyloidosis (AL) is caused by deposition of clonal LCs produced by an underlying plasma cell neoplasm. The clonotypic LC sequences are unique to each patient, and they cannot be reliably detected by either immunoassays or standard proteomic workflows that target the constant regions of LCs. We addressed this issue by developing a novel sequence template-based workflow to detect LC variable (LCV) region peptides directly from AL amyloid deposits. The workflow was implemented in a CAP/CLIA compliant clinical laboratory dedicated to proteomic subtyping of amyloid deposits extracted from either formalin-fixed paraffin-embedded tissues or subcutaneous fat aspirates. We evaluated the performance of the workflow on a validation cohort of 30 AL patients, whose amyloidogenic clone was identified using a novel proteogenomics method, and 30 controls. The recall and negative predictive values of the workflow, when identifying the gene family of the AL clone, were 93 and 98%, respectively. Application of the workflow on a clinical cohort of 500 AL amyloidosis samples highlighted a bias in the LCV gene families used by the AL clones. We also detected similarity between AL clones deposited in multiple organs of systemic AL patients. In summary, AL proteomic data sets are rich in LCV region peptides of potential clinical significance that are recoverable with advanced bioinformatics.

PatternCNV: a versatile tool for detecting copy number changes from exome sequencing data.

MOTIVATION: Exome sequencing (exome-seq) data, which are typically used for calling exonic mutations, have also been utilized in detecting DNA copy number variations (CNVs). Despite the existence of several CNV detection tools, there is still a great need for a sensitive and an accurate CNV-calling algorithm with built-in QC steps, and does not require a paired reference for each sample. RESULTS: We developed a novel method named PatternCNV, which (i) accounts for the read coverage variations between exons while leveraging the consistencies of this variability across different samples; (ii) reduces alignment BAM files to WIG format and therefore greatly accelerates computation; (iii) incorporates multiple QC measures designed to identify outlier samples and batch effects; and (iv) provides a variety of visualization options including chromosome, gene and exon-level views of CNVs, along with a tabular summarization of the exon-level CNVs. Compared with other CNV-calling algorithms using data from a lymphoma exome-seq study, PatternCNV has higher sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: The software for PatternCNV is implemented using Perl and R, and can be used in Mac or Linux environments. Software and user manual are available at http://bioinformaticstools.mayo.edu/research/patterncnv/, and R package at https://github.com/topsoil/patternCNV/.

Specification of human regional brain lineages using orthogonal gradients of WNT and SHH in organoids.

The repertory of neurons generated by progenitor cells depends on their location along antero-posterior and dorso-ventral axes of the neural tube. To understand if recreating those axes was sufficient to specify human brain neuronal diversity, we designed a mesofluidic device termed Duo-MAPS to expose induced pluripotent stem cells (iPSC) to concomitant orthogonal gradients of a posteriorizing and a ventralizing morphogen, activating WNT and SHH signaling, respectively. Comparison of single cell transcriptomes with fetal human brain revealed that Duo-MAPS-patterned organoids generated the major neuronal lineages of the forebrain, midbrain, and hindbrain. Morphogens crosstalk translated into early patterns of gene expression programs predicting the generation of specific brain lineages. Human iPSC lines from six different genetic backgrounds showed substantial differences in response to morphogens, suggesting that interindividual genomic and epigenomic variations could impact brain lineages formation. Morphogen gradients promise to be a key approach to model the brain in its entirety.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments.

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

Transcriptomic classification of diffuse large B-cell lymphoma identifies a high-risk activated B-cell-like subpopulation with targetable MYC dysregulation.

Immunochemotherapy has been the mainstay of treatment for newly diagnosed diffuse large B-cell lymphoma (ndDLBCL) yet is inadequate for many patients. In this work, we perform unsupervised clustering on transcriptomic features from a large cohort of ndDLBCL patients and identify seven clusters, one called A7 with poor prognosis, and develop a classifier to identify these clusters in independent ndDLBCL cohorts. This high-risk cluster is enriched for activated B-cell cell-of-origin, low immune infiltration, high MYC expression, and copy number aberrations. We compare and contrast our methodology with recent DLBCL classifiers to contextualize our clusters and show improved prognostic utility. Finally, using pre-clinical models, we demonstrate a mechanistic rationale for IKZF1/3 degraders such as lenalidomide to overcome the low immune infiltration phenotype of A7 by inducing T-cell trafficking into tumors and upregulating MHC I and II on tumor cells, and demonstrate that TCF4 is an important regulator of MYC-related biology in A7.

Multiomic analysis identifies a high-risk signature that predicts early clinical failure in DLBCL.

Recent genetic and molecular classification of DLBCL has advanced our knowledge of disease biology, yet were not designed to predict early events and guide anticipatory selection of novel therapies. To address this unmet need, we used an integrative multiomic approach to identify a signature at diagnosis that will identify DLBCL at high risk of early clinical failure. Tumor biopsies from 444 newly diagnosed DLBCL were analyzed by WES and RNAseq. A combination of weighted gene correlation network analysis and differential gene expression analysis was used to identify a signature associated with high risk of early clinical failure independent of IPI and COO. Further analysis revealed the signature was associated with metabolic reprogramming and identified cases with a depleted immune microenvironment. Finally, WES data was integrated into the signature and we found that inclusion of ARID1A mutations resulted in identification of 45% of cases with an early clinical failure which was validated in external DLBCL cohorts. This novel and integrative approach is the first to identify a signature at diagnosis, in a real-world cohort of DLBCL, that identifies patients at high risk for early clinical failure and may have significant implications for design of therapeutic options.

Molecular Clusters and Tumor-Immune Drivers of IgM Monoclonal Gammopathies.

PURPOSE: IgM monoclonal gammopathy of undetermined significance (MGUS) and Waldenström macroglobulinemia (WM) represent a disease spectrum with highly varied therapeutic management, ranging from observation to chemoimmunotherapy. The current classification relies solely on clinical features and does not explain the heterogeneity that exists within each of these conditions. Further investigation is warranted to shed light on the biology that may account for the clinical differences. EXPERIMENTAL DESIGN: We used bone marrow (BM) clonal CD19+ and/or CD138+ sorted cells, matched BM supernatant, and peripheral blood serum from 32 patients (7 MGUS, 25 WM) to perform the first multi-omics approach including whole-exome sequencing, RNA sequencing, proteomics, metabolomics, and mass cytometry. RESULTS: We identified three clusters with distinct pathway activation, immune content, metabolomic, and clinical features. Cluster 1 included only patients with WM and was characterized by transcriptional silencing of genes involved in cell cycle and immune response, enrichment of mitochondrial metabolism, infiltration of senescent T effector memory cells, and aggressive clinical behavior. Genetic/structural alterations of TNFAIP3 were distinct events of this cluster. Cluster 2 comprised both MGUS and WM patients with upregulation of inflammatory response, senescence and glycolysis signatures, increased activated T follicular helper and T regulatory cells, and indolent clinical behavior. Cluster 3 also included both MGUS and WM patients and exhibited intermediate features, including proliferative and inflammatory signaling, as well as glycolysis and mitochondrial metabolism. CONCLUSIONS: We have identified three distinct molecular clusters, suggesting a potential biologic classification that may have therapeutic implications.

Multiomic Analysis Identifies a High-Risk Metabolic and TME Depleted Signature that Predicts Early Clinical Failure in DLBCL.

PURPOSE: 60-70% of newly diagnosed diffuse large B-cell lymphoma (DLBCL) patients avoid events within 24 months of diagnosis (EFS24) and the remainder have poor outcomes. Recent genetic and molecular classification of DLBCL has advanced our knowledge of disease biology, yet were not designed to predict early events and guide anticipatory selection of novel therapies. To address this unmet need, we used an integrative multiomic approach to identify a signature at diagnosis that will identify DLBCL at high risk of early clinical failure. PATIENTS AND METHODS: Tumor biopsies from 444 newly diagnosed DLBCL were analyzed by WES and RNAseq. A combination of weighted gene correlation network analysis and differential gene expression analysis followed by integration with clinical and genomic data was used to identify a multiomic signature associated with high risk of early clinical failure. RESULTS: Current DLBCL classifiers are unable to discriminate cases who fail EFS24. We identified a high risk RNA signature that had a hazard ratio (HR, 18.46 [95% CI 6.51-52.31] P < .001) in a univariate model, which did not attenuate after adjustment for age, IPI and COO (HR, 20.8 [95% CI, 7.14-61.09] P < .001). Further analysis revealed the signature was associated with metabolic reprogramming and a depleted immune microenvironment. Finally, WES data was integrated into the signature and we found that inclusion of ARID1A mutations resulted in identification of 45% of cases with an early clinical failure which was validated in external DLBCL cohorts. CONCLUSION: This novel and integrative approach is the first to identify a signature at diagnosis that will identify DLBCL at high risk for early clinical failure and may have significant implications for design of therapeutic options.

Frequent POLE-driven hypermutation in ovarian endometrioid cancer revealed by mutational signatures in RNA sequencing.

BACKGROUND: DNA polymerase epsilon (POLE) is encoded by the POLE gene, and POLE-driven tumors are characterized by high mutational rates. POLE-driven tumors are relatively common in endometrial and colorectal cancer, and their presence is increasingly recognized in ovarian cancer (OC) of endometrioid type. POLE-driven cases possess an abundance of TCT > TAT and TCG > TTG somatic mutations characterized by mutational signature 10 from the Catalog of Somatic Mutations in Cancer (COSMIC). By quantifying the contribution of COSMIC mutational signature 10 in RNA sequencing (RNA-seq) we set out to identify POLE-driven tumors in a set of unselected Mayo Clinic OC. METHODS: Mutational profiles were calculated using expressed single-nucleotide variants (eSNV) in the Mayo Clinic OC tumors (n = 195), The Cancer Genome Atlas (TCGA) OC tumors (n = 419), and the Genotype-Tissue Expression (GTEx) normal ovarian tissues (n = 84). Non-negative Matrix Factorization (NMF) of the mutational profiles inferred the contribution per sample of four distinct mutational signatures, one of which corresponds to COSMIC mutational signature 10. RESULTS: In the Mayo Clinic OC cohort we identified six tumors with a predicted contribution from COSMIC mutational signature 10 of over five mutations per megabase. These six cases harbored known POLE hotspot mutations (P286R, S297F, V411L, and A456P) and were of endometrioid histotype (P = 5e-04). These six tumors had an early onset (average age of patients at onset, 48.33 years) when compared to non-POLE endometrioid OC cohort (average age at onset, 60.13 years; P = .008). Samples from TCGA and GTEx had a low COSMIC signature 10 contribution (median 0.16 mutations per megabase; maximum 1.78 mutations per megabase) and carried no POLE hotspot mutations. CONCLUSIONS: From the largest cohort of RNA-seq from endometrioid OC to date (n = 53), we identified six hypermutated samples likely driven by POLE (frequency, 11%). Our result suggests the clinical need to screen for POLE driver mutations in endometrioid OC, which can guide enrollment in immunotherapy clinical trials.

Early developmental asymmetries in cell lineage trees in living individuals.

Mosaic mutations can be used to track cell lineages in humans. We used cell cloning to analyze embryonic cell lineages in two living individuals and a postmortem human specimen. Of 10 reconstructed postzygotic divisions, none resulted in balanced contributions of daughter lineages to tissues. In both living individuals, one of two lineages from the first cleavage was dominant across tissues, with 90% frequency in blood. We propose that the efficiency of DNA repair contributes to lineage imbalance. Allocation of lineages in postmortem brain correlated with anterior-posterior axis, associating lineage history with cell fate choices in embryos. We establish a minimally invasive framework for defining cell lineages in any living individual, which paves the way for studying their relevance in health and disease.