Dissecting the treatment-naive ecosystem of human melanoma brain metastasis.

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

Spatial multi-omic map of human myocardial infarction.

Myocardial infarction is a leading cause of death worldwide1. Although advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality2. Here we generate an integrative high-resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from patients with myocardial infarction and controls. Multi-modal data integration enabled us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell states of major cell types and analysed them in their spatial context, evaluating their dependency on other cell types. Our data elucidate the molecular principles of human myocardial tissue organization, recapitulating a gradual cardiomyocyte and myeloid continuum following ischaemic injury. In sum, our study provides an integrative molecular map of human myocardial infarction, represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.

Multimodal CRISPR perturbations of GWAS loci associated with coronary artery disease in vascular endothelial cells.

Genome-wide association studies have identified >250 genetic variants associated with coronary artery disease (CAD), but the causal variants, genes and molecular mechanisms remain unknown at most loci. We performed pooled CRISPR screens to test the impact of sequences at or near CAD-associated genetic variants on vascular endothelial cell functions. Using CRISPR knockout, inhibition and activation, we targeted 1998 variants at 83 CAD loci to assess their effect on three adhesion proteins (E-selectin, ICAM1, VCAM1) and three key endothelial functions (nitric oxide and reactive oxygen species production, calcium signalling). At a false discovery rate ≤10%, we identified significant CRISPR perturbations near 42 variants located within 26 CAD loci. We used base editing to validate a putative causal variant in the promoter of the FES gene. Although a few of the loci include genes previously characterized in endothelial cells (e.g. AIDA, ARHGEF26, ADAMTS7), most are implicated in endothelial dysfunction for the first time. Detailed characterization of one of these new loci implicated the RNA helicase DHX38 in vascular endothelial cell senescence. While promising, our results also highlighted several limitations in using CRISPR perturbations to functionally dissect GWAS loci, including an unknown false negative rate and potential off-target effects.

Cells and Molecules Underpinning Cannabis-Related Variations in Cortical Thickness during Adolescence.

During adolescence, cannabis experimentation is common, and its association with interindividual variations in brain maturation well studied. Cellular and molecular underpinnings of these system-level relationships are, however, unclear. We thus conducted a three-step study. First, we exposed adolescent male mice to Δ-9-tetrahydrocannabinol (THC) or a synthetic cannabinoid WIN 55,212-2 (WIN) and assessed differentially expressed genes (DEGs), spine numbers, and dendritic complexity in their frontal cortex. Second, in human (male) adolescents, we examined group differences in cortical thickness in 34 brain regions, using magnetic resonance imaging, between those who experimented with cannabis before age 16 (n = 140) and those who did not (n = 327). Finally, we correlated spatially these group differences with gene expression of human homologs of mouse-identified DEGs. The spatial expression of 13 THC-related human homologs of DEGs correlated with cannabis-related variations in cortical thickness, and virtual histology revealed coexpression patterns of these 13 genes with cell-specific markers of astrocytes, microglia, and a type of pyramidal cells enriched in dendrite-regulating genes. Similarly, the spatial expression of 18 WIN-related human homologs of DEGs correlated with group differences in cortical thickness and showed coexpression patterns with the same three cell types. Gene ontology analysis indicated that 37 THC-related human homologs are enriched in neuron projection development, while 33 WIN-related homologs are enriched in processes associated with learning and memory. In mice, we observed spine loss and lower dendritic complexity in pyramidal cells of THC-exposed animals (vs controls). Experimentation with cannabis during adolescence may influence cortical thickness by impacting glutamatergic synapses and dendritic arborization.

A 4-lineage Statistical Suite to Evaluate the Support of Large-Scale Retrotransposon Insertion Data to Reconstruct Evolutionary Trees.

Retrophylogenomics makes use of genome-wide retrotransposon presence/absence insertion patterns to resolve questions in phylogeny and population genetics. In the genomics era, evaluating high-throughput data requires the associated development of appropriately powerful statistical tools. The currently used KKSC 3-lineage statistical test for estimating the significance of retrophylogenomic data is limited by the number of possible tree topologies it can assess in one step. To improve on this, we have extended the analysis to simultaneously compare four lineages, enabling us to evaluate ten distinct presence/absence insertion patterns for 26 possible tree topologies plus 129 trees with different incidences of hybridization or introgression. The new tool provides statistics for cases involving multiple ancestral hybridizations/introgressions, ancestral incomplete lineage sorting, bifurcation, and polytomy. The test is embedded in a user-friendly web R application (http://retrogenomics.uni-muenster.de:3838/hammlet/) and is available for use by the scientific community. [ancestral hybridization/introgression; ancestral incomplete lineage sorting (ILS); empirical distribution; KKSC-statistics; 4-lineage (4-LIN) insertion polymorphism; polytomy; retrophylogenomics.].

Acan downregulation in parvalbumin GABAergic cells reduces spontaneous recovery of fear memories.

While persistence of fear memories is essential for survival, a failure to inhibit fear in response to harmless stimuli is a feature of anxiety disorders. Extinction training only temporarily suppresses fear memory recovery in adults, but it is highly effective in juvenile rodents. Maturation of GABAergic circuits, in particular of parvalbumin-positive (PV+) cells, restricts plasticity in the adult brain, thus reducing PV+ cell maturation could promote the suppression of fear memories following extinction training in adults. Epigenetic modifications such as histone acetylation control gene accessibility for transcription and help couple synaptic activity to changes in gene expression. Histone deacetylase 2 (Hdac2), in particular, restrains both structural and functional synaptic plasticity. However, whether and how Hdac2 controls the maturation of postnatal PV+ cells is not well understood. Here, we show that PV+- cell specific Hdac2 deletion limits spontaneous fear memory recovery in adult mice, while enhancing PV+ cell bouton remodeling and reducing perineuronal net aggregation around PV+ cells in prefrontal cortex and basolateral amygdala. Prefrontal cortex PV+ cells lacking Hdac2, show reduced expression of Acan, a critical perineuronal net component, which is rescued by Hdac2 re-expression. Pharmacological inhibition of Hdac2 before extinction training is sufficient to reduce both spontaneous fear memory recovery and Acan expression in wild-type adult mice, while these effects are occluded in PV+-cell specific Hdac2 conditional knockout mice. Finally, a brief knock-down of Acan expression mediated by intravenous siRNA delivery before extinction training but after fear memory acquisition is sufficient to reduce spontaneous fear recovery in wild-type mice. Altogether, these data suggest that controlled manipulation of PV+ cells by targeting Hdac2 activity, or the expression of its downstream effector Acan, promotes the long-term efficacy of extinction training in adults.

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Family Based Whole Exome Sequencing Reveals the Multifaceted Role of Notch Signaling in Congenital Heart Disease.

Left-ventricular outflow tract obstructions (LVOTO) encompass a wide spectrum of phenotypically heterogeneous heart malformations which frequently cluster in families. We performed family based whole-exome and targeted re-sequencing on 182 individuals from 51 families with multiple affected members. Central to our approach is the family unit which serves as a reference to identify causal genotype-phenotype correlations. Screening a multitude of 10 overlapping phenotypes revealed disease associated and co-segregating variants in 12 families. These rare or novel protein altering mutations cluster predominantly in genes (NOTCH1, ARHGAP31, MAML1, SMARCA4, JARID2, JAG1) along the Notch signaling cascade. This is in line with a significant enrichment (Wilcoxon, p< 0.05) of variants with a higher pathogenicity in the Notch signaling pathway in patients compared to controls. The significant enrichment of novel protein truncating and missense mutations in NOTCH1 highlights the allelic and phenotypic heterogeneity in our pediatric cohort. We identified novel co-segregating pathogenic mutations in NOTCH1 associated with left and right-sided cardiac malformations in three independent families with a total of 15 affected individuals. In summary, our results suggest that a small but highly pathogenic fraction of family specific mutations along the Notch cascade are a common cause of LVOTO.

Cell type mapping reveals tissue niches and interactions in subcortical multiple sclerosis lesions.

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system. Inflammation is gradually compartmentalized and restricted to specific tissue niches such as the lesion rim. However, the precise cell type composition of such niches, their interactions and changes between chronic active and inactive stages are incompletely understood. We used single-nucleus and spatial transcriptomics from subcortical MS and corresponding control tissues to map cell types and associated pathways to lesion and nonlesion areas. We identified niches such as perivascular spaces, the inflamed lesion rim or the lesion core that are associated with the glial scar and a cilia-forming astrocyte subtype. Focusing on the inflamed rim of chronic active lesions, we uncovered cell-cell communication events between myeloid, endothelial and glial cell types. Our results provide insight into the cellular composition, multicellular programs and intercellular communication in tissue niches along the conversion from a homeostatic to a dysfunctional state underlying lesion progression in MS.

Triglyceride-derived fatty acids reduce autophagy in a model of retinal angiomatous proliferation.

Dyslipidemia and autophagy have been implicated in the pathogenesis of blinding neovascular age-related macular degeneration (NV-AMD). VLDL receptor (VLDLR), expressed in photoreceptors with a high metabolic rate, facilitates the uptake of triglyceride-derived fatty acids. Since fatty acid uptake is reduced in Vldlr-/- tissues, more remain in circulation, and the retina is fuel deficient, driving the formation in mice of neovascular lesions reminiscent of retinal angiomatous proliferation (RAP), a subtype of NV-AMD. Nutrient scarcity and energy failure are classically mitigated by increasing autophagy. We found that excess circulating lipids restrained retinal autophagy, which contributed to pathological angiogenesis in the Vldlr-/- RAP model. Triglyceride-derived fatty acid sensed by free fatty acid receptor 1 (FFAR1) restricted autophagy and oxidative metabolism in photoreceptors. FFAR1 suppressed transcription factor EB (TFEB), a master regulator of autophagy and lipid metabolism. Reduced TFEB, in turn, decreased sirtuin-3 expression and mitochondrial respiration. Metabolomic signatures of mouse RAP-like retinas were consistent with a role in promoting angiogenesis. This signature was also found in human NV-AMD vitreous. Restoring photoreceptor autophagy in Vldlr-/- retinas, either pharmacologically or by deleting Ffar1, enhanced metabolic efficiency and suppressed pathological angiogenesis. Dysregulated autophagy by circulating lipids might therefore contribute to the energy failure of photoreceptors driving neovascular eye diseases, and FFAR1 may be a target for intervention.

Loss of ADAMTS19 causes progressive non-syndromic heart valve disease.

Valvular heart disease is observed in approximately 2% of the general population1. Although the initial observation is often localized (for example, to the aortic or mitral valve), disease manifestations are regularly observed in the other valves and patients frequently require surgery. Despite the high frequency of heart valve disease, only a handful of genes have so far been identified as the monogenic causes of disease2-7. Here we identify two consanguineous families, each with two affected family members presenting with progressive heart valve disease early in life. Whole-exome sequencing revealed homozygous, truncating nonsense alleles in ADAMTS19 in all four affected individuals. Homozygous knockout mice for Adamts19 show aortic valve dysfunction, recapitulating aspects of the human phenotype. Expression analysis using a lacZ reporter and single-cell RNA sequencing highlight Adamts19 as a novel marker for valvular interstitial cells; inference of gene regulatory networks in valvular interstitial cells positions Adamts19 in a highly discriminatory network driven by the transcription factor lymphoid enhancer-binding factor 1 downstream of the Wnt signaling pathway. Upregulation of endocardial Krüppel-like factor 2 in Adamts19 knockout mice precedes hemodynamic perturbation, showing that a tight balance in the Wnt-Adamts19-Klf2 axis is required for proper valve maturation and maintenance.

Validation of Genome-Wide Polygenic Risk Scores for Coronary Artery Disease in French Canadians.

BACKGROUND: Coronary artery disease (CAD) represents one of the leading causes of morbidity and mortality worldwide. Given the healthcare risks and societal impacts associated with CAD, their clinical management would benefit from improved prevention and prediction tools. Polygenic risk scores (PRS) based on an individual's genome sequence are emerging as potentially powerful biomarkers to predict the risk to develop CAD. Two recently derived genome-wide PRS have shown high specificity and sensitivity to identify CAD cases in European-ancestry participants from the UK Biobank. However, validation of the PRS predictive power and transferability in other populations is now required to support their clinical utility. METHODS: We calculated both PRS (GPSCAD and metaGRSCAD) in French-Canadian individuals from 3 cohorts totaling 3639 prevalent CAD cases and 7382 controls and tested their power to predict prevalent, incident, and recurrent CAD. We also estimated the impact of the founder French-Canadian familial hypercholesterolemia deletion ( LDLR delta >15 kb deletion) on CAD risk in one of these cohorts and used this estimate to calibrate the impact of the PRS. RESULTS: Our results confirm the ability of both PRS to predict prevalent CAD comparable to the original reports (area under the curve=0.72-0.89). Furthermore, the PRS identified about 6% to 7% of individuals at CAD risk similar to carriers of the LDLR delta >15 kb mutation, consistent with previous estimates. However, the PRS did not perform as well in predicting an incident or recurrent CAD (area under the curve=0.56-0.60), maybe because of confounding because 76% of the participants were on statin treatment. This result suggests that additional work is warranted to better understand how ascertainment biases and study design impact PRS for CAD. CONCLUSIONS: Collectively, our results confirm that novel, genome-wide PRS is able to predict CAD in French Canadians; with further improvements, this is likely to pave the way towards more targeted strategies to predict and prevent CAD-related adverse events.

Copy number variation analysis in bicuspid aortic valve-related aortopathy identifies TBX20 as a contributing gene.

Bicuspid aortic valve (BAV) is the most common congenital heart defect (CHD), affecting 1-2% of the population. BAV is associated with thoracic aortic aneurysms (TAAs). Deleterious copy number variations (CNVs) were found previously in up to 10% of CHD cases. This study aimed at unravelling the contribution of deleterious deletions or duplications in 95 unrelated BAV/TAA patients. Seven unique or rare CNVs were validated, harbouring protein-coding genes with a role in the cardiovascular system. Based on the presence of overlapping CNVs in patients with cardiovascular phenotypes in the DECIPHER database, the identification of similar CNVs in whole-exome sequencing data of 67 BAV/TAA patients and suggested topological domain involvement from Hi-C data, supportive evidence was obtained for two genes (DGCR6 and TBX20) of the seven initially validated CNVs. A rare variant burden analysis using next-generation sequencing data from 637 BAV/TAA patients was performed for these two candidate genes. This revealed a suggestive genetic role for TBX20 in BAV/TAA aetiology, further reinforced by segregation of a rare TBX20 variant with the phenotype within a BAV/TAA family. To conclude, our results do not confirm a significant contribution for deleterious CNVs in BAV/TAA as only one potentially pathogenic CNV (1.05%) was identified. We cannot exclude the possibility that BAV/TAA is occasionally attributed to causal CNVs though, or that certain CNVs act as genetic risk factors by creating a sensitised background for BAV/TAA. Finally, accumulative evidence for TBX20 involvement in BAV/TAA aetiology underlines the importance of this transcription factor in cardiovascular disease.

ROBO4 variants predispose individuals to bicuspid aortic valve and thoracic aortic aneurysm.

Bicuspid aortic valve (BAV) is a common congenital heart defect (population incidence, 1-2%)1-3 that frequently presents with ascending aortic aneurysm (AscAA)4. BAV/AscAA shows autosomal dominant inheritance with incomplete penetrance and male predominance. Causative gene mutations (for example, NOTCH1, SMAD6) are known for ≤1% of nonsyndromic BAV cases with and without AscAA5-8, impeding mechanistic insight and development of therapeutic strategies. Here, we report the identification of variants in ROBO4 (which encodes a factor known to contribute to endothelial performance) that segregate with disease in two families. Targeted sequencing of ROBO4 showed enrichment for rare variants in BAV/AscAA probands compared with controls. Targeted silencing of ROBO4 or mutant ROBO4 expression in endothelial cell lines results in impaired barrier function and a synthetic repertoire suggestive of endothelial-to-mesenchymal transition. This is consistent with BAV/AscAA-associated findings in patients and in animal models deficient for ROBO4. These data identify a novel endothelial etiology for this common human disease phenotype.

Corrigendum: Candidate Gene Resequencing in a Large Bicuspid Aortic Valve-Associated Thoracic Aortic Aneurysm Cohort: SMAD6 as an Important Contributor.

[This corrects the article on p. 400 in vol. 8, PMID: 28659821.].

Candidate Gene Resequencing in a Large Bicuspid Aortic Valve-Associated Thoracic Aortic Aneurysm Cohort: SMAD6 as an Important Contributor.

Bicuspid aortic valve (BAV) is the most common congenital heart defect. Although many BAV patients remain asymptomatic, at least 20% develop thoracic aortic aneurysm (TAA). Historically, BAV-related TAA was considered as a hemodynamic consequence of the valve defect. Multiple lines of evidence currently suggest that genetic determinants contribute to the pathogenesis of both BAV and TAA in affected individuals. Despite high heritability, only very few genes have been linked to BAV or BAV/TAA, such as NOTCH1, SMAD6, and MAT2A. Moreover, they only explain a minority of patients. Other candidate genes have been suggested based on the presence of BAV in knockout mouse models (e.g., GATA5, NOS3) or in syndromic (e.g., TGFBR1/2, TGFB2/3) or non-syndromic (e.g., ACTA2) TAA forms. We hypothesized that rare genetic variants in these genes may be enriched in patients presenting with both BAV and TAA. We performed targeted resequencing of 22 candidate genes using Haloplex target enrichment in a strictly defined BAV/TAA cohort (n = 441; BAV in addition to an aortic root or ascendens diameter ≥ 4.0 cm in adults, or a Z-score ≥ 3 in children) and in a collection of healthy controls with normal echocardiographic evaluation (n = 183). After additional burden analysis against the Exome Aggregation Consortium database, the strongest candidate susceptibility gene was SMAD6 (p = 0.002), with 2.5% (n = 11) of BAV/TAA patients harboring causal variants, including two nonsense, one in-frame deletion and two frameshift mutations. All six missense mutations were located in the functionally important MH1 and MH2 domains. In conclusion, we report a significant contribution of SMAD6 mutations to the etiology of the BAV/TAA phenotype.

Aortic Dilatation Associated With a De Novo Mutation in the SOX18 Gene: Expanding the Clinical Spectrum of Hypotrichosis-Lymphedema-Telangiectasia Syndrome.

BACKGROUND: We report a 13-year-old female patient followed since birth for multiple rare congenital defects, including hypotrichosis, telangiectasia, and severe dilatation of the ascending aorta. METHODS: Comprehensive phenotype assessment throughout childhood included repeated echocardiographic measurements, evaluation of renal function, and immunohistochemical analysis of skin biopsy samples. Whole-exome sequencing was performed for the patient and both unaffected parents. RESULTS: We identified a novel de novo mutation in the transcription factor SOX18 (c.481C>T:p.Gln161*) in the patient, which was absent in all unaffected family members. Echocardiography revealed early onset and progressive dilatation of the ascending aorta. Skin biopsy results confirmed the defects of the blood vasculature in the presence of intact lymphatic vessels. Assessment of renal function did not show any signs of renal problems or renal failure in the patient. CONCLUSIONS: The genetic finding of a pathogenic SOX18 mutation enabled the diagnosis of the rare hypotrichosis-lymphedema-telangiectasia syndrome in our patient. The identification of a novel stop gain mutation in the SOX18 gene in association with dilatation of the aorta highlights the importance of this gene during the development of the circulatory system. Our study highlights the importance of whole-exome sequencing in the rapid identification of genes and gene mutations involved in rare conditions and thus expanding the knowledge and spectrum of clinical manifestations associated with them.

Mutations in SGOL1 cause a novel cohesinopathy affecting heart and gut rhythm.

The pacemaking activity of specialized tissues in the heart and gut results in lifelong rhythmic contractions. Here we describe a new syndrome characterized by Chronic Atrial and Intestinal Dysrhythmia, termed CAID syndrome, in 16 French Canadians and 1 Swede. We show that a single shared homozygous founder mutation in SGOL1, a component of the cohesin complex, causes CAID syndrome. Cultured dermal fibroblasts from affected individuals showed accelerated cell cycle progression, a higher rate of senescence and enhanced activation of TGF-ß signaling. Karyotypes showed the typical railroad appearance of a centromeric cohesion defect. Tissues derived from affected individuals displayed pathological changes in both the enteric nervous system and smooth muscle. Morpholino-induced knockdown of sgol1 in zebrafish recapitulated the abnormalities seen in humans with CAID syndrome. Our findings identify CAID syndrome as a novel generalized dysrhythmia, suggesting a new role for SGOL1 and the cohesin complex in mediating the integrity of human cardiac and gut rhythm.