Resistance to mesenchymal reprogramming sustains clonal propagation in metastatic breast cancer

The acquisition of mesenchymal traits is considered a hallmark of breast cancer progression. However, the functional relevance of epithelial-to-mesenchymal transition (EMT) remains controversial and context dependent. Here, we isolate epithelial and mesenchymal populations from human breast cancer metastatic biopsies and assess their functional potential in vivo. Strikingly, progressively decreasing epithelial cell adhesion molecule (EPCAM) levels correlate with declining disease propagation. Mechanistically, we find that persistent EPCAM expression marks epithelial clones that resist EMT induction and propagate competitively. In contrast, loss of EPCAM defines clones arrested in a mesenchymal state, with concomitant suppression of tumorigenicity and metastatic potential. This dichotomy results from distinct clonal trajectories impacting global epigenetic programs that are determined by the interplay between human ZEB1 and its target GRHL2. Collectively, our results indicate that susceptibility to irreversible EMT restrains clonal propagation, whereas resistance to mesenchymal reprogramming sustains disease spread in multiple models of human metastatic breast cancer, including patient-derived cells in vivo.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics

BACKGROUND

Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics.

METHODS

We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage.

RESULTS

We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions.

CONCLUSION

Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.

A macrophage-hepatocyte glucocorticoid receptor axis coordinates fasting ketogenesis

Fasting metabolism and immunity are tightly linked; however, it is largely unknown how immune cells contribute to metabolic homeostasis during fasting in healthy subjects. Here, we combined cell-type-resolved genomics and computational approaches to map crosstalk between hepatocytes and liver macrophages during fasting. We identified the glucocorticoid receptor (GR) as a key driver of fasting-induced reprogramming of the macrophage secretome including fasting-suppressed cytokines and showed that lack of macrophage GR impaired induction of ketogenesis during fasting as well as endotoxemia. Mechanistically, macrophage GR suppressed the expression of tumor necrosis factor (TNF) and promoted nuclear translocation of hepatocyte GR to activate a fat oxidation/ketogenesis-related gene program, cooperatively induced by GR and peroxisome proliferator-activated receptor alpha (PPARα) in hepatocytes. Together, our results demonstrate how resident liver macrophages directly influence ketogenesis in hepatocytes, thereby also outlining a strategy by which the immune system can set the metabolic tone during inflammatory disease and infection.

Cell cycle defects underlie childhood-onset cardiomyopathy associated with Noonan syndrome

Childhood-onset myocardial hypertrophy and cardiomyopathic changes are associated with significant morbidity and mortality in early life, particularly in patients with Noonan syndrome, a multisystemic genetic disorder caused by autosomal dominant mutations in genes of the Ras-MAPK pathway. Although the cardiomyopathy associated with Noonan syndrome (NS-CM) shares certain cardiac features with the hypertrophic cardiomyopathy caused by mutations in sarcomeric proteins (HCM), such as pathological myocardial remodeling, ventricular dysfunction, and increased risk for malignant arrhythmias, the clinical course of NS-CM significantly differs from HCM. This suggests a distinct pathophysiology that remains to be elucidated. Here, through analysis of sarcomeric myosin conformational states, histopathology, and gene expression in left ventricular myocardial tissue from NS-CM, HCM, and normal hearts complemented with disease modeling in cardiomyocytes differentiated from patient-derived PTPN11 N308S/+ induced pluripotent stem cells, we demonstrate distinct disease phenotypes between NS-CM and HCM and uncover cell cycle defects as a potential driver of NS-CM.

SP8 Promotes an Aggressive Phenotype in Hepatoblastoma via FGF8 Activation

Hepatoblastoma (HB) is the most common malignant liver tumor in childhood and it generally has a good prognosis. However, if associated with aggressive metastatic disease, outcome is still poor. The molecular mechanisms leading to metastatic spread in HB patients are still unknown. By combining RNA-sequencing and a genome-wide methylome analysis, we identified the transcription factor SP8 and the growth factor FGF8 among the most strongly upregulated genes in metastatic HB cases, with a concomitant robust demethylation of the respective promoter regions. Of note, high expression of both candidates was associated with the aggressive C2 subtype of the 16-gene signature and poor survival. Chromatin immunoprecipitation revealed a direct transcriptional regulation of FGF8 through binding of SP8 to the FGF8 promoter. Gain- and loss-of-function experiments proved promoting effects of SP8 on motility, self-renewal, migration, and the invasive potential of HB cells. Moreover, stable overexpression of SP8 in Hep3B cells resulted in the acquisition of a mesenchymal phenotype and a strong upregulation of epithelial-mesenchymal transition-associated genes. Using KRAB-mediated CRISPR-dCas9 interference directed against FGF8, we could show that FGF8 is essential for the SP8-mediated aggressive tumor behavior. Treatment of HB cell lines with the pan SP family inhibitor mithramycin A resulted in a significant inhibition of their clonogenic growth. In summary, we identified SP8 and FGF8 as key players in aggressive traits of HB and propose SP8 inhibiting drugs as a new effective treatment strategy especially for metastatic tumors.

MAU2 and NIPBL Variants Impair the Heterodimerization of the Cohesin Loader Subunits and Cause Cornelia de Lange Syndrome

The NIPBL/MAU2 heterodimer loads cohesin onto chromatin. Mutations in NIPBL account for most cases of the rare developmental disorder Cornelia de Lange syndrome (CdLS). Here we report a MAU2 variant causing CdLS, a deletion of seven amino acids that impairs the interaction between MAU2 and the NIPBL N terminus. Investigating this interaction, we discovered that MAU2 and the NIPBL N terminus are largely dispensable for normal cohesin and NIPBL function in cells with a NIPBL early truncating mutation. Despite a predicted fatal outcome of an out-of-frame single nucleotide duplication in NIPBL, engineered in two different cell lines, alternative translation initiation yields a form of NIPBL missing N-terminal residues. This form cannot interact with MAU2, but binds DNA and mediates cohesin loading. Altogether, our work reveals that cohesin loading can occur independently of functional NIPBL/MAU2 complexes and highlights a novel mechanism protective against out-of-frame mutations that is potentially relevant for other genetic conditions.

Bi-allelic variants in RNF170 are associated with hereditary spastic paraplegia

Alterations of Ca2+ homeostasis have been implicated in a wide range of neurodegenerative diseases. Ca2+ efflux from the endoplasmic reticulum into the cytoplasm is controlled by binding of inositol 1,4,5-trisphosphate to its receptor. Activated inositol 1,4,5-trisphosphate receptors are then rapidly degraded by the endoplasmic reticulum-associated degradation pathway. Mutations in genes encoding the neuronal isoform of the inositol 1,4,5-trisphosphate receptor (ITPR1) and genes involved in inositol 1,4,5-trisphosphate receptor degradation (ERLIN1, ERLIN2) are known to cause hereditary spastic paraplegia (HSP) and cerebellar ataxia. We provide evidence that mutations in the ubiquitin E3 ligase gene RNF170, which targets inositol 1,4,5-trisphosphate receptors for degradation, are the likely cause of autosomal recessive HSP in four unrelated families and functionally evaluate the consequences of mutations in patient fibroblasts, mutant SH-SY5Y cells and by gene knockdown in zebrafish. Our findings highlight inositol 1,4,5-trisphosphate signaling as a candidate key pathway for hereditary spastic paraplegias and cerebellar ataxias and thus prioritize this pathway for therapeutic interventions.

Epstein-Barr virus reprograms human B lymphocytes immediately in the prelatent phase of infection

Epstein-Barr virus (EBV) is a human tumor virus and a model of herpesviral latency. The virus efficiently infects resting human B lymphocytes and induces their continuous proliferation in vitro, which mimics certain aspects of EBV's oncogenic potential in vivo. How lymphoblastoid cell lines (LCLs) evolve from the infected lymphocytes is uncertain. We conducted a systematic time-resolved longitudinal study of cellular functions and transcriptional profiles of newly infected naïve primary B lymphocytes. EBV reprograms the cells comprehensively and globally. Rapid and extensive transcriptional changes occur within 24 h and precede any metabolic and phenotypic changes. Within 72 h, the virus activates the cells, changes their phenotypes with respect to cell size, RNA, and protein content, and induces metabolic pathways to cope with the increased demand for energy, supporting an efficient cell cycle entry on day 3 postinfection. The transcriptional program that EBV initiates consists of 3 waves of clearly discernable clusters of cellular genes that peak on day 2, 3, or 4 and regulate RNA synthesis, metabolic pathways, and cell division, respectively. Upon onset of cell doublings on day 4, the cellular transcriptome appears to be completely reprogrammed to support the proliferating cells, but 3 additional clusters of EBV-regulated genes fine-tune cell signaling, migration, and immune response pathways, eventually. Our study reveals that more than 11,000 genes are regulated upon EBV infection as naïve B cells exit quiescence to enter a germinal center-like differentiation program, which culminates in immortalized, proliferating cells that partially resemble plasmablasts and early plasma cells.

Intermuscular adipose tissue directly modulates skeletal muscle insulin sensitivity in humans

Intermuscular adipose tissue (IMAT) is negatively related to insulin sensitivity, but a causal role of IMAT in the development of insulin resistance is unknown. IMAT was sampled in humans to test for the ability to induce insulin resistance in vitro and characterize gene expression to uncover how IMAT may promote skeletal muscle insulin resistance. Human primary muscle cells were incubated with conditioned media from IMAT, visceral (VAT), or subcutaneous adipose tissue (SAT) to evaluate changes in insulin sensitivity. RNAseq analysis was performed on IMAT with gene expression compared with skeletal muscle and SAT, and relationships to insulin sensitivity were determined in men and women spanning a wide range of insulin sensitivity measured by hyperinsulinemic-euglycemic clamp. Conditioned media from IMAT and VAT decreased insulin sensitivity similarly compared with SAT. Multidimensional scaling analysis revealed distinct gene expression patterns in IMAT compared with SAT and muscle. Pathway analysis revealed that IMAT expression of genes in insulin signaling, oxidative phosphorylation, and peroxisomal metabolism related positively to donor insulin sensitivity, whereas expression of macrophage markers, inflammatory cytokines, and secreted extracellular matrix proteins were negatively related to insulin sensitivity. Perilipin 5 gene expression suggested greater IMAT lipolysis in insulin-resistant individuals. Combined, these data show that factors secreted from IMAT modulate muscle insulin sensitivity, possibly via secretion of inflammatory cytokines and extracellular matrix proteins, and by increasing local FFA concentration in humans. These data suggest IMAT may be an important regulator of skeletal muscle insulin sensitivity and could be a novel therapeutic target for skeletal muscle insulin resistance.

Mutations in PPCS, Encoding Phosphopantothenoylcysteine Synthetase, Cause Autosomal-Recessive Dilated Cardiomyopathy

Coenzyme A (CoA) is an essential metabolic cofactor used by around 4% of cellular enzymes. Its role is to carry and transfer acetyl and acyl groups to other molecules. Cells can synthesize CoA de novo from vitamin B5 (pantothenate) through five consecutive enzymatic steps. Phosphopantothenoylcysteine synthetase (PPCS) catalyzes the second step of the pathway during which phosphopantothenate reacts with ATP and cysteine to form phosphopantothenoylcysteine. Inborn errors of CoA biosynthesis have been implicated in neurodegeneration with brain iron accumulation (NBIA), a group of rare neurological disorders characterized by accumulation of iron in the basal ganglia and progressive neurodegeneration. Exome sequencing in five individuals from two unrelated families presenting with dilated cardiomyopathy revealed biallelic mutations in PPCS, linking CoA synthesis with a cardiac phenotype. Studies in yeast and fruit flies confirmed the pathogenicity of identified mutations. Biochemical analysis revealed a decrease in CoA levels in fibroblasts of all affected individuals. CoA biosynthesis can occur with pantethine as a source independent from PPCS, suggesting pantethine as targeted treatment for the affected individuals still alive.

Overexpression of UHRF1 promotes silencing of tumor suppressor genes and predicts outcome in hepatoblastoma

Background

Hepatoblastoma (HB) is the most common liver tumor of childhood and occurs predominantly within the first 3 years of life. In accordance to its early manifestation, HB has been described to display an extremely low mutation rate. As substitute, epigenetic modifiers seem to play an exceptional role in its tumorigenesis, which holds promise to develop targeted therapies and establish biomarkers for patient risk stratification.

Results

We examined the role of a newly described protein complex consisting of three epigenetic regulators, namely E3 ubiquitin-like containing PHD and RING finger domain 1 (UHRF1), ubiquitin-specific-processing protease 7 (USP7), and DNA methyltransferase 1 (DNMT1), in HB. We found the complex to be located on the promoter regions of the pivotal HB-associated tumor suppressor genes (TSGs) HHIP, IGFBP3, and SFRP1 in HB cells, thereby leading to strong repression through DNA methylation and histone modifications. Consequently, knockdown of UHRF1 led to DNA demethylation and loss of the repressive H3K9me2 histone mark at the TSG loci with their subsequent transcriptional reactivation. The observed growth impairment of HB cells upon UHRF1 knockdown could be attributed to reduced expression of genes involved in cell cycle progression, negative regulation of cell death, LIN28B signaling, and the adverse 16-gene signature, as revealed by global RNA sequencing. Clinically, overexpression of UHRF1 in primary tumor tissues was significantly associated with poor survival and the prognostic high-risk 16-gene signature.

Conclusion

These findings suggest that UHRF1 is critical for aberrant TSG silencing and sustained growth signaling in HB and that UHRF1 overexpression levels might serve as a prognostic biomarker and potential molecular target for HB patients.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0462-7) contains supplementary material, which is available to authorized users.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

Long-Term Cold Adaptation Does Not Require FGF21 or UCP1

Brown adipose tissue (BAT)-dependent thermogenesis and its suggested augmenting hormone, FGF21, are potential therapeutic targets in current obesity and diabetes research. Here, we studied the role of UCP1 and FGF21 for metabolic homeostasis in the cold and dissected underlying molecular mechanisms using UCP1-FGF21 double-knockout mice. We report that neither UCP1 nor FGF21, nor even compensatory increases of FGF21 serum levels in UCP1 knockout mice, are required for defense of body temperature or for maintenance of energy metabolism and body weight. Remarkably, cold-induced browning of inguinal white adipose tissue (iWAT) is FGF21 independent. Global RNA sequencing reveals major changes in response to UCP1- but not FGF21-ablation in BAT, iWAT, and muscle. Markers of mitochondrial failure and inflammation are observed in BAT, but in particular the enhanced metabolic reprogramming in iWAT supports the thermogenic role of UCP1 and excludes an important thermogenic role of endogenous FGF21 in normal cold acclimation.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Coexisting variants in OSTM1 and MANEAL cause a complex neurodegenerative disorder with NBIA-like brain abnormalities

Coexistence of different hereditary diseases is a known phenomenon in populations with a high consanguinity rate. The resulting clinical phenotypes are extremely challenging for physicians involved in the care of these patients. Here we describe a 6-year-old boy with co-occurrence of a homozygous splice defect in OSTM1, causing infantile malignant osteopetrosis, and a loss-of-function variant in MANEAL, which has not been associated with human disease so far. The child suffered from severe infantile-onset neurodegeneration that could not be stopped by bone marrow transplantation. Magnetic resonance imaging demonstrated global brain atrophy and showed hypointensities of globus pallidus, corpora mamillaria, and cerebral peduncles, which were comparable to findings in neurodegeneration with brain iron accumulation disorders. LC-MS/MS analysis of urine and cerebrospinal fluid samples revealed a distinct metabolic profile with accumulation of mannose tetrasaccharide molecules, suggestive of an oligosaccharide storage disease. Our results demonstrate that exome sequencing is a very effective tool in dissecting complex neurological diseases. Moreover, we suggest that MANEAL is an interesting candidate gene that should be considered in the context of neurological disorders with brain iron accumulation and/or indications of an oligosaccharide storage disease.

Genetic diagnosis of Mendelian disorders via RNA sequencing

Across a variety of Mendelian disorders, ∼50-75% of patients do not receive a genetic diagnosis by exome sequencing indicating disease-causing variants in non-coding regions. Although genome sequencing in principle reveals all genetic variants, their sizeable number and poorer annotation make prioritization challenging. Here, we demonstrate the power of transcriptome sequencing to molecularly diagnose 10% (5 of 48) of mitochondriopathy patients and identify candidate genes for the remainder. We find a median of one aberrantly expressed gene, five aberrant splicing events and six mono-allelically expressed rare variants in patient-derived fibroblasts and establish disease-causing roles for each kind. Private exons often arise from cryptic splice sites providing an important clue for variant prioritization. One such event is found in the complex I assembly factor TIMMDC1 establishing a novel disease-associated gene. In conclusion, our study expands the diagnostic tools for detecting non-exonic variants and provides examples of intronic loss-of-function variants with pathological relevance.

Genetic Landscape of Sporadic Unilateral Adrenocortical Adenomas Without PRKACA p.Leu206Arg Mutation

CONTEXT

Adrenocortical adenomas (ACAs) are among the most frequent human neoplasias. Genetic alterations affecting the cAMP/protein kinase A signaling pathway are common in cortisol-producing ACAs, whereas activating mutations in the gene encoding β-catenin (CTNNB1) have been reported in a subset of both benign and malignant adrenocortical tumors. However, the molecular pathogenesis of most ACAs is still largely unclear.

OBJECTIVE

The aim of the study was to define the genetic landscape of sporadic unilateral ACAs.

DESIGN AND SETTING

Next-generation whole-exome sequencing was performed on fresh-frozen tumor samples and corresponding normal tissue samples.

PATIENTS

Ninety-nine patients with ACAs (74 cortisol-producing and 25 endocrine inactive) negative for p.Leu206Arg PRKACA mutation.

MAIN OUTCOME MEASURES

Identification of known and/or new genetic alterations potentially involved in adrenocortical tumorigenesis and autonomous hormone secretion, genotype-phenotype correlation.

RESULTS

A total of 706 somatic protein-altering mutations were detected in 88 of 99 tumors (median, six per tumor). We identified several mutations in genes of the cAMP/protein kinase A pathway, including three novel mutations in PRKACA, associated with female sex and Cushing's syndrome. We also found genetic alterations in different genes involved in the Wnt/β-catenin pathway, associated with larger tumors and endocrine inactivity, and notably, in many genes of the Ca(2+)-signaling pathway. Finally, by comparison of our genetic data with those available in the literature, we describe a comprehensive genetic landscape of unilateral ACAs.

CONCLUSIONS

This study provides the largest sequencing effort on ACAs to date. We thereby identified somatic alterations affecting known and novel pathways potentially involved in adrenal tumorigenesis.

Recurrent EZH1 mutations are a second hit in autonomous thyroid adenomas

Autonomous thyroid adenomas (ATAs) are a frequent cause of hyperthyroidism. Mutations in the genes encoding the TSH receptor (TSHR) or the Gs protein α subunit (GNAS) are found in approximately 70% of ATAs. The involvement of other genes and the pathogenesis of the remaining cases are presently unknown. Here, we performed whole-exome sequencing in 19 ATAs that were paired with normal DNA samples and identified a recurrent hot-spot mutation (c.1712A>G; p.Gln571Arg) in the enhancer of zeste homolog 1 (EZH1) gene, which codes for a catalytic subunit of the polycomb complex. Targeted screening in an independent cohort confirmed that this mutation occurs with high frequency (27%) in ATAs. EZH1 mutations were strongly associated with known (TSHR, GNAS) or presumed (adenylate cyclase 9 [ADCY9]) alterations in cAMP pathway genes. Furthermore, functional studies revealed that the p.Gln571Arg EZH1 mutation caused increased histone H3 trimethylation and increased proliferation of thyroid cells. In summary, this study revealed that a hot-spot mutation in EZH1 is the second most frequent genetic alteration in ATAs. The association between EZH1 and TSHR mutations suggests a 2-hit model for the pathogenesis of these tumors, whereby constitutive activation of the cAMP pathway and EZH1 mutations cooperate to induce the hyperproliferation of thyroid cells.

PRKACA Somatic Mutations Are Rare Findings in Aldosterone-Producing Adenomas

CONTEXT

Somatic mutations have been found causative for endocrine autonomy in aldosterone-producing adenomas (APAs). Whereas mutations of PRKACA (catalytic subunit of protein kinase A) have been identified in cortisol-producing adenomas, the presence of PRKACA variants in APAs is unknown, especially in those that display cosecretion of cortisol.

OBJECTIVE

The objective of the study was to investigate PRKACA somatic variants identified in APA cases.

DESIGN

Identification of PRKACA somatic variants in APAs by whole-exome sequencing followed by in vitro analysis of the enzymatic activity of PRKACA variants and functional characterization by double immunofluorescence of CYP11B2 and CYP11B1 expression in the corresponding tumor tissues.

SETTING AND PATIENTS

APA tissues were collected from 122 patients who underwent unilateral adrenalectomy for primary aldosteronism between 2005 and 2015 at a single institution.

RESULTS

PRKACA somatic mutations were identified in two APA cases (1.6%). One APA carried a newly identified p.His88Asp variant, whereas in a second case, a p.Leu206Arg mutation was found, previously described only in cortisol-producing adenomas with overt Cushing's syndrome. Functional analysis showed that the p.His88Asp variant was not associated with gain of function. Although CYP11B2 was strongly expressed in the p.His88Asp-mutated APA, the p.Leu206Arg carrying APA predominantly expressed CYP11B1. Accordingly, biochemical Cushing's syndrome was present only in the patient with the p.Leu206Arg mutation. After adrenalectomy, both patients improved with a reduced number of antihypertensive medications and normalized serum potassium levels.

CONCLUSIONS

We describe for the first time PRKACA mutations as rare findings associated with unilateral primary aldosteronism. As cortisol cosecretion occurs in a subgroup of APAs, other molecular mechanisms are likely to exist.

Riboflavin-Responsive and -Non-responsive Mutations in FAD Synthase Cause Multiple Acyl-CoA Dehydrogenase and Combined Respiratory-Chain Deficiency

Multiple acyl-CoA dehydrogenase deficiencies (MADDs) are a heterogeneous group of metabolic disorders with combined respiratory-chain deficiency and a neuromuscular phenotype. Despite recent advances in understanding the genetic basis of MADD, a number of cases remain unexplained. Here, we report clinically relevant variants in FLAD1, which encodes FAD synthase (FADS), as the cause of MADD and respiratory-chain dysfunction in nine individuals recruited from metabolic centers in six countries. In most individuals, we identified biallelic frameshift variants in the molybdopterin binding (MPTb) domain, located upstream of the FADS domain. Inasmuch as FADS is essential for cellular supply of FAD cofactors, the finding of biallelic frameshift variants was unexpected. Using RNA sequencing analysis combined with protein mass spectrometry, we discovered FLAD1 isoforms, which only encode the FADS domain. The existence of these isoforms might explain why affected individuals with biallelic FLAD1 frameshift variants still harbor substantial FADS activity. Another group of individuals with a milder phenotype responsive to riboflavin were shown to have single amino acid changes in the FADS domain. When produced in E. coli, these mutant FADS proteins resulted in impaired but detectable FADS activity; for one of the variant proteins, the addition of FAD significantly improved protein stability, arguing for a chaperone-like action similar to what has been reported in other riboflavin-responsive inborn errors of metabolism. In conclusion, our studies identify FLAD1 variants as a cause of potentially treatable inborn errors of metabolism manifesting with MADD and shed light on the mechanisms by which FADS ensures cellular FAD homeostasis.

Landscape of somatic mutations in sporadic GH-secreting pituitary adenomas

CONTEXT

Alterations in the cAMP signaling pathway are common in hormonally active endocrine tumors. Somatic mutations at GNAS are causative in 30-40% of GH-secreting adenomas. Recently, mutations affecting the USP8 and PRKACA gene have been reported in ACTH-secreting pituitary adenomas and cortisol-secreting adrenocortical adenomas respectively. However, the pathogenesis of many GH-secreting adenomas remains unclear.

AIM

Comprehensive genetic characterization of sporadic GH-secreting adenomas and identification of new driver mutations.

DESIGN

Screening for somatic mutations was performed in 67 GH-secreting adenomas by targeted sequencing for GNAS, PRKACA, and USP8 mutations (n=31) and next-generation exome sequencing (n=36).

RESULTS

By targeted sequencing, known activating mutations in GNAS were detected in five cases (16.1%), while no somatic mutations were observed in both PRKACA and USP8. Whole-exome sequencing identified 132 protein-altering somatic mutations in 31/36 tumors with a median of three mutations per sample (range: 1-13). The only recurrent mutations have been observed in GNAS (31.4% of cases). However, seven genes involved in cAMP signaling pathway were affected in 14 of 36 samples and eight samples harbored variants in genes involved in the calcium signaling or metabolism. At the enrichment analysis, several altered genes resulted to be associated with developmental processes. No significant correlation between genetic alterations and the clinical data was observed.

CONCLUSION

This study provides a comprehensive analysis of somatic mutations in a large series of GH-secreting adenomas. No novel recurrent genetic alterations have been observed, but the data suggest that beside cAMP pathway, calcium signaling might be involved in the pathogenesis of these tumors.

Time- and compartment-resolved proteome profiling of the extracellular niche in lung injury and repair

The extracellular matrix (ECM) is a key regulator of tissue morphogenesis and repair. However, its composition and architecture are not well characterized. Here, we monitor remodeling of the extracellular niche in tissue repair in the bleomycin-induced lung injury mouse model. Mass spectrometry quantified 8,366 proteins from total tissue and bronchoalveolar lavage fluid (BALF) over the course of 8 weeks, surveying tissue composition from the onset of inflammation and fibrosis to its full recovery. Combined analysis of proteome, secretome, and transcriptome highlighted post-transcriptional events during tissue fibrogenesis and defined the composition of airway epithelial lining fluid. To comprehensively characterize the ECM, we developed a quantitative detergent solubility profiling (QDSP) method, which identified Emilin-2 and collagen-XXVIII as novel constituents of the provisional repair matrix. QDSP revealed which secreted proteins interact with the ECM, and showed drastically altered association of morphogens to the insoluble matrix upon injury. Thus, our proteomic systems biology study assigns proteins to tissue compartments and uncovers their dynamic regulation upon lung injury and repair, potentially contributing to the development of anti-fibrotic strategies.

Biallelic Mutations in NBAS Cause Recurrent Acute Liver Failure with Onset in Infancy

Acute liver failure (ALF) in infancy and childhood is a life-threatening emergency. Few conditions are known to cause recurrent acute liver failure (RALF), and in about 50% of cases, the underlying molecular cause remains unresolved. Exome sequencing in five unrelated individuals with fever-dependent RALF revealed biallelic mutations in NBAS. Subsequent Sanger sequencing of NBAS in 15 additional unrelated individuals with RALF or ALF identified compound heterozygous mutations in an additional six individuals from five families. Immunoblot analysis of mutant fibroblasts showed reduced protein levels of NBAS and its proposed interaction partner p31, both involved in retrograde transport between endoplasmic reticulum and Golgi. We recommend NBAS analysis in individuals with acute infantile liver failure, especially if triggered by fever.

Frequency and clinical correlates of somatic Ying Yang 1 mutations in sporadic insulinomas

CONTEXT

Insulinomas represent pancreatic neuroendocrine neoplasms that cause severe morbidity attributed to their often pronounced endocrine activity. Apart from hereditary forms such as multiple endocrine neoplasia type 1 (MEN-1), genetic causes for sporadic insulinoma development had remained obscure until recently. Applying next-generation sequencing methods, disease-causing genetic alterations have been identified in various endocrine tumors.

OBJECTIVE AND DESIGN

Paired tumor and blood DNA from eight patients with sporadic insulinomas (five females and two malignant tumors) were analyzed by whole-exome sequencing. After this initial analysis, Ying Yang 1 (YY1) mutation status was assessed in a larger cohort of 39 additional insulinomas (including eight malignant and one liver metastasis) from three German hospitals by targeted sequencing. The mutation status was correlated with various clinical parameters.

RESULTS

A range of one to 12 somatic genetic variants were identified by exome sequencing. A recurrent somatic Thr372Arg YY1 point mutation was detected in two patients of the initial cohort and four patients of the second cohort (total, six of 47; 13%). The presence of the mutation was associated with a trend toward higher age (63.5 y; IQR, 48.0-74.0 vs 45.0 y; IQR, 33.0-63.0; P = .05), and all affected patients were females (six of six; P = .04). All other clinical parameters, including the presence of malignancy and metastatic spread, tumor localization, and hypoglycemic episodes were not different between YY1-mutated and nonmutated tumor carriers.

CONCLUSIONS

The somatic Thr372Arg YY1 mutation is a relevant finding in female patients with sporadic insulinomas. The prevalence of this mutation in this Caucasian population is considerably lower compared to that of a recently described Asian cohort.

Compound heterozygosity of low-frequency promoter deletions and rare loss-of-function mutations in TXNL4A causes Burn-McKeown syndrome

Mutations in components of the major spliceosome have been described in disorders with craniofacial anomalies, e.g., Nager syndrome and mandibulofacial dysostosis type Guion-Almeida. The U5 spliceosomal complex of eight highly conserved proteins is critical for pre-mRNA splicing. We identified biallelic mutations in TXNL4A, a member of this complex, in individuals with Burn-McKeown syndrome (BMKS). This rare condition is characterized by bilateral choanal atresia, hearing loss, cleft lip and/or palate, and other craniofacial dysmorphisms. Mutations were found in 9 of 11 affected families. In 8 families, affected individuals carried a rare loss-of-function mutation (nonsense, frameshift, or microdeletion) on one allele and a low-frequency 34 bp deletion (allele frequency 0.76%) in the core promoter region on the other allele. In a single highly consanguineous family, formerly diagnosed as oculo-oto-facial dysplasia, the four affected individuals were homozygous for a 34 bp promoter deletion, which differed from the promoter deletion in the other families. Reporter gene and in vivo assays showed that the promoter deletions led to reduced expression of TXNL4A. Depletion of TXNL4A (Dib1) in yeast demonstrated reduced assembly of the tri-snRNP complex. Our results indicate that BMKS is an autosomal-recessive condition, which is frequently caused by compound heterozygosity of low-frequency promoter deletions in combination with very rare loss-of-function mutations.

The genomic landscape of hepatoblastoma and their progenies with HCC-like features

BACKGROUND & AIMS

Hepatoblastoma (HB) is the most common childhood liver cancer and occasionally presents with histological and clinical features reminiscent of hepatocellular carcinoma (HCC). Identification of molecular mechanisms that drive the neoplastic continuation towards more aggressive HCC phenotypes may help to guide the new stage of targeted therapies.

METHODS

We performed comprehensive studies on genetic and chromosomal alterations as well as candidate gene function and their clinical relevance.

RESULTS

Whole-exome sequencing identified HB as a genetically very simple tumour (2.9 mutations per tumour) with recurrent mutations in ß-catenin (CTNNB1) (12/15 cases) and the transcription factor NFE2L2 (2/15 cases). Their HCC-like progenies share the common CTNNB1 mutation, but additionally exhibit a significantly increased mutation number and chromosomal instability due to deletions of the genome guardians RAD17 and TP53, accompanied by telomerase reverse-transcriptase (TERT) promoter mutations. Targeted genotyping of 33 primary tumours and cell lines revealed CTNNB1, NFE2L2, and TERT mutations in 72.5%, 9.8%, and 5.9% of cases, respectively. All NFE2L2 mutations affected residues of the NFE2L2 protein that are recognized by the KEAP1/CUL3 complex for proteasomal degradation. Consequently, cells transfected with mutant NFE2L2 were insensitive to KEAP1-mediated downregulation of NFE2L2 signalling. Clinically, overexpression of the NFE2L2 target gene NQO1 in tumours was significantly associated with metastasis, vascular invasion, the adverse prognostic C2 gene signature, as well as poor outcome.

CONCLUSIONS

Our study demonstrates the importance of CTNNB1 mutations and NFE2L2-KEAP1 pathway activation in HB development and defines loss of genomic stability and TERT promoter mutations as prominent characteristics of aggressive HB with HCC features.

Abstract LB-182: Constitutive activation of PRKACA in adrenal Cushing's syndrome

Corticotropin-independent Cushing's syndrome may be caused by tumors or hyperplasia of the adrenal cortex. Until now genetic alterations explain only a small fraction of cases. The observation that a subset of adrenal adenomas is characterized by abnormal PKA activity despite the absence of mutations in candidate genes suggested as yet unknown alterations in the cAMP/PKA signaling cascade in these tumors. The aim of this study was the analysis of the genetic basis of Cushing's syndrome in order to reveal the gene/s responsible for the disease. Exome sequencing was performed in ten cortisol-producing adenomas and recurrent mutations in candidate genes were evaluated in additional 171 patients with adrenocortical tumors. Genome-wide copy number analysis was performed in 35 patients with cortisol secreting bilateral hyperplasias. The effects of these genetic defects were studied both clinically and in vitro. Exome sequencing in 8/10 adenomas revealed somatic mutations in the PRKACA gene, which encodes the main catalytic subunit of cyclic AMP-dependent protein kinase (PKA) (c.617A>C in seven and c.595_596insCAC in one). Overall, PRKACA somatic mutations were identified in a total of 22/59 (37%) adenomas from patients with overt Cushing's syndrome while these mutations were not detectable in patients with subclinical hypercortisolism (n=40) or in other adrenal tumors (n=82). Among 35 patients with cortisol producing hyperplasias, 5 (with two patients as first degree relatives) carried germline copy number gain of the chromosome 19 region including the PRKACA gene. In vitro studies demonstrated impaired inhibition of the mutant PRKACA by the PKA regulatory subunit, while cells from patients with germline chromosomal gains showed increased protein levels; in both cases, PKA activity was increased. The present study shows that more than one third of cortisol-producing adenomas associated with overt Cushing syndrome harbor unique somatic mutations of the main cAMP-dependent kinase catalytic subunit, PRKACA resulting in constitutive PKA activation. While in these patients the mutation is present only in tumor cells, germline duplication of the PRKACA gene was identified in a group of patients with bilateral adrenal hyperplasias. This is the first report of genetic alterations of the catalytic subunit of PKA linked to human disease: Germline PRKACA duplications with bilateral adrenal hyperplasias and somatic PRKACA mutations with unilateral cortisol producing adrenal adenomas.

Citation Format: Fabio R. Faucz, Felix Beuschlei, Martin Fassnacht, Guilaume Assie, Davide Calebiro, Constantine Stratakis, Andrea Osswald, Cristina L. Ronchi, Thomas Wieland, Silviu Sbiera, Katrin Schaak, Anett Schmittfull, Thomas Schwarzmayr, Olivia Barreau, Delphine Vezzosi, Marthe Rizk-Rabbin, Ulrike Zabel, Eva Szarek, Paraskevi Salpea, Antonella Forlino, Annalisa Vetro, Orsetta Zuffardi, Caroline Kisker, Susanne Diener, Thomas Meitinger, Martin J. Lohse, Martin Reincke, Jerome Bertherat, Tim M. Strom, Bruno Allolio. Constitutive activation of PRKACA in adrenal Cushing's syndrome. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr LB-182. doi:10.1158/1538-7445.AM2014-LB-182

Exome sequencing identifies compound heterozygous mutations in C12orf57 in two siblings with severe intellectual disability, hypoplasia of the corpus callosum, chorioretinal coloboma, and intractable seizures

In patients with genetically heterogeneous disorders such as intellectual disability or epilepsy, exome sequencing is a powerful tool to elucidate the underlying genetic cause. Homozygous and compound heterozygous mutations in C12orf57 have recently been described to cause an autosomal recessive syndromic form of intellectual disability, including agenesis/hypoplasia of the corpus callosum, optic coloboma, and intractable seizures. Here, we report on two siblings from nonconsanguineous parents harboring two compound heterozygous loss-of-function mutations in C12orf57 identified by exome sequencing, including a novel nonsense mutation, and review the patients described in the literature.

Constitutive activation of PKA catalytic subunit in adrenal Cushing's syndrome

BACKGROUND

Corticotropin-independent Cushing's syndrome is caused by tumors or hyperplasia of the adrenal cortex. The molecular pathogenesis of cortisol-producing adrenal adenomas is not well understood.

METHODS

We performed exome sequencing of tumor-tissue specimens from 10 patients with cortisol-producing adrenal adenomas and evaluated recurrent mutations in candidate genes in an additional 171 patients with adrenocortical tumors. We also performed genomewide copy-number analysis in 35 patients with cortisol-secreting bilateral adrenal hyperplasias. We studied the effects of these genetic defects both clinically and in vitro.

RESULTS

Exome sequencing revealed somatic mutations in PRKACA, which encodes the catalytic subunit of cyclic AMP-dependent protein kinase (protein kinase A [PKA]), in 8 of 10 adenomas (c.617A→C in 7 and c.595_596insCAC in 1). Overall, PRKACA somatic mutations were identified in 22 of 59 unilateral adenomas (37%) from patients with overt Cushing's syndrome; these mutations were not detectable in 40 patients with subclinical hypercortisolism or in 82 patients with other adrenal tumors. Among 35 patients with cortisol-producing hyperplasias, 5 (including 2 first-degree relatives) carried a germline copy-number gain (duplication) of the genomic region on chromosome 19 that includes PRKACA. In vitro studies showed impaired inhibition of both PKA catalytic subunit mutants by the PKA regulatory subunit, whereas cells from patients with germline chromosomal gains showed increased protein levels of the PKA catalytic subunit; in both instances, basal PKA activity was increased.

CONCLUSIONS

Genetic alterations of the catalytic subunit of PKA were found to be associated with human disease. Germline duplications of this gene resulted in bilateral adrenal hyperplasias, whereas somatic PRKACA mutations resulted in unilateral cortisol-producing adrenal adenomas. (Funded by the European Commission Seventh Framework Program and others.).

Exome sequence reveals mutations in CoA synthase as a cause of neurodegeneration with brain iron accumulation

Neurodegeneration with brain iron accumulation (NBIA) comprises a clinically and genetically heterogeneous group of disorders with progressive extrapyramidal signs and neurological deterioration, characterized by iron accumulation in the basal ganglia. Exome sequencing revealed the presence of recessive missense mutations in COASY, encoding coenzyme A (CoA) synthase in one NBIA-affected subject. A second unrelated individual carrying mutations in COASY was identified by Sanger sequence analysis. CoA synthase is a bifunctional enzyme catalyzing the final steps of CoA biosynthesis by coupling phosphopantetheine with ATP to form dephospho-CoA and its subsequent phosphorylation to generate CoA. We demonstrate alterations in RNA and protein expression levels of CoA synthase, as well as CoA amount, in fibroblasts derived from the two clinical cases and in yeast. This is the second inborn error of coenzyme A biosynthesis to be implicated in NBIA.

Transcriptome and genome sequencing uncovers functional variation in humans

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of mRNA and miRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project – the first uniformly processed RNA-seq data from multiple human populations with high-quality genome sequences. We discovered extremely widespread genetic variation affecting regulation of the majority of genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on cellular mechanisms of regulatory and loss-of-function variation, and allowed us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Exome sequencing reveals de novo WDR45 mutations causing a phenotypically distinct, X-linked dominant form of NBIA

Neurodegeneration with brain iron accumulation (NBIA) is a group of genetic disorders characterized by abnormal iron deposition in the basal ganglia. We report that de novo mutations in WDR45, a gene located at Xp11.23 and encoding a beta-propeller scaffold protein with a putative role in autophagy, cause a distinctive NBIA phenotype. The clinical features include early-onset global developmental delay and further neurological deterioration (parkinsonism, dystonia, and dementia developing by early adulthood). Brain MRI revealed evidence of iron deposition in the substantia nigra and globus pallidus. Males and females are phenotypically similar, an observation that might be explained by somatic mosaicism in surviving males and germline or somatic mutations in females, as well as skewing of X chromosome inactivation. This clinically recognizable disorder is among the more common forms of NBIA, and we suggest that it be named accordingly as beta-propeller protein-associated neurodegeneration.

Range of genetic mutations associated with severe non-syndromic sporadic intellectual disability: an exome sequencing study

BACKGROUND

The genetic cause of intellectual disability in most patients is unclear because of the absence of morphological clues, information about the position of such genes, and suitable screening methods. Our aim was to identify de-novo variants in individuals with sporadic non-syndromic intellectual disability.

METHODS

In this study, we enrolled children with intellectual disability and their parents from ten centres in Germany and Switzerland. We compared exome sequences between patients and their parents to identify de-novo variants. 20 children and their parents from the KORA Augsburg Diabetes Family Study were investigated as controls.

FINDINGS

We enrolled 51 participants from the German Mental Retardation Network. 45 (88%) participants in the case group and 14 (70%) in the control group had de-novo variants. We identified 87 de-novo variants in the case group, with an exomic mutation rate of 1·71 per individual per generation. In the control group we identified 24 de-novo variants, which is 1·2 events per individual per generation. More participants in the case group had loss-of-function variants than in the control group (20/51 vs 2/20; p=0·022), suggesting their contribution to disease development. 16 patients carried de-novo variants in known intellectual disability genes with three recurrently mutated genes (STXBP1, SYNGAP1, and SCN2A). We deemed at least six loss-of-function mutations in six novel genes to be disease causing. We also identified several missense alterations with potential pathogenicity.

INTERPRETATION

After exclusion of copy-number variants, de-novo point mutations and small indels are associated with severe, sporadic non-syndromic intellectual disability, accounting for 45-55% of patients with high locus heterogeneity. Autosomal recessive inheritance seems to contribute little in the outbred population investigated. The large number of de-novo variants in known intellectual disability genes is only partially attributable to known non-specific phenotypes. Several patients did not meet the expected syndromic manifestation, suggesting a strong bias in present clinical syndrome descriptions.

FUNDING

German Ministry of Education and Research, European Commission 7th Framework Program, and Swiss National Science Foundation.