Regularized estimation in sparse high-dimensional multivariate regression, with application to a DNA methylation study

In this article, we consider variable selection for correlated high dimensional DNA methylation markers as multivariate outcomes. A novel weighted square-root LASSO procedure is proposed to estimate the regression coefficient matrix. A key feature of this method is tuning-insensitivity, which greatly simplifies the computation by obviating cross validation for penalty parameter selection. A precision matrix obtained via the constrained ℓ1 minimization method is used to account for the within-subject correlation among multivariate outcomes. Oracle inequalities of the regularized estimators are derived. The performance of our proposed method is illustrated via extensive simulation studies. We apply our method to study the relation between smoking and high dimensional DNA methylation markers in the Normative Aging Study (NAS).

Whole-exome sequencing is a valuable diagnostic tool for inherited peripheral neuropathies: Outcomes from a cohort of 50 families

The inherited peripheral neuropathies (IPNs) are characterized by marked clinical and genetic heterogeneity and include relatively frequent presentations such as Charcot-Marie-Tooth disease and hereditary motor neuropathy, as well as more rare conditions where peripheral neuropathy is associated with additional features. There are over 250 genes known to cause IPN-related disorders but it is estimated that in approximately 50% of affected individuals a molecular diagnosis is not achieved. In this study, we examine the diagnostic utility of whole-exome sequencing (WES) in a cohort of 50 families with 1 or more affected individuals with a molecularly undiagnosed IPN with or without additional features. Pathogenic or likely pathogenic variants in genes known to cause IPN were identified in 24% (12/50) of the families. A further 22% (11/50) of families carried sequence variants in IPN genes in which the significance remains unclear. An additional 12% (6/50) of families had variants in novel IPN candidate genes, 3 of which have been published thus far as novel discoveries (KIF1A, TBCK, and MCM3AP). This study highlights the use of WES in the molecular diagnostic approach of highly heterogeneous disorders, such as IPNs, places it in context of other published neuropathy cohorts, while further highlighting associated benefits for discovery.

Noncompaction cardiomyopathy in an infant with Walker-Warburg syndrome

Walker-Warburg syndrome (WWS) is a rare autosomal recessive, congenital muscular dystrophy that is associated with brain and eye anomalies. Several genes encoding proteins involved in α-dystroglycan glycosylation have been implicated in the aetiology of WWS. We describe a patient with nonclassical features of WWS presenting with heart failure related to noncompaction cardiomyopathy resulting in death at 4 months of age. Muscle biopsy revealed absent α-dystroglycan on immunostaining and genetic testing confirmed the diagnosis with two previously described POMT2 mutations. This is the first reported case of WWS syndrome associated with noncompaction cardiomyopathy.

Impact of diabetes in the Friedreich ataxia clinical outcome measures study

OBJECTIVE

Friedreich ataxia (FA) is a progressive neuromuscular disorder caused by GAA triplet repeat expansions or point mutations in the FXN gene. FA is associated with increased risk of diabetes mellitus (DM). This study assessed the age-specific prevalence of FA-associated DM and its impact on neurologic outcomes.

RESEARCH DESIGN AND METHODS

Participants were 811 individuals with FA from 12 international sites in a prospective natural history study (FA Clinical Outcome Measures Study, FACOMS). Physical function was assessed, using validated instruments. Multivariable regression analyses examined the independent association of DM with outcomes.

RESULTS

Mean age of participants was 30.1 years (SD 15.3, range: 7-82), 50% were female, and 94% were non-Hispanic white. 9% (42/459) of adults and 3% (10/352) of children had DM. Individuals with FA-associated DM were older (P < 0.001), had longer GAA repeat length on the least affected FXN allele (P = 0.037), and more severe FA (P = 0.0001). Of individuals with DM, 65% (34/52) were taking insulin. Even after accounting statistically for both age and GAA repeat length, DM was independently associated with greater FA symptom burden (P = 0.010), reduced capacity to perform activities of daily living (P = 0.021), and a decrease of 0.33 SDs on a composite performance measure (95% CI: -0.56-0.11, P = 0.004); the relative impact of DM was most apparent in younger individuals.

CONCLUSIONS

DM-associated FA has an independent adverse impact on well-being in affected individuals, particularly at younger ages. In future, evidence-based approaches for identification and management of FA-related DM may improve both health and function.

Ultra-high dimensional variable selection with application to normative aging study: DNA methylation and metabolic syndrome

Background

Metabolic syndrome has become a major public health challenge worldwide. The association between metabolic syndrome and DNA methylation is of great research interest.

Results

We constructed a binomial model to investigate the association between a metabolic syndrome index and DNA methylation in the Normative Aging Study. We applied the Iterative Sure Independence Screening (ISIS) method with elastic net penalty to DNA methylation levels at 484,548 CpG markers from 659 human subjects, and demonstrated that the screening step in ISIS can significantly improve the performance of the elastic net.

Conclusion

The proposed method identifies four CpGs which can be mapped to two biologically relevant and functional genes. Identification of significant CpG markers may potentially have practical implications for disease prevention and treatment.

Interferon-Stimulated Gene Expression as a Preferred Biomarker for Disease Activity in Aicardi-Goutières Syndrome

Aicardi-Goutières syndrome (AGS) is an early-onset, genetic disease characterized by recurrent fever, multifocal lesions of the brain, and systemic autoimmunity. We report on 3 AGS patients, 2 siblings with an RNASEH2A gene mutation and 1 patient with a SAMHD1 gene mutation. Serial analysis of peripheral blood from all 3 AGS patients showed consistently elevated expression of the interferon-stimulated genes (ISGs): ISG15, RSAD2, and IFI27, not observed in unaffected family members. Enumeration of circulating white blood cells and platelets and examination of C-reactive protein showed no significant deviation from the normal range for Patient 2 with the RNASEH2A mutation and Patient 3 with the SAMHD1 mutation, even when Patient 2 had magnetic resonance imaging abnormalities and ongoing febrile episodes. Erythrocyte sedimentation rates fluctuated within the normal range for Patient 2, with some elevation, yet, were in the normal range during the second febrile episode when there were accompanying neurological abnormalities. These preliminary data suggest that ISG expression may be a more specific indicator of disease activity in comparison to standard inflammatory markers.

Congenital myopathy with "corona" fibres, selective muscle atrophy, and craniosynostosis associated with novel recessive mutations in SCN4A

We describe two brothers with lower facial weakness, highly arched palate, scaphocephaly due to synostosis of the sagittal and metopic sutures, axial hypotonia, proximal muscle weakness, and mild scoliosis. The muscle MRI of the younger sibling revealed a selective pattern of atrophy of the gluteus maximus, adductor magnus and soleus muscles. Muscle biopsy of the younger sibling revealed myofibres with internalized nuclei, myofibrillar disarray, and "corona" fibres. Both affected siblings were found to be compound heterozygous for c.3425G>A (p.Arg1142Gln) and c.1123T>C (p.Cys375Arg) mutations in SCN4A on exome sequencing, and the parents were confirmed carriers of one of the mutations. Electrophysiological characterization of the mutations revealed the Cys375Arg confers full and Arg1142Gln mild partial loss-of-function. Loss of function of the Na_v1.4 channel leads to a decrement of the action potential and subsequent reduction of muscle contraction. The unusual muscle biopsy features suggest a more complex pathomechanism, and broaden the phenotype associated with SCN4A mutations.

XRCC1 mutation is associated with PARP1 hyperactivation and cerebellar ataxia

XRCC1 is a molecular scaffold protein that assembles multi-protein complexes involved in DNA single-strand break repair. Here we show that biallelic mutations in the human XRCC1 gene are associated with ocular motor apraxia, axonal neuropathy, and progressive cerebellar ataxia. Cells from a patient with mutations in XRCC1 exhibited not only reduced rates of single-strand break repair but also elevated levels of protein ADP-ribosylation. This latter phenotype is recapitulated in a related syndrome caused by mutations in the XRCC1 partner protein PNKP and implicates hyperactivation of poly(ADP-ribose) polymerase/s as a cause of cerebellar ataxia. Indeed, remarkably, genetic deletion of Parp1 rescued normal cerebellar ADP-ribose levels and reduced the loss of cerebellar neurons and ataxia in Xrcc1-defective mice, identifying a molecular mechanism by which endogenous single-strand breaks trigger neuropathology. Collectively, these data establish the importance of XRCC1 protein complexes for normal neurological function and identify PARP1 as a therapeutic target in DNA strand break repair-defective disease.

Clinical and genetic study of hereditary spastic paraplegia in Canada

Objective:

To describe the clinical, genetic, and epidemiologic features of hereditary spastic paraplegia (HSP) in Canada and to determine which clinical, radiologic, and genetic factors determine functional outcomes for patients with HSP.

Methods:

We conducted a multicenter observational study of patients who met clinical criteria for the diagnosis of HSP in the provinces of Alberta, Ontario, and Quebec from 2012 to 2015. Characteristics of the participants were analyzed using descriptive statistics. The main outcome measure for a subset of the cohort (n = 48) was the Spastic Paraplegia Rating Scale. We also used the SPATAX-EUROSPA disability stage (disability score) to assess disability (n = 65).

Results:

A total of 526 patients were identified with HSP across the country, and 150 patients had a confirmed genetic diagnosis. Mutations were identified in 15 different genes; the most common were SPAST (SPG4, 48%), ATL1 (SPG3A, 16%), SPG11 (8%), SPG7 (7%), and KIAA0196 (SPG8, 5%). The diagnosis of SPG4 was associated with older age at symptom onset (p = 0.0017). SPG4 and SPG3A were less associated with learning disabilities compared to other subtypes of HSP, and SPG11 was strongly associated with progressive cognitive deficits (odds ratio 87.75, 95% confidence interval 14.04–548.24, p < 0.0001). SPG3A was associated with better functional outcomes compared to other HSP subtypes (p = 0.04) on multivariate analysis. The strongest predictor of significant disability was abnormal brain MRI (p = 0.014).

Conclusions:

The most important predictors of disability in our HSP cohort were SPG11 mutations and abnormal brain MRI. Accurate molecular characterization of well-phenotyped cohorts and international collaboration are essential to establish the natural history of these rare neurodegenerative disorders.

Validation of the finding of hypertrophy of the clava in infantile neuroaxonal dystrophy/PLA2G6 by biometric analysis

INTRODUCTION

Infantile neuroaxonal dystrophy (INAD), an autosomal recessive neurodegenerative disorder due to PLA2G6 mutation, is classified both as a PLA2G6-associated neurodegeneration (PLAN) disorder and as one of the neurodegeneration with brain iron accumulation (NBIA) disorders. Age of onset and clinical presentation in INAD is variable. Typically described imaging features of cerebellar atrophy, cerebellar cortex bright FLAIR signal, and globus pallidus iron deposition are variable or late findings. We characterize clinical and neuroimaging phenotypes in nine children with confirmed PLA2G6 mutations and show a useful imaging feature, clava hypertrophy, which may aid in earlier identification of patients. Measurements of the clava confirm actual enlargement, rather than apparent enlargement due to volume loss of the other brain stem structures.

METHODS

A retrospective clinical and MRI review was performed. Brain stem measurements were performed and compared with age-matched controls.

RESULTS

We identified nine patients, all with novel PLA2G6 gene mutations. MRI, available in eight, showed clava hypertrophy, regardless of age or the absence of other more typically described neuroimaging findings. Brain autopsy in our cohort confirmed prominent spheroid bodies in the clava nuclei.

CONCLUSION

Clava hypertrophy is an important early imaging feature which may aid in indentification of children who would benefit from specific testing for PLA2G6 mutations.

Comorbid Medical Conditions in Friedreich Ataxia: Association With Inflammatory Bowel Disease and Growth Hormone Deficiency

Friedreich ataxia is a progressive degenerative disease with neurologic and cardiac involvement. This study characterizes comorbid medical conditions in a large cohort of patients with Friedreich ataxia. Patient diagnoses were collected in a large natural history study of 641 subjects. Prevalence of diagnoses in the cohort with Friedreich ataxia was compared with prevalence in the population without Friedreich ataxia. Ten patients (1.6%) had inflammatory bowel disease, 3.5 times more common in this cohort of individuals with Friedreich ataxia than in the general population. Four subjects were growth hormone deficient, reflecting a prevalence in Friedreich ataxia that is 28 times greater than the general population. The present study identifies specific diagnoses not traditionally associated with Friedreich ataxia that are found at higher frequency in this disease. These associations could represent coincidence, shared genetic background, or potentially interactive disease mechanisms with Friedreich ataxia.

Progression of Friedreich ataxia: quantitative characterization over 5 years

OBJECTIVE

Friedreich ataxia (FRDA) is a progressive neurodegenerative disorder of adults and children. This study analyzed neurological outcomes and changes to identify predictors of progression and generate power calculations for clinical trials.

METHODS

Eight hundred and twelve subjects in a natural history study were evaluated annually across 12 sites using the Friedreich Ataxia Rating Scale (FARS), 9-Hole Peg Test, Timed 25-Foot Walk, visual acuity tests, self-reported surveys and disability scales. Cross-sectional outcomes were assessed from recent visits, and longitudinal changes were gaged over 5 years from baseline.

RESULTS

Cross-sectional outcomes correlated with measures of disease severity. Age, genetic severity (guanine-adenine-adenine [GAA] repeat length), and testing site predicted performance. Serial progression was relatively linear using FARS and composite measures of performance, while individual performance outcomes were nonlinear over time. Age strongly predicted change from baseline until removing the effects of baseline FARS scores, when GAA becomes a more important factor. Progression is fastest in younger subjects and subjects with longer GAA repeats. Improved coefficients of variation show that progression results are more reproducible over longer assessment durations.

INTERPRETATION

While age predicted progression speed in simple analyses and may provide an effective way to stratify cohorts, separating the effects of age and genetic severity is difficult. Controlling for baseline severity, GAA is the major determinant of progression rate in FRDA. Clinical trials will benefit from enrollment of younger subjects, and sample size requirements will shrink with longer assessment periods. These findings should prove useful in devising gene therapy trials in the near future.

Indexing Effects of Copy Number Variation on Genes Involved in Developmental Delay

A challenge in clinical genomics is to predict whether copy number variation (CNV) affecting a gene or multiple genes will manifest as disease. Increasing recognition of gene dosage effects in neurodevelopmental disorders prompted us to develop a computational approach based on critical-exon (highly expressed in brain, highly conserved) examination for potential etiologic effects. Using a large CNV dataset, our updated analyses revealed significant (P < 1.64 × 10⁻¹⁵) enrichment of critical-exons within rare CNVs in cases compared to controls. Separately, we used a weighted gene co-expression network analysis (WGCNA) to construct an unbiased protein module from prenatal and adult tissues and found it significantly enriched for critical exons in prenatal (P < 1.15 × 10⁻⁵⁰, OR = 2.11) and adult (P < 6.03 × 10⁻¹⁸, OR = 1.55) tissues. WGCNA yielded 1,206 proteins for which we prioritized the corresponding genes as likely to have a role in neurodevelopmental disorders. We compared the gene lists obtained from critical-exon and WGCNA analysis and found 438 candidate genes associated with CNVs annotated as pathogenic, or as variants of uncertain significance (VOUS), from among 10,619 developmental delay cases. We identified genes containing CNVs previously considered to be VOUS to be new candidate genes for neurodevelopmental disorders (GIT1, MVB12B and PPP1R9A) demonstrating the utility of this strategy to index the clinical effects of CNVs.

Estimating and testing high-dimensional mediation effects in epigenetic studies

MOTIVATION

High-dimensional DNA methylation markers may mediate pathways linking environmental exposures with health outcomes. However, there is a lack of analytical methods to identify significant mediators for high-dimensional mediation analysis.

RESULTS

Based on sure independent screening and minimax concave penalty techniques, we use a joint significance test for mediation effect. We demonstrate its practical performance using Monte Carlo simulation studies and apply this method to investigate the extent to which DNA methylation markers mediate the causal pathway from smoking to reduced lung function in the Normative Aging Study. We identify 2 CpGs with significant mediation effects.

AVAILABILITY AND IMPLEMENTATION

R package, source code, and simulation study are available at https://github.com/YinanZheng/HIMA CONTACT: lei.liu@northwestern.edu.