Parafibromin Abnormalities in Ossifying Fibroma.

Ossifying fibromas are very rare tumors that are sometimes seen as part of the hyperparathyroidism-jaw tumor syndrome (HPT-JT), which is caused by inactivating mutations of the HRPT2/CDC73 tumor suppressor gene. CDC73 mutations have been identified in a subset of sporadic cases but aberrant expression of the encoded protein, parafibromin, has not been demonstrated in ossifying fibroma. We sought to determine if loss of parafibromin regularly contributes to the development of sporadic, nonsyndromic ossifying fibroma. We examined a series of 9 ossifying fibromas, including ossifying, cemento-ossifying, and juvenile active variants, for parafibromin protein expression by immunohistochemistry and for CDC73 sequence abnormalities by Sanger sequencing and/or targeted AmpliSeq panel sequencing. Four ossifying fibromas showed a complete absence of nuclear parafibromin expression; loss of parafibromin expression was coupled with aberrant cytoplasmic parafibromin expression in 1 case. CDC73 mutations were detected in 2 cases with aberrant parafibromin expression. These results provide novel evidence, at the level of protein expression, that loss of the parathyroid CDC73/parafibromin tumor suppressor may play a role in the pathogenesis of a subset of ossifying fibromas.

A Novel Framework for Characterizing Genomic Haplotype Diversity in the Human Immunoglobulin Heavy Chain Locus.

An incomplete ascertainment of genetic variation within the highly polymorphic immunoglobulin heavy chain locus (IGH) has hindered our ability to define genetic factors that influence antibody-mediated processes. Due to locus complexity, standard high-throughput approaches have failed to accurately and comprehensively capture IGH polymorphism. As a result, the locus has only been fully characterized two times, severely limiting our knowledge of human IGH diversity. Here, we combine targeted long-read sequencing with a novel bioinformatics tool, IGenotyper, to fully characterize IGH variation in a haplotype-specific manner. We apply this approach to eight human samples, including a haploid cell line and two mother-father-child trios, and demonstrate the ability to generate high-quality assemblies (>98% complete and >99% accurate), genotypes, and gene annotations, identifying 2 novel structural variants and 15 novel IGH alleles. We show multiplexing allows for scaling of the approach without impacting data quality, and that our genotype call sets are more accurate than short-read (>35% increase in true positives and >97% decrease in false-positives) and array/imputation-based datasets. This framework establishes a desperately needed foundation for leveraging IG genomic data to study population-level variation in antibody-mediated immunity, critical for bettering our understanding of disease risk, and responses to vaccines and therapeutics.

Medically actionable pathogenic variants in a population of 13,131 healthy elderly individuals.

PURPOSE: To measure the prevalence of medically actionable pathogenic variants (PVs) among a population of healthy elderly individuals. METHODS: We used targeted sequencing to detect pathogenic or likely pathogenic variants in 55 genes associated with autosomal dominant medically actionable conditions, among a population of 13,131 individuals aged 70 or older (mean age 75 years) enrolled in the ASPirin in Reducing Events in the Elderly (ASPREE) trial. Participants had no previous diagnosis or current symptoms of cardiovascular disease, physical disability or dementia, and no current diagnosis of life-threatening cancer. Variant curation followed American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) standards. RESULTS: One in 75 (1.3%) healthy elderly individuals carried a PV. This was lower than rates reported from population-based studies, which have ranged from 1.8% to 3.4%. We detected 20 PV carriers for Lynch syndrome (MSH6/MLH1/MSH2/PMS2) and 13 for familial hypercholesterolemia (LDLR/APOB/PCSK9). Among 7056 female participants, we detected 15 BRCA1/BRCA2 PV carriers (1 in 470 females). We detected 86 carriers of PVs in lower-penetrance genes associated with inherited cardiac disorders. CONCLUSION: Medically actionable PVs are carried in a healthy elderly population. Our findings raise questions about the actionability of lower-penetrance genes, especially when PVs are detected in the absence of symptoms and/or family history of disease.

Epigenetic preconditioning with decitabine sensitizes glioblastoma to temozolomide via induction of MLH1.

INTRODUCTION: To improve the standard treatment paradigm for glioblastoma (GBM), efforts have been made to explore the efficacy of epigenetic agents as chemosensitizers. Recent data suggest possible synergy between decitabine (DAC), a DNA hypomethylating agent, and temozolomide (TMZ) in GBM, but the mechanism remains unclear. The objective of this study was to determine the effects of DAC on TMZ sensitization in a consecutively derived set of primary GBM cultures, with a focus on mismatch repair (MMR) proteins. METHODS: Half maximal inhibitory concentrations (IC50) of TMZ were calculated in eleven consecutive patient-derived GBM cell lines before and after preconditioning with DAC. MMR protein expression changes were determined by quantitative immunoblots and qPCR arrays. Single-molecule real-time (SMRT) sequencing of bisulfite (BS)-converted PCR amplicons of the MLH1 promoter was performed to determine methylation status. RESULTS: TMZ IC50 significantly changed in 6 of 11 GBM lines of varying MGMT promoter methylation status in response to DAC preconditioning. Knockdown of MLH1 after preconditioning reversed TMZ sensitization. SMRT-BS sequencing of the MLH1 promoter region revealed higher levels of baseline methylation at proximal CpGs in desensitized lines compared to sensitized lines. CONCLUSIONS: DAC enhances TMZ cytotoxicity in a subset of GBM cell lines, comprising lines both MGMT methylated and unmethylated tumors. This effect may be driven by levels of MLH1 via E2F1 transcription factor binding. Using unbiased long-range next-generation bisulfite-sequencing, we identified a region of the proximal MLH1 promoter with differential methylation patterns that has potential utility as a clinical biomarker for TMZ sensitization.

Genetic associations with clozapine-induced myocarditis in patients with schizophrenia.

Clozapine is the most effective antipsychotic drug for schizophrenia, yet it can cause life-threatening adverse drug reactions, including myocarditis. The aim of this study was to determine whether schizophrenia patients with clozapine-induced myocarditis have a genetic predisposition compared with clozapine-tolerant controls. We measured different types of genetic variation, including genome-wide single-nucleotide polymorphisms (SNPs), coding variants that alter protein expression, and variable forms of human leucocyte antigen (HLA) genes, alongside traditional clinical risk factors in 42 cases and 67 controls. We calculated a polygenic risk score (PRS) based on variation at 96 different genetic sites, to estimate the genetic liability to clozapine-induced myocarditis. Our genome-wide association analysis identified four SNPs suggestive of increased myocarditis risk (P < 1 × 10-6), with odds ratios ranging 5.5-13.7. The SNP with the lowest P value was rs74675399 (chr19p13.3, P = 1.21 × 10-7; OR = 6.36), located in the GNA15 gene, previously associated with heart failure. The HLA-C*07:01 allele was identified as potentially predisposing to clozapine-induced myocarditis (OR = 2.89, 95% CI: 1.11-7.53), consistent with a previous report of association of the same allele with clozapine-induced agranulocytosis. Another seven HLA alleles, including HLA-B*07:02 (OR = 0.25, 95% CI: 0.05-1.2) were found to be putatively protective. Long-read DNA sequencing provided increased resolution of HLA typing and validated the HLA associations. The PRS explained 66% of liability (P value = 9.7 × 10-5). Combining clinical and genetic factors together increased the proportion of variability accounted for (r2 0.73, P = 9.8 × 10-9). However, due to the limited sample size, individual genetic associations were not statistically significant after correction for multiple testing. We report novel candidate genetic associations with clozapine-induced myocarditis, which may have potential clinical utility, but larger cohorts are required for replication.

Evolution and Diversification of FRUITFULL Genes in Solanaceae.

Ecologically and economically important fleshy edible fruits have evolved from dry fruit numerous times during angiosperm diversification. However, the molecular mechanisms that underlie these shifts are unknown. In the Solanaceae there has been a major shift to fleshy fruits in the subfamily Solanoideae. Evidence suggests that an ortholog of FRUITFULL (FUL), a transcription factor that regulates cell proliferation and limits the dehiscence zone in the silique of Arabidopsis, plays a similar role in dry-fruited Solanaceae. However, studies have shown that FUL orthologs have taken on new functions in fleshy fruit development, including regulating elements of tomato ripening such as pigment accumulation. FUL belongs to the core eudicot euFUL clade of the angiosperm AP1/FUL gene lineage. The euFUL genes fall into two paralogous clades, euFULI and euFULII. While most core eudicots have one gene in each clade, Solanaceae have two: FUL1 and FUL2 in the former, and MBP10 and MBP20 in the latter. We characterized the evolution of the euFUL genes to identify changes that might be correlated with the origin of fleshy fruit in Solanaceae. Our analyses revealed that the Solanaceae FUL1 and FUL2 clades probably originated through an early whole genome multiplication event. By contrast, the data suggest that the MBP10 and MBP20 clades are the result of a later tandem duplication event. MBP10 is expressed at weak to moderate levels, and its atypical short first intron lacks putative transcription factor binding sites, indicating possible pseudogenization. Consistent with this, our analyses show that MBP10 is evolving at a faster rate compared to MBP20. Our analyses found that Solanaceae euFUL gene duplications, evolutionary rates, and changes in protein residues and expression patterns are not correlated with the shift in fruit type. This suggests deeper analyses are needed to identify the mechanism underlying the change in FUL ortholog function.

Mapping and characterizing N6-methyladenine in eukaryotic genomes using single-molecule real-time sequencing.

N6-Methyladenine (m6dA) has been discovered as a novel form of DNA methylation prevalent in eukaryotes; however, methods for high-resolution mapping of m6dA events are still lacking. Single-molecule real-time (SMRT) sequencing has enabled the detection of m6dA events at single-nucleotide resolution in prokaryotic genomes, but its application to detecting m6dA in eukaryotic genomes has not been rigorously examined. Herein, we identified unique characteristics of eukaryotic m6dA methylomes that fundamentally differ from those of prokaryotes. Based on these differences, we describe the first approach for mapping m6dA events using SMRT sequencing specifically designed for the study of eukaryotic genomes and provide appropriate strategies for designing experiments and carrying out sequencing in future studies. We apply the novel approach to study two eukaryotic genomes. For green algae, we construct the first complete genome-wide map of m6dA at single-nucleotide and single-molecule resolution. For human lymphoblastoid cells (hLCLs), it was necessary to integrate SMRT sequencing data with independent sequencing data. The joint analyses suggest putative m6dA events are enriched in the promoters of young full-length LINE-1 elements (L1s), but call for validation by additional methods. These analyses demonstrate a general method for rigorous mapping and characterization of m6dA events in eukaryotic genomes.

Treatment-associated TP53 DNA-binding domain missense mutations in the pathogenesis of secondary gliosarcoma.

BACKGROUND: Gliosarcoma is a rare variant of glioblastoma (GBM) that exhibits frequent mutations in TP53 and can develop in a secondary fashion after chemoradiation of a primary GBM. Whether temozolomide (TMZ)-induced mutagenesis of the TP53 DNA-binding domain (DBD) can drive the pathogenesis of gliosarcoma is unclear. METHODS: We identified a case of a primary GBM that rapidly progressed into secondary gliosarcoma shortly after chemoradiation was initiated. Bulk tumor was collected and gliomasphere cultures derived from both the pre- and post-treatment tumors. We performed targeted DNA sequencing and transcriptome analyses of the specimens to understand their phylogenetic relationship and identify differentially expressed gene pathways. Gliomaspheres from the primary GBM were treated with TMZ and then analyzed to compare patterns of mutagenesis in vivo and ex vivo. RESULTS: The pre- and post-treatment tumors shared EGFR, CDKN2A, and PTEN mutations, but only the secondary gliosarcoma exhibited TP53 DBD missense mutations. Two mutations, R110C, and R175H, were identified, each in distinct clones. Both were base transitions characteristic of TMZ mutagenesis. Gene expression analysis identified increased JAK-STAT signaling in the gliosarcoma, together with reduced expression of microRNAs known to regulate epithelial-mesenchymal transition. Ex vivo treatment of the GBM spheres with TMZ generated numerous variants in cancer driver genes, including TP53 and CDH1, which were mutated in the post-treatment tumor. CONCLUSIONS: TMZ-induced TP53 gain-of-function mutations can have a driving role in secondary gliosarcoma pathogenesis. Analysis of variants identified in ex vivo TMZ-treated gliomaspheres may have utility in predicting GBM evolutionary trajectories in vivo during standard chemoradiation.

Genomic profiling reveals mutational landscape in parathyroid carcinomas.

Parathyroid carcinoma (PC) is an extremely rare malignancy lacking effective therapeutic intervention. We generated and analyzed whole-exome sequencing data from 17 patients to identify somatic and germline genetic alterations. A panel of selected genes was sequenced in a 7-tumor expansion cohort. We show that 47% (8 of 17) of the tumors harbor somatic mutations in the CDC73 tumor suppressor, with germline inactivating variants in 4 of the 8 patients. The PI3K/AKT/mTOR pathway was altered in 21% of the 24 cases, revealing a major oncogenic pathway in PC. We observed CCND1 amplification in 29% of the 17 patients, and a previously unreported recurrent mutation in putative kinase ADCK1. We identified the first sporadic PCs with somatic mutations in the Wnt canonical pathway, complementing previously described epigenetic mechanisms mediating Wnt activation. This is the largest genomic sequencing study of PC, and represents major progress toward a full molecular characterization of this rare malignancy to inform improved and individualized treatments.

Variation in the flowering time orthologs BrFLC and BrSOC1 in a natural population of Brassica rapa.

Understanding the genetic basis of natural phenotypic variation is of great importance, particularly since selection can act on this variation to cause evolution. We examined expression and allelic variation in candidate flowering time loci in Brassica rapa plants derived from a natural population and showing a broad range in the timing of first flowering. The loci of interest were orthologs of the Arabidopsis genes FLC and SOC1 (BrFLC and BrSOC1, respectively), which in Arabidopsis play a central role in the flowering time regulatory network, with FLC repressing and SOC1 promoting flowering. In B. rapa, there are four copies of FLC and three of SOC1. Plants were grown in controlled conditions in the lab. Comparisons were made between plants that flowered the earliest and latest, with the difference in average flowering time between these groups ∼30 days. As expected, we found that total expression of BrSOC1 paralogs was significantly greater in early than in late flowering plants. Paralog-specific primers showed that expression was greater in early flowering plants in the BrSOC1 paralogs Br004928, Br00393 and Br009324, although the difference was not significant in Br009324. Thus expression of at least 2 of the 3 BrSOC1 orthologs is consistent with their predicted role in flowering time in this natural population. Sequences of the promoter regions of the BrSOC1 orthologs were variable, but there was no association between allelic variation at these loci and flowering time variation. For the BrFLC orthologs, expression varied over time, but did not differ between the early and late flowering plants. The coding regions, promoter regions and introns of these genes were generally invariant. Thus the BrFLC orthologs do not appear to influence flowering time in this population. Overall, the results suggest that even for a trait like flowering time that is controlled by a very well described genetic regulatory network, understanding the underlying genetic basis of natural variation in such a quantitative trait is challenging.