Single-nucleus DNA sequencing reveals hidden somatic loss-of-heterozygosity in Cerebral Cavernous Malformations

Cerebral Cavernous Malformations (CCMs) are vascular malformations of the central nervous system which can lead to moderate to severe neurological phenotypes in patients. A majority of CCM lesions are driven by a cancer-like three-hit mutational mechanism, including a somatic, activating mutation in the oncogene PIK3CA, as well as biallelic loss-of-function mutations in a CCM gene. However, standard sequencing approaches often fail to yield a full complement of pathogenic mutations in many CCMs. We suggest this reality reflects the limited sensitivity to identify low-frequency variants and the presence of mutations undetectable with bulk short-read sequencing. Here we report a single-nucleus DNA-sequencing approach that leverages the underlying biology of CCMs to identify lesions with somatic loss-of-heterozygosity, a class of such hidden mutations. We identify an alternative genetic mechanism for CCM pathogenesis and establish a method that can be repurposed to investigate the genetic underpinning of other disorders with multiple somatic mutations.

Gonomics: uniting high performance and readability for genomics with Go

SUMMARY

Many existing software libraries for genomics require researchers to pick between competing considerations: the performance of compiled languages and the accessibility of interpreted languages. Go, a modern compiled language, provides an opportunity to address this conflict. We introduce Gonomics, an open-source collection of command line programs and bioinformatic libraries implemented in Go that unites readability and performance for genomic analyses. Gonomics contains packages to read, write, and manipulate a wide array of file formats (e.g. FASTA, FASTQ, BED, BEDPE, SAM, BAM, and VCF), and can convert and interface between these formats. Furthermore, our modular library structure provides a flexible platform for researchers developing their own software tools to address specific questions. These commands can be combined and incorporated into complex pipelines to meet the growing need for high-performance bioinformatic resources.

AVAILABILITY AND IMPLEMENTATION

Gonomics is implemented in the Go programming language. Source code, installation instructions, and documentation are freely available at https://github.com/vertgenlab/gonomics.

Adaptive sequence divergence forged new neurodevelopmental enhancers in humans

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

A novel somatic mutation in GNAQ in a capillary malformation provides insight into molecular pathogenesis

Sturge-Weber syndrome (SWS) is a sporadic, congenital, neuro-cutaneous disorder characterized by a mosaic, capillary malformation. SWS and non-syndromic capillary malformations are both caused by a somatic activating mutation in GNAQ encoding the G protein subunit alpha-q protein. The missense mutation R183Q is the sole GNAQ mutation identified thus far in 90% of SWS-associated or isolated capillary malformations. In this study, we sequenced skin biopsies of capillary malformations from 9 patients. We identified the R183Q mutation in nearly all samples, but one sample exhibited a Q209R mutation. This new mutation occurs at the same residue as the constitutively-activating Q209L mutation, commonly seen in tumors. However, Q209R is a rare variant in this gene. To compare the effect of the Q209R mutation on downstream signaling, we performed reporter assays with a GNAQ-responsive reporter co-transfected with either GNAQ WT, R183Q, Q209L, Q209R, or C9X (representing a null allele). Q209L showed the highest reporter activation, with R183Q and Q209R showing significantly lower activation. To determine whether these mutations had similar or different downstream consequences we performed RNA-seq analysis in microvascular endothelial cells (HMEC-1) electroporated with the same GNAQ variants. The R183 and Q209 missense variants caused extensive dysregulation of a broad range of transcripts compared to the WT or null allele, confirming that these are all activating mutations. However, the missense variants exhibited very few differentially expressed genes (DEGs) when compared to each other. These data suggest that these activating GNAQ mutations differ in magnitude of activation but have similar downstream effects.

Developmental venous anomalies are a genetic primer for cerebral cavernous malformations

Cerebral cavernous malformations (CCM) are a neurovascular anomaly that may occur sporadically, or be inherited due to autosomal dominant mutations in KRIT1 , CCM2 , or PDCD10 . Individual lesions are caused by somatic mutations which have been identified in KRIT1, CCM2, PDCD10, MAP3K3, and PIK3CA . However, the interactions between mutations, and their relative contributions to sporadic versus familial cases remain unclear. We show that mutations in KRIT1, CCM2, PDCD10, and MAP3K3 are mutually exclusive, but may co-occur with mutations in PIK3CA. We also find that MAP3K3 mutations may cause sporadic, but not familial CCM. Furthermore, we find identical PIK3CA mutations in CCMs and adjacent developmental venous anomalies (DVA), a common vascular malformation frequently found in the vicinity of sporadic CCMs. However, somatic mutations in MAP3K3 are found only in the CCM. This suggests that sporadic CCMs are derived from cells of the DVA which have acquired an additional mutation in MAP3K3 .

Cerebral Cavernous Malformation: From Mechanism to Therapy

Cerebral cavernous malformations are acquired vascular anomalies that constitute a common cause of central nervous system hemorrhage and stroke. The past 2 decades have seen a remarkable increase in our understanding of the pathogenesis of this vascular disease. This new knowledge spans genetic causes of sporadic and familial forms of the disease, molecular signaling changes in vascular endothelial cells that underlie the disease, unexpectedly strong environmental effects on disease pathogenesis, and drivers of disease end points such as hemorrhage. These novel insights are the integrated product of human clinical studies, human genetic studies, studies in mouse and zebrafish genetic models, and basic molecular and cellular studies. This review addresses the genetic and molecular underpinnings of cerebral cavernous malformation disease, the mechanisms that lead to lesion hemorrhage, and emerging biomarkers and therapies for clinical treatment of cerebral cavernous malformation disease. It may also serve as an example for how focused basic and clinical investigation and emerging technologies can rapidly unravel a complex disease mechanism.

PIK3CA and CCM mutations fuel cavernomas through a cancer-like mechanism

Vascular malformations are thought to be monogenic disorders that result in dysregulated growth of blood vessels. In the brain, cerebral cavernous malformations (CCMs) arise owing to inactivation of the endothelial CCM protein complex, which is required to dampen the activity of the kinase MEKK31-4. Environmental factors can explain differences in the natural history of CCMs between individuals5, but why single CCMs often exhibit sudden, rapid growth, culminating in strokes or seizures, is unknown. Here we show that growth of CCMs requires increased signalling through the phosphatidylinositol-3-kinase (PI3K)-mTOR pathway as well as loss of function of the CCM complex. We identify somatic gain-of-function mutations in PIK3CA and loss-of-function mutations in the CCM complex in the same cells in a majority of human CCMs. Using mouse models, we show that growth of CCMs requires both PI3K gain of function and CCM loss of function in endothelial cells, and that both CCM loss of function and increased expression of the transcription factor KLF4 (a downstream effector of MEKK3) augment mTOR signalling in endothelial cells. Consistent with these findings, the mTORC1 inhibitor rapamycin effectively blocks the formation of CCMs in mouse models. We establish a three-hit mechanism analogous to cancer, in which aggressive vascular malformations arise through the loss of vascular 'suppressor genes' that constrain vessel growth and gain of a vascular 'oncogene' that stimulates excess vessel growth. These findings suggest that aggressive CCMs could be treated using clinically approved mTORC1 inhibitors.

Somatic Mutations in Vascular Malformations of Hereditary Hemorrhagic Telangiectasia Result in Bi-allelic Loss of ENG or ACVRL1

Hereditary hemorrhagic telangiectasia (HHT) is a Mendelian disease characterized by vascular malformations (VMs) including visceral arteriovenous malformations and mucosal telangiectasia. HHT is caused by loss-of-function (LoF) mutations in one of three genes, ENG, ACVRL1, or SMAD4, and is inherited as an autosomal-dominant condition. Intriguingly, the constitutional mutation causing HHT is present throughout the body, yet the multiple VMs in individuals with HHT occur focally, rather than manifesting as a systemic vascular defect. This disconnect between genotype and phenotype suggests that a local event is necessary for the development of VMs. We investigated the hypothesis that local somatic mutations seed the formation HHT-related telangiectasia in a genetic two-hit mechanism. We identified low-frequency somatic mutations in 9/19 telangiectasia through the use of next-generation sequencing. We established phase for seven of nine samples, which confirms that the germline and somatic mutations in all seven samples exist in trans configuration; this is consistent with a genetic two-hit mechanism. These combined data suggest that bi-allelic loss of ENG or ACVRL1 may be a required event in the development of telangiectasia, and that rather than haploinsufficiency, VMs in HHT are caused by a Knudsonian two-hit mechanism.

Cerebral Cavernous Malformations Develop Through Clonal Expansion of Mutant Endothelial Cells

RATIONALE

Vascular malformations arise in vessels throughout the entire body. Causative genetic mutations have been identified for many of these diseases; however, little is known about the mutant cell lineage within these malformations.

OBJECTIVE

We utilize an inducible mouse model of cerebral cavernous malformations (CCMs) coupled with a multicolor fluorescent reporter to visualize the contribution of mutant endothelial cells (ECs) to the malformation.

METHODS AND RESULTS

We combined a Ccm3 mouse model with the confetti fluorescent reporter to simultaneously delete Ccm3 and label the mutant EC with 1 of 4 possible colors. We acquired Z-series confocal images from serial brain sections and created 3-dimensional reconstructions of entire CCMs to visualize mutant ECs during CCM development. We observed a pronounced pattern of CCMs lined with mutant ECs labeled with a single confetti color (n=42). The close 3-dimensional distribution, as determined by the nearest neighbor analysis, of the clonally dominant ECs within the CCM was statistically different than the background confetti labeling of ECs in non-CCM control brain slices as well as a computer simulation ( P<0.001). Many of the small (<100 μm diameter) CCMs consisted, almost exclusively, of the clonally dominant mutant ECs labeled with the same confetti color, whereas the large (>100 μm diameter) CCMs contained both the clonally dominant mutant cells and wild-type ECs. We propose of model of CCM development in which an EC acquires a second somatic mutation, undergoes clonal expansion to initiate CCM formation, and then incorporates neighboring wild-type ECs to increase the size of the malformation.

CONCLUSIONS

This is the first study to visualize, with single-cell resolution, the clonal expansion of mutant ECs within CCMs. The incorporation of wild-type ECs into the growing malformation presents another series of cellular events whose elucidation would enhance our understanding of CCMs and may provide novel therapeutic opportunities.