YBX1-interacting small RNAs and RUNX2 can be blocked in primary bone cancer using CADD522.

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGlyTCC cleavage into 5' end tRF-GlyTCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-GlyTCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-GlyTCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-GlyTCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

Somatic variant calling from single-cell DNA sequencing data.

Single-cell sequencing has gained popularity in recent years. Despite its numerous applications, single-cell DNA sequencing data is highly error-prone due to technical biases arising from uneven sequencing coverage, allelic dropout, and amplification error. With these artifacts, the identification of somatic genomic variants becomes a challenging task, and over the years, several methods have been developed explicitly for this type of data. Single-cell variant callers implement distinct strategies, make different use of the data, and typically result in many discordant calls when applied to real data. Here, we review current approaches for single-cell variant calling, emphasizing single nucleotide variants. We highlight their potential benefits and shortcomings to help users choose a suitable tool for their data at hand.

Characterization of Single Nucleotide Variants of OPN3 Gene in Melanocytic Nevi and Melanoma.

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.

Polygenic Risk Scores for Atherosclerotic Cardiovascular Disease in the Asia-Pacific Region.

Approximately one-half of the phenotypic susceptibility to atherosclerotic cardiovascular disease (ASCVD) has a genetic basis. Although individual allelic variants generally impart a small effect on risk for ASCVD, an emerging body of data has shown that the aggregation and weighting of many of these genetic variations into "scores" can further discriminate an individual's risk beyond traditional risk factors alone. Consistent with the theory of population genetics, such polygenic risk scores (PRS) appear to be ethnicity specific because their elements comprise single-nucleotide variants that are always ethnicity specific. The currently available PRS are derived predominantly from European ancestry and thus predictably perform less well among non-European participants, a fact that has implications for their use in the Asia-Pacific region. This paper describes the current state of knowledge of PRS, the available data that support their use in this region, and highlights the needs moving forward to safely and effectively implement them in clinical care in the Asia-Pacific region.

New structural variations responsible for Charcot-Marie-Tooth disease: The first two large KIF5A deletions detected by CovCopCan software.

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.

Identification of GPC3 mutation and upregulation in a multidrug resistant osteosarcoma and its spheroids as therapeutic target.

BACKGROUND: Drug resistance and the lack of molecular therapeutic target are the main challenges in the management of osteosarcomas (OSs). Identification of novel genetic alteration(s) related with OS recurrence and chemotherapeutic resistance would be of scientific and clinical significance. METHODS: To identify potential genetic alterations related with OS recurrence and chemotherapeutic resistance, the biopsies of a 20-year-old male osteosarcoma patient were collected at primary site (p-OS) and from its metastatic tumor (m-OS) formed after 5 months of adjuvant chemotherapy. Both OS specimens were subjected to cancer-targeted next generation sequencing (NGS) and their cell suspensions were cultured under three-dimensional condition to establish spheroid therapeutic model. Transcript-oriented Sanger sequencing for GPC3, the detected mutated gene, was performed on RNA samples of p-OS and m-OS tissues and spheroids. The effects of anti-GPC3 antibody and its combination with cisplatin on m-OS spheroids were elucidated. RESULTS: NGS revealed 4 mutations (GPC3, SOX10, MDM4 and MAPK8) and 6 amplifications (MDM2, CDK4, CCND3, RUNX2, GLI1 and FRS2) in p-OS, and 3 mutations (GPC3, SOX10 and EGF) and 10 amplifications (CDK4, CCND3, MDM2, RUNX2, GLI1, FRS2, CARD11, RAC1, SLC16A7 and PMS2) in m-OS. Among those alterations, the mutation abundance of GPC3 was the highest (56.49%) in p-OS and showed 1.54 times increase in m-OS. GPC3 transcript-oriented Sanger sequencing confirmed the mutation at 1046 in Exon 4, and immunohistochemical staining showed increased GPC3 production in m-OS tissues and its spheroids. EdU cell proliferation and Calcein/PI cell viability assays revealed that of the anti-OS first line drugs (doxorubicin, cisplatin, methotrexate, ifosfamide and carboplatin), 10 µM carboplatin exerted the best inhibitory effects on the p-OS but not the m-OS spheroids. 2 µg/mL anti-GPC3 antibody effectively committed m-OS spheroids to death by itself (76.43%) or in combination with cisplatin (92.93%). CONCLUSION: This study demonstrates increased abundance and up-regulated expression of mutant GPC3 in metastatic osteosarcoma and its spheroids with multidrug resistance. As GPC3-targeting therapy has been used to treat hepatocellular carcinomas and it is also effective to OS PDSs, GPC3 would be a novel prognostic parameter and therapeutic target of osteosarcomas.

The molecular feature of macrophages in tumor immune microenvironment of glioma patients.

BACKGROUND: Gliomas are one of the most common types of primary tumors in central nervous system. Previous studies have found that macrophages actively participate in tumor growth. METHODS: Weighted gene co-expression network analysis was used to identify meaningful macrophage-related gene genes for clustering. Pamr, SVM, and neural network were applied for validating clustering results. Somatic mutation and methylation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of MScores. The expression of macrophage-specific genes were evaluated in tumor microenvironment based on single cell sequencing analysis. A total of 2365 samples from 15 glioma datasets and 5842 pan-cancer samples were used for external validation of MScore. RESULTS: Macrophages were identified to be negatively associated with the survival of glioma patients. Twenty-six macrophage-specific DEGs obtained by elastic regression and PCA were highly expressed in macrophages at single-cell level. The prognostic value of MScores in glioma was validated by the active proinflammatory and metabolic profile of infiltrating microenvironment and response to immunotherapies of samples with this signature. MScores managed to stratify patient survival probabilities in 15 external glioma datasets and pan-cancer datasets, which predicted worse survival outcome. Sequencing data and immunohistochemistry of Xiangya glioma cohort confirmed the prognostic value of MScores. A prognostic model based on MScores demonstrated high accuracy rate. CONCLUSION: Our findings strongly support a modulatory role of macrophages, especially M2 macrophages in glioma progression and warrants further experimental studies.

Clinical Interpretation and Management of Genetic Variants.

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.

Machine Learning Identifies Clinical and Genetic Factors Associated With Anthracycline Cardiotoxicity in Pediatric Cancer Survivors.

BACKGROUND: Despite known clinical risk factors, predicting anthracycline cardiotoxicity remains challenging. OBJECTIVES: This study sought to develop a clinical and genetic risk prediction model for anthracycline cardiotoxicity in childhood cancer survivors. METHODS: We performed exome sequencing in 289 childhood cancer survivors at least 3 years from anthracycline exposure. In a nested case-control design, 183 case patients with reduced left ventricular ejection fraction despite low-dose doxorubicin (≤250 mg/m2), and 106 control patients with preserved left ventricular ejection fraction despite doxorubicin >250 mg/m2 were selected as extreme phenotypes. Rare/low-frequency variants were collapsed to identify genes differentially enriched for variants between case patients and control patients. The expression levels of 5 top-ranked genes were evaluated in human induced pluripotent stem cell-derived cardiomyocytes, and variant enrichment was confirmed in a replication cohort. Using random forest, a risk prediction model that included genetic and clinical predictors was developed. RESULTS: Thirty-one genes were differentially enriched for variants between case patients and control patients (p < 0.001). Only 42.6% case patients harbored a variant in these genes compared to 89.6% control patients (odds ratio: 0.09; 95% confidence interval: 0.04 to 0.17; p = 3.98 × 10-15). A risk prediction model for cardiotoxicity that included clinical and genetic factors had a higher prediction accuracy and lower misclassification rate compared to the clinical-only model. In vitro inhibition of gene-associated pathways (PI3KR2, ZNF827) provided protection from cardiotoxicity in cardiomyocytes. CONCLUSIONS: Our study identified variants in cardiac injury pathway genes that protect against cardiotoxicity and informed the development of a prediction model for delayed anthracycline cardiotoxicity, and it also provided new targets in autophagy genes for the development of cardio-protective drugs. (Preventing Cardiac Sequelae in Pediatric Cancer Survivors [PCS2]; NCT01805778).

Finding new cancer epigenetic and genetic biomarkers from cell-free DNA by combining SALP-seq and machine learning.

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.

FOXO3 on the Road to Longevity: Lessons From SNPs and Chromatin Hubs.

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.

Characterization of a Unique Form of Arrhythmic Cardiomyopathy Caused by Recessive Mutation in LEMD2.

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.

Strong spontaneous tumor neoantigen responses induced by a natural human carcinogen.

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.

Discovery of Genetic Variants of the Kinases That Activate Tenofovir in a Compartment-specific Manner.

Tenofovir (TFV) is used in combination with other antiretroviral drugs for human immunodeficiency virus (HIV) treatment and prevention. TFV requires two phosphorylation steps to become pharmacologically active; however, the kinases that activate TFV in cells and tissues susceptible to HIV infection have yet to be identified. Peripheral blood mononuclear cells (PBMC), vaginal, and colorectal tissues were transfected with siRNA targeting nucleotide kinases, incubated with TFV, and TFV-monophosphate (TFV-MP) and TFV-diphosphate (TFV-DP) were measured using mass spectrometry-liquid chromatography. Adenylate kinase 2 (AK2) performed the first TFV phosphorylation step in PBMC, vaginal, and colorectal tissues. Interestingly, both pyruvate kinase isozymes, muscle (PKM) or liver and red blood cell (PKLR), were able to phosphorylate TFV-MP to TFV-DP in PBMC and vaginal tissue, while creatine kinase, muscle (CKM) catalyzed this conversion in colorectal tissue. In addition, next-generation sequencing of the Microbicide Trials Network MTN-001 clinical samples detected 71 previously unreported genetic variants in the genes encoding these kinases. In conclusion, our results demonstrate that TFV is activated in a compartment-specific manner. Further, genetic variants have been identified that could negatively impact TFV activation, thereby compromising TFV efficacy in HIV treatment and prevention.

Exploring genetic polymorphism in innate immune genes in Indian cattle (Bos indicus) and buffalo (Bubalus bubalis) using next generation sequencing technology.

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.

Decoding the Molecular and Mutational Ambiguities of Gastroenteropancreatic Neuroendocrine Neoplasm Pathobiology.

Gastroenteropancreatic neuroendocrine neoplasms (GEP-NEN), considered a heterogeneous neoplasia, exhibit ill-defined pathobiology and protean symptomatology and are ubiquitous in location. They are difficult to diagnose, challenging to manage, and outcome depends on cell type, secretory product, histopathologic grading, and organ of origin. A morphologic and molecular genomic review of these lesions highlights tumor characteristics that can be used clinically, such as somatostatin-receptor expression, and confirms features that set them outside the standard neoplasia paradigm. Their unique pathobiology is useful for developing diagnostics using somatostatin-receptor targeted imaging or uptake of radiolabeled amino acids specific to secretory products or metabolism. Therapy has evolved via targeting of protein kinase B signaling or somatostatin receptors with drugs or isotopes (peptide-receptor radiotherapy). With DNA sequencing, rarely identified activating mutations confirm that tumor suppressor genes are relevant. Genomic approaches focusing on cancer-associated genes and signaling pathways likely will remain uninformative. Their uniquely dissimilar molecular profiles mean individual tumors are unlikely to be easily or uniformly targeted by therapeutics currently linked to standard cancer genetic paradigms. The prevalence of menin mutations in pancreatic NEN and P27KIP1 mutations in small intestinal NEN represents initial steps to identifying a regulatory commonality in GEP-NEN. Transcriptional profiling and network-based analyses may define the cellular toolkit. Multianalyte diagnostic tools facilitate more accurate molecular pathologic delineations of NEN for assessing prognosis and identifying strategies for individualized patient treatment. GEP-NEN remain unique, poorly understood entities, and insight into their pathobiology and molecular mechanisms of growth and metastasis will help identify the diagnostic and therapeutic weaknesses of this neoplasia.