The eye and brain are composed of elaborately organized tissues, development of which is supported by spatiotemporally precise expression of a number of transcription factors and developmental regulators. Here we report the molecular and genetic characterization of Integrator complex subunit 15 (INTS15). INTS15 was identified in search for the causative gene(s) for an autosomal-dominant eye disease with variable individual manifestation found in a large pedigree. While homozygous Ints15 knockout mice are embryonic lethal, mutant mice lacking a small C-terminal region of Ints15 show ocular malformations similar to the human patients. INTS15 is highly expressed in the eye and brain during embryogenesis and stably interacts with the Integrator complex to support small nuclear RNA 3' end processing. Its knockdown resulted in missplicing of a large number of genes, probably as a secondary consequence, and substantially affected genes associated with eye and brain development. Moreover, studies using human iPS cells-derived neural progenitor cells showed that INTS15 is critical for axonal outgrowth in retinal ganglion cells. This study suggests a new link between general transcription machinery and a highly specific hereditary disease.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

In humans, mutations in calmodulin cause cardiac arrhythmia. These mutations disrupt the ability of calmodulin to sense calcium concentrations and correctly regulate two central calcium channels, together obstructing heart rhythm. This correlation is well established, but also surprising since calmodulin is expressed in all tissues and interacts with hundreds of proteins. Until now, most studies have focused on cardiac cell function and regulation of specific cardiac targets, and thus, potential other effects of these mutations have largely been unexplored. Here, we introduce the nematode Caenorhabditis elegans as an in vivo model to study effects of three human calmodulin mutations with different impairment on calcium binding. We find that arrhythmic effects of the calmodulin mutations N54I and D96V can be recapitulated in disruption of two rhythmic behaviors, pharynx pumping and defecation motor program. Interestingly, we also find that these mutations affect neuronal function, but in different ways. Whereas D96V sensitizes signaling at the neuromuscular junction, N54I has a protective effect. The mutation N98S did not affect rhythmic behavior, but impaired chemosensing. Therefore, pathogenic calmodulin mutations act through different mechanisms in rhythmic behavior and neuronal function in C. elegans, emphasizing the strength of using live multicellular models. Finally, our results support the hypothesis that human calmodulin mutations could also contribute to neurological diseases.

BACKGROUND: To understand the shared genetic basis between colorectal cancer (CRC) and other cancers and identify potential pleiotropic loci for compensating the missing genetic heritability of CRC. METHODS: We conducted a systematic genome-wide pleiotropy scan to appraise associations between cancer-related genetic variants and CRC risk among European populations. Single nucleotide polymorphism (SNP)-set analysis was performed using data from the UK Biobank and the Study of Colorectal Cancer in Scotland (10 039 CRC cases and 30 277 controls) to evaluate the overlapped genetic regions for susceptibility of CRC and other cancers. The variant-level pleiotropic associations between CRC and other cancers were examined by CRC genome-wide association study meta-analysis and the pleiotropic analysis under composite null hypothesis (PLACO) pleiotropy test. Gene-based, co-expression and pathway enrichment analyses were performed to explore potential shared biological pathways. The interaction between novel genetic variants and common environmental factors was further examined for their effects on CRC. RESULTS: Genome-wide pleiotropic analysis identified three novel SNPs (rs2230469, rs9277378 and rs143190905) and three mapped genes (PIP4K2A, HLA-DPB1 and RTEL1) to be associated with CRC. These genetic variants were significant expressions quantitative trait loci in colon tissue, influencing the expression of their mapped genes. Significant interactions of PIP4K2A and HLA-DPB1 with environmental factors, including smoking and alcohol drinking, were observed. All mapped genes and their co-expressed genes were significantly enriched in pathways involved in carcinogenesis. CONCLUSION: Our findings provide an important insight into the shared genetic basis between CRC and other cancers. We revealed several novel CRC susceptibility loci to help understand the genetic architecture of CRC.

Recessive variants in the oxidoreductase PYROXD1 are reported to cause a myopathy in 22 affected individuals from 15 families. Here, we describe two female probands from unrelated families presenting with features of a congenital connective tissue disorder including osteopenia, blue sclera, soft skin, joint hypermobility and neuromuscular junction dysfunction in addition to known features of PYROXD1 myopathy including respiratory difficulties, weakness, hypotonia and oromotor dysfunction. Proband AII:1 is compound heterozygous for the recurrent PYROXD1 variant Chr12(GRCh38):g.21452130A>G;NM_024854.5:c.464A>G;p.(N155S) and Chr12(GRCh38):g.21462019_21462022del;NM_024854.5:c.892_895del;p.(V298Mfs*4) and proband BII:1 is compound heterozygous for Chr12(GRCh38):g.21468739-21468741del;NM_024854.5:c.1488_1490del;p.(E496del) and Chr12(GRCh38):g.21467619del;NM_024854.5:c.1254+1del. RNA studies demonstrate c.892_895del;p.(V298Mfs*4) is targeted by nonsense mediated decay and c.1254+1delG elicits in-frame skipping of exon-11. Western blot from cultured fibroblasts shows reduced PYROXD1 protein levels in both probands. Testing urine from BII:1 and six individuals with PYROXD1 myopathy showed elevated levels of deoxypyridinoline, a mature collagen crosslink, correlating with PYROXD1-disorder severity. Urine and serum amino acid testing of the same individuals revealed no reportable changes. In contrast to PYROXD1 knock-out, we find no evidence for disrupted tRNA ligase activity, as measured via XBP1 splicing, in fibroblasts expressing PYROXD1 variants. In summary, we expand the clinical spectrum of PYROXD1-related disorders to include an overlapping connective tissue and myopathy presentation, identify three novel, pathogenic PYROXD1 variants, and provide preliminary evidence that elevated urine DPD crosslinks may provide a clinical biomarker for PYROXD1 disorders. Our results advocate consideration of PYROXD1 variants in the differential diagnosis for undiagnosed individuals presenting with a connective tissue disorder and myopathy.

Patient mutations have been identified throughout dynamin-related protein 1 (Drp1), the key protein mediator of mitochondrial fission. These changes generally impact young children and often result in severe neurological defects and, in some instances, death. Until now, the underlying functional defect leading to patient phenotypes has been largely speculative. We therefore analyzed six disease-associated mutations throughout the GTPase and middle domains (MD) of Drp1. The MD plays a role in Drp1 oligomerization, and three mutations in this region were predictably impaired in self-assembly. However, another mutant in this region (F370C) retained oligomerization capability on pre-curved membranes despite being assembly-limited in solution. Instead, this mutation impaired membrane remodeling of liposomes, which highlights the importance of Drp1 in generating local membrane curvature before fission. Two GTPase domain mutations were also observed in different patients. The G32A mutation was impaired in GTP hydrolysis both in solution and in the presence of lipid but remains capable of self-assembly on these lipid templates. The G223V mutation also exhibited decreased GTPase activity and was able to assemble on pre-curved lipid templates; however, this change impaired membrane remodeling of unilamellar liposomes similar to F370C. This demonstrates that the Drp1 GTPase domain also contributes to self-assembly interactions that drive membrane curvature. Overall, the functional defects caused by mutations in Drp1 are highly variable even for mutations that reside within the same functional domain. This study provides a framework for characterizing additional Drp1 mutations to provide a comprehensive understanding of functional sites within this essential protein.

Facial dysmorphology is a hallmark of 22q11.2 deletion syndrome (22q11DS). Nearly all affected individuals have facial features characteristic of the syndrome: a vertically long face with broad nasal bridge, narrow palpebral fissures and mild micrognathia, sometimes accompanied by facial skeletal and oropharyngeal anomalies. Despite the frequency of craniofacial dysmorphology due to 22q11.2 deletion, there is still incomplete understanding of the contribution of individual 22q11 genes to craniofacial and oropharyngeal development. We asked whether homozygous or heterozygous loss of function of single 22q11 genes compromises craniofacial and/or oropharyngeal morphogenesis related to these 22q11DS phenotypes. We found that Ranbp1, a 22q11DS gene that mediates nucleocytoplasmic protein trafficking, is a dosage-dependent modulator of craniofacial development. Ranbp1-/- embryos have variably penetrant facial phenotypes, including altered facial morphology and cleft palate. This 22q11DS-related dysmorphology is particularly evident in the midline of the facial skeleton, as evidenced by a robustly quantifiable dysmorphology of the vomer, an unpaired facial midline bone. 22q11DS-related oropharyngeal phenotypes reflect Ranbp1 function in both the cranial neural crest and cranial ectoderm based upon tissue-selective Ranbp1 deletion. Analyses of genetic interaction show that Ranbp1 mutation disrupts BMP signaling-dependent midline gene expression and BMP-mediated craniofacial and cranial skeletal morphogenesis. Finally, midline defects that parallel those in Ranbp1 mutant mice are observed at similar frequencies in the LgDel 22q112DS mouse model. Apparently, Ranbp1 is a modulator of craniofacial development, and in the context of broader 22q11 deletion, Ranbp1 mutant phenotypes mirror key aspects of 22q11DS midline facial anomalies.

Barth syndrome is an X-linked disorder caused by loss-of-function mutations in Tafazzin (TAZ), an acyltransferase that catalyzes remodeling of cardiolipin, a signature phospholipid of the inner mitochondrial membrane. Patients develop cardiac and skeletal muscle weakness, growth delay and neutropenia, although phenotypic expression varies considerably between patients. Taz knockout mice recapitulate many of the hallmark features of the disease. We used mouse genetics to test the hypothesis that genetic modifiers alter the phenotypic manifestations of Taz inactivation. We crossed TazKO/X females in the C57BL6/J inbred strain to males from eight inbred strains and evaluated the phenotypes of first-generation (F1) TazKO/Y progeny, compared to TazWT/Y littermates. We observed that genetic background strongly impacted phenotypic expression. C57BL6/J and CAST/EiJ[F1] TazKO/Y mice developed severe cardiomyopathy, whereas A/J[F1] TazKO/Y mice had normal heart function. C57BL6/J and WSB/EiJ[F1] TazKO/Y mice had severely reduced treadmill endurance, whereas endurance was normal in A/J[F1] and CAST/EiJ[F1] TazKO/Y mice. In all genetic backgrounds, cardiolipin showed similar abnormalities in knockout mice, and transcriptomic and metabolomic investigations identified signatures of mitochondrial uncoupling and activation of the integrated stress response. TazKO/Y cardiac mitochondria were small, clustered and had reduced cristae density in knockouts in severely affected genetic backgrounds but were relatively preserved in the permissive A/J[F1] strain. Gene expression and mitophagy measurements were consistent with reduced mitophagy in knockout mice in genetic backgrounds intolerant of Taz mutation. Our data demonstrate that genetic modifiers powerfully modulate phenotypic expression of Taz loss-of-function and act downstream of cardiolipin, possibly by altering mitochondrial quality control.

Genome-wide association studies (GWAS) have successfully identified associations for cervical cancer, but the underlying mechanisms of cervical biology and pathology remain uncharacterised. Our GWAS meta-analyses fill this gap, as we characterise the genetic architecture of cervical phenotypes, including cervical ectropion, cervicitis, cervical dysplasia, as well as up to 9229 cases and 490 304 controls for cervical cancer from diverse ancestries. Leveraging the latest computational methods and gene expression data, we refine the association signals for cervical cancer and propose potential causal variants and genes at each locus. We prioritise PAX8/PAX8-AS1, LINC00339, CDC42, CLPTM1L, HLA-DRB1 and GSDMB as the most likely candidate genes for cervical cancer signals, providing insights into cervical cancer pathogenesis and supporting the involvement of reproductive tract development, immune response and cellular proliferation/apoptosis. We construct a genetic risk score (GRS) that is associated with cervical cancer [hazard ratios (HR) = 3.1 (1.7-5.6) for the top 15% vs lowest 15% of individuals], and with other HPV- and immune-system-related diagnoses in a phenome-wide association study analysis. Our results propose valuable leads for further functional studies and present a GRS for cervical cancer that allows additional risk stratification and could potentially be used to personalise the conventional screening strategies for groups more susceptible to cervical cancer.

Zinc is an essential trace mineral. Dietary zinc deficiency results in stunted growth, skin lesions, hypogonadism and frequent infections in humans. Mice genetically lacking Slc30a7 suffer from mild zinc deficiency and are prone to development of prostate cancer and insulin resistance. Disease-causing variants or mutations in the human SLC30A7 (ZNT7) gene have not been previously reported. Here, we describe two-boy siblings from a French family with stunted growth, testicular hypoplasia and bone marrow failure. Exome sequencing revealed compound heterozygous variants in ZNT7 consisting of NM_133496.5:c.21dup; p.Asp8ArgfsTer3 and c.842 + 15 T > C inherited from their unaffected mother and father, respectively. The c.21dup variant led to a premature stop codon generated in exon 1 of the ZNT7 coding sequence. RNA-seq analysis demonstrated that the c.842 + 15 T > C variant resulted in a leaky mRNA splicing event generating a premature stop codon right after the splicing donor site of exon 8. Moreover, the expression of ZNT7 protein was remarkably reduced by 80-96% in the affected brothers compared to the control cells. These findings strongly suggest that biallelic variants in SLC30A7 should be considered as a cause of growth retardation, testicular hypoplasia and syndromic bone marrow failure.

SURF1 deficiency (OMIM # 220110) causes Leigh syndrome (LS, OMIM # 256000), a mitochondrial disorder typified by stress-induced metabolic strokes, neurodevelopmental regression and progressive multisystem dysfunction. Here, we describe two novel surf1-/- zebrafish knockout models generated by CRISPR/Cas9 technology. While gross larval morphology, fertility, and survival into adulthood appeared unaffected, surf1-/- mutants manifested adult-onset ocular anomalies and decreased swimming activity, as well as classical biochemical hallmarks of human SURF1 disease, including reduced complex IV expression and enzymatic activity and increased tissue lactate. surf1-/- larvae also demonstrated oxidative stress and stressor hypersensitivity to the complex IV inhibitor, azide, which exacerbated their complex IV deficiency, reduced supercomplex formation, and induced acute neurodegeneration typical of LS including brain death, impaired neuromuscular responses, reduced swimming activity, and absent heartrate. Remarkably, prophylactic treatment of surf1-/- larvae with either cysteamine bitartrate or N-acetylcysteine, but not other antioxidants, significantly improved animal resiliency to stressor-induced brain death, swimming and neuromuscular dysfunction, and loss of heartbeat. Mechanistic analyses demonstrated cysteamine bitartrate pretreatment did not improve complex IV deficiency, ATP deficiency, or increased tissue lactate but did reduce oxidative stress and restore glutathione balance in surf1-/- animals. Overall, two novel surf1-/- zebrafish models recapitulate the gross neurodegenerative and biochemical hallmarks of LS, including azide stressor hypersensitivity that was associated with glutathione deficiency and ameliorated by cysteamine bitartrate or N-acetylcysteine therapy.

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes that ligate tRNA molecules to cognate amino acids. Heterozygosity for missense variants or small in-frame deletions in six ARS genes causes dominant axonal peripheral neuropathy. These pathogenic variants reduce enzyme activity without significantly decreasing protein levels and reside in genes encoding homo-dimeric enzymes. These observations raise the possibility that neuropathy-associated ARS variants exert a dominant-negative effect, reducing overall ARS activity below a threshold required for peripheral nerve function. To test such variants for dominant-negative properties, we developed a humanized yeast assay to co-express pathogenic human alanyl-tRNA synthetase (AARS1) mutations with wild-type human AARS1. We show that multiple loss-of-function AARS1 mutations impair yeast growth through an interaction with wild-type AARS1, but that reducing this interaction rescues yeast growth. This suggests that neuropathy-associated AARS1 variants exert a dominant-negative effect, which supports a common, loss-of-function mechanism for ARS-mediated dominant peripheral neuropathy.

Pathogenic heterozygous variants in SCN2A, which encodes the neuronal sodium channel NaV1.2, cause different types of epilepsy or intellectual disability (ID)/autism without seizures. Previous studies using mouse models or heterologous systems suggest that NaV1.2 channel gain-of-function typically causes epilepsy, whereas loss-of-function leads to ID/autism. How altered channel biophysics translate into patient neurons remains unknown. Here, we investigated iPSC-derived early-stage cortical neurons from ID patients harboring diverse pathogenic SCN2A variants [p.(Leu611Valfs*35); p.(Arg937Cys); p.(Trp1716*)] and compared them with neurons from an epileptic encephalopathy (EE) patient [p.(Glu1803Gly)] and controls. ID neurons consistently expressed lower NaV1.2 protein levels. In neurons with the frameshift variant, NaV1.2 mRNA and protein levels were reduced by ~ 50%, suggesting nonsense-mediated decay and haploinsufficiency. In other ID neurons, only protein levels were reduced implying NaV1.2 instability. Electrophysiological analysis revealed decreased sodium current density and impaired action potential (AP) firing in ID neurons, consistent with reduced NaV1.2 levels. In contrast, epilepsy neurons displayed no change in NaV1.2 levels or sodium current density, but impaired sodium channel inactivation. Single-cell transcriptomics identified dysregulation of distinct molecular pathways including inhibition of oxidative phosphorylation in neurons with SCN2A haploinsufficiency and activation of calcium signaling and neurotransmission in epilepsy neurons. Together, our patient iPSC-derived neurons reveal characteristic sodium channel dysfunction consistent with biophysical changes previously observed in heterologous systems. Additionally, our model links the channel dysfunction in ID to reduced NaV1.2 levels and uncovers impaired AP firing in early-stage neurons. The altered molecular pathways may reflect a homeostatic response to NaV1.2 dysfunction and can guide further investigations.

OBJECTIVE: In Friedreich's ataxia (FRDA), the most affected tissues are not accessible to sampling and available transcriptomic findings originate from blood-derived cells and animal models. Herein, we aimed at dissecting for the first time the pathophysiology of FRDA by means of RNA-sequencing in an affected tissue sampled in vivo. METHODS: Skeletal muscle biopsies were collected from seven FRDA patients before and after treatment with recombinant human Erythropoietin (rhuEPO) within a clinical trial. Total RNA extraction, 3'-mRNA library preparation and sequencing were performed according to standard procedures. We tested for differential gene expression with DESeq2 and performed gene set enrichment analysis with respect to control subjects. RESULTS: FRDA transcriptomes showed 1873 genes differentially expressed from controls. Two main signatures emerged: (1) a global downregulation of the mitochondrial transcriptome as well as of ribosome/translational machinery and (2) an upregulation of genes related to transcription and chromatin regulation, especially of repressor terms. Downregulation of the mitochondrial transcriptome was more profound than previously shown in other cellular systems. Furthermore, we observed in FRDA patients a marked upregulation of leptin, the master regulator of energy homeostasis. RhuEPO treatment further enhanced leptin expression. INTERPRETATION: Our findings reflect a double hit in the pathophysiology of FRDA: a transcriptional/translational issue and a profound mitochondrial failure downstream. Leptin upregulation in the skeletal muscle in FRDA may represent a compensatory mechanism of mitochondrial dysfunction, which is amenable to pharmacological boosting. Skeletal muscle transcriptomics is a valuable biomarker to monitor therapeutic interventions in FRDA.

There is increasing evidence that exosome-mediated transmission of microRNAs helps to connect tumor-associated macrophages and cancer cells, including lung adenocarcinoma (LUAD) cells. The objective of this study is to identify the role of miR-3153 in LUAD progression and M2 macrophage polarization and explore its regulatory mechanism. The relevant molecular mechanisms were analyzed and validated through mechanistic assays. In vitro functional assays followed by in vivo experiments were implemented to evaluate the role of exosomes in mediating M2 macrophage polarization and LUAD progression. LUAD cells transmitted miR-3153 through exosomes. Heterogeneous nuclear ribonucleoprotein A2B1 promoted miR-3153 biosynthesis and exosomal sorting. Exosomal miR-3153 targeted zinc finger protein 91 to suppress the ubiquitination and degradation of misshapen-like kinase 1, thereby activating the c-Jun N-terminal kinase (JNK) signaling pathway and inducing M2 macrophage polarization. M2 macrophage polarization induced by LUAD cell-derived exosomes promoted the malignant process of LUAD cells. Transmission of exosomal miR-3153 by LUAD cells activates the JNK signaling pathway and induces M2 macrophage polarization, thus promoting the progression of LUAD.

The susceptibility single nucleotide polymorphisms (SNPs) obtained by genome-wide association studies leave some thorny questions, such as prioritization, false positives and unknown pathogenesis. Previous studies suggested that genetic variation may perturb the RNA secondary structure, influence protein recruitment and binding and ultimately affect splicing processes. Therefore, exploring the perturbation of SNPs to structure-function correlations may provide an effective bridge toward understanding the genetic contribution to diseases. Here, aiming to decipher the regulatory mechanism of myopia susceptibility variants, we systematically evaluated the roles of SNP-induced structural changes during splicing. In addition, 7.53% of myopia-related SNPs exhibited significant global structural changes, 19.53% presented noteworthy local structural disturbance and there were wide-ranging structural perturbations in the splice-related motifs. We established a comprehensive evaluation system for structural disturbance in the splicing-related motifs and gave the priority ranking for the SNPs at RNA structural level. These high-priority SNPs were revealed to widely disturb the molecular interaction properties between splicing-related proteins and pre-mRNAs by HDOCK. Moreover, mini-gene assays confirmed that structural perturbation could influence splicing efficiency through structural remodelling. This study deepens our understanding of the potential molecular regulatory mechanisms of susceptible SNPs in myopia and contributes to personalized diagnosis, personalized medicine, disease-risk prediction and functional verification study by guiding the prioritization of the susceptibility SNPs.

Cachexia occurrence and development are associated with loss of white adipose tissues, which may be involved with cancer-derived exosomes. This study attempted to characterize the functional mechanisms of breast cancer (BC) cell-derived exosome-loaded microRNA (miR)-155 in cancer cachexia-related fat loss. Exosomes were incubated with preadipocytes and cellular lipid droplet accumulation was observed using Oil Red O staining. Western blotting evaluated the cellular levels of lipogenesis marker peroxisome proliferator activated receptor gamma (PPARγ) and adiponectin, C1Q and collagen domain containing (AdipoQ). Differentiated adipocytes were incubated with exosomes, and phosphate hormone sensitive lipase (P-HSL), adipose triglyceride lipase (ATGL) and glycerol were detected in adipocytes, in addition to uncoupling protein 1 (UCP1) and leptin levels. A mouse model of cancer cachexia was established where cancer exosomes were injected intravenously. The changes in body weight and tumor-free body weights were recorded and serum glycerol levels and lipid accumulation in adipose tissues were determined. Also, the relationship between miR-155 and UBQLN1 was predicted and verified. BC exosome treatment reduced PPARγ and AdipoQ protein levels, promoted the levels of P-HSL and ATGL proteins, facilitated glycerol release, increased UCP1 expression and lowered leptin expression in adipocytes. Exosomal miR-155 inhibited lipogenesis in preadipocytes and boosted the browning of white adipose tissues. miR-155 downregulation alleviated cancer exosome-induced browning of white adipose tissues and fat loss. Mechanistically, miR-155 targeted UBQLN1, and UBQLN1 upregulation reversed the impacts of cancer exosomes. miR-155 loaded by BC cell-derived exosomes significantly affects white adipose browning and inhibition of cancer-derived exosomes.

Kabuki syndrome (KS) is a rare, multisystem disorder with a variable clinical phenotype. The majority of KS is caused by dominant loss-of-function mutations in KMT2D (lysine methyltransferase 2D). KMT2D mediates chromatin accessibility by adding methyl groups to lysine residue 4 of histone 3, which plays a critical role in cell differentiation and homeostasis. The molecular underpinnings of KS remain elusive partly because of a lack of histone modification data from human samples. Consequently, we profiled and characterized alterations in histone modification and gene transcription in peripheral blood mononuclear cells (PBMCs) from 33 patients with KMT2D mutations and 36 unaffected healthy controls. Our analysis identified unique enhancer signatures in H3K4me1 and H3K4me2 in KS compared with controls. Reduced enhancer signals were present for promoter-distal sites of immune-related genes for which co-binding of PBMC-specific transcription factors was predicted; 31% of super-enhancers of normal blood cells overlapped with disrupted enhancers in KS, supporting an association of reduced enhancer activity of immune-related genes with immune deficiency phenotypes. In contrast, increased enhancer signals were observed for promoter-proximal regions of metabolic genes enriched with EGR1 and E2F2 motifs, whose transcriptional levels were significantly increased in KS. Additionally, we identified ~100 de novo enhancers in genes, such as in MYO1F and AGAP2. Together, our results underscore the effect of KMT2D haploinsufficiency on dysregulation of enhancer states and gene transcription and provide a framework for the identification of therapeutic targets and biomarkers in preparation for clinical trial readiness.

SOX7 is a transcription factor-encoding gene located in a region on chromosome 8p23.1 that is recurrently deleted in individuals with ventricular septal defects (VSDs). We have previously shown that Sox7-/- embryos die of heart failure around E11.5. Here, we demonstrate that these embryos have hypocellular endocardial cushions with severely reduced numbers of mesenchymal cells. Ablation of Sox7 in the endocardium also resulted in hypocellular endocardial cushions, and we observed VSDs in rare E15.5 Sox7flox/-;Tie2-Cre and Sox7flox/flox;Tie2-Cre embryos that survived to E15.5. In atrioventricular explant studies, we showed that SOX7 deficiency leads to a severe reduction in endocardial-to-mesenchymal transition (EndMT). RNA-seq studies performed on E9.5 Sox7-/- heart tubes revealed severely reduced Wnt4 transcript levels. Wnt4 is expressed in the endocardium and promotes EndMT by acting in a paracrine manner to increase the expression of Bmp2 in the myocardium. Both WNT4 and BMP2 have been previously implicated in the development of VSDs in individuals with 46,XX sex reversal with dysgenesis of kidney, adrenals and lungs (SERKAL) syndrome and in individuals with short stature, facial dysmorphism and skeletal anomalies with or without cardiac anomalies 1 (SSFSC1) syndrome, respectively. We now show that Sox7 and Wnt4 interact genetically in the development of VSDs through their additive effects on endocardial cushion development with Sox7+/-;Wnt4+/- double heterozygous embryos having hypocellular endocardial cushions and perimembranous and muscular VSDs not seen in their Sox7+/- and Wnt4+/- littermates. These results provide additional evidence that SOX7, WNT4 and BMP2 function in the same pathway during mammalian septal development and that their deficiency can contribute to the development of VSDs in humans.

Canine RPGRIP1-cone-rod dystrophy (CRD), a model for human inherited retinal diseases (IRDs), was originally identified as autosomal recessive early-onset blindness. However, later studies revealed extensive phenotypic variability among RPGRIP1 mutants. This led to the identification of a homozygous MAP9 variant as a modifier associated with early-onset disease. Based on further phenotypic variation affecting cone photoreceptor function, we report mapping of L3 as an additional modifier locus, within a 4.1-Mb locus on canine chromosome 30. We establish the natural disease history of RPGRIP1-CRD based on up to 9-year long-term functional and structural retinal data from 58 dogs including 44 RPGRIP1 mutants grouped according to the modifier status. RPGRIP1 mutants affected by both MAP9 and L3 modifiers exhibited the most severe phenotypes with rapid disease progression. MAP9 alone was found to act as an overall accelerator of rod and cone diseases, while L3 had a cone-specific effect. Ultrastructural analysis of photoreceptors revealed varying degrees of rod and cone damage, while the connecting cilia appeared structurally preserved in all groups. We conclude that RPGRIP1-CRD is an oligogenic disease with at least three loci contributing to the pathogenesis. While the RPGRIP1 variant is required for developing the disease, MAP9 and L3 modifiers exacerbate the phenotype, individually and cumulatively. Oligogenic canine RPGRIP1-CRD illustrates the impact of multiple genetic modifiers on disease phenotype and thus has the potential to reveal new targets for broad-spectrum therapies for oligogenic or polygenic forms of human IRDs.