Rare and de novo coding variants in chromodomain genes in Chiari I malformation.

Chiari I malformation (CM1), the displacement of the cerebellum through the foramen magnum into the spinal canal, is one of the most common pediatric neurological conditions. Individuals with CM1 can present with neurological symptoms, including severe headaches and sensory or motor deficits, often as a consequence of brainstem compression or syringomyelia (SM). We conducted whole-exome sequencing (WES) on 668 CM1 probands and 232 family members and performed gene-burden and de novo enrichment analyses. A significant enrichment of rare and de novo non-synonymous variants in chromodomain (CHD) genes was observed among individuals with CM1 (combined p = 2.4 × 10-10), including 3 de novo loss-of-function variants in CHD8 (LOF enrichment p = 1.9 × 10-10) and a significant burden of rare transmitted variants in CHD3 (p = 1.8 × 10-6). Overall, individuals with CM1 were found to have significantly increased head circumference (p = 2.6 × 10-9), with many harboring CHD rare variants having macrocephaly. Finally, haploinsufficiency for chd8 in zebrafish led to macrocephaly and posterior hindbrain displacement reminiscent of CM1. These results implicate chromodomain genes and excessive brain growth in CM1 pathogenesis.

Clinical and metabolomic risk factors associated with rapid renal function decline in sickle cell disease.

Sickle cell disease (SCD) nephropathy and lower estimated glomerular filtration rate (eGFR) are risk factors for early mortality. Furthermore, rate of eGFR decline predicts progression to end-stage renal disease in many clinical settings. However, factors predicting renal function decline in SCD are poorly documented. Using clinical, laboratory, genetic, and metabolomic data, we evaluated predictors of renal function decline in a longitudinal cohort of 288 adults (mean age 33.0 years). In 193 subjects with 5-year follow-up data, mean rate of eGFR decline was 2.35 mL/min/1.73 m2 /year, nearly twice that of African American adults overall. Hyperfiltration was prevalent at baseline (61.1%), and 36.8% of subjects experienced rapid eGFR decline (≥3 mL/min/1.73 m2 /year). Severe Hb genotype; proteinuria; higher platelet and reticulocyte counts, and systolic BP; and lower Hb level and BMI were associated with rapid decline. A risk scoring system was created using these 7 variables and was highly predictive of rapid eGFR decline, with odds of rapid decline increasing 1.635-fold for every point increment (P < 0.0001). Rapid eGFR decline was also associated with higher organ system severity score and peak creatinine. Additionally, two metabolites (asymmetric dimethylarginine and quinolinic acid) were associated with rapid decline. Further investigation into longitudinal SCD nephropathy (SCDN) trajectory, early markers of SCDN, and tools for risk stratification should inform interventional studies targeted to slowing GFR decline and improving SCD outcomes.

In vivo Modeling Implicates APOL1 in Nephropathy: Evidence for Dominant Negative Effects and Epistasis under Anemic Stress.

African Americans have a disproportionate risk for developing nephropathy. This disparity has been attributed to coding variants (G1 and G2) in apolipoprotein L1 (APOL1); however, there is little functional evidence supporting the role of this protein in renal function. Here, we combined genetics and in vivo modeling to examine the role of apol1 in glomerular development and pronephric filtration and to test the pathogenic potential of APOL1 G1 and G2. Translational suppression or CRISPR/Cas9 genome editing of apol1 in zebrafish embryos results in podocyte loss and glomerular filtration defects. Complementation of apol1 morphants with wild-type human APOL1 mRNA rescues these defects. However, the APOL1 G1 risk allele does not ameliorate defects caused by apol1 suppression and the pathogenicity is conferred by the cis effect of both individual variants of the G1 risk haplotype (I384M/S342G). In vivo complementation studies of the G2 risk allele also indicate that the variant is deleterious to protein function. Moreover, APOL1 G2, but not G1, expression alone promotes developmental kidney defects, suggesting a possible dominant-negative effect of the altered protein. In sickle cell disease (SCD) patients, we reported previously a genetic interaction between APOL1 and MYH9. Testing this interaction in vivo by co-suppressing both transcripts yielded no additive effects. However, upon genetic or chemical induction of anemia, we observed a significantly exacerbated nephropathy phenotype. Furthermore, concordant with the genetic interaction observed in SCD patients, APOL1 G2 reduces myh9 expression in vivo, suggesting a possible interaction between the altered APOL1 and myh9. Our data indicate a critical role for APOL1 in renal function that is compromised by nephropathy-risk encoding variants. Moreover, our interaction studies indicate that the MYH9 locus is also relevant to the phenotype in a stressed microenvironment and suggest that consideration of the context-dependent functions of both proteins will be required to develop therapeutic paradigms.

Joint eQTL assessment of whole blood and dura mater tissue from individuals with Chiari type I malformation.

BACKGROUND: Expression quantitative trait loci (eQTL) play an important role in the regulation of gene expression. Gene expression levels and eQTLs are expected to vary from tissue to tissue, and therefore multi-tissue analyses are necessary to fully understand complex genetic conditions in humans. Dura mater tissue likely interacts with cranial bone growth and thus may play a role in the etiology of Chiari Type I Malformation (CMI) and related conditions, but it is often inaccessible and its gene expression has not been well studied. A genetic basis to CMI has been established; however, the specific genetic risk factors are not well characterized. RESULTS: We present an assessment of eQTLs for whole blood and dura mater tissue from individuals with CMI. A joint-tissue analysis identified 239 eQTLs in either dura or blood, with 79% of these eQTLs shared by both tissues. Several identified eQTLs were novel and these implicate genes involved in bone development (IPO8, XYLT1, and PRKAR1A), and ribosomal pathways related to marrow and bone dysfunction, as potential candidates in the development of CMI. CONCLUSIONS: Despite strong overall heterogeneity in expression levels between blood and dura, the majority of cis-eQTLs are shared by both tissues. The power to detect shared eQTLs was improved by using an integrative statistical approach. The identified tissue-specific and shared eQTLs provide new insight into the genetic basis for CMI and related conditions.

Gene-centric association study of acute chest syndrome and painful crisis in sickle cell disease patients.

Patients with sickle cell disease (SCD) present with a wide range of clinical complications. Understanding this clinical heterogeneity offers the prospects to tailor the right treatments to the right patients and also guide the development of novel therapies. Several environmental (eg, nutrition) and nonenvironmental (eg, fetal hemoglobin levels, α-thalassemia status) factors are known to modify SCD severity. To find new genetic modifiers of SCD severity, we performed a gene-centric association study in 1514 African American participants from the Cooperative Study of Sickle Cell Disease (CSSCD) for acute chest syndrome (ACS) and painful crisis. From the initial results, we selected 36 single nucleotide polymorphism (SNPs) and genotyped them for replication in 387 independent patients from the CSSCD, 318 SCD patients recruited at Georgia Health Sciences University, and 449 patients from the Duke SCD cohort. In the combined analysis, an association between ACS and rs6141803 reached array-wide significance (P = 4.1 × 10(-7)). This SNP is located 8.2 kilobases upstream of COMMD7, a gene highly expressed in the lung that interacts with nuclear factor-κB signaling. Our results provide new leads to gaining a better understanding of clinical variability in SCD, a "simple" monogenic disease.

Genetic evaluation and application of posterior cranial fossa traits as endophenotypes for Chiari type I malformation.

Chiari Type I Malformation (CMI) is characterized by herniation of the cerebellar tonsils through the base of the skull. Although cerebellar tonsillar herniation (CTH) is hypothesized to result from an underdeveloped posterior cranial fossa (PF), patients are frequently diagnosed by the extent of CTH without cranial morphometric assessment. We recently completed the largest CMI whole genome qualitative linkage screen to date. Despite an initial lack of statistical evidence, stratified analyses using clinical criteria to reduce heterogeneity resulted in a striking increase in evidence for linkage. The present study focused on the use of cranial base morphometrics to further dissect this heterogeneity and increase power to identify disease genes. We characterized the genetic contribution for a series of PF traits and evaluated the use of heritable, disease-relevant PF traits in ordered subset analysis (OSA). Consistent with a genetic hypothesis for CMI, much of the PF morphology was found to be heritable and multiple genomic regions were strongly implicated from OSA, including regions on Chromosomes 1 (LOD = 3.07, p = 3 × 10(-3) ) and 22 (LOD = 3.45, p = 6 × 10(-5) ) containing several candidates warranting further investigation. This study underscores the genetic heterogeneity of CMI and the utility of PF traits in CMI genetic studies.

Missing genetic risk in neural tube defects: can exome sequencing yield an insight?

BACKGROUND: Neural tube defects (NTD) have a strong genetic component, with up to 70% of variance in human prevalence determined by heritable factors. Although the identification of causal DNA variants by sequencing candidate genes from functionally relevant pathways and model organisms has provided some success, alternative approaches are demanded. METHODS: Next generation sequencing platforms are facilitating the production of massive amounts of sequencing data, primarily from the protein coding regions of the genome, at a faster rate and cheaper cost than has previously been possible. These platforms are permitting the identification of variants (de novo, rare, and common) that are drivers of NYTD etiology, and the cost of the approach allows for the screening of increased numbers of affected and unaffected individuals from NTD families and in simplex cases. CONCLUSION: The next generation sequencing platforms represent a powerful tool in the armory of the genetics researcher to identify the causal genetic basis of NTDs.

APOL1-mediated monovalent cation transport contributes to APOL1-mediated podocytopathy in kidney disease.

Two coding variants of apolipoprotein L1 (APOL1), called G1 and G2, explain much of the excess risk of kidney disease in African Americans. While various cytotoxic phenotypes have been reported in experimental models, the proximal mechanism by which G1 and G2 cause kidney disease is poorly understood. Here, we leveraged 3 experimental models and a recently reported small molecule blocker of APOL1 protein, VX-147, to identify the upstream mechanism of G1-induced cytotoxicity. In HEK293 cells, we demonstrated that G1-mediated Na+ import/K+ efflux triggered activation of GPCR/IP3-mediated calcium release from the ER, impaired mitochondrial ATP production, and impaired translation, which were all reversed by VX-147. In human urine-derived podocyte-like epithelial cells (HUPECs), we demonstrated that G1 caused cytotoxicity that was again reversible by VX-147. Finally, in podocytes isolated from APOL1 G1 transgenic mice, we showed that IFN-γ-mediated induction of G1 caused K+ efflux, activation of GPCR/IP3 signaling, and inhibition of translation, podocyte injury, and proteinuria, all reversed by VX-147. Together, these results establish APOL1-mediated Na+/K+ transport as the proximal driver of APOL1-mediated kidney disease.

Genetic association analyses of nitric oxide synthase genes and neural tube defects vary by phenotype.

Neural tube defects (NTDs) are caused by improper neural tube closure during the early stages of embryonic development. NTDs are hypothesized to have a complex genetic origin and numerous candidate genes have been proposed. The nitric oxide synthase 3 (NOS3) G594T polymorphism has been implicated in risk for spina bifida, and interactions between that single nucleotide polymorphism (SNP) and the methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism have also been observed. To evaluate other genetic variation in the NO pathway in the development of NTDs, we examined all three NOS genes: NOS1, NOS2, and NOS3. Using 3109 Caucasian samples in 745 families, we evaluated association in the overall dataset and within specific phenotypic subsets. Haplotype tagging SNPs in the NOS genes were tested for genetic association with NTD subtypes, both for main effects as well as for the presence of interactions with the MTHFR C677T polymorphism. Nominal main effect associations were found with all subtypes, across all three NOS genes, and interactions were observed between SNPs in all three NOS genes and MTHFR C677T. Unlike the previous report, the most significant associations in our dataset were with cranial subtypes and the AG genotype of rs4795067 in NOS2 (p = 0.0014) and the interaction between the rs9658490 G allele in NOS1 and MTHFR 677TT genotype (p = 0.0014). Our data extend the previous findings by implicating a role for all three NOS genes, independently and through interactions with MTHFR, in risk not only for spina bifida, but all NTD subtypes.

Genome-wide meta-analysis identifies new candidate genes for sickle cell disease nephropathy.

Sickle cell disease nephropathy (SCDN), a common SCD complication, is strongly associated with mortality. Polygenic risk scores calculated from recent transethnic meta-analyses of urinary albumin-to-creatinine ratio and estimated glomerular filtration rate (eGFR) trended toward association with proteinuria and eGFR in SCD but the model fit was poor (R2 < 0.01), suggesting that there are likely unique genetic risk factors for SCDN. Therefore, we performed genome-wide association studies (GWAS) for 2 critical manifestations of SCDN, proteinuria and decreased eGFR, in 2 well-characterized adult SCD cohorts, representing, to the best of our knowledge, the largest SCDN sample to date. Meta-analysis identified 6 genome-wide significant associations (false discovery rate, q ≤ 0.05): 3 for proteinuria (CRYL1, VWF, and ADAMTS7) and 3 for eGFR (LRP1B, linc02288, and FPGT-TNNI3K/TNNI3K). These associations are independent of APOL1 risk and represent novel SCDN loci, many with evidence for regulatory function. Moreover, GWAS SNPs in CRYL1, VWF, ADAMTS7, and linc02288 are associated with gene expression in kidney and pathways important to both renal function and SCD biology, supporting the hypothesis that SCDN pathophysiology is distinct from other forms of kidney disease. Together, these findings provide new targets for functional follow-up that could be tested prospectively and potentially used to identify patients with SCD who are at risk, before onset of kidney dysfunction.

JAK inhibitor blocks COVID-19 cytokine-induced JAK/STAT/APOL1 signaling in glomerular cells and podocytopathy in human kidney organoids.

COVID-19 infection causes collapse of glomerular capillaries and loss of podocytes, culminating in a severe kidney disease called COVID-19-associated nephropathy (COVAN). The underlying mechanism of COVAN is unknown. We hypothesized that cytokines induced by COVID-19 trigger expression of pathogenic APOL1 via JAK/STAT signaling, resulting in podocyte loss and COVAN phenotype. Here, based on 9 biopsy-proven COVAN cases, we demonstrated for the first time, to the best of our knowledge, that APOL1 protein was abundantly expressed in podocytes and glomerular endothelial cells (GECs) of COVAN kidneys but not in controls. Moreover, a majority of patients with COVAN carried 2 APOL1 risk alleles. We show that recombinant cytokines induced by SARS-CoV-2 acted synergistically to drive APOL1 expression through the JAK/STAT pathway in primary human podocytes, GECs, and kidney micro-organoids derived from a carrier of 2 APOL1 risk alleles, but expression was blocked by a JAK1/2 inhibitor, baricitinib. We demonstrate that cytokine-induced JAK/STAT/APOL1 signaling reduced the viability of kidney organoid podocytes but was rescued by baricitinib. Together, our results support the conclusion that COVID-19-induced cytokines are sufficient to drive COVAN-associated podocytopathy via JAK/STAT/APOL1 signaling and that JAK inhibitors could block this pathogenic process. These findings suggest JAK inhibitors may have therapeutic benefits for managing cytokine-induced, APOL1-mediated podocytopathy.

MYH9 and APOL1 are both associated with sickle cell disease nephropathy.

Renal failure occurs in 5-18% of sickle cell disease (SCD) patients and is associated with early mortality. At-risk SCD patients cannot be identified prior to the appearance of proteinuria and the pathobiology is not well understood. The myosin, heavy chain 9, non-muscle (MYH9) and apolipoprotein L1 (APOL1) genes have been associated with risk for focal segmental glomerulosclerosis and end-stage renal disease in African Americans. We genotyped 26 single nucleotide polymorphisms (SNPs) in MYH9 and 2 SNPs in APOL1 (representing the G1 and G2 tags) in 521 unrelated adult (18-83 years) SCD patients screened for proteinuria. Using logistic regression, SNPs were evaluated for association with proteinuria. Seven SNPs in MYH9 and one in APOL1 remained significantly associated with proteinuria after multiple testing correction (P < 0·0025). An MYH9 risk haplotype (P = 0·001) and the APOL1 G1/G2 recessive model (P < 0·0001) were strongly associated with proteinuria, even when accounting for the other. Glomerular filtration rate was negatively correlated with proteinuria (P < 0·0001), and was significantly predicted by an interaction between MYH9 and APOL1 in age-adjusted analyses. Our data provide insight into the pathobiology of renal dysfunction in SCD, suggesting that MYH9 and APOL1 are both associated with risk.

Rare functional genetic variants in COL7A1, COL6A5, COL1A2 and COL5A2 frequently occur in Chiari Malformation Type 1.

Chiari Malformation Type 1 (CM-1) is characterized by herniation of the cerebellar tonsils below the foramen magnum and the presence of headaches and other neurologic symptoms. Cranial bone constriction is suspected to be the most common biologic mechanism leading to CM-1. However, other mechanisms may also contribute, particularly in the presence of connective tissue disorders (CTDs), such as Ehlers Danlos Syndrome (EDS). Accumulating data suggest CM-1 with connective tissue disorders (CTD+) may have a different patho-mechanism and different genetic risk factors than CM-1 without CTDs (CTD-). To identify CM-1 genetic risk variants, we performed whole exome sequencing on a single large, multiplex family from Spain and targeted sequencing on a cohort of 186 unrelated adult, Caucasian females with CM-1. Targeted sequencing captured the coding regions of 21 CM-1 and EDS candidate genes, including two genes identified in the Spanish family. Using gene burden analysis, we compared the frequency of rare, functional variants detected in CM-1 cases versus publically available ethnically-matched controls from gnomAD. A secondary analysis compared the presence of rare variants in these genes between CTD+ and CTD- CM-1 cases. In the Spanish family, rare variants co-segregated with CM-1 in COL6A5, ADGRB3 and DST. A variant in COL7A1 was present in affected and unaffected family members. In the targeted sequencing analysis, rare variants in six genes (COL7A1, COL5A2, COL6A5, COL1A2, VEGFB, FLT1) were significantly more frequent in CM-1 cases compared to public controls. In total, 47% of CM-1 cases presented with rare variants in at least one of the four significant collagen genes and 10% of cases harbored variants in multiple significant collagen genes. Moreover, 26% of CM-1 cases presented with rare variants in the COL6A5 gene. We also identified two genes (COL7A1, COL3A1) for which the burden of rare variants differed significantly between CTD+ and CTD- CM-1 cases. A higher percentage of CTD+ patients had variants in COL7A1 compared to CTD+ patients, while CTD+ patients had fewer rare variants in COL3A1 than did CTD- patients. In summary, rare variants in several collagen genes are particularly frequent in CM-1 cases and those in COL6A5 co-segregated with CM-1 in a Spanish multiplex family. COL6A5 has been previously associated with musculoskeletal phenotypes, but this is the first association with CM-1. Our findings underscore the contribution of rare genetic variants in collagen genes to CM-1, and suggest that CM-1 in the presence and absence of CTD symptoms is driven by different genes.

Sex dependent glial-specific changes in the chromatin accessibility landscape in late-onset Alzheimer's disease brains.

BACKGROUND: In the post-GWAS era, there is an unmet need to decode the underpinning genetic etiologies of late-onset Alzheimer's disease (LOAD) and translate the associations to causation. METHODS: We conducted ATAC-seq profiling using NeuN sorted-nuclei from 40 frozen brain tissues to determine LOAD-specific changes in chromatin accessibility landscape in a cell-type specific manner. RESULTS: We identified 211 LOAD-specific differential chromatin accessibility sites in neuronal-nuclei, four of which overlapped with LOAD-GWAS regions (±100 kb of SNP). While the non-neuronal nuclei did not show LOAD-specific differences, stratification by sex identified 842 LOAD-specific chromatin accessibility sites in females. Seven of these sex-dependent sites in the non-neuronal samples overlapped LOAD-GWAS regions including APOE. LOAD loci were functionally validated using single-nuclei RNA-seq datasets. CONCLUSIONS: Using brain sorted-nuclei enabled the identification of sex-dependent cell type-specific LOAD alterations in chromatin structure. These findings enhance the interpretation of LOAD-GWAS discoveries, provide potential pathomechanisms, and suggest novel LOAD-loci.

PCM1 is necessary for focal ciliary integrity and is a candidate for severe schizophrenia.

The neuronal primary cilium and centriolar satellites have functions in neurogenesis, but little is known about their roles in the postnatal brain. We show that ablation of pericentriolar material 1 in the mouse leads to progressive ciliary, anatomical, psychomotor, and cognitive abnormalities. RNAseq reveals changes in amine- and G-protein coupled receptor pathways. The physiological relevance of this phenotype is supported by decreased available dopamine D2 receptor (D2R) levels and the failure of antipsychotic drugs to rescue adult behavioral defects. Immunoprecipitations show an association with Pcm1 and D2Rs. Finally, we sequence PCM1 in two human cohorts with severe schizophrenia. Systematic modeling of all discovered rare alleles by zebrafish in vivo complementation reveals an enrichment for pathogenic alleles. Our data emphasize a role for the pericentriolar material in the postnatal brain, with progressive degenerative ciliary and behavioral phenotypes; and they support a contributory role for PCM1 in some individuals diagnosed with schizophrenia.

RNA sequencing of isolated cell populations expressing human APOL1 G2 risk variant reveals molecular correlates of sickle cell nephropathy in zebrafish podocytes.

Kidney failure occurs in 5-13% of individuals with sickle cell disease and is associated with early mortality. Two APOL1 alleles (G1 and G2) have been identified as risk factors for sickle cell disease nephropathy. Both risk alleles are prevalent in individuals with recent African ancestry and have been associated with nephropathic complications in other diseases. Despite the association of G1 and G2 with kidney dysfunction, the mechanisms by which these variants contribute to increased risk remain poorly understood. Previous work in zebrafish models suggest that the G2 risk allele functions as a dominant negative, whereas the G1 allele is a functional null. To understand better the cellular pathology attributed to APOL1 G2, we investigated the in vivo effects of the G2 risk variant on distinct cell types using RNA sequencing. We surveyed APOL1 G2 associated transcriptomic alterations in podocytes and vascular endothelial cells isolated from zebrafish larvae expressing cell-type specific reporters. Our analysis identified many transcripts (n = 7,523) showing differential expression between APOL1 G0 (human wild-type) and APOL1 G2 exposed podocytes. Conversely, relatively few transcripts (n = 107) were differentially expressed when comparing APOL1 G0 and APOL1 G2 exposed endothelial cells. Pathway analysis of differentially expressed transcripts in podocytes showed enrichment for autophagy associated terms such as "Lysosome" and "Phagosome", implicating these pathways in APOL1 G2 associated kidney dysfunction. This work provides insight into the molecular pathology of APOL1 G2 nephropathy which may offer new therapeutic strategies for multiple disease contexts such as sickle cell nephropathy.