Search

Sensitive PCR detection of viral nucleic acids plays a critical role in infectious disease research, diagnosis and monitoring. In the context of SARS-CoV-2 detection, recent reports indicate that digital PCR-based tests are significantly more sensitive than traditional qPCR tests. Numerous factors can influence digital PCR reaction sensitivity. In this review, using a model for human HIV infection and the Raindance ddPCR platform as an example, we describe technical aspects that contribute to sensitive viral signal detection in DNA and RNA from tissue samples, which often harbor viral reservoirs and serve as better predictors of disease outcome and indicators of treatment efficacy.

Failure of dermal protection or repair mechanisms might lead to visibly aged skin. The study aimed to identify genetic associations with perceived age. A genome-wide association study was undertaken in 423,992 adult participants of UK Biobank, using questionnaire data on perceived age and genetic data imputed to the Haplotype Reference Consortium imputation panel. Seventy-four to our knowledge previously unreported and independent associated genetic loci were identified (P

Graphical abstract, The effective non-invasive diagnosis and prognosis are critical for cancer treatment. The plasma cell-free DNA (cfDNA) provides a good material for cancer liquid biopsy and its worth in this field is increasingly explored. Here we describe a new pipeline for effectively finding new cfDNA-based biomarkers for cancers by combining SALP-seq and machine learning. Using the pipeline, 30 cfDNA samples from 26 esophageal cancer (ESCA) patients and 4 healthy people were analyzed as an example. As a result, 103 epigenetic markers (including 54 genome-wide and 49 promoter markers) and 37 genetic markers were identified for this cancer. These markers provide new biomarkers for ESCA diagnosis, prognosis and therapy. Importantly, these markers, especially epigenetic markers, not only shed important new insights on the regulatory mechanisms of this cancer, but also could be used to classify the cfDNA samples. We therefore developed a new pipeline for effectively finding new cfDNA-based biomarkers for cancers by combining SALP-seq and machine learning. In this study, we also discovered new clinical worth of cfDNA distinct from other reported characters.

Analysis of nucleotide variants is a cornerstone of cancer medicine. Although only 2% of the genomic sequence is protein coding, mutations occurring in these regions have the potential to influence protein structure or modification status and may have severe impact on disease aetiology. Proteogenomics enables the analysis of sample-specific nonsynonymous nucleotide variants with regard to their effect at the proteome and phosphoproteome levels. Here, we developed a proof-of-concept proteogenomics workflow and applied it to the malignant melanoma cell line A375. Initially, we studied the resistance to serine/threonine-protein kinase B-raf (BRAF) inhibitor (BRAFi) vemurafenib in A375 cells. This allowed identification of several oncogenic nonsynonymous nucleotide variants, including a gain-of-function variant on aurora kinase A (AURKA) at F31I. We also detected significant changes in abundance among (phospho)proteins, which led to reactivation of the MAPK signaling pathway in BRAFi-resistant A375 cells. Upon reconstruction of the multiomic integrated signaling networks, we predicted drug therapies with the potential to disrupt BRAFi resistance mechanism in A375 cells. Notably, we showed that AURKA inhibition is effective and specific against BRAFi-resistant A375 cells. Subsequently, we investigated amino acid variants that interfere with protein posttranslational modification (PTM) status and potentially influence A375 cell signaling irrespective of BRAFi resistance. Mass spectrometry (MS) measurements confirmed variant-driven PTM changes in 12 proteins. Among them was the runt-related transcription factor 1 (RUNX1) displaying a variant on a known phosphorylation site S(Ph)276L. We confirmed the loss of phosphorylation site by MS and demonstrated the impact of this variant on RUNX1 interactome., Graphical Abstract, Highlights • Proteogenomics workflow to study the malignant melanoma cell line A375. • Key signaling pathways are perturbed in BRAFi-sensitive and -resistant A375. • The signaling network in BRAFi-resistant cells can be targeted by several drugs. • Multiple amino acid variants directly affect protein phosphorylation status. • Loss of a phosphorylation site on RUNX1 alters its interactome., In Brief Proteogenomics is a powerful tool to study the mode of action of disease-associated mutations at the genome, transcriptome, proteome, and PTM level. Here, we applied a proteogenomics workflow to study the malignant melanoma cell line A375. Such workflow, used here as a proof of concept on A375 cells, may be applicable to other cancer types, cell lines, or even patient-derived samples.

Highlights • FGFR2 rearrangement frequency in Chinese cases is only 5.45%. • The first genomic aberration profiling of FGFRs in Chinese cases is presented. • Cases with FGFRs rearrangement have distinct clinical phenotype. • Patients without FGFR2 rearrangement may harbor other FGFR alterations., Genomic aberrations (GAs) in fibroblast growth factor receptors (FGFRs) are involved in the pathogenesis of intrahepatic cholangiocarcinoma (ICC), and clinical trials have shown efficacy of FGFR inhibitors in treating ICC patients with FGFR GAs such as FGFR2 rearrangement. To clarify the FGFRs GA profile and corresponding clinicopathological features in Chinese patients with ICC, a total of 257 cases were identified. Fourteen cases (5.45%) were positive for FGFR2 rearrangement. Further analysis on the 110 FGFR2 rearrangement negative cases showed that 13 patients present additional FGFRs GAs, including FGFR3 rearrangement (2.73%), and FGFRs mutations. When compared with patients without FGFRs GAs, those with FGFR2 or FGFR3 rearrangement presented more under the age of 58 years, female sex, HBsAb positivity, CD10 expression, and PD-L1 expression. The clinical characteristics between patients with FGFRs mutation and those without FGFRs GAs were similar, with the exception that cases with FGFRs mutation have more hepatolithiasis. We concluded that FGFR rearrangement is associated with unique clinical phenotypes in ICC.

Visual Abstract, Highlights • The homozygous c.38T>G mutation in the LEMD2 gene causes arrhythmic cardiomyopathy with bilateral juvenile cataract in the Hutterite population. • The cardiac phenotype is characterized by localized inferior and inferolateral fibrosis of the left ventricle and mild impairment of left ventricular systolic function but severe ventricular arrhythmias leading to sudden cardiac death. • Affected heart tissue and fibroblasts exhibit abnormally shaped nuclei with condensed peripheral heterochromatin. • Functional assays on affected fibroblasts show decreased proliferation capacity, cellular senescence, and a prolonged G1 phase, suggesting premature aging and cellular senescence in proliferating cells., Summary Nuclear envelope proteins have been shown to play an important role in the pathogenesis of inherited dilated cardiomyopathy. Here, we present a remarkable cardiac phenotype caused by a homozygous LEMD2 mutation in patients of the Hutterite population with juvenile cataract. Mutation carriers develop arrhythmic cardiomyopathy with mild impairment of left ventricular systolic function but severe ventricular arrhythmias leading to sudden cardiac death. Affected cardiac tissue from a deceased patient and fibroblasts exhibit elongated nuclei with abnormal condensed heterochromatin at the periphery. The patient fibroblasts demonstrate cellular senescence and reduced proliferation capacity, which may suggest an involvement of LEM domain containing protein 2 in chromatin remodeling processes and premature aging.

Health span is driven by a precise interplay between genes and the environment. Cell response to environmental cues is mediated by signaling cascades and genetic variants that affect gene expression by regulating chromatin plasticity. Indeed, they can promote the interaction of promoters with regulatory elements by forming active chromatin hubs. FOXO3 encodes a transcription factor with a strong impact on aging and age-related phenotypes, as it regulates stress response, therefore affecting lifespan. A significant association has been shown between human longevity and several FOXO3 variants located in intron 2. This haplotype block forms a putative aging chromatin hub in which FOXO3 has a central role, as it modulates the physical connection and activity of neighboring genes involved in age-related processes. Here we describe the role of FOXO3 and its single-nucleotide polymorphisms (SNPs) in healthy aging, with a focus on the enhancer region encompassing the SNP rs2802292, which upregulates FOXO3 expression and can promote the activity of the aging hub in response to different stress stimuli. FOXO3 protective effect on lifespan may be due to the accessibility of this region to transcription factors promoting its expression. This could in part explain the differences in FOXO3 association with longevity between genders, as its activity in females may be modulated by estrogens through estrogen receptor response elements located in the rs2802292-encompassing region. Altogether, the molecular mechanisms described here may help establish whether the rs2802292 SNP can be taken advantage of in predictive medicine and define the potential of targeting FOXO3 for age-related diseases., Graphical Abstract Unlabelled Image

Graphical abstract, Next-generation sequencing (NGS) allows the detection of mutations in inherited genetic diseases, like the Charcot-Marie-Tooth disease (CMT) which is the most common hereditary peripheral neuropathy. The majority of mutations detected by NGS are single nucleotide variants (SNVs) or small indels, while structural variants (SVs) are often underdiagnosed. PMP22 was the first gene described as being involved in CMT via a SV of duplication type. To date, more than 90 genes are known to be involved in CMT, with mainly SNVs and short indels described. Herein targeted NGS and the CovCopCan bioinformatic tool were used in two unrelated families, both presenting with typical CMT symptoms with pyramidal involvement. We have discovered two large SVs in KIF5A, a gene known to cause axonal forms of CMT (CMT2) in which no SVs have yet been described. In the first family, the patient presented with a large deletion of 12 kb in KIF5A from Chr12:57,956,278 to Chr12:57,968,335 including exons 2–15, that could lead to mutation c.(130-943_c.1717-533del), p.(Gly44_Leu572del). In the second family, two cases presented with a large deletion of 3 kb in KIF5A from Chr12:57,974,133 to Chr12:57,977,210 including exons 24–28, that could lead to mutation c.(2539-605_*36 + 211del), p.(Leu847_Ser1032delins33). In addition, bioinformatic sequence analysis revealed that a NAHR (Non-Allelic-Homologous-Recombination) mechanism, such as those in the PMP22 duplication, could be responsible for one of the KIF5A SVs and could potentially be present in a number of other patients. This study reveals that large KIF5A deletions can cause CMT2 and highlights the importance of analyzing not only the SNVs but also the SVs during diagnosis of neuropathies.

The ion pump Na+,K+-ATPase is a critical determinant of neuronal excitability; however, its role in the etiology of diseases of the central nervous system (CNS) is largely unknown. We describe here the molecular phenotype of a Trp931Arg mutation of the Na+,K+-ATPase catalytic α1 subunit in an infant diagnosed with therapy-resistant lethal epilepsy. In addition to the pathological CNS phenotype, we also detected renal wasting of Mg2+. We found that membrane expression of the mutant α1 protein was low, and ion pumping activity was lost. Arginine insertion into membrane proteins can generate water-filled pores in the plasma membrane, and our molecular dynamic (MD) simulations of the principle states of Na+, K+-ATPase transport demonstrated massive water inflow into mutant α1 and destabilization of the ion-binding sites. MD simulations also indicated that a water pathway was created between the mutant arginine residue and the cytoplasm, and analysis of oocytes expressing mutant α1 detected a nonspecific cation current. Finally, neurons expressing mutant α1 were observed to be depolarized compared with neurons expressing wild-type protein, compatible with a lowered threshold for epileptic seizures. The results imply that Na+,K+-ATPase should be considered a neuronal locus minoris resistentia in diseases associated with epilepsy and with loss of plasma membrane integrity.

In this study, we examined single nucleotide variants (SNVs) of the OPN3 gene in malignant melanoma and melanocytic nevi. A total of 20 variants of SNVs were detected. Of these variants, five nonsynonymous mutations of OPN3 were identified, including c.T152C, c.T401C, c.G547A, c.G768A, and c.G992A. Three prediction tools, MutationTaster2, Polymorphism Phenotyping version 2, and PROVEAN (Protein Variation Effect Analyzer), which predict possible impact of an amino acid substitution, suggested that the mutations could be deleterious. Nine SNVs occurred in 3′ untranslated regions, whereas two were observed in 5′ untranslated regions. In all cases, four intronic variants were identified. In addition, we identified nine 3′ untranslated region SNVs in OPN3; one of them (OPN3[NM_014322:c.∗83C>T]) is predicted to disrupt a conserved microRNA (has-miR-376c-3p) target site, located in position 86–93 of OPN3 3′ untranslated region. Our findings suggest that there is a strong possibility that OPN3 SNVs play a role in the pathogenesis of melanocytic tumors via prediction of functional phenotype.

Background The usefulness of circulating tumor DNA (ctDNA) in detecting mutations and monitoring treatment response has not been well studied beyond a few actionable biomarkers in non-small cell lung cancer (NSCLC). Research Question How does the usefulness of ctDNA analysis compare with that of solid tumor biopsy analysis in patients with NSCLC? Methods We retrospectively evaluated 370 adult patients with NSCLC treated at the City of Hope between November 2015 and August 2019 to assess the usefulness of ctDNA in mutation identification, survival, concordance with matched tissue samples in 32 genes, and tumor evolution. Results A total of 1,688 somatic mutations were detected in 473 ctDNA samples from 370 patients with NSCLC. Of the 473 samples, 177 showed at least one actionable mutation with currently available Food and Drug Administration-approved NSCLC therapies. MET and CDK6 amplifications co-occurred with BRAF amplifications (false discovery rate [FDR], < 0.01), and gene-level mutations were mutually exclusive in KRAS and EGFR (FDR, 0.0009). Low cumulative percent ctDNA levels were associated with longer progression-free survival (hazard ratio [HR], 0.56; 95% CI, 0.37-0.85; P = .006). Overall survival was shorter in patients harboring BRAF mutations (HR, 2.35; 95% CI, 1.24-4.6; P = .009), PIK3CA mutations (HR, 2.77; 95% CI, 1.56-4.9; P < .001) and KRAS mutations (HR, 2.32; 95% CI, 1.30-4.1; P = .004). Gene-level concordance was 93.8%, whereas the positive concordance rate was 41.6%. More mutations in targetable genes were found in ctDNA than in tissue biopsy samples. Treatment response and tumor evolution over time were detected in repeated ctDNA samples. Interpretation Although ctDNA analysis exhibited similar usefulness to tissue biopsy analysis, more mutations in targetable genes were missed in tissue biopsy analyses. Therefore, the evaluation of ctDNA in conjunction with tissue biopsy samples may help to detect additional targetable mutations to improve clinical outcomes in advanced NSCLC.

Highlights • The human genome contains approximately 4 million variants, whose population frequencies vary according to the ethnic backgrounds. • Genetic diversity of humans in part determines interindividual variability in susceptibility to diseases, response to therapy, and the clinical outcomes. • Genetic variants exert a gradient of biological and clinical effect sizes. In general, variants with the largest effect sizes are responsible for the single-gene disorders, whereas those with moderate and modest effect sizes are responsible for oligogenic and polygenic diseases, respectively. • A phenotype is the consequence of nonlinear stochastic interactions among multiple genetic and nongenetic determinants. • Discerning pathogenicity of the genetic variants, identified through genetic testing, in the clinical phenotype is challenging and requires complementary expertise in human molecular genetics and clinical medicine., Summary Genetic variants are major determinants of susceptibility to disease, response to therapy, and clinical outcomes. Advances in the short-read sequencing technologies, despite some shortcomings, have enabled identification of the vast majority of the genetic variants in each genome. The major challenge is in identifying the pathogenic variants in cardiovascular diseases. The yield of the genetic testing has been limited because of technological shortcomings and our incomplete understanding of the genetic basis of cardiovascular disorders. To advance the field, a shift to long-read sequencing platforms is necessary. In addition, to discern the pathogenic variants, genetic diseases should be considered as a continuum and the genetic variants as probabilistic factors with a gradient of effect sizes. Moreover, disease-specific physician-scientists with expertise in the clinical medicine and molecular genetics are best equipped to discern functional and clinical significance of the genetic variants. The changes would be expected to enhance clinical utilities of the genetic discoveries., Central Illustration

Background: The tumor microenvironment can be classified into immunologically active “inflamed” tumors and inactive “non-inflamed” tumors based on the infiltration of cytotoxic immune cells. Previous studies on liver cancer have reported a superior prognosis for inflamed tumors compared to non-inflamed tumors. However, liver cancer is highly heterogeneous immunologically and genetically, and a finer classification of the liver cancer microenvironment may improve our understanding of its immunological diversity and response to immune therapy. Methods: We characterized the immune gene signatures of 234 primary liver cancers, mainly virus-related, from a Japanese population using RNA-Seq of tumors and matched non-tumorous hepatitis livers. We then compared them with the somatic alterations detected using the whole-genome sequencing. Findings: Liver cancers expressed lower levels of immune marker genes than non-tumorous hepatitis livers, indicating immunosuppression in the tumor microenvironment. Several immunosuppression mechanisms functioned actively and mutually exclusively, resulting in four immune subclasses of liver cancer: tumor-associated macrophage (TAM), CTNNB1, cytolytic activity (CYT), and regulatory T cell (Treg). The CYT and Treg subclasses represented inflamed tumors, while the TAM and CTNNB1 subclasses represented non-inflamed tumors. The TAM subclass, which comprised 31% of liver cancers, showed a poor survival, expressed elevated levels of extracellular matrix genes, and was associated with somatic mutations of chromatin regulator ARID2. The results of cell line experiments suggested a functional link between ARID2 and chemokine production by liver cancer cells. Interpretation: Primary liver cancer was classified into four subclasses based on mutually exclusive mechanisms for immunosuppression. This classification indicate the importance of immunosuppression mechanisms, such as TAM and Treg, as therapeutic targets for liver cancer. Funding: The Japan Agency for Medical Research and Development (AMED). Keywords: Liver cancer, Tumor microenvironment, Tumor-associated macrophage, Regulatory T cell

Background & Aims RNase H2 is a holoenzyme, composed of 3 subunits (ribonuclease H2 subunits A, B, and C), that cleaves RNA:DNA hybrids and removes mis-incorporated ribonucleotides from genomic DNA through ribonucleotide excision repair. Ribonucleotide incorporation by eukaryotic DNA polymerases occurs during every round of genome duplication and produces the most frequent type of naturally occurring DNA lesion. We investigated whether intestinal epithelial proliferation requires RNase H2 function and whether RNase H2 activity is disrupted during intestinal carcinogenesis. Methods We generated mice with epithelial-specific deletion of ribonuclease H2 subunit B (H2bΔIEC) and mice that also had deletion of tumor-suppressor protein p53 (H2b/p53ΔIEC); we compared phenotypes with those of littermate H2bfl/fl or H2b/p53fl/fl (control) mice at young and old ages. Intestinal tissues were collected and analyzed by histology. We isolated epithelial cells, generated intestinal organoids, and performed RNA sequence analyses. Mutation signatures of spontaneous tumors from H2b/p53ΔIEC mice were characterized by exome sequencing. We collected colorectal tumor specimens from 467 patients, measured levels of ribonuclease H2 subunit B, and associated these with patient survival times and transcriptome data. Results The H2bΔIEC mice had DNA damage to intestinal epithelial cells and proliferative exhaustion of the intestinal stem cell compartment compared with controls and H2b/p53ΔIEC mice. However, H2b/p53ΔIEC mice spontaneously developed small intestine and colon carcinomas. DNA from these tumors contained T>G base substitutions at GTG trinucleotides. Analyses of transcriptomes of human colorectal tumors associated lower levels of RNase H2 with shorter survival times. Conclusions In analyses of mice with disruption of the ribonuclease H2 subunit B gene and colorectal tumors from patients, we provide evidence that RNase H2 functions as a colorectal tumor suppressor. H2b/p53ΔIEC mice can be used to study the roles of RNase H2 in tissue-specific carcinogenesis., Graphical abstract

Immune checkpoint inhibitors and adoptive lymphocyte transfer–based therapies have shown great therapeutic potential in cancers with high tumor mutational burden (TMB), such as melanoma, but not in cancers with low TMB, such as mutant epidermal growth factor receptor (EGFR)–driven lung adenocarcinoma. Precision immunotherapy is an unmet need for most cancers, particularly for cancers that respond inadequately to immune checkpoint inhibitors. Here, we employed large-scale MS-based proteogenomic profiling to identify potential immunogenic human leukocyte antigen (HLA) class I-presented peptides in melanoma and EGFR-mutant lung adenocarcinoma. Similar numbers of peptides were identified from both tumor types. Cell line and patient-specific databases (DBs) were constructed using variants identified from whole-exome sequencing. A de novo search algorithm was used to interrogate the HLA class I immunopeptidome MS data. We identified 12 variant peptides and several classes of tumor-associated antigen-derived peptides. We constructed a cancer germ line (CG) antigen DB with 285 antigens. This allowed us to identify 40 class I-presented CG antigen–derived peptides. The class I immunopeptidome comprised more than 1000 post-translationally modified (PTM) peptides representing 58 different PTMs, underscoring the critical role PTMs may play in HLA binding. Finally, leveraging de novo search algorithm and an annotated long noncoding RNA (lncRNA) DB, we developed a novel lncRNA-encoded peptide discovery pipeline to identify 44 lncRNA-derived peptides that are presented by class I. We validated tandem MS spectra of select variant, CG antigen, and lncRNA-derived peptides using synthetic peptides and performed HLA class I-binding assays to demonstrate binding to class I proteins. In summary, we provide direct evidence of HLA class I presentation of a large number of variant and tumor-associated peptides in both low and high TMB cancer. These results can potentially be useful for precision immunotherapies, such as vaccine or adoptive cell therapies in melanoma and EGFR-mutant lung cancers., Graphical Abstract, Highlights • Proteogenomics identified ∼35,000 class I-presented peptides. • CG antigen and PTM peptides identified in melanoma and lung cancer. • De novo search identified variant and lncRNA-derived peptides. • A new strategy to identify class I-presented lncRNA-derived peptides developed., In Brief Cancer immunotherapy is ineffective in low TMB EGFR-mutant lung adenocarcinoma. Qi et al. performed a comprehensive proteogenomic profiling of HLA class I-presented immunopeptides in high TMB melanoma and low TMB EGFR-mutant lung cancer. Similar numbers of immunopeptides were identified from both. Variant, CG antigen, PTM, and lncRNA-derived peptides were identified. A novel strategy to identify lncRNA-derived peptides was developed. The direct identification of class I-presented immunopeptides will potentially accelerate precision immunotherapy for low TMB tumors.

Mycobacterium tuberculosis (M. tuberculosis) is considered innately resistant to β-lactam antibiotics. However, there is evidence that susceptibility to β-lactam antibiotics in combination with β–lactamase inhibitors is variable among clinical isolates, and these may present therapeutic options for drug-resistant cases. Here we report our investigation of susceptibility to β-lactam/β–lactamase inhibitor combinations among clinical isolates of M. tuberculosis, and the use of comparative genomics to understand the observed heterogeneity in susceptibility. Eighty-nine South African clinical isolates of varying first and second-line drug susceptibility patterns and two reference strains of M. tuberculosis underwent minimum inhibitory concentration (MIC) determination to two β-lactams: amoxicillin and meropenem, both alone and in combination with clavulanate, a β–lactamase inhibitor. 41/91 (45%) of tested isolates were found to be hypersusceptible to amoxicillin/clavulanate relative to reference strains, including 14/24 (58%) of multiple drug-resistant (MDR) and 22/38 (58%) of extensively drug-resistant (XDR) isolates. Genome-wide polymorphisms identified using whole-genome sequencing were used in a phylogenetically-aware linear mixed model to identify polymorphisms associated with amoxicillin/clavulanate susceptibility. Susceptibility to amoxicillin/clavulanate was over-represented among isolates within a specific clade (LAM4), in particular among XDR strains. Twelve sets of polymorphisms were identified as putative markers of amoxicillin/clavulanate susceptibility, five of which were confined solely to LAM4. Within the LAM4 clade, ‘paradoxical hypersusceptibility’ to amoxicillin/clavulanate has evolved in parallel to first and second-line drug resistance. Given the high prevalence of LAM4 among XDR TB in South Africa, our data support an expanded role for β-lactam/β-lactamase inhibitor combinations for treatment of drug-resistant M. tuberculosis., Highlights • Paradoxical hypersusceptibility is observed drug susceptibility despite innate resistance in the wild type state. • Many MDR and XDR M. tuberculosis strains are susceptible to amoxicillin/clavulanate. • Whole-genome sequencing identified mutations associated with paradoxical hypersusceptibility. • An expanded role for β-lactams in drug-resistant M. tuberculosis is supported. The global increase in drug-resistant tuberculosis has prompted a search for alternative therapies, including repurposing existing antibiotics. β-lactam antibiotics are safe drugs, however, they have previously been thought to be of limited use for tuberculosis due to innate resistance to this drug class. In this study, the authors found many drug-resistant tuberculosis isolates from South Africa to be susceptible to a β-lactam and β-lactamase combination, amoxicillin/clavulanate. With the use of comparative genomics, multiple genetic mutations were identified to be associated with this hypersusceptible phenotype. These findings support an expanded role of β-lactam/β-lactamase inhibitor combinations for treatment of drug-resistant TB.

Background Neurofibromatosis type I (NF1) is caused by heterozygous loss-of-function variants in the NF1 gene encoding neurofibromin which serves as a tumor suppressor that inhibits RAS signaling and regulates cell proliferation and differentiation. While, the only well-established functional domain in the NF1 protein is the GAP-related domain (GRD), most of the identified non-truncating disease-causing variants are located outside of this domain, supporting the existence of other important disease-associated domains. Identifying these domains may reveal novel functions of NF1. Methods By implementing inferential statistics combined with machine-learning methods, we developed a novel NF1-specific functional prediction model that focuses on nonsynonymous single nucleotide variants (SNVs). The model enables annotating all possible NF1 nonsynonymous variants, thus mapping the range of pathogenic non-truncating variants at the codon level across the NF1 gene. Findings The generated model demonstrates high absolute prediction value for missense and splice-site variations (area under the ROC curve of 0.96) outperforming 14 other established models. By reviewing the entire dataset of nonsynonymous variants, two novel domains (Armadillo type fold 1 and 2) were identified as being associated with pathogenicity (OR 1.86; CI 1.04 to 3.34 and OR 2.08; CI 1.08 to 4.04, respectively; P

PURPOSE:To compare the efficacy of whole genome sequencing (WGS) with targeted next-generation sequencing (NGS) in the diagnosis of inherited retinal disease (IRD).DESIGN:Case series.PARTICIPANTS:A total of 562 patients diagnosed with IRD.METHODS:We performed a direct comparative analysis of current molecular diagnostics with WGS. We retrospectively reviewed the findings from a diagnostic NGS DNA test for 562 patients with IRD. A subset of 46 of 562 patients (encompassing potential clinical outcomes of diagnostic analysis) also underwent WGS, and we compared mutation detection rates and molecular diagnostic yields. In addition, we compared the sensitivity and specificity of the 2 techniques to identify known single nucleotide variants (SNVs) using 6 control samples with publically available genotype data.MAIN OUTCOME MEASURES:Diagnostic yield of genomic testing.RESULTS:Across known disease-causing genes, targeted NGS and WGS achieved similar levels of sensitivity and specificity for SNV detection. However, WGS also identified 14 clinically relevant genetic variants through WGS that had not been identified by NGS diagnostic testing for the 46 individuals with IRD. These variants included large deletions and variants in noncoding regions of the genome. Identification of these variants confirmed a molecular diagnosis of IRD for 11 of the 33 individuals referred for WGS who had not obtained a molecular diagnosis through targeted NGS testing. Weighted estimates, accounting for population structure, suggest that WGS methods could result in an overall 29% (95% confidence interval, 15-45) uplift in diagnostic yield.CONCLUSIONS:We show that WGS methods can detect disease-causing genetic variants missed by current NGS diagnostic methodologies for IRD and thereby demonstrate the clinical utility and additional value of WGS.Copyright © 2016 American Academy of Ophthalmology. Published by Elsevier Inc. All rights reserved.

Activation of innate immunity initiates various cascades of reactions that largely contribute to defense against physical, microbial or chemical damage, prompt for damage repair and removal of causative organisms as well as restoration of tissue homeostasis. Genetic polymorphism in innate immune genes plays prominent role in disease resistance capabilities in various breeds of cattle and buffalo. Here we studied single nucleotide variations (SNP/SNV) and haplotype structure in innate immune genes viz CHGA, CHGB, CHGC, NRAMP1, NRAMP2, DEFB1, BNBD4, BNBD5, TAP and LAP in Gir cattle and Murrah buffalo. Targeted sequencing of exonic regions of these genes was performed by Ion Torrent PGM sequencing platform. The sequence reads obtained corresponding to coding regions of these genes were mapped to reference genome of cattle BosTau7 by BWA program using genome analysis tool kit (GATK). Further variant analysis by Unified Genotyper revealed 54 and 224 SNPs in Gir and Murrah respectively and also 32 SNVs was identified. Among these SNPs 43, 36, 11,32,81,21 and 22 variations were in CHGA, CHGB, CHGC, NRAMP1, NRAMP2, DEFB1 and TAP genes respectively. Among these identified 278 SNPs, 24 were found to be reported in the dbSNP database. Variant analysis was followed by structure formation of haplotypes based on multiple SNPs using SAS software revealed a large number of haplotypes. The SNP discovery in innate immune genes in cattle and buffalo breeds of India would advance our understanding of role of these genes in determining the disease resistance/susceptibility in Indian breeds. The identified SNPs and haplotype data would also provide a wealth of sequence information for conservation studies, selective breeding and designing future strategies for identifying disease associations involving samples from distinct populations., Highlights • Sequencing of exonic regions of 10 selected innate immune genes, 895.60 Mb data was generated. • Variant analysis using GATK pipeline revealed 278 SNPs, 32 SNVs and from these 24 were found to be reported. • 43, 36 and 11 SNPs were observed in chromagranin genes viz. CHGA, CHGB and CHGC respectively. • 32 and 81 SNPs were identified in NRAMP genes viz. NRAMP1 and NRAMP2 respectively. • 21 and 22 SNPs were observed in β-defensin genes viz. DEFB1 and TAP respectively.

A key to improving cancer immunotherapy will be the identification of tumor-specific "neoantigens" that arise from mutations and augment the resultant host immune response. In this study we identified single nucleotide variants (SNVs) by RNA sequencing of asbestos-induced murine mesothelioma cell lines AB1 and AB1-HA. Using the NetMHCpan 2.8 algorithm, the theoretical binding affinity of predicted peptides arising from high-confidence, exonic, non-synonymous SNVs was determined for the BALB/c strain. The immunoreactivity to 20 candidate mutation-carrying peptides of increased affinity and the corresponding wild-type peptides was determined using interferon-γ ELISPOT assays and lymphoid organs of non-manipulated tumor-bearing mice. A strong endogenous immune response was demonstrated to one of the candidate neoantigens, Uqcrc2; this response was detected in the draining lymph node and spleen. Antigen reactive cells were not detected in non-tumor bearing mice. The magnitude of the response to the Uqcrc2 neoantigen was similar to that of the strong influenza hemagglutinin antigen, a model tumor neoantigen. This work confirms that the approach of RNAseq plus peptide prediction and ELISPOT testing is sufficient to identify natural tumor neoantigens.