Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization.

BACKGROUND: As genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype. RESULTS: In this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation . CONCLUSIONS: GaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

Physician and informal care use explained by the Pediatric Quality of Life Inventory (PedsQL) in children with suspected genetic disorders.

PURPOSE: To examine associations between Pediatric Quality of Life Inventory (PedsQL) 4.0 Generic Core Scales and PedsQL Infant Scales with formal health care resource utilization (HCRU) and informal caregiver burden. METHODS: We studied a pediatric cohort of 837 patients (median age: 8.4 years) with suspected genetic disorders enrolled January 2019 through July 2021 in the NYCKidSeq program for diagnostic sequencing. Using linked ~ nine-month longitudinal survey and physician claims data collected through May 2022, we modeled the association between baseline PedsQL scores and post-baseline HCRU (median follow-up: 21.1 months) and informal care. We also assessed the longitudinal change in PedsQL scores with physician services using linear mixed-effects models. RESULTS: Lower PedsQL total and physical health scores were independently associated with increases in 18-month physician services, encounters, and weekly informal care. Comparing low vs. median total scores, increases were 10.6 services (95% CI: 1.0-24.6), 3.3 encounters (95% CI: 0.5-6.8), and $668 (95% CI: $350-965), respectively. For the psychosocial domain, higher scores were associated with decreased informal care. Based on adjusted linear mixed-effects modeling, every additional ten physician services was associated with diminished improvement in longitudinal PedsQL total score trajectories by 1.1 point (95% confidence interval: 0.6-1.6) on average. Similar trends were observed in the physical and psychosocial domains. CONCLUSION: PedsQL scores were independently associated with higher utilization of physician services and informal care. Moreover, longitudinal trajectories of PedsQL scores became less favorable with increased physician services. Adding PedsQL survey instruments to conventional measures for improved risk stratification should be evaluated in further research.

The Pediatric Quality of Life Inventory (PedsQL) is widely used to measure health-related quality of life in pediatric patients; however, few studies have examined whether the PedsQL is indicative of longitudinal outcomes of morbidity and health care needs. This study captures associations between PedsQL scores with utilization of physician and informal care in children with suspected genetic disorders. We demonstrate that lower PedsQL total and physical health scores are independently associated with greater utilization of physician services and informal care. Moreover, longitudinal trajectories of PedsQL scores become less favorable with increased physician services. Results can inform future applications of PedsQL instruments.

Physician services and costs after disclosure of diagnostic sequencing results in the NYCKidSeq program.

PURPOSE: To better understand the effects of returning diagnostic sequencing results on clinical actions and economic outcomes for pediatric patients with suspected genetic disorders. METHODS: Longitudinal physician claims data after diagnostic sequencing were obtained for patients aged 0 to 21 years with neurologic, cardiac, and immunologic disorders with suspected genetic etiology. We assessed specialist consultation rates prompted by primary diagnostic results, as well as marginal effects on overall 18-month physician services and costs. RESULTS: We included data on 857 patients (median age: 9.6 years) with a median follow-up of 17.3 months after disclosure of diagnostic sequencing results. The likelihood of having ≥1 recommendation for specialist consultation in 155 patients with positive findings was high (72%) vs 23% in 443 patients with uncertain findings and 21% in 259 patients with negative findings (P < .001). Follow-through consultation occurred in 30%. Increases in 18-month physician services and costs following a positive finding diminished after multivariable adjustment. Also, no significant differences between those with uncertain and negative findings were demonstrated. CONCLUSION: Our study did not provide evidence for significant increases in downstream physician services and costs after returning positive or uncertain diagnostic sequencing findings. More large-scale longitudinal studies are needed to confirm these findings.

Acute Lymphoblastic Leukemia with Myeloid Mutations Is a High-Risk Disease Associated with Clonal Hematopoiesis.

Myeloid neoplasms arise from preexisting clonal hematopoiesis (CH); however, the role of CH in the pathogenesis of acute lymphoblastic leukemia (ALL) is unknown. We found that 18% of adult ALL cases harbored TP53, and 16% had myeloid CH-associated gene mutations. ALL with myeloid mutations (MyM) had distinct genetic and clinical characteristics, associated with inferior survival. By using single-cell proteogenomic analysis, we demonstrated that myeloid mutations were present years before the diagnosis of ALL, and a subset of these clones expanded over time to manifest as dominant clones in ALL. Single-cell RNA sequencing revealed upregulation of genes associated with cell survival and resistance to apoptosis in B-ALL with MyM, which responds better to newer immunotherapeutic approaches. These findings define ALL with MyM as a high-risk disease that can arise from antecedent CH and offer new mechanistic insights to develop better therapeutic and preventative strategies. SIGNIFICANCE: CH is a precursor lesion for lymphoblastic leukemogenesis. ALL with MyM has distinct genetic and clinical characteristics, associated with adverse survival outcomes after chemotherapy. CH can precede ALL years before diagnosis, and ALL with MyM is enriched with activated T cells that respond to immunotherapies such as blinatumomab. See related commentary by Iacobucci, p. 142.

The NYCKidSeq randomized controlled trial: Impact of GUÍA digitally enhanced genetic results disclosure in diverse families.

Digital solutions are needed to support rapid increases in the application of genetic/genomic tests (GTs) in diverse clinical settings and patient populations. We developed GUÍA, a bilingual digital application that facilitates disclosure of GT results. The NYCKidSeq randomized controlled trial enrolled diverse children with neurologic, cardiac, and immunologic conditions who underwent GTs. The trial evaluated GUÍA's impact on understanding the GT results by randomizing families to results disclosure genetic counseling with GUÍA (intervention) or standard of care (SOC). Parents/legal guardians (participants) completed surveys at baseline, post-results disclosure, and 6 months later. Survey measures assessed the primary study outcomes of participants' perceived understanding of and confidence in explaining their child's GT results and the secondary outcome of objective understanding. The analysis included 551 diverse participants, 270 in the GUÍA arm and 281 in SOC. Participants in the GUÍA arm had significantly higher perceived understanding post-results (OR = 2.8, CI[1.004, 7.617], p = 0.049) and maintained higher objective understanding over time (OR = 1.1, CI[1.004, 1.127], p = 0.038) compared to SOC. There was no impact on perceived confidence. Hispanic/Latino(a) individuals in the GUÍA arm maintained higher perceived understanding (OR = 3.9, CI[1.603, 9.254], p = 0.003), confidence (OR = 2.7, CI[1.021, 7.277], p = 0.046), and objective understanding (OR = 1.1, CI[1.009, 1.212], p = 0.032) compared to SOC. This trial demonstrates that GUÍA positively impacts understanding of GT results in diverse parents of children with suspected genetic conditions and builds a case for utilizing GUÍA to deliver complex results. Continued development and evaluation of digital applications in diverse populations are critical for equitably scaling GT offerings in specialty clinics.

Early-life stress and ovarian hormones alter transcriptional regulation in the nucleus accumbens resulting in sex-specific responses to cocaine.

Early-life stress and ovarian hormones contribute to increased female vulnerability to cocaine addiction. Here, we reveal molecular substrates in the reward area, the nucleus accumbens, through which these female-specific factors affect immediate and conditioning responses to cocaine. We find shared involvement of X chromosome inactivation-related and estrogen signaling-related gene regulation in enhanced conditioning responses following early-life stress and during the low-estrogenic state in females. Low-estrogenic females respond to acute cocaine by opening neuronal chromatin enriched for the sites of ΔFosB, a transcription factor implicated in chronic cocaine response and addiction. Conversely, high-estrogenic females respond to cocaine by preferential chromatin closing, providing a mechanism for limiting cocaine-driven chromatin and synaptic plasticity. We find that physiological estrogen withdrawal, early-life stress, and absence of one X chromosome all nullify the protective effect of a high-estrogenic state on cocaine conditioning in females. Our findings offer a molecular framework to enable understanding of sex-specific neuronal mechanisms underlying cocaine use disorder.

Regulatory Landscape Enrichment Analysis (RLEA) using gaiaAssociation.

Motivation: To understand whether sets of genomic loci are enriched at the regulatory loci of one or more cell types, we developed the gaiaAssociation package to perform Regulatory Landscape Enrichment Analysis (RLEA). RLEA is a novel analytical process that tests for enrichment of sets of loci in cell type-specific open chromatin regions (OCRs) in the genome. Results: We demonstrate that the application of RLEA to genome-wide association study (GWAS) data reveals cell types likely to be mediating the phenotype studied, and clusters OCRs based on their shared regulatory profiles. GaiaAssociation is Python code that is freely available for use in functional genomics studies. Availability and Implementation: Gaia Association is available on PyPi (https://pypi.org/project/gaiaAssociation/0.6.0/#description) for pip download and use on the command line or as an inline Python package. Gaia Association can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation.

Urine Proteomics Link Complement Activation with Interstitial Fibrosis/Tubular Atrophy in Lupus Nephritis Patients.

BACKGROUND: Intrarenal complement activation has been implicated in the pathogenesis of tubulointerstitial fibrosis in lupus nephritis (LN) based on prior animal studies. The assembly of the membrane attack complex (MAC) by complement C5b to C9 on the cell membrane leads to cytotoxic pores and cell lysis, while CD59 inhibits MAC formation by preventing C9 from joining the complex. We hypothesize that complement activation and imbalance between complement activation and inhibition, as defined by increased production of individual complement components and uncontrolled MAC activation relative to CD59 inhibition, are associated with interstitial fibrosis and tubular atrophy (IFTA) in LN and correlate with the key mediators of kidney fibrosis- transforming growth factor receptors beta (TGFRß), platelet-derived growth factor beta (PDGFß) and platelet-derived growth factor receptor beta (PDGFRß). METHODS: We included urine samples from 46 adults and pediatric biopsy-proven lupus nephritis patients who underwent clinically indicated kidney biopsies between 2010 and 2019. We compared individual urinary complement components and the urinary C9-to-CD59 ratio between LN patients with moderate/severe IFTA and none/mild IFTA. IFTA was defined as none/mild (<25% of interstitium affected) versus moderate/severe (≥ 25% of interstitium affected). Proteomics analysis was performed using mass spectrometry (Orbitrap Fusion Lumos, Thermo Scientific) and processed by the Proteome Discoverer. Urinary complement proteins enriched in LN patients with moderate/severe IFTA were correlated with serum creatinine, TGFßR1, TGFßR2, PDGFß, and PDGFRß. RESULTS: Of the 46 LN patients included in the study, 41 (89.1%) were women, 20 (43.5%) self-identified as Hispanic or Latino, and 26 (56.5%) self-identified as Black or African American. Ten of the 46 (21.7%) LN patients had moderate/severe IFTA on kidney biopsy. LN patients with moderate/severe IFTA had an increased urinary C9-to-CD59 ratio [median 0.91 (0.83-1.05) vs 0.81 (0.76-0.91), p=0.01]. Urinary C3 and CFI levels in LN patients with moderate/severe IFTA were higher compared to those with none/mild IFTA [C3 median (IQR) 24.4(23.5-25.5) vs. 20.2 (18.5-22.2), p= 0.02], [CFI medium (IQR) 28.8 (21.8-30.6) vs. 20.4 (18.5-22.9), p=0.01]. Complement C9, CD59, C3 and CFI correlated with TGFßR1, PDGFß, and PDGFRß, while C9, CD59 and C3 correlated with TGFßR2. CONCLUSION: This study is one of the first to compare the urinary complement profile in LN patients with moderate/severe IFTA and none/mild IFTA in human tissues. This study identified C3, CFI, and C9-to-CD59 ratio as potential markers of tubulointerstitial fibrosis in LN.

Prenatal vitamin D deficiency alters immune cell proportions of young adult offspring through alteration of long-term stem cell fates.

Vitamin D deficiency is a common deficiency worldwide, particularly among women of reproductive age. During pregnancy, it increases the risk of immune-related diseases in offspring later in life. However, exactly how the body remembers exposure to an adverse environment during development is poorly understood. Herein, we explore the effects of prenatal vitamin D deficiency on immune cell proportions in offspring using vitamin D deficient mice established by dietary manipulation. We show that prenatal vitamin D deficiency alters immune cell proportions in offspring by changing the transcriptional properties of genes downstream of vitamin D receptor signaling in hematopoietic stem and progenitor cells of both the fetus and adults. Further investigations of the associations between maternal vitamin D levels and cord blood immune cell profiles from 75 healthy pregnant women and their term babies also confirm that maternal vitamin D levels significantly affect immune cell proportions in the babies. Thus, lack of prenatal vitamin D, particularly at the time of hematopoietic stem cell migration from the liver to the bone marrow, has long-lasting effects on immune cell proportions. This highlights the importance of providing vitamin D supplementation at specific stages of pregnancy.

The NYCKidSeq randomized controlled trial: Impact of GUÍA digitally enhanced genetic counseling in racially and ethnically diverse families.

Background: Digital solutions are needed to support rapid increases in the application of genetic and genomic tests (GT) in diverse clinical settings and patient populations. We developed GUÍA, a bi-lingual web-based platform that facilitates disclosure of GT results. The NYCKidSeq randomized controlled trial evaluated GUÍA's impact on understanding of GT results. Methods: NYCKidSeq enrolled diverse children with neurologic, cardiac, and immunologic conditions who underwent GT. Families were randomized to genetic counseling with GUÍA (intervention) or standard of care (SOC) genetic counseling for results disclosure. Parents/legal guardians (participants) completed surveys at baseline, post-results disclosure, and 6-months later. Survey measures assessed the primary study outcomes of perceived understanding of and confidence in explaining their child's GT results and the secondary outcome of objective understanding. We used regression models to evaluate the association between the intervention and the study outcomes. Results: The analysis included 551 participants, 270 in the GUÍA arm and 281 in SOC. Participants' mean age was 41.1 years and 88.6% were mothers. Most participants were Hispanic/Latino(a) (46.3%), White/European American (24.5%), or Black/African American (15.8%). Participants in the GUÍA arm had significantly higher perceived understanding post-results (OR=2.8, CI[1.004,7.617], P=0.049) and maintained higher objective understanding over time (OR=1.1, CI[1.004, 1.127], P=0.038) compared to those in the SOC arm. There was no impact on perceived confidence. Hispanic/Latino(a) individuals in the GUÍA arm maintained higher perceived understanding (OR=3.9, CI[1.6, 9.3], P=0.003), confidence (OR=2.7, CI[1.021, 7.277], P=0.046), and objective understanding (OR=1.1, CI[1.009, 1.212], P=0.032) compared to SOC . Conclusions: This trial demonstrates that GUÍA positively impacts understanding of GT results in diverse parents of children with suspected genetic conditions. These findings build a case for utilizing GUÍA to deliver complex and often ambiguous genetic results. Continued development and evaluation of digital applications in diverse populations are critical for equitably scaling GT offerings in specialty clinics. Trial Registration: Clinicaltrials.gov identifier NCT03738098.

Identification of copy number variants with genome sequencing: Clinical experiences from the NYCKidSeq program.

Copy number variations (CNVs) play a significant role in human disease. While chromosomal microarray has traditionally been the first-tier test for CNV detection, use of genome sequencing (GS) is increasing. We report the frequency of CNVs detected with GS in a diverse pediatric cohort from the NYCKidSeq program and highlight specific examples of its clinical impact. A total of 1052 children (0-21 years) with neurodevelopmental, cardiac, and/or immunodeficiency phenotypes received GS. Phenotype-driven analysis was used, resulting in 183 (17.4%) participants with a diagnostic result. CNVs accounted for 20.2% of participants with a diagnostic result (37/183) and ranged from 0.5 kb to 16 Mb. Of participants with a diagnostic result (n = 183) and phenotypes in more than one category, 5/17 (29.4%) were solved by a CNV finding, suggesting a high prevalence of diagnostic CNVs in participants with complex phenotypes. Thirteen participants with a diagnostic CNV (35.1%) had previously uninformative genetic testing, of which nine included a chromosomal microarray. This study demonstrates the benefits of GS for reliable detection of CNVs in a pediatric cohort with variable phenotypes.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

PURPOSE: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. METHODS: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. RESULTS: A total of 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses (P < .001). Yield was greater for GS vs TGPs in Hispanic/Latino(a) (17.2% vs 9.5%, P < .001) and White/European American (19.8% vs 7.9%, P < .001) but not in Black/African American (11.5% vs 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. CONCLUSION: GS may yield up to twice as many diagnoses in pediatric patients compared with TGP testing but not yet across all population groups.

An optimized approach for multiplexing single-nuclear ATAC-seq using oligonucleotide-conjugated antibodies.

BACKGROUND: Single-cell technologies to analyze transcription and chromatin structure have been widely used in many research areas to reveal the functions and molecular properties of cells at single-cell resolution. Sample multiplexing techniques are valuable when performing single-cell analysis, reducing technical variation and permitting cost efficiencies. Several commercially available methods have been used in many scRNA-seq studies. On the other hand, while several methods have been published, multiplexing techniques for single nuclear assay for transposase-accessible chromatin (snATAC)-seq assays remain under development. We developed a simple nucleus hashing method using oligonucleotide-conjugated antibodies recognizing nuclear pore complex proteins, NuHash, to perform snATAC-seq library preparations by multiplexing. RESULTS: We performed multiplexing snATAC-seq analyses on a mixture of human and mouse cell samples (two samples, 2-plex, and four samples, 4-plex) using NuHash. The analyses on nuclei with at least 10,000 read counts showed that the demultiplexing accuracy of NuHash was high, and only ten out of 9144 nuclei (2-plex) and 150 of 12,208 nuclei (4-plex) had discordant classifications between NuHash demultiplexing and discrimination using reference genome alignments. The differential open chromatin region (OCR) analysis between female and male samples revealed that male-specific OCRs were enriched in chromosome Y (four out of nine). We also found that five female-specific OCRs (20 OCRs) were on chromosome X. A comparative analysis between snATAC-seq and deeply sequenced bulk ATAC-seq on the same samples revealed that the bulk ATAC-seq signal intensity was positively correlated with the number of cell clusters detected in snATAC-seq. Moreover, when we categorized snATAC-seq peaks based on the number of cell clusters in which the peak was present, we observed different distributions over different genomic features between the groups. This result suggests that the peak intensities of bulk ATAC-seq can be used to identify different types of functional loci. CONCLUSIONS: Our multiplexing method using oligo-conjugated anti-nuclear pore complex proteins, NuHash, permits high-accuracy demultiplexing of samples. The NuHash protocol is straightforward, works on frozen samples, and requires no modifications for snATAC-seq library preparation.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

Purpose: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. Methods: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. Results: 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses ( P < .001). Yield was greater for GS vs . TGPs in Hispanic/Latino(a) (17.2% vs . 9.5%, P < .001) and White/European American (19.8% vs . 7.9%, P < .001), but not in Black/African American (11.5% vs . 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs . White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. Conclusion: GS may yield up to twice as many diagnoses in pediatric patients compared to TGP testing, but not yet across all population groups.

The TeleKidSeq pilot study: incorporating telehealth into clinical care of children from diverse backgrounds undergoing whole genome sequencing.

BACKGROUND: The COVID-19 pandemic forced healthcare institutions and many clinical research programs to adopt telehealth modalities in order to mitigate viral spread. With the expanded use of telehealth, there is the potential to increase access to genomic medicine to medically underserved populations, yet little is known about how best to communicate genomic results via telehealth while also ensuring equitable access. NYCKidSeq, a multi-institutional clinical genomics research program in New York City, launched the TeleKidSeq pilot study to assess alternative forms of genomic communication and telehealth service delivery models with families from medically underserved populations. METHODS: We aim to enroll 496 participants between 0 and 21 years old to receive clinical genome sequencing. These individuals have a neurologic, cardiovascular, and/or immunologic disease. Participants will be English- or Spanish-speaking and predominantly from underrepresented groups who receive care in the New York metropolitan area. Prior to enrollment, participants will be randomized to either genetic counseling via videoconferencing with screen-sharing or genetic counseling via videoconferencing without screen-sharing. Using surveys administered at baseline, results disclosure, and 6-months post-results disclosure, we will evaluate the impact of the use of screen-sharing on participant understanding, satisfaction, and uptake of medical recommendations, as well as the psychological and socioeconomic implications of obtaining genome sequencing. Clinical utility, cost, and diagnostic yield of genome sequencing will also be assessed. DISCUSSION: The TeleKidSeq pilot study will contribute to innovations in communicating genomic test results to diverse populations through telehealth technology. In conjunction with NYCKidSeq, this work will inform best practices for the implementation of genomic medicine in diverse, English- and Spanish-speaking populations.

Using epigenomics to understand cellular responses to environmental influences in diseases.

It is a generally accepted model that environmental influences can exert their effects, at least in part, by changing the molecular regulators of transcription that are described as epigenetic. As there is biochemical evidence that some epigenetic regulators of transcription can maintain their states long term and through cell division, an epigenetic model encompasses the idea of maintenance of the effect of an exposure long after it is no longer present. The evidence supporting this model is mostly from the observation of alterations of molecular regulators of transcription following exposures. With the understanding that the interpretation of these associations is more complex than originally recognised, this model may be oversimplistic; therefore, adopting novel perspectives and experimental approaches when examining how environmental exposures are linked to phenotypes may prove worthwhile. In this review, we have chosen to use the example of nonalcoholic fatty liver disease (NAFLD), a common, complex human disease with strong environmental and genetic influences. We describe how epigenomic approaches combined with emerging functional genetic and single-cell genomic techniques are poised to generate new insights into the pathogenesis of environmentally influenced human disease phenotypes exemplified by NAFLD.

Detection of mosaic variants using genome sequencing in a large pediatric cohort.

The increased use of next-generation sequencing has expanded our understanding of the involvement and prevalence of mosaicism in genetic disorders. We describe a total of eleven cases: nine in which mosaic variants detected by genome sequencing (GS) and/or targeted gene panels (TGPs) were considered to be causative for the proband's phenotype, and two of apparent parental mosaicism. Variants were identified in the following genes: PHACTR1, SCN8A, KCNT1, CDKL5, NEXMIF, CUX1, TSC2, GABRB2, and SMARCB1. In addition, we identified one large duplication including three genes, UBE3A, GABRB3, and MAGEL2, and one large deletion including deletion of ARFGAP1, EEF1A2, CHRNA4, and KCNQ2. All patients were enrolled in the NYCKidSeq study, a research program studying the communication of genomic information in clinical care, as well as the clinical utility and diagnostic yield of GS for children with suspected genetic disorders in diverse populations in New York City. We observed variability in the correlation between reported variant allele fraction and the severity of the patient's phenotype, although we were not able to determine the mosaicism percentage in clinically relevant tissue(s). Although our study was not sufficiently powered to assess differences in mosaicism detection between the two testing modalities, we saw a trend toward better detection by GS as compared with TGP testing. This case series supports the importance of mosaicism in childhood-onset genetic conditions and informs guidelines for laboratory and clinical interpretation of mosaic variants detected by GS.

Expansion and mechanistic insights into de novo DEAF1 variants in DEAF1-associated neurodevelopmental disorders.

De novo deleterious and heritable biallelic mutations in the DNA binding domain (DBD) of the transcription factor deformed epidermal autoregulatory factor 1 (DEAF1) result in a phenotypic spectrum of disorders termed DEAF1-associated neurodevelopmental disorders (DAND). RNA-sequencing using hippocampal RNA from mice with conditional deletion of Deaf1 in the central nervous system indicate that loss of Deaf1 activity results in the altered expression of genes involved in neuronal function, dendritic spine maintenance, development, and activity, with reduced dendritic spines in hippocampal regions. Since DEAF1 is not a dosage-sensitive gene, we assessed the dominant negative activity of previously identified de novo variants and a heritable recessive DEAF1 variant on selected DEAF1-regulated genes in 2 different cell models. While no altered gene expression was observed in cells over-expressing the recessive heritable variant, the gene expression profiles of cells over-expressing de novo variants resulted in similar gene expression changes as observed in CRISPR-Cas9-mediated DEAF1-deleted cells. Altered expression of DEAF1-regulated genes was rescued by exogenous expression of WT-DEAF1 but not by de novo variants in cells lacking endogenous DEAF1. De novo heterozygous variants within the DBD of DEAF1 were identified in 10 individuals with a phenotypic spectrum including autism spectrum disorder, developmental delays, sleep disturbance, high pain tolerance, and mild dysmorphic features. Functional assays demonstrate these variants alter DEAF1 transcriptional activity. Taken together, this study expands the clinical phenotypic spectrum of individuals with DAND, furthers our understanding of potential roles of DEAF1 on neuronal function, and demonstrates dominant negative activity of identified de novo variants.

Egr1 is a sex-specific regulator of neuronal chromatin, synaptic plasticity, and behaviour.

Sex differences are found in brain structure and function across species, and across brain disorders in humans1-3. The major source of brain sex differences is differential secretion of steroid hormones from the gonads across the lifespan4. Specifically, ovarian hormones oestrogens and progesterone are known to dynamically change structure and function of the adult female brain, having a major impact on psychiatric risk5-7. However, due to limited molecular studies in female rodents8, very little is still known about molecular drivers of female-specific brain and behavioural plasticity. Here we show that overexpressing Egr1, a candidate oestrous cycle-dependent transcription factor9, induces sex-specific changes in ventral hippocampal neuronal chromatin, gene expression, and synaptic plasticity, along with hippocampus-dependent behaviours. Importantly, Egr1 overexpression mimics the high-oestrogenic phase of the oestrous cycle, and affects behaviours in ovarian hormone-depleted females but not in males. We demonstrate that Egr1 opens neuronal chromatin directly across the sexes, although with limited genomic overlap. Our study not only reveals the first sex-specific chromatin regulator in the brain, but also provides functional evidence that this sex-specific gene regulation drives neuronal gene expression, synaptic plasticity, and anxiety- and depression-related behaviour. Our study exemplifies an innovative sex-based approach to studying neuronal gene regulation1 in order to understand sex-specific synaptic and behavioural plasticity and inform novel brain disease treatments.