ABSTRACT
Convenient, rapid, and accurate detection of nitroaromatic organic toxins and harmful substances is of great significance in research. In the present study, two-dimensional layered rare-earth hydroxides (LYH) were used as ion-exchange matrix materials, and the anionic fluorescent dye molecules (HPTS) were successfully introduced into the LYH structures in situ via a simple and effective "plug-and-play" strategy, which gave the compounds ultra-sensitive fluorescence sensing detection of nitrobenzene, p-nitrotoluene and p-nitrophenol (Fluorescence response time < 1 sec, and the LOD for nitrobenzene, p-nitrophenol and p-nitrotoluene reached an impressive 349 ppb, 22 ppb and 98 ppb, respectively). Combined with theoretical calculations, we elucidated in detail the fluorescence quenching response mechanism of the LYH-HPTS towards nitroaromatic. Additionally, we also constructed fluorescent paper sensor, which effectively transformed the LYH-HPTS from theoretical detection to device application. The LYH-HPTS material is not only simple to synthesize, cost-effective and stable, but also has the features of fast response, excellent sensitivity and selectivity, and good reproducibility, which provides a new approach for the rapid and accurate detection of nitroaromatic.
ABSTRACT
Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (â¼30 tissues × â¼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.
Subject(s)
Epigenome , Quantitative Trait Loci , Genome-Wide Association Study , Genomics , Phenotype , Polymorphism, Single NucleotideABSTRACT
BACKGROUND: Germline variants in fumarate hydratase (FH) are associated with autosomal dominant (AD) hereditary leiomyomatosis and renal cell cancer (HLRCC) and autosomal recessive (AR) fumarase deficiency (FMRD). The prevalence and cancer penetrance across different FH variants remain unclear. METHODS: A database containing 120,061 records from individuals undergoing cancer germline testing was obtained. FH variants were classified into 3 categories: AD HLRCC variants, AR FMRD variants, and variants of unknown significance (VUSs). Individuals with variants from these categories were compared with those with negative genetic testing. RESULTS: FH variants were detected in 1.3% of individuals (AD HLRCC, 0.3%; AR FMRD, 0.4%; VUS, 0.6%). The rate of AD HLRCC variants discovered among reportedly asymptomatic individuals without a clear indication for HLRCC testing was 1 in 2668 (0.04%). In comparison with those with negative genetic testing, the renal cell carcinoma (RCC) prevalence was elevated with AD HLRCC variants (17.0% vs 4.5%; P < .01) and VUSs (6.4% vs 4.5%; P = .02) but not with AR FMRD variants. CONCLUSIONS: The prevalence of HLRCC discovered incidentally on germline testing is similar to recent population carrier estimates, and this suggests that this is a relatively common cancer syndrome. Compared with those with negative genetic testing, those with VUSs had an elevated risk of RCC, whereas those with AR FMRD variants did not.
Subject(s)
Carcinoma, Renal Cell , Fumarate Hydratase , Kidney Neoplasms , Leiomyomatosis , Neoplastic Syndromes, Hereditary , Skin Neoplasms , Uterine Neoplasms , Carcinoma, Renal Cell/epidemiology , Carcinoma, Renal Cell/genetics , Female , Fumarate Hydratase/genetics , Germ Cells , Germ-Line Mutation , Humans , Kidney Neoplasms/epidemiology , Kidney Neoplasms/genetics , Leiomyomatosis/epidemiology , Leiomyomatosis/genetics , Leiomyomatosis/pathology , Neoplastic Syndromes, Hereditary/epidemiology , Neoplastic Syndromes, Hereditary/genetics , Neoplastic Syndromes, Hereditary/pathology , Prevalence , Skin Neoplasms/epidemiology , Skin Neoplasms/genetics , Skin Neoplasms/pathology , Uterine Neoplasms/epidemiology , Uterine Neoplasms/genetics , Uterine Neoplasms/pathologyABSTRACT
BACKGROUND: Inflammatory breast cancer (IBC) has a highly invasive and metastatic phenotype. However, little is known about its genetic drivers. To address this, we report the largest cohort of whole-genome sequencing (WGS) of IBC cases. METHODS: We performed WGS of 20 IBC samples and paired normal blood DNA to identify genomic alterations. For comparison, we used 23 matched non-IBC samples from the Cancer Genome Atlas Program (TCGA). We also validated our findings using WGS data from the International Cancer Genome Consortium (ICGC) and the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium. We examined a wide selection of genomic features to search for differences between IBC and conventional breast cancer. These include (i) somatic and germline single-nucleotide variants (SNVs), in both coding and non-coding regions; (ii) the mutational signature and the clonal architecture derived from these SNVs; (iii) copy number and structural variants (CNVs and SVs); and (iv) non-human sequence in the tumors (i.e., exogenous sequences of bacterial origin). RESULTS: Overall, IBC has similar genomic characteristics to non-IBC, including specific alterations, overall mutational load and signature, and tumor heterogeneity. In particular, we observed similar mutation frequencies between IBC and non-IBC, for each gene and most cancer-related pathways. Moreover, we found no exogenous sequences of infectious agents specific to IBC samples. Even though we could not find any strongly statistically distinguishing genomic features between the two groups, we did find some suggestive differences in IBC: (i) The MAST2 gene was more frequently mutated (20% IBC vs. 0% non-IBC). (ii) The TGF ß pathway was more frequently disrupted by germline SNVs (50% vs. 13%). (iii) Different copy number profiles were observed in several genomic regions harboring cancer genes. (iv) Complex SVs were more frequent. (v) The clonal architecture was simpler, suggesting more homogenous tumor-evolutionary lineages. CONCLUSIONS: Whole-genome sequencing of IBC manifests a similar genomic architecture to non-IBC. We found no unique genomic alterations shared in just IBCs; however, subtle genomic differences were observed including germline alterations in TGFß pathway genes and somatic mutations in the MAST2 kinase that could represent potential therapeutic targets.
Subject(s)
Genome, Human , Inflammatory Breast Neoplasms/genetics , Mutation/genetics , Whole Genome Sequencing , Clone Cells , DNA Copy Number Variations/genetics , Evolution, Molecular , Humans , Inflammatory Breast Neoplasms/microbiology , Inflammatory Breast Neoplasms/pathology , Molecular Sequence Annotation , Phenotype , Signal Transduction/geneticsABSTRACT
PURPOSE: Molecular profiling of cancer is increasingly common as part of routine care in oncology, and germline and somatic profiling may provide insights and actionable targets for men with metastatic prostate cancer. However, all reported cases are of deidentified individuals without full medical and genomic data available in the public domain. PATIENT AND METHODS: We present a case of whole-genome tumor and germline sequencing in a patient with advanced prostate cancer, who has agreed to make his genomic and clinical data publicly available. RESULTS: We describe an 84-year-old Caucasian male with a Gleason 10 oligometastastic hormone-sensitive prostate cancer. Whole-genome sequencing provided insights into his tumor's underlying mutational processes and the development of an SPOP mutation. It also revealed an androgen-receptor dependency of his cancer which was reflected in his durable response to radiation and hormonal therapy. Potentially actionable genomic lesions in the tumor were identified through a personalized medicine approach for potential future therapy, but at the moment, he remains in remission, illustrating the hormonal sensitivity of his SPOP-driven prostate cancer. We also placed this patient in the context of a large prostate-cancer cohort from the PCAWG (Pan-cancer Analysis of Whole Genomes) group. In this comparison, the patient's cancer appears typical in terms of the number and type of somatic mutations, but it has a somewhat larger contribution from the mutational process associated with aging. CONCLUSION: We combined the expertise of medical oncology and genomics approaches to develop a molecular tumor board to integrate the care and study of this patient, who continues to have an outstanding response to his combined modality treatment. This identifiable case potentially helps overcome barriers to clinical and genomic data sharing.
Subject(s)
Biomarkers, Tumor/genetics , Genome, Human , Genomics/methods , Precision Medicine , Prostatic Neoplasms/pathology , Whole Genome Sequencing/methods , Aged, 80 and over , Humans , Male , Molecular Medicine , Prognosis , Prostatic Neoplasms/geneticsABSTRACT
Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20128-w.
ABSTRACT
The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) curated consensus somatic mutation calls using whole exome sequencing (WES) and whole genome sequencing (WGS), respectively. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2,658 cancers across 38 tumour types, we compare WES and WGS side-by-side from 746 TCGA samples, finding that ~80% of mutations overlap in covered exonic regions. We estimate that low variant allele fraction (VAF < 15%) and clonal heterogeneity contribute up to 68% of private WGS mutations and 71% of private WES mutations. We observe that ~30% of private WGS mutations trace to mutations identified by a single variant caller in WES consensus efforts. WGS captures both ~50% more variation in exonic regions and un-observed mutations in loci with variable GC-content. Together, our analysis highlights technological divergences between two reproducible somatic variant detection efforts.
Subject(s)
Genome, Human/genetics , Mutation , Neoplasms/genetics , Base Composition , DNA, Intergenic , Databases, Genetic , Exome/genetics , Exons , Humans , Retrospective Studies , Exome Sequencing , Whole Genome SequencingABSTRACT
Proteolysis is a major posttranslational regulator of biology inside and outside of cells. Broad identification of optimal cleavage sites and natural substrates of proteases is critical for drug discovery and to understand protease biology. Here, we present a method that employs two genetically encoded substrate phage display libraries coupled with next generation sequencing (SPD-NGS) that allows up to 10,000-fold deeper sequence coverage of the typical six- to eight-residue protease cleavage sites compared to state-of-the-art synthetic peptide libraries or proteomics. We applied SPD-NGS to two classes of proteases, the intracellular caspases, and the ectodomains of the sheddases, ADAMs 10 and 17. The first library (Lib 10AA) allowed us to identify 104 to 105 unique cleavage sites over a 1,000-fold dynamic range of NGS counts and produced consensus and optimal cleavage motifs based position-specific scoring matrices. A second SPD-NGS library (Lib hP), which displayed virtually the entire human proteome tiled in contiguous 49 amino acid sequences with 25 amino acid overlaps, enabled us to identify candidate human proteome sequences. We identified up to 104 natural linear cut sites, depending on the protease, and captured most of the examples previously identified by proteomics and predicted 10- to 100-fold more. Structural bioinformatics was used to facilitate the identification of candidate natural protein substrates. SPD-NGS is rapid, reproducible, simple to perform and analyze, inexpensive, and renewable, with unprecedented depth of coverage for substrate sequences, and is an important tool for protease biologists interested in protease specificity for specific assays and inhibitors and to facilitate identification of natural protein substrates.
Subject(s)
Caspase 3/metabolism , Proteome , Caspase 3/genetics , Gene Expression Regulation, Enzymologic , Humans , Peptide Library , Substrate SpecificityABSTRACT
A central issue in drug risk-benefit assessment is identifying frequencies of side effects in humans. Currently, frequencies are experimentally determined in randomised controlled clinical trials. We present a machine learning framework for computationally predicting frequencies of drug side effects. Our matrix decomposition algorithm learns latent signatures of drugs and side effects that are both reproducible and biologically interpretable. We show the usefulness of our approach on 759 structurally and therapeutically diverse drugs and 994 side effects from all human physiological systems. Our approach can be applied to any drug for which a small number of side effect frequencies have been identified, in order to predict the frequencies of further, yet unidentified, side effects. We show that our model is informative of the biology underlying drug activity: individual components of the drug signatures are related to the distinct anatomical categories of the drugs and to the specific drug routes of administration.
Subject(s)
Drug-Related Side Effects and Adverse Reactions , Machine Learning , Algorithms , Computational Biology/methods , Databases, Pharmaceutical , Humans , Pharmaceutical Preparations/administration & dosage , ProbabilityABSTRACT
Multiple mutational processes drive carcinogenesis, leaving characteristic signatures in tumor genomes. Determining the active signatures from a full repertoire of potential ones helps elucidate mechanisms of cancer development. This involves optimally decomposing the counts of cancer mutations, tabulated according to their trinucleotide context, into a linear combination of known signatures. Here, we develop sigLASSO (a software tool at github.com/gersteinlab/siglasso) to carry out this optimization efficiently. sigLASSO has four key aspects: (1) It jointly optimizes the likelihood of sampling and signature fitting, by explicitly factoring multinomial sampling into the objective function. This is particularly important when mutation counts are low and sampling variance is high (e.g., in exome sequencing). (2) sigLASSO uses L1 regularization to parsimoniously assign signatures, leading to sparse and interpretable solutions. (3) It fine-tunes model complexity, informed by data scale and biological priors. (4) Consequently, sigLASSO can assess model uncertainty and abstain from making assignments in low-confidence contexts.
Subject(s)
Computational Biology/methods , Mutation , Neoplasms/genetics , Genome, Human , Humans , Probability , SoftwareABSTRACT
ENCODE comprises thousands of functional genomics datasets, and the encyclopedia covers hundreds of cell types, providing a universal annotation for genome interpretation. However, for particular applications, it may be advantageous to use a customized annotation. Here, we develop such a custom annotation by leveraging advanced assays, such as eCLIP, Hi-C, and whole-genome STARR-seq on a number of data-rich ENCODE cell types. A key aspect of this annotation is comprehensive and experimentally derived networks of both transcription factors and RNA-binding proteins (TFs and RBPs). Cancer, a disease of system-wide dysregulation, is an ideal application for such a network-based annotation. Specifically, for cancer-associated cell types, we put regulators into hierarchies and measure their network change (rewiring) during oncogenesis. We also extensively survey TF-RBP crosstalk, highlighting how SUB1, a previously uncharacterized RBP, drives aberrant tumor expression and amplifies the effect of MYC, a well-known oncogenic TF. Furthermore, we show how our annotation allows us to place oncogenic transformations in the context of a broad cell space; here, many normal-to-tumor transitions move towards a stem-like state, while oncogene knockdowns show an opposing trend. Finally, we organize the resource into a coherent workflow to prioritize key elements and variants, in addition to regulators. We showcase the application of this prioritization to somatic burdening, cancer differential expression and GWAS. Targeted validations of the prioritized regulators, elements and variants using siRNA knockdowns, CRISPR-based editing, and luciferase assays demonstrate the value of the ENCODE resource.
Subject(s)
Databases, Genetic , Genomics , Neoplasms/genetics , Cell Line, Tumor , Cell Transformation, Neoplastic/genetics , Gene Regulatory Networks , Humans , Mutation/genetics , Reproducibility of Results , Transcription Factors/metabolismABSTRACT
BACKGROUND: Hereditary leiomyomatosis and renal cancer (HLRCC) is a cancer syndrome associated with a germline mutation in fumarate hydratase (FH). The syndrome is associated with cutaneous and uterine leiomyomas, and some patients develop a lethal form of kidney cancer. This study provides estimates for the FH carrier frequency and kidney cancer penetrance. METHODS: Data sets containing sequencing data for the FH gene were used: the 1000 Genomes Project (1000GP) and the Exome Aggregation Consortium (ExAC). Alterations in the FH gene were characterized on the basis of different variant risk tiers: 1) ClinVar annotated variants, 2) loss-of-function alterations, and 3) highly impactful missense alterations. The cumulative incidence of FH alterations overall and by different world populations was evaluated in 1000GP and ExAC. A lifetime penetrance of HLRCC kidney cancer risk was generated with 3 estimates of the annual incidence. RESULTS: The overall allele frequencies of tier 1 to 3 FH alterations in the ExAC and 1000GP data sets were 2.54 × 10-3 (1 in 393) and 1.20 × 10-3 (1 in 835), respectively. There were differences in the allele frequencies of FH alterations between world populations. Based on various estimates of the percentage of kidney cancers with FH alterations, the lifetime kidney cancer penetrance for carrier estimate 3 in ExAC was 1.7% to 5.8%. CONCLUSIONS: FH alterations are common and are carried by approximately 1 in 1000 individuals according to the more conservative estimates. The lifetime kidney cancer penetrance appears lower than previously estimated. Although databases are not population cohorts, they provide a useful quantitative estimate of rare variants with low penetrance.
Subject(s)
Fumarate Hydratase/genetics , Genetic Predisposition to Disease , Kidney Neoplasms/genetics , Leiomyomatosis/genetics , Neoplastic Syndromes, Hereditary/genetics , Skin Neoplasms/genetics , Uterine Neoplasms/genetics , Adult , Exome/genetics , Female , Gene Frequency , Germ-Line Mutation/genetics , Heterozygote , Humans , Kidney/metabolism , Kidney/pathology , Kidney Neoplasms/epidemiology , Kidney Neoplasms/etiology , Kidney Neoplasms/pathology , Leiomyomatosis/complications , Leiomyomatosis/epidemiology , Leiomyomatosis/pathology , Male , Middle Aged , Neoplastic Syndromes, Hereditary/complications , Neoplastic Syndromes, Hereditary/epidemiology , Neoplastic Syndromes, Hereditary/pathology , Risk Factors , Skin Neoplasms/complications , Skin Neoplasms/epidemiology , Skin Neoplasms/pathology , Uterine Neoplasms/complications , Uterine Neoplasms/epidemiology , Uterine Neoplasms/pathologyABSTRACT
The dichotomous model of "drivers" and "passengers" in cancer posits that only a few mutations in a tumor strongly affect its progression, with the remaining ones being inconsequential. Here, we leveraged the comprehensive variant dataset from the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) project to demonstrate that-in addition to the dichotomy of high- and low-impact variants-there is a third group of medium-impact putative passengers. Moreover, we also found that molecular impact correlates with subclonal architecture (i.e., early versus late mutations), and different signatures encode for mutations with divergent impact. Furthermore, we adapted an additive-effects model from complex-trait studies to show that the aggregated effect of putative passengers, including undetected weak drivers, provides significant additional power (â¼12% additive variance) for predicting cancerous phenotypes, beyond PCAWG-identified driver mutations. Finally, this framework allowed us to estimate the frequency of potential weak-driver mutations in PCAWG samples lacking any well-characterized driver alterations.
Subject(s)
Genome, Human/genetics , Genomics/methods , Mutation/genetics , Neoplasms/genetics , DNA Mutational Analysis/methods , Disease Progression , Humans , Neoplasms/pathology , Whole Genome SequencingABSTRACT
Data science allows the extraction of practical insights from large-scale data. Here, we contextualize it as an umbrella term, encompassing several disparate subdomains. We focus on how genomics fits as a specific application subdomain, in terms of well-known 3 V data and 4 M process frameworks (volume-velocity-variety and measurement-mining-modeling-manipulation, respectively). We further analyze the technical and cultural "exports" and "imports" between genomics and other data-science subdomains (e.g., astronomy). Finally, we discuss how data value, privacy, and ownership are pressing issues for data science applications, in general, and are especially relevant to genomics, due to the persistent nature of DNA.
Subject(s)
Data Science , GenomicsABSTRACT
Retroduplications come from reverse transcription of mRNAs and their insertion back into the genome. Here, we performed comprehensive discovery and analysis of retroduplications in a large cohort of 2,535 individuals from 26 human populations, as part of 1000 Genomes Phase 3. We developed an integrated approach to discover novel retroduplications combining high-coverage exome and low-coverage whole-genome sequencing data, utilizing information from both exon-exon junctions and discordant paired-end reads. We found 503 parent genes having novel retroduplications absent from the reference genome. Based solely on retroduplication variation, we built phylogenetic trees of human populations; these represent superpopulation structure well and indicate that variable retroduplications are effective population markers. We further identified 43 retroduplication parent genes differentiating superpopulations. This group contains several interesting insertion events, including a SLMO2 retroduplication and insertion into CAV3, which has a potential disease association. We also found retroduplications to be associated with a variety of genomic features: (1) Insertion sites were correlated with regular nucleosome positioning. (2) They, predictably, tend to avoid conserved functional regions, such as exons, but, somewhat surprisingly, also avoid introns. (3) Retroduplications tend to be co-inserted with young L1 elements, indicating recent retrotranspositional activity, and (4) they have a weak tendency to originate from highly expressed parent genes. Our investigation provides insight into the functional impact and association with genomic elements of retroduplications. We anticipate our approach and analytical methodology to have application in a more clinical context, where exome sequencing data is abundant and the discovery of retroduplications can potentially improve the accuracy of SNP calling.
Subject(s)
Gene Duplication/genetics , Genetic Variation/genetics , Genome, Human/genetics , High-Throughput Nucleotide Sequencing/methods , Retroelements/genetics , Sequence Analysis, DNA/methods , DNA Transposable Elements/genetics , Exome/genetics , Humans , Species SpecificityABSTRACT
To date, studies on papillary renal-cell carcinoma (pRCC) have largely focused on coding alterations in traditional drivers, particularly the tyrosine-kinase, Met. However, for a significant fraction of tumors, researchers have been unable to determine a clear molecular etiology. To address this, we perform the first whole-genome analysis of pRCC. Elaborating on previous results on MET, we find a germline SNP (rs11762213) in this gene predicting prognosis. Surprisingly, we detect no enrichment for small structural variants disrupting MET. Next, we scrutinize noncoding mutations, discovering potentially impactful ones associated with MET. Many of these are in an intron connected to a known, oncogenic alternative-splicing event; moreover, we find methylation dysregulation nearby, leading to a cryptic promoter activation. We also notice an elevation of mutations in the long noncoding RNA NEAT1, and these mutations are associated with increased expression and unfavorable outcome. Finally, to address the origin of pRCC heterogeneity, we carry out whole-genome analyses of mutational processes. First, we investigate genome-wide mutational patterns, finding they are governed mostly by methylation-associated C-to-T transitions. We also observe significantly more mutations in open chromatin and early-replicating regions in tumors with chromatin-modifier alterations. Finally, we reconstruct cancer-evolutionary trees, which have markedly different topologies and suggested evolutionary trajectories for the different subtypes of pRCC.
Subject(s)
Carcinoma, Renal Cell/genetics , Kidney Neoplasms/genetics , Proto-Oncogene Proteins c-met/genetics , RNA, Long Noncoding/genetics , Carcinoma, Renal Cell/pathology , Chromatin/genetics , DNA Methylation/genetics , Female , Gene Expression Regulation, Neoplastic , Genome, Human , Humans , Introns/genetics , Kidney Neoplasms/pathology , Male , Middle Aged , Mutation , Polymorphism, Single NucleotideABSTRACT
The rapidly growing volume of data being produced by next-generation sequencing initiatives is enabling more in-depth analyses of conservation than previously possible. Deep sequencing is uncovering disease loci and regions under selective constraint, despite the fact that intuitive biophysical reasons for such constraint are sometimes absent. Allostery may often provide the missing explanatory link. We use models of protein conformational change to identify allosteric residues by finding essential surface pockets and information-flow bottlenecks, and we develop a software tool that enables users to perform this analysis on their own proteins of interest. Though fundamentally 3D-structural in nature, our analysis is computationally fast, thereby allowing us to run it across the PDB and to evaluate general properties of predicted allosteric residues. We find that these tend to be conserved over diverse evolutionary time scales. Finally, we highlight examples of allosteric residues that help explain poorly understood disease-associated variants.
Subject(s)
Allosteric Site , Conserved Sequence , Sequence Analysis, Protein/methods , Software , Animals , Humans , Proteome/chemistryABSTRACT
As the cost of sequencing continues to decrease and the amount of sequence data generated grows, new paradigms for data storage and analysis are increasingly important. The relative scaling behavior of these evolving technologies will impact genomics research moving forward.