A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

The diabetes susceptibility gene Clec16a regulates mitophagy.

Clec16a has been identified as a disease susceptibility gene for type 1 diabetes, multiple sclerosis, and adrenal dysfunction, but its function is unknown. Here we report that Clec16a is a membrane-associated endosomal protein that interacts with E3 ubiquitin ligase Nrdp1. Loss of Clec16a leads to an increase in the Nrdp1 target Parkin, a master regulator of mitophagy. Islets from mice with pancreas-specific deletion of Clec16a have abnormal mitochondria with reduced oxygen consumption and ATP concentration, both of which are required for normal ß cell function. Indeed, pancreatic Clec16a is required for normal glucose-stimulated insulin release. Moreover, patients harboring a diabetogenic SNP in the Clec16a gene have reduced islet Clec16a expression and reduced insulin secretion. Thus, Clec16a controls ß cell function and prevents diabetes by controlling mitophagy. This pathway could be targeted for prevention and control of diabetes and may extend to the pathogenesis of other Clec16a- and Parkin-associated diseases.

Managing differential performance of polygenic risk scores across groups: Real-world experience of the eMERGE Network.

The differential performance of polygenic risk scores (PRSs) by group is one of the major ethical barriers to their clinical use. It is also one of the main practical challenges for any implementation effort. The social repercussions of how people are grouped in PRS research must be considered in communications with research participants, including return of results. Here, we outline the decisions faced and choices made by a large multi-site clinical implementation study returning PRSs to diverse participants in handling this issue of differential performance. Our approach to managing the complexities associated with the differential performance of PRSs serves as a case study that can help future implementers of PRSs to plot an anticipatory course in response to this issue.

Dependency-aware deep generative models for multitasking analysis of spatial omics data.

Spatially resolved transcriptomics (SRT) technologies have significantly advanced biomedical research, but their data analysis remains challenging due to the discrete nature of the data and the high levels of noise, compounded by complex spatial dependencies. Here, we propose spaVAE, a dependency-aware, deep generative spatial variational autoencoder model that probabilistically characterizes count data while capturing spatial correlations. spaVAE introduces a hybrid embedding combining a Gaussian process prior with a Gaussian prior to explicitly capture spatial correlations among spots. It then optimizes the parameters of deep neural networks to approximate the distributions underlying the SRT data. With the approximated distributions, spaVAE can contribute to several analytical tasks that are essential for SRT data analysis, including dimensionality reduction, visualization, clustering, batch integration, denoising, differential expression, spatial interpolation, resolution enhancement and identification of spatially variable genes. Moreover, we have extended spaVAE to spaPeakVAE and spaMultiVAE to characterize spatial ATAC-seq (assay for transposase-accessible chromatin using sequencing) data and spatial multi-omics data, respectively.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets.

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Prospective, multi-site study of healthcare utilization after actionable monogenic findings from clinical sequencing.

As large-scale genomic screening becomes increasingly prevalent, understanding the influence of actionable results on healthcare utilization is key to estimating the potential long-term clinical impact. The eMERGE network sequenced individuals for actionable genes in multiple genetic conditions and returned results to individuals, providers, and the electronic health record. Differences in recommended health services (laboratory, imaging, and procedural testing) delivered within 12 months of return were compared among individuals with pathogenic or likely pathogenic (P/LP) findings to matched individuals with negative findings before and after return of results. Of 16,218 adults, 477 unselected individuals were found to have a monogenic risk for arrhythmia (n = 95), breast cancer (n = 96), cardiomyopathy (n = 95), colorectal cancer (n = 105), or familial hypercholesterolemia (n = 86). Individuals with P/LP results more frequently received services after return (43.8%) compared to before return (25.6%) of results and compared to individuals with negative findings (24.9%; p < 0.0001). The annual cost of qualifying healthcare services increased from an average of $162 before return to $343 after return of results among the P/LP group (p < 0.0001); differences in the negative group were non-significant. The mean difference-in-differences was $149 (p < 0.0001), which describes the increased cost within the P/LP group corrected for cost changes in the negative group. When stratified by individual conditions, significant cost differences were observed for arrhythmia, breast cancer, and cardiomyopathy. In conclusion, less than half of individuals received billed health services after monogenic return, which modestly increased healthcare costs for payors in the year following return.

Complex hierarchical structures in single-cell genomics data unveiled by deep hyperbolic manifold learning.

With the advances in single-cell sequencing techniques, numerous analytical methods have been developed for delineating cell development. However, most are based on Euclidean space, which would distort the complex hierarchical structure of cell differentiation. Recently, methods acting on hyperbolic space have been proposed to visualize hierarchical structures in single-cell RNA-seq (scRNA-seq) data and have been proven to be superior to methods acting on Euclidean space. However, these methods have fundamental limitations and are not optimized for the highly sparse single-cell count data. To address these limitations, we propose scDHMap, a model-based deep learning approach to visualize the complex hierarchical structures of scRNA-seq data in low-dimensional hyperbolic space. The evaluations on extensive simulation and real experiments show that scDHMap outperforms existing dimensionality-reduction methods in various common analytical tasks as needed for scRNA-seq data, including revealing trajectory branches, batch correction, and denoising the count matrix with high dropout rates. In addition, we extend scDHMap to visualize single-cell ATAC-seq data.

Single-Cell RNA Sequencing Technology Landscape in 2023.

Single-cell RNA sequencing (scRNA-seq) has revolutionized our understanding of cellular heterogeneity and the dynamics of gene expression, bearing profound significance in stem cell research. Depending on the starting materials used for analysis, scRNA-seq encompasses scRNA-seq and single-nucleus RNA sequencing (snRNA-seq). scRNA-seq excels in capturing cellular heterogeneity and characterizing rare cell populations within complex tissues, while snRNA-seq is advantageous in situations where intact cell dissociation is challenging or undesirable (eg, epigenomic studies). A number of scRNA-seq technologies have been developed as of late, including but not limited to droplet-based, plate-based, hydrogel-based, and spatial transcriptomics. The number of cells, sequencing depth, and sequencing length in scRNA-seq can vary across different studies. Addressing current technical challenges will drive the future of scRNA-seq, leading to more comprehensive and precise insights into cellular biology and disease mechanisms informing therapeutic interventions.

Elevated Levels of the Cytokine LIGHT in Pediatric Crohn's Disease.

LIGHT (homologous to lymphotoxins, exhibits inducible expression, and competes with HSV glycoprotein D for herpes virus entry mediator, a receptor expressed by T lymphocytes), encoded by the TNFSF14 gene, is a cytokine belonging to the TNF superfamily. On binding to its receptors, herpes virus entry mediator and lymphotoxin ß receptor, it activates inflammatory responses. We conducted this study to determine whether plasma LIGHT levels are elevated in Crohn's disease (CD) in a pediatric population with the aim of nominating this cytokine as a therapeutic target. We used a single-molecule immunoassay to determine the circulating levels of free LIGHT in plasma from pediatric patients with CD in our biobank (n = 183), a panel of healthy pediatric (n = 9) or adult (n = 22) reference samples, and pediatric biobank controls (n = 19). We performed correlational analyses between LIGHT levels and the clinical characteristics of the CD cohort, including age, Montreal classification, family history, medical/surgical therapy, and routine blood test parameters. LIGHT levels were greatly elevated in CD, with an average of 305 versus 32.4 pg/ml for controls from the biobank (p < 0.0001). The outside reference samples showed levels of 57 pg/ml in pediatric controls and 55 pg/ml in adults (p < 0.0001). We found a statistically significant correlation between white blood cell count and free LIGHT (p < 0.046). We conclude that free, soluble LIGHT is increased 5- to 10-fold in pediatric CD across an array of disease subtypes and characteristics.

NASP gene contributes to autism by epigenetic dysregulation of neural and immune pathways.

BACKGROUND: Epigenetics makes substantial contribution to the aetiology of autism spectrum disorder (ASD) and may harbour a unique opportunity to prevent the development of ASD. We aimed to identify novel epigenetic genes involved in ASD aetiology. METHODS: Trio-based whole exome sequencing was conducted on ASD families. Genome editing technique was used to knock out the candidate causal gene in a relevant cell line. ATAC-seq, ChIP-seq and RNA-seq were performed to investigate the functional impact of knockout (KO) or mutation in the candidate gene. RESULTS: We identified a novel candidate gene NASP (nuclear autoantigenic sperm protein) for epigenetic dysregulation in ASD in a Chinese nuclear family including one proband with autism and comorbid atopic disease. The de novo likely gene disruptive variant tNASP(Q289X) subjects the expression of tNASP to nonsense-mediated decay. tNASP KO increases chromatin accessibility, promotes the active promoter state of genes enriched in synaptic signalling and leads to upregulated expression of genes in the neural signalling and immune signalling pathways. Compared with wild-type tNASP, tNASP(Q289X) enhances chromatin accessibility of the genes with enriched expression in the brain. RNA-seq revealed that genes involved in neural and immune signalling are affected by the tNASP mutation, consistent with the phenotypic impact and molecular effects of nasp-1 mutations in Caenorhabditis elegans. Two additional patients with ASD were found carrying deletion or deleterious mutation in the NASP gene. CONCLUSION: We identified novel epigenetic mechanisms mediated by tNASP which may contribute to the pathogenesis of ASD and its immune comorbidity.

CSF-1 maintains pathogenic but not homeostatic myeloid cells in the central nervous system during autoimmune neuroinflammation.

The receptor for colony stimulating factor 1 (CSF-1R) is important for the survival and function of myeloid cells that mediate pathology during experimental autoimmune encephalomyelitis (EAE), an animal model of multiple sclerosis (MS). CSF-1 and IL-34, the ligands of CSF-1R, have similar bioactivities but distinct tissue and context-dependent expression patterns, suggesting that they have different roles. This could be the case in EAE, given that CSF-1 expression is up-regulated in the CNS, while IL-34 remains constitutively expressed. We found that targeting CSF-1 with neutralizing antibody halted ongoing EAE, with efficacy superior to CSF-1R inhibitor BLZ945, whereas IL-34 neutralization had no effect, suggesting that pathogenic myeloid cells were maintained by CSF-1. Both anti­CSF-1 and BLZ945 treatment greatly reduced the number of monocyte-derived cells and microglia in the CNS. However, anti­CSF-1 selectively depleted inflammatory microglia and monocytes in inflamed CNS areas, whereas BLZ945 depleted virtually all myeloid cells, including quiescent microglia, throughout the CNS. Anti­CSF-1 treatment reduced the size of demyelinated lesions and microglial activation in the gray matter. Lastly, we found that bone marrow­derived immune cells were the major mediators of CSF-1R­dependent pathology, while microglia played a lesser role. Our findings suggest that targeting CSF-1 could be effective in ameliorating MS pathology, while preserving the homeostatic functions of myeloid cells, thereby minimizing risks associated with ablation of CSF-1R­dependent cells.

Target genes regulated by CLEC16A intronic region associated with common variable immunodeficiency.

BACKGROUND: CLEC16A intron 19 has been identified as a candidate locus for common variable immunodeficiency (CVID). OBJECTIVES: This study sought to elucidate the molecular mechanism by which variants at the CLEC16A intronic locus may contribute to the pathogenesis of CVID. METHODS: The investigators performed fine-mapping of the CLEC16A locus in a CVID cohort, then deleted the candidate functional SNP in T-cell lines by the CRISPR-Cas9 technique and conducted RNA-sequencing to identify target gene(s). The interactions between the CLEC16A locus and its target genes were identified using circular chromosome conformation capture. The transcription factor complexes mediating the chromatin interactions were determined by proteomic approach. The molecular pathways regulated by the CLEC16A locus were examined by RNA-sequencing and reverse phase protein array. RESULTS: This study showed that the CLEC16A locus is an enhancer regulating expression of multiple target genes including a distant gene ATF7IP2 through chromatin interactions. Distinct transcription factor complexes mediate the chromatin interactions in an allele-specific manner. Disruption of the CLEC16A locus affects the AKT signaling pathway, as well as the molecular response of CD4+ T cells to immune stimulation. CONCLUSIONS: Through multiomics and targeted experimental approaches, this study elucidated the underlying target genes and signaling pathways involved in the genetic association of CLEC16A with CVID, and highlighted plausible molecular targets for developing novel therapeutics.

The common variable immunodeficiency IgM repertoire narrowly recognizes erythrocyte and platelet glycans.

BACKGROUND: Autoimmune cytopenias (AICs) regularly occur in profoundly IgG-deficient patients with common variable immunodeficiency (CVID). The isotypes, antigenic targets, and origin(s) of their disease-causing autoantibodies are unclear. OBJECTIVE: We sought to determine reactivity, clonality, and provenance of AIC-associated IgM autoantibodies in patients with CVID. METHODS: We used glycan arrays, patient erythrocytes, and platelets to determine targets of CVID IgM autoantibodies. Glycan-binding profiles were used to identify autoreactive clones across B-cell subsets, specifically circulating marginal zone (MZ) B cells, for sorting and IGH sequencing. The locations, transcriptomes, and responses of tonsillar MZ B cells to different TH- cell subsets were determined by confocal microscopy, RNA-sequencing, and cocultures, respectively. RESULTS: Autoreactive IgM coated erythrocytes and platelets from many CVID patients with AICs (CVID+AIC). On glycan arrays, CVID+AIC plasma IgM narrowly recognized erythrocytic i antigens and platelet i-related antigens and failed to bind hundreds of pathogen- and tumor-associated carbohydrates. Polyclonal i antigen-recognizing B-cell receptors were highly enriched among CVID+AIC circulating MZ B cells. Within tonsillar tissues, MZ B cells secreted copious IgM when activated by the combination of IL-10 and IL-21 or when cultured with IL-10/IL-21-secreting FOXP3-CD25hi T follicular helper (Tfh) cells. In lymph nodes from immunocompetent controls, MZ B cells, plentiful FOXP3+ regulatory T cells, and rare FOXP3-CD25+ cells that represented likely CD25hi Tfh cells all localized outside of germinal centers. In CVID+AIC lymph nodes, cellular positions were similar but CD25hi Tfh cells greatly outnumbered regulatory cells. CONCLUSIONS: Our findings indicate that glycan-reactive IgM autoantibodies produced outside of germinal centers may contribute to the autoimmune pathogenesis of CVID.

Mutation burden analysis of six common mental disorders in African Americans by whole genome sequencing.

Mental disorders present a global health concern and have limited treatment options. In today's medical practice, medications such as antidepressants are prescribed not only for depression but also for conditions such as anxiety and attention deficit hyperactivity disorder (ADHD). Therefore, identifying gene targets for specific disorders is important and offers improved precision. In this study, we performed a genetic analysis of six common mental disorders-ADHD, anxiety, depression, delays in mental development, intellectual disabilities (IDs) and speech/language disorder-in the ethnic minority of African Americans (AAs) using whole genome sequencing (WGS). WGS data were generated from blood-derived DNA from 4178 AA individuals, including 1384 patients with the diagnosis of at least one mental disorder. Mutation burden analysis was applied based on rare and deleterious mutations in the AA population between cases and controls, and further analyzed in the context of patients with single mental disorder diagnosis. Certain genes uncovered demonstrated significant P-values in mutation burden analysis. In addition, exclusive recurrences in specific type of disorder were scanned through gene-drug interaction databases to assess for availability of potential medications. We uncovered 15 genes harboring deleterious mutations, including 3-Hydroxy-3-Methylglutaryl-CoA Reductase (HMGCR) and Uronyl 2-Sulfotransferase (UST) for ADHD; Farnesyltransferase, CAAX Box, Beta (FNTB) for anxiety; Xin Actin Binding Repeat Containing 2 (XIRP2), Natriuretic Peptide C (NPPC), Serine/Threonine Kinase 33 (STK33), Pannexin 1 (PANX1) and Neurotensin (NTS) for depression; RUNX Family Transcription Factor 3 (RUNX3), Tachykinin Receptor 1 (TACR1) and NADH:Ubiquinone Oxidoreductase Core Subunit S7 (NDUFS7) for delays in mental development; Hepsin (HPN) for ID and Collagen Type VI Alpha 3 Chain (COL6A3), Damage Specific DNA Binding Protein 1 (DDB1) and NADH:Ubiquinone Oxidoreductase Subunit A11 (NDUFA11) for speech/language disorder. Taken together, we have established critical insights into the development of new precision medicine approaches for mental disorders in AAs.

Bi-allelic variants in SPATA5L1 lead to intellectual disability, spastic-dystonic cerebral palsy, epilepsy, and hearing loss.

Spermatogenesis-associated 5 like 1 (SPATA5L1) represents an orphan gene encoding a protein of unknown function. We report 28 bi-allelic variants in SPATA5L1 associated with sensorineural hearing loss in 47 individuals from 28 (26 unrelated) families. In addition, 25/47 affected individuals (53%) presented with microcephaly, developmental delay/intellectual disability, cerebral palsy, and/or epilepsy. Modeling indicated damaging effect of variants on the protein, largely via destabilizing effects on protein domains. Brain imaging revealed diminished cerebral volume, thin corpus callosum, and periventricular leukomalacia, and quantitative volumetry demonstrated significantly diminished white matter volumes in several individuals. Immunofluorescent imaging in rat hippocampal neurons revealed localization of Spata5l1 in neuronal and glial cell nuclei and more prominent expression in neurons. In the rodent inner ear, Spata5l1 is expressed in the neurosensory hair cells and inner ear supporting cells. Transcriptomic analysis performed with fibroblasts from affected individuals was able to distinguish affected from controls by principal components. Analysis of differentially expressed genes and networks suggested a role for SPATA5L1 in cell surface adhesion receptor function, intracellular focal adhesions, and DNA replication and mitosis. Collectively, our results indicate that bi-allelic SPATA5L1 variants lead to a human disease characterized by sensorineural hearing loss (SNHL) with or without a nonprogressive mixed neurodevelopmental phenotype.

Monogenic etiologies of persistent human papillomavirus infections: A comprehensive systematic review.

PURPOSE: Persistent human papillomavirus infection (PHPVI) causes cutaneous, anogenital, and mucosal warts. Cutaneous warts include common warts, Treeman syndrome, and epidermodysplasia verruciformis, among others. Although more reports of monogenic predisposition to PHPVI have been published with the development of genomic technologies, genetic testing is rarely incorporated into clinical assessments. To encourage broader molecular testing, we compiled a list of the various monogenic etiologies of PHPVI. METHODS: We conducted a systematic literature review to determine the genetic, immunological, and clinical characteristics of patients with PHPVI. RESULTS: The inclusion criteria were met by 261 of 40,687 articles. In 842 patients, 83 PHPVI-associated genes were identified, including 42, 6, and 35 genes with strong, moderate, and weak evidence for causality, respectively. Autosomal recessive inheritance predominated (69%). PHPVI onset age was 10.8 ± 8.6 years, with an interquartile range of 5 to 14 years. GATA2,IL2RG,DOCK8, CXCR4, TMC6, TMC8, and CIB1 are the most frequently reported PHPVI-associated genes with strong causality. Most genes (74 out of 83) belong to a catalog of 485 inborn errors of immunity-related genes, and 40 genes (54%) are represented in the nonsyndromic and syndromic combined immunodeficiency categories. CONCLUSION: PHPVI has at least 83 monogenic etiologies and a genetic diagnosis is essential for effective management.

De novo variants in RNF213 are associated with a clinical spectrum ranging from Leigh syndrome to early-onset stroke.

PURPOSE: RNF213, encoding a giant E3 ubiquitin ligase, has been recognized for its role as a key susceptibility gene for moyamoya disease. Case reports have also implicated specific variants in RNF213 with an early-onset form of moyamoya disease with full penetrance. We aimed to expand the phenotypic spectrum of monogenic RNF213-related disease and to evaluate genotype-phenotype correlations. METHODS: Patients were identified through reanalysis of exome sequencing data of an unselected cohort of unsolved pediatric cases and through GeneMatcher or ClinVar. Functional characterization was done by proteomics analysis and oxidative phosphorylation enzyme activities using patient-derived fibroblasts. RESULTS: We identified 14 individuals from 13 unrelated families with (de novo) missense variants in RNF213 clustering within or around the Really Interesting New Gene (RING) domain. Individuals presented either with early-onset stroke (n = 11) or with Leigh syndrome (n = 3). No genotype-phenotype correlation could be established. Proteomics using patient-derived fibroblasts revealed no significant differences between clinical subgroups. 3D modeling revealed a clustering of missense variants in the tertiary structure of RNF213 potentially affecting zinc-binding suggesting a gain-of-function or dominant negative effect. CONCLUSION: De novo missense variants in RNF213 clustering in the E3 RING or other regions affecting zinc-binding lead to an early-onset syndrome characterized by stroke or Leigh syndrome.

Revealing Novel Genomic Insights and Therapeutic Targets for Juvenile Idiopathic Arthritis through Omics.

BACKGROUND: The genetic architecture of juvenile idiopathic arthritis (JIA) remains only partially comprehended. There is a clear imperative for continued endeavors to uncover insights into the underlying causes of JIA. METHODS: This study encompassed a comprehensive spectrum of endeavors, including conducting a JIA GWAS meta-analysis that incorporated data from 4,550 JIA cases and 18 446 controls. We employed in silico and genome-editing approaches to prioritize target genes. To investigate pleiotropic effects, we conducted phenome-wide association studies. Cell-type enrichment analyses were performed by integrating bulk and single-cell sequencing data. Finally, we delved into potential druggable targets for JIA. RESULTS: Fourteen genome-wide significant non-HLA loci were identified including four novel loci, each exhibiting pleiotropic associations with other autoimmune diseases or musculoskeletal traits. We uncovered strong genetic correlation between JIA and bone mineral density (BMD) traits at 52 genomic regions, including three GWAS loci for JIA. Candidate genes with immune functions were captured by in silico analyses at each novel locus, with additional findings identified through our experimental approach. Cell-type enrichment analysis revealed 21 specific immune cell types crucial for affected organs in JIA, indicating their potential contribution to the disease. Finally, 24 known or candidate druggable target genes were prioritized. CONCLUSIONS: Our identification of four novel JIA associated genes, CD247, RHOH, COLEC10 and IRF8, broadens novel potential drug repositioning opportunities. We established a new genetic link between COLEC10, TNFRSF11B and JIA/BMD. Additionally, the identification of RHOH underscores its role in positive thymocyte selection, thereby illuminating a critical facet of JIA's underlying biological mechanisms.

Shared molecular mechanisms and transdiagnostic potential of neurodevelopmental disorders and immune disorders.

The co-occurrence and familial clustering of neurodevelopmental disorders and immune disorders suggest shared genetic risk factors. Based on genome-wide association summary statistics from five neurodevelopmental disorders and four immune disorders, we conducted genome-wide, local genetic correlation and polygenic overlap analysis. We further performed a cross-trait GWAS meta-analysis. Pleotropic loci shared between the two categories of diseases were mapped to candidate genes using multiple algorithms and approaches. Significant genetic correlations were observed between neurodevelopmental disorders and immune disorders, including both positive and negative correlations. Neurodevelopmental disorders exhibited higher polygenicity compared to immune disorders. Around 50%-90% of genetic variants of the immune disorders were shared with neurodevelopmental disorders. The cross-trait meta-analysis revealed 154 genome-wide significant loci, including 8 novel pleiotropic loci. Significant associations were observed for 30 loci with both types of diseases. Pathway analysis on the candidate genes at these loci revealed common pathways shared by the two types of diseases, including neural signaling, inflammatory response, and PI3K-Akt signaling pathway. In addition, 26 of the 30 lead SNPs were associated with blood cell traits. Neurodevelopmental disorders exhibit complex polygenic architecture, with a subset of individuals being at a heightened genetic risk for both neurodevelopmental and immune disorders. The identification of pleiotropic loci has important implications for exploring opportunities for drug repurposing, enabling more accurate patient stratification, and advancing genomics-informed precision in the medical field of neurodevelopmental disorders.

Microcystic lymphatic malformations in Turner syndrome are due to somatic mosaicism of PIK3CA.

Turner syndrome (45,X) is caused by a complete or partial absence of a single X chromosome. Vascular malformations occur due to abnormal development of blood and/or lymphatic vessels. They arise from either somatic or germline pathogenic variants in the genes regulating growth and apoptosis of vascular channels. Aortic abnormalities are a common, known vascular anomaly of Turner syndrome. However, previous studies have described other vascular malformations as a rare feature of Turner syndrome and suggested that vascular abnormalities in individuals with Turner syndrome may be more generalized. In this study, we describe two individuals with co-occurrence of Turner syndrome and vascular malformations with a lymphatic component. In these individuals, genetic testing of the lesional tissue revealed a somatic pathogenic variant in PIK3CA-a known and common cause of lymphatic malformations. Based on this finding, we conclude that the vascular malformations presented here and likely those previously in the literature are not a rare part of the clinical spectrum of Turner syndrome, but rather a separate clinical entity that may or may not co-occur in individuals with Turner syndrome.