Dissecting genetics of spectrum of epilepsies with eyelid myoclonia by exome sequencing.

OBJECTIVE: Epilepsy with eyelid myoclonia (EEM) spectrum is a generalized form of epilepsy characterized by eyelid myoclonia with or without absences, eye closure-induced seizures with electroencephalographic paroxysms, and photosensitivity. Based on the specific clinical features, age at onset, and familial occurrence, a genetic cause has been postulated. Pathogenic variants in CHD2, SYNGAP1, NEXMIF, RORB, and GABRA1 have been reported in individuals with photosensitivity and eyelid myoclonia, but whether other genes are also involved, or a single gene is uniquely linked with EEM, or its subtypes, is not yet known. We aimed to dissect the genetic etiology of EEM. METHODS: We studied a cohort of 105 individuals by using whole exome sequencing. Individuals were divided into two groups: EEM- (isolated EEM) and EEM+ (EEM accompanied by intellectual disability [ID] or any other neurodevelopmental/psychiatric disorder). RESULTS: We identified nine variants classified as pathogenic/likely pathogenic in the entire cohort (8.57%); among these, eight (five in CHD2, one in NEXMIF, one in SYNGAP1, and one in TRIM8) were found in the EEM+ subcohort (28.57%). Only one variant (IFIH1) was found in the EEM- subcohort (1.29%); however, because the phenotype of the proband did not fit with published data, additional evidence is needed before considering IFIH1 variants and EEM- an established association. Burden analysis did not identify any single burdened gene or gene set. SIGNIFICANCE: Our results suggest that for EEM, as for many other epilepsies, the identification of a genetic cause is more likely with comorbid ID and/or other neurodevelopmental disorders. Pathogenic variants were mostly found in CHD2, and the association of CHD2 with EEM+ can now be considered a reasonable gene-disease association. We provide further evidence to strengthen the association of EEM+ with NEXMIF and SYNGAP1. Possible new associations between EEM+ and TRIM8, and EEM- and IFIH1, are also reported. Although we provide robust evidence for gene variants associated with EEM+, the core genetic etiology of EEM- remains to be elucidated.

SRSF1 haploinsufficiency is responsible for a syndromic developmental disorder associated with intellectual disability.

SRSF1 (also known as ASF/SF2) is a non-small nuclear ribonucleoprotein (non-snRNP) that belongs to the arginine/serine (R/S) domain family. It recognizes and binds to mRNA, regulating both constitutive and alternative splicing. The complete loss of this proto-oncogene in mice is embryonically lethal. Through international data sharing, we identified 17 individuals (10 females and 7 males) with a neurodevelopmental disorder (NDD) with heterozygous germline SRSF1 variants, mostly de novo, including three frameshift variants, three nonsense variants, seven missense variants, and two microdeletions within region 17q22 encompassing SRSF1. Only in one family, the de novo origin could not be established. All individuals featured a recurrent phenotype including developmental delay and intellectual disability (DD/ID), hypotonia, neurobehavioral problems, with variable skeletal (66.7%) and cardiac (46%) anomalies. To investigate the functional consequences of SRSF1 variants, we performed in silico structural modeling, developed an in vivo splicing assay in Drosophila, and carried out episignature analysis in blood-derived DNA from affected individuals. We found that all loss-of-function and 5 out of 7 missense variants were pathogenic, leading to a loss of SRSF1 splicing activity in Drosophila, correlating with a detectable and specific DNA methylation episignature. In addition, our orthogonal in silico, in vivo, and epigenetics analyses enabled the separation of clearly pathogenic missense variants from those with uncertain significance. Overall, these results indicated that haploinsufficiency of SRSF1 is responsible for a syndromic NDD with ID due to a partial loss of SRSF1-mediated splicing activity.

Variant-specific changes in RAC3 function disrupt corticogenesis in neurodevelopmental phenotypes.

Variants in RAC3, encoding a small GTPase RAC3 which is critical for the regulation of actin cytoskeleton and intracellular signal transduction, are associated with a rare neurodevelopmental disorder with structural brain anomalies and facial dysmorphism. We investigated a cohort of 10 unrelated participants presenting with global psychomotor delay, hypotonia, behavioural disturbances, stereotyped movements, dysmorphic features, seizures and musculoskeletal abnormalities. MRI of brain revealed a complex pattern of variable brain malformations, including callosal abnormalities, white matter thinning, grey matter heterotopia, polymicrogyria/dysgyria, brainstem anomalies and cerebellar dysplasia. These patients harboured eight distinct de novo RAC3 variants, including six novel variants (NM_005052.3): c.34G > C p.G12R, c.179G > A p.G60D, c.186_188delGGA p.E62del, c.187G > A p.D63N, c.191A > G p.Y64C and c.348G > C p.K116N. We then examined the pathophysiological significance of these novel and previously reported pathogenic variants p.P29L, p.P34R, p.A59G, p.Q61L and p.E62K. In vitro analyses revealed that all tested RAC3 variants were biochemically and biologically active to variable extent, and exhibited a spectrum of different affinities to downstream effectors including p21-activated kinase 1. We then focused on the four variants p.Q61L, p.E62del, p.D63N and p.Y64C in the Switch II region, which is essential for the biochemical activity of small GTPases and also a variation hot spot common to other Rho family genes, RAC1 and CDC42. Acute expression of the four variants in embryonic mouse brain using in utero electroporation caused defects in cortical neuron morphology and migration ending up with cluster formation during corticogenesis. Notably, defective migration by p.E62del, p.D63N and p.Y64C were rescued by a dominant negative version of p21-activated kinase 1. Our results indicate that RAC3 variants result in morphological and functional defects in cortical neurons during brain development through variant-specific mechanisms, eventually leading to heterogeneous neurodevelopmental phenotypes.

Aldo-Keto Reductase 1C1 (AKR1C1) as the First Mutated Gene in a Family with Nonsyndromic Primary Lipedema.

Lipedema is an often underdiagnosed chronic disorder that affects subcutaneous adipose tissue almost exclusively in women, which leads to disproportionate fat accumulation in the lower and upper body extremities. Common comorbidities include anxiety, depression, and pain. The correlation between mood disorder and subcutaneous fat deposition suggests the involvement of steroids metabolism and neurohormones signaling, however no clear association has been established so far. In this study, we report on a family with three patients affected by sex-limited autosomal dominant nonsyndromic lipedema. They had been screened by whole exome sequencing (WES) which led to the discovery of a missense variant p.(Leu213Gln) in AKR1C1, the gene encoding for an aldo-keto reductase catalyzing the reduction of progesterone to its inactive form, 20-α-hydroxyprogesterone. Comparative molecular dynamics simulations of the wild-type vs. variant enzyme, corroborated by a thorough structural and functional bioinformatic analysis, suggest a partial loss-of-function of the variant. This would result in a slower and less efficient reduction of progesterone to hydroxyprogesterone and an increased subcutaneous fat deposition in variant carriers. Overall, our results suggest that AKR1C1 is the first candidate gene associated with nonsyndromic lipedema.

An Alu-mediated duplication in NMNAT1, involved in NAD biosynthesis, causes a novel syndrome, SHILCA, affecting multiple tissues and organs.

We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients' fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.

Sinus pericranii, skull defects, and structural brain anomalies in TRAF7-related disorder.

BACKGROUND: Several somatic mutations in TRAF7 have been reported in cancers, whereas a few germline heterozygous mutations have been recently linked to a neurodevelopmental disorder, characterized by craniofacial dysmorphisms, congenital heart defects, and digital anomalies. CASES: We report two subjects harboring de novo heterozygous missense variants in TRAF7, namely the recurrent 1964G>A(p.Arg655Gln) and the novel missense c.1204C>G(p.Leu402Val) variants. In addition to the typical hallmarks of the TRAF7-related disorder, both subjects presented with a recognizable "pear-shaped" skull due to multiple craniosynostosis, sinus pericranii, skull base/cranio-cervical junction anomalies, dysgyria, and inferior cerebellar vermis hypoplasia. CONCLUSIONS: Hence, we expand the genotypic and phenotypic spectrum of this neurodevelopmental disorder, discussing possible implications for clinical management of subjects with germline TRAF7 mutations.

Expansion of the phenotype of lateral meningocele syndrome.

Lateral meningocele syndrome (LMS) is due to specific pathogenic variants in the last exon of NOTCH3 gene. Besides the lateral meningoceles, this condition presents with dysmorphic features, short stature, congenital heart defects, and feeding difficulties. Here, we report a girl with neurosensorial hearing loss, severe gastroesophageal reflux disease, congenital heart defects, multiple renal cysts, kyphosis and left-convex scoliosis, dysmorphic features, and mild developmental delay. Exome sequencing detected the previously unreported de novo loss-of-function variant in exon 33 of NOTCH3 p.(Lys2137fs). Following the identification of the gene defect, MRI of the brain and spine revealed temporal encephaloceles, inner ears anomalies, multiple spinal lateral meningoceles, and intra- and extra-dural arachnoid spinal cysts. This case illustrates the power of reverse phenotyping to establish clinical diagnosis and expands the spectrum of clinical manifestations related to LMS to include inner ear abnormalities and multi-cystic kidney disease.

Identifying Fabry patients in dialysis population: prevalence of GLA mutations by renal clinic screening, 1995-2019.

BACKGROUND: Fabry disease (FD) is a rare X-linked genetic disorder of glycosphingolipid catabolism caused by mutations in the GLA gene. Its heterogeneous presentation, the paucity of specific early markers, and the absence of a genotype-phenotype correlation are associated with a delayed or missed diagnosis. The true prevalence of FD remains so far unknown. METHODS: A systematic search of FD screening studies in dialysis patients published from January 1995 until January 2019 was performed to reanalyze the prevalence of GLA mutations in this population after assigning their correct phenotype. RESULTS: Twenty five screening studies involving 39,621 dialysis patients were included. Of them, 116 [91 males (0.23%) and 25 females (0.06%)] were positive to the GLA sequencing analysis. 56 (48.2%) had benign variant, 52 (44.8%) a pathogenic GLA mutation (39 classic and 13 late onset mutations) and 8 (6.9%) a mutation of uncertain significance. The overall prevalence of GLA variants was 0.24% [CI 95%, 0.17-0.32] while the overall prevalence recalculated on basis of only pathogenetic mutations was 0.14% [CI 95%, 0.08-0.20]. This difference was significant (P = 0.048). CONCLUSIONS: Although the real prevalence of classic FD is low, the screening in the high-risk renal population remains of primary interest as an early diagnosis is fundamental for a timely specific therapy; moreover, the identification of index cases could allow patients' relatives to be investigated and promptly treated.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes.

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

Three de novo DDX3X variants associated with distinctive brain developmental abnormalities and brain tumor in intellectually disabled females.

De novo DDX3X variants account for 1-3% of syndromic intellectual disability (ID) in females and have been occasionally reported in males. Furthermore, somatic DDX3X variants occur in several aggressive cancers, including medulloblastoma. We report three unrelated females with severe ID, dysmorphic features, and a common brain malformative pattern characterized by malformations of cortical development, callosal dysgenesis, basal ganglia anomalies, and midbrain-hindbrain malformations. A pilocytic astrocytoma was incidentally diagnosed in Patient 1 and trigonocephaly was found in Patient 2. With the use of family based whole exome sequencing (WES), we identified three distinct de novo variants in DDX3X. These findings expand the phenotypic spectrum of DDX3X-related disorders, demonstrating unique neuroradiological features resembling those of the tubulinopathies, and support a role for DDX3X in neuronal development. Our observations further suggest a possible link between germline DDX3X variants and cancer development.

Gait disturbance and lower limb pain in a patient with PIK3CA-related disorder.

Post-zygotic activating mutations in PIK3CA and other genes encoding members of PI3K-AKT-mTOR pathway have been found in various overgrowth syndromes that have been grouped together as PIK3CA-related overgrowth spectrum (PROS). We report a female patient with gait disturbance, leg pain, isolated macrodactyly of the foot, and mild intellectual disability. Imaging of the lower limb showed a lipoblastoma of the right thigh. A mosaic gain-of-function mutation in the catalytic domain of PIK3CA (c.3140 A > G; p.His1047Arg) was detected in the adipose tissue and in skin cultured fibroblasts from the macrodactyly but not in blood. The leg pain and the severe walking disturbance improved slightly over time and serial MRI of the lower limbs suggested that the size of the lipoblastoma relative to the lower limb muscles or to the whole lower limb was unchanged as consequence of limb growth. This case report illustrates that pain and gait disturbance can be features of PROS and highlights the need of better knowledge about the natural history of the disease.

The PPARγ2 Pro12Ala variant is protective against progression of nephropathy in people with type 2 diabetes.

OBJECTIVE: Cross-sectional studies suggest the association between diabetic nephropathy and the PPARγ2 Pro12Ala polymorphism of the peroxisome proliferator-activated receptor γ2 (PPARγ2). Prospective data are limited to microalbuminuria and no information on renal function is available to date. The present study evaluates the association between the Pro12Ala polymorphism of PPARγ2 and the progression of albuminuria and decay in glomerular filtration rate (GFR) in type 2 diabetes. PATIENTS AND MEASUREMENTS: We studied 256 patients with an average 5-year follow-up. Among others, urinary albumin excretion rate (UAER) was measured on spot sample, GFR was estimated with the CKD-EPI Equation. RESULTS: Baseline UAER and GFR were similar for carriers or non-carriers of the polymorphism. At follow-up no significant changes from baseline were observed for UAER or eGFR in carriers of the Pro12Ala polymorphism whereas a significant increase in UAER [17 (11.3-37.9) versus 24.5 (13.8-49.9) µg/mg, p < 0.006)] and a significant reduction in the eGFR (82.8 ± 14.5 versus 80.3 ± 17.3 ml/min/1.73, m(2) p = 0.02), were observed in non carriers of the Pro12Ala polymorphism. Progression of nephropathy - defined according to a combined end point of UAER and eGFR- i.e. doubling of baseline UAER to at least 100 µg/mg, or new onset microalbuminuria, or progression from micro to macroalbuminuria, or 25% reduction of eGFR, or annualized eGFR decline >3 ml/min/year - was significantly less frequent in Ala carriers than non carriers (11.4% vs 35.8%; p < 0.01); HR adjusted for baseline age, AER, eGFR, HbA1c, diabetes duration and blood pressure was 0.32 (0.12-0.80). CONCLUSIONS: This study found that among patients with type 2 diabetes, the PPARγ2 Pro12Ala polymorphism is protective against progression of nephropathy and decay of renal function independent of major confounders.

Schizophrenia and vitamin D related genes could have been subject to latitude-driven adaptation.

BACKGROUND: Many natural phenomena are directly or indirectly related to latitude. Living at different latitudes, indeed, has its consequences with being exposed to different climates, diets, light/dark cycles, etc. In humans, one of the best known examples of genetic traits following a latitudinal gradient is skin pigmentation. Nevertheless, also several diseases show latitudinal clinals such as hypertension, cancer, dismetabolic conditions, schizophrenia, Parkinson's disease and many more. RESULTS: We investigated, for the first time on a wide genomic scale, the latitude-driven adaptation phenomena. In particular, we selected a set of genes showing signs of latitude-dependent population differentiation. The biological characterization of these genes showed enrichment for neural-related processes. In light of this, we investigated whether genes associated to neuropsychiatric diseases were enriched by Latitude-Related Genes (LRGs). We found a strong enrichment of LRGs in the set of genes associated to schizophrenia. In an attempt to try to explain this possible link between latitude and schizophrenia, we investigated their associations with vitamin D. We found in a set of vitamin D related genes a significant enrichment of both LRGs and of genes involved in schizophrenia. CONCLUSIONS: Our results suggest a latitude-driven adaptation for both schizophrenia and vitamin D related genes. In addition we confirm, at a molecular level, the link between schizophrenia and vitamin D. Finally, we discuss a model in which schizophrenia is, at least partly, a maladaptive by-product of latitude dependent adaptive changes in vitamin D metabolism.

Uncoupling protein 2 G(-866)A polymorphism: a new gene polymorphism associated with C-reactive protein in type 2 diabetic patients.

BACKGROUND: This study evaluated the relationship between the G(-866)A polymorphism of the uncoupling protein 2 (UCP2) gene and high-sensitivity C reactive protein (hs-CRP) plasma levels in diabetic patients. METHODS: We studied 383 unrelated people with type 2 diabetes aged 40-70 years. Anthropometry, fasting lipids, glucose, HbA1c, and hs-CRP were measured. Participants were genotyped for the G (-866)A polymorphism of the uncoupling protein 2 gene. RESULTS: Hs-CRP (mg/L) increased progressively across the three genotype groups AA, AG, or GG, being respectively 3.0 ± 3.2, 3.6 ± 5.0, and 4.8 ± 5.3 (p for trend = 0.03). Since hs-CRP values were not significantly different between AA and AG genotype, these two groups were pooled for further analyses. Compared to participants with the AA/AG genotypes, homozygotes for the G allele (GG genotype) had significantly higher hs-CRP levels (4.8 ± 5.3 vs 3.5 ± 4.7 mg/L, p = 0.01) and a larger proportion (53.9% vs 46.1%, p = 0.013) of elevated hs-CRP (> 2 mg/L). This was not explained by major confounders such as age, gender, BMI, waist circumference, HbA1c, smoking, or medications use which were comparable in the two genotype groups. CONCLUSIONS: The study shows for the first time, in type 2 diabetic patients, a significant association of hs-CRP levels with the G(-866)A polymorphism of UCP2 beyond the effect of major confounders.

A novel approach to simulate gene-environment interactions in complex diseases.

BACKGROUND: Complex diseases are multifactorial traits caused by both genetic and environmental factors. They represent the major part of human diseases and include those with largest prevalence and mortality (cancer, heart disease, obesity, etc.). Despite a large amount of information that has been collected about both genetic and environmental risk factors, there are few examples of studies on their interactions in epidemiological literature. One reason can be the incomplete knowledge of the power of statistical methods designed to search for risk factors and their interactions in these data sets. An improvement in this direction would lead to a better understanding and description of gene-environment interactions. To this aim, a possible strategy is to challenge the different statistical methods against data sets where the underlying phenomenon is completely known and fully controllable, for example simulated ones. RESULTS: We present a mathematical approach that models gene-environment interactions. By this method it is possible to generate simulated populations having gene-environment interactions of any form, involving any number of genetic and environmental factors and also allowing non-linear interactions as epistasis. In particular, we implemented a simple version of this model in a Gene-Environment iNteraction Simulator (GENS), a tool designed to simulate case-control data sets where a one gene-one environment interaction influences the disease risk. The main aim has been to allow the input of population characteristics by using standard epidemiological measures and to implement constraints to make the simulator behaviour biologically meaningful. CONCLUSIONS: By the multi-logistic model implemented in GENS it is possible to simulate case-control samples of complex disease where gene-environment interactions influence the disease risk. The user has full control of the main characteristics of the simulated population and a Monte Carlo process allows random variability. A knowledge-based approach reduces the complexity of the mathematical model by using reasonable biological constraints and makes the simulation more understandable in biological terms. Simulated data sets can be used for the assessment of novel statistical methods or for the evaluation of the statistical power when designing a study.

PPAR-gamma agonist Azelaoyl PAF increases frataxin protein and mRNA expression: new implications for the Friedreich's ataxia therapy.

Friedreich's ataxia is a neurodegenerative disease due to frataxin deficiency, and thus, drugs increasing the frataxin amount are excellent candidates for therapy. By screening Gene Expression Omnibus profiles, we identified records showing a frataxin response to the peroxisome proliferator-activated receptors gamma (PPAR-gamma) agonist rosiglitazone. We decided to investigate the effect of the PPAR-gamma agonist Azelaoyl PAF on the frataxin protein and mRNA expression profile. We treated human neuroblastoma cells SKNBE and primary fibroblasts from skin biopsies from Friedreich's ataxia (FRDA) patients and healthy controls with the PPAR-gamma agonist Azelaoyl PAF. We show in this paper for the first time that Azelaoyl PAF significantly increases the intracellular frataxin levels by twofold in the neuroblastoma cell line SKNBE and fibroblasts from FRDA patients and that Azelaoyl PAF increases frataxin protein through a transcriptional mechanism. PPAR-gamma agonist Azelaoyl PAF increases both messenger RNA and protein levels of frataxin. We hypothesize that PPAR-gamma agonists could play a role in the treatment of FRDA, and our results offer the logical bases to further investigate the usefulness of this group of agents for the treatment of the FRDA.