Analysis of non-motor symptoms in amyotrophic lateral sclerosis

OBJECTIVE

We investigated non-motor symptoms in ALS using sequential questionnaires; here we report the findings of the second questionnaire.

METHODS

A social media platform (Twitter, now known as X) was used to publicize the questionnaires. Data were downloaded from SurveyMonkey and analyzed by descriptive statistics, comparison of means, and regression models.

RESULTS

There were 182 people with ALS and 57 controls. The most important non-motor symptoms were cold limbs (60.4% cases, 14% controls, p = 9.67 x 10-10) and appetite loss (29.7% cases, 5.3% controls, p = 1.6 x 10-4). The weaker limb was most likely to feel cold (p = 9.67 x 10-10), and symptoms were more apparent in the evening and night. Appetite loss was reported as due to feeling full and the time taken to eat. People with ALS experienced medium-intensity pain, more usually shock-like pain than burning or cold-like pain, although the most prevalent type of pain was non-differentiated.

CONCLUSIONS

Non-motor symptoms are an important feature of ALS. Further investigation is needed to understand their physiological basis and whether they represent phenotypic differences useful for subtyping ALS.

Statistical inference on qualitative differences in the magnitude of an effect

Qualitative interactions occur when a treatment effect or measure of association varies in sign by sub-population. Of particular interest in many biomedical settings are absence/presence qualitative interactions, which occur when an effect is present in one sub-population but absent in another. Absence/presence interactions arise in emerging applications in precision medicine, where the objective is to identify a set of predictive biomarkers that have prognostic value for clinical outcomes in some sub-population but not others. They also arise naturally in gene regulatory network inference, where the goal is to identify differences in networks corresponding to diseased and healthy individuals, or to different subtypes of disease; such differences lead to identification of network-based biomarkers for diseases. In this paper, we argue that while the absence/presence hypothesis is important, developing a statistical test for this hypothesis is an intractable problem. To overcome this challenge, we approximate the problem in a novel inference framework. In particular, we propose to make inferences about absence/presence interactions by quantifying the relative difference in effect size, reasoning that when the relative difference is large, an absence/presence interaction occurs. The proposed methodology is illustrated through a simulation study as well as an analysis of breast cancer data from the Cancer Genome Atlas.

Plasma Metabolomics Reveals Distinct Biological and Diagnostic Signatures for Melioidosis

Rationale: The global burden of sepsis is greatest in low-resource settings. Melioidosis, infection with the gram-negative bacterium Burkholderia pseudomallei, is a frequent cause of fatal sepsis in endemic tropical regions such as Southeast Asia. Objectives: To investigate whether plasma metabolomics would identify biological pathways specific to melioidosis and yield clinically meaningful biomarkers. Methods: Using a comprehensive approach, differential enrichment of plasma metabolites and pathways was systematically evaluated in individuals selected from a prospective cohort of patients hospitalized in rural Thailand with infection. Statistical and bioinformatics methods were used to distinguish metabolomic features and processes specific to patients with melioidosis and between fatal and nonfatal cases. Measurements and Main Results: Metabolomic profiling and pathway enrichment analysis of plasma samples from patients with melioidosis (n = 175) and nonmelioidosis infections (n = 75) revealed a distinct immuno-metabolic state among patients with melioidosis, as suggested by excessive tryptophan catabolism in the kynurenine pathway and significantly increased levels of sphingomyelins and ceramide species. We derived a 12-metabolite classifier to distinguish melioidosis from other infections, yielding an area under the receiver operating characteristic curve of 0.87 in a second validation set of patients. Melioidosis nonsurvivors (n = 94) had a significantly disturbed metabolome compared with survivors (n = 81), with increased leucine, isoleucine, and valine metabolism, and elevated circulating free fatty acids and acylcarnitines. A limited eight-metabolite panel showed promise as an early prognosticator of mortality in melioidosis. Conclusions: Melioidosis induces a distinct metabolomic state that can be examined to distinguish underlying pathophysiological mechanisms associated with death. A 12-metabolite signature accurately differentiates melioidosis from other infections and may have diagnostic applications.

Non-motor symptoms in amyotrophic lateral sclerosis

OBJECTIVE

While motor symptoms are well-known in ALS, non-motor symptoms are often under-reported and may have a significant impact on quality of life. In this study, we aimed to examine the nature and extent of non-motor symptoms in ALS.

METHODS

A 20-item questionnaire was developed covering the domains of autonomic function, sleep, pain, gastrointestinal disturbance, and emotional lability, posted online and shared on social media platforms to target people with ALS and controls.

RESULTS

A total of 1018 responses were received, of which 927 were complete from 506 people with ALS and 421 unaffected individuals. Cold limbs (p 1.66 × 10-36), painful limbs (p 6.14 × 10-28), and urinary urgency (p 4.70 × 10-23) were associated with ALS. People with ALS were more likely to report autonomic symptoms, pain, and psychiatric symptoms than controls (autonomic symptoms B = 0.043, p 6.10 × 10-5, pain domain B = 0.18, p 3.72 × 10-11 and psychiatric domain B = 0.173, p 1.32 × 10-4).

CONCLUSIONS

Non-motor symptoms in ALS are common. The identification and management of non-motor symptoms should be integrated into routine clinical care for people with ALS. Further research is warranted to investigate the relationship between non-motor symptoms and disease progression, as well as to develop targeted interventions to improve the quality of life for people with ALS.

GAD65Abs Are Not Associated With Beta-Cell Dysfunction in Patients With T2D in the GRADE Study

Context

Autoantibodies directed against the 65-kilodalton isoform of glutamic acid decarboxylase (GAD65Abs) are markers of autoimmune type 1 diabetes (T1D) but are also present in patients with Latent Autoimmune Diabetes of Adults and autoimmune neuromuscular diseases, and also in healthy individuals. Phenotypic differences between these conditions are reflected in epitope-specific GAD65Abs and anti-idiotypic antibodies (anti-Id) against GAD65Abs. We previously reported that 7.8% of T2D patients in the GRADE study have GAD65Abs but found that GAD65Ab positivity was not correlated with beta-cell function, glycated hemoglobin (HbA1c), or fasting glucose levels.

Context

In this study, we aimed to better characterize islet autoantibodies in this T2D cohort. This is an ancillary study to NCT01794143.

Methods

We stringently defined GAD65Ab positivity with a competition assay, analyzed GAD65Ab-specific epitopes, and measured GAD65Ab-specific anti-Id in serum.

Results

Competition assays confirmed that 5.9% of the patients were GAD65Ab positive, but beta-cell function was not associated with GAD65Ab positivity, GAD65Ab epitope specificity or GAD65Ab-specific anti-Id. GAD65-related autoantibody responses in GRADE T2D patients resemble profiles in healthy individuals (low GAD65Ab titers, presence of a single autoantibody, lack of a distinct epitope pattern, and presence of anti-Id to diabetes-associated GAD65Ab). In this T2D cohort, GAD65Ab positivity is likely unrelated to the pathogenesis of beta-cell dysfunction.

Conclusion

Evidence for islet autoimmunity in the pathophysiology of T2D beta-cell dysfunction is growing, but T1D-associated autoantibodies may not accurately reflect the nature of their autoimmune process.

GENERALIZED MATRIX DECOMPOSITION REGRESSION: ESTIMATION AND INFERENCE FOR TWO-WAY STRUCTURED DATA

Motivated by emerging applications in ecology, microbiology, and neuroscience, this paper studies high-dimensional regression with two-way structured data. To estimate the high-dimensional coefficient vector, we propose the generalized matrix decomposition regression (GMDR) to efficiently leverage auxiliary information on row and column structures. GMDR extends the principal component regression (PCR) to two-way structured data, but unlike PCR, GMDR selects the components that are most predictive of the outcome, leading to more accurate prediction. For inference on regression coefficients of individual variables, we propose the generalized matrix decomposition inference (GMDI), a general high-dimensional inferential framework for a large family of estimators that include the proposed GMDR estimator. GMDI provides more flexibility for incorporating relevant auxiliary row and column structures. As a result, GMDI does not require the true regression coefficients to be sparse, but constrains the coordinate system representing the regression coefficients according to the column structure. GMDI also allows dependent and heteroscedastic observations. We study the theoretical properties of GMDI in terms of both the type-I error rate and power and demonstrate the effectiveness of GMDR and GMDI in simulation studies and an application to human microbiome data.

Identification of circulating proteins associated with general cognitive function among middle-aged and older adults

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p < 4.9E-5, n up to 7289). Of these, 45 were replicated using SomaScan data, and three were replicated using Olink data at Bonferroni-corrected significance. Enrichment analysis linked the proteins associated with general cognitive function to cell signaling pathways and synapse architecture. Mendelian randomization analysis implicated higher levels of NECTIN2, a protein mediating viral entry into neuronal cells, with higher Alzheimer's disease (AD) risk (p = 2.5E-26). Levels of 14 other protein measures were implicated as consequences of AD susceptibility (p < 2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets.

Determinants of mosaic chromosomal alteration fitness

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R 2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified TCL1A , NRIP1 , and TERT locus variants as modulators of mCA clonal expansion rate.

The functional impact of rare variation across the regulatory cascade

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Human microglia show unique transcriptional changes in Alzheimer's disease

Microglia, the innate immune cells of the brain, influence Alzheimer's disease (AD) progression and are potential therapeutic targets. However, microglia exhibit diverse functions, the regulation of which is not fully understood, complicating therapeutics development. To better define the transcriptomic phenotypes and gene regulatory networks associated with AD, we enriched for microglia nuclei from 12 AD and 10 control human dorsolateral prefrontal cortices (7 males and 15 females, all aged >60 years) before single-nucleus RNA sequencing. Here we describe both established and previously unrecognized microglial molecular phenotypes, the inferred gene networks driving observed transcriptomic change, and apply trajectory analysis to reveal the putative relationships between microglial phenotypes. We identify microglial phenotypes more prevalent in AD cases compared with controls. Further, we describe the heterogeneity in microglia subclusters expressing homeostatic markers. Our study demonstrates that deep profiling of microglia in human AD brain can provide insight into microglial transcriptional changes associated with AD.

Non-motor symptoms in amyotrophic lateral sclerosis: lessons from Parkinson's disease

Amyotrophic lateral sclerosis and Parkinson's disease are neurodegenerative diseases of the motor system which are now recognized also to affect non-motor pathways. Non-motor symptoms have been acknowledged as important determinants of quality of life in Parkinson's disease, and there is increasing interest in understanding the extent and role of non-motor symptoms in amyotrophic lateral sclerosis. We therefore reviewed what is known about non-motor symptoms in amyotrophic lateral sclerosis, using lessons from Parkinson's disease.

Plasma Proteomic Associations With Incident Ischemic Stroke in Older Adults: The Cardiovascular Health Study

BACKGROUND

Plasma proteomics may elucidate novel insights into the pathophysiology of ischemic stroke (IS), identify biomarkers of IS risk, and guide development of nascent prevention strategies. We evaluated the relationship between the plasma proteome and IS risk in the population-based Cardiovascular Health Study (CHS).

METHODS

Eligible CHS participants were free of prevalent stroke and underwent quantification of 1298 plasma proteins using the aptamer-based SOMAScan assay platform from the 1992-1993 study visit. Multivariable Cox proportional hazards regression was used to evaluate associations between a 1-standard deviation increase in the log-2 transformed estimated plasma protein concentrations and incident IS, adjusting for demographics, IS risk factors, and estimated glomerular filtration rate. For proteins independently associated with incident IS, a secondary stratified analysis evaluated associations in subgroups defined by sex and race. Exploratory analyses evaluated plasma proteomic associations with cardioembolic and non-cardioembolic IS as well as proteins associated with IS risk in participants with left atrial dysfunction but without atrial fibrillation.

RESULTS

Of 2983 eligible participants, the mean age was 74.3 (± 4.8) years, 61.2% were women, and 15.4% were Black. Over a median follow-up of 12.6 years, 450 participants experienced an incident IS. N-terminal pro-brain natriuretic peptide (NTproBNP, adjusted HR 1.37, 95% CI 1.23-1.53, P=2.08x10^-08) and macrophage metalloelastase (MMP12, adjusted HR 1.30, 95% CI 1.16-1.45, P=4.55x10^-06) were independently associated with IS risk. These two associations were similar in men and women and in Black and non-Black participants. In exploratory analyses, NTproBNP was independently associated with incident cardioembolic IS, E-selectin with incident non-cardioembolic IS, and secreted frizzled-related protein 1 with IS risk in participants with left atrial dysfunction.

CONCLUSIONS

In a cohort of older adults, NTproBNP and MMP12 were independently associated with IS risk. We identified plasma proteomic determinants of incident cardioembolic and non-cardioembolic IS and found a novel protein associated with IS risk in those with left atrial dysfunction.

Alteration of oral microbiome composition in children living with pesticide-exposed farm workers

Our prior work shows that azinphos-methyl pesticide exposure is associated with altered oral microbiomes in exposed farmworkers. Here we extend this analysis to show the same association pattern is also evident in their children. Oral buccal swab samples were analyzed at two time points, the apple thinning season in spring-summer 2005 for 78 children and 101 adults and the non-spray season in winter 2006 for 62 children and 82 adults. The pesticide exposure for the children were defined by the farmworker occupation of the cohabitating household adult and the blood azinphos-methyl detection of the cohabitating adult. Oral buccal swab 16S rRNA sequencing determined taxonomic microbiota proportional composition from concurrent samples from both adults and children. Analysis of the identified bacteria showed significant proportional changes for 12 of 23 common oral microbiome genera in association with azinphos-methyl detection and farmworker occupation. The most common significantly altered genera had reductions in the abundance of Streptococcus, suggesting an anti-microbial effect of the pesticide. Principal component analysis of the microbiome identified two primary clusters, with association of principal component 1 to azinphos-methyl blood detection and farmworker occupational status of the household. The children's buccal microbiota composition clustered with their household adult in ∼95% of the households. Household adult farmworker occupation and household pesticide exposure is associated with significant alterations in their children's oral microbiome composition. This suggests that parental occupational exposure and pesticide take-home exposure pathways elicit alteration of their children's microbiomes.

Metabolomic Signatures Associated With Pulmonary Arterial Hypertension Outcomes

BACKGROUND

Pulmonary arterial hypertension (PAH) is a complex disease characterized by progressive right ventricular (RV) failure leading to significant morbidity and mortality. Investigating metabolic features and pathways associated with RV dilation, mortality, and measures of disease severity can provide insight into molecular mechanisms, identify subphenotypes, and suggest potential therapeutic targets.

METHODS

We collected data from a prospective cohort of PAH participants and performed untargeted metabolomic profiling on 1045 metabolites from circulating blood. Analyses were intended to identify metabolomic differences across a range of common metrics in PAH (eg, dilated versus nondilated RV). Partial least squares discriminant analysis was first applied to assess the distinguishability of relevant outcomes. Significantly altered metabolites were then identified using linear regression, and Cox regression models (as appropriate for the specific outcome) with adjustments for age, sex, body mass index, and PAH cause. Models exploring RV maladaptation were further adjusted for pulmonary vascular resistance. Pathway enrichment analysis was performed to identify significantly dysregulated processes.

RESULTS

A total of 117 participants with PAH were included. Partial least squares discriminant analysis showed cluster differentiation between participants with dilated versus nondilated RVs, survivors versus nonsurvivors, and across a range of NT-proBNP (N-terminal pro-B-type natriuretic peptide) levels, REVEAL 2.0 composite scores, and 6-minute-walk distances. Polyamine and histidine pathways were associated with differences in RV dilation, mortality, NT-proBNP, REVEAL score, and 6-minute walk distance. Acylcarnitine pathways were associated with NT-proBNP, REVEAL score, and 6-minute walk distance. Sphingomyelin pathways were associated with RV dilation and NT-proBNP after adjustment for pulmonary vascular resistance.

CONCLUSIONS

Distinct plasma metabolomic profiles are associated with RV dilation, mortality, and measures of disease severity in PAH. Polyamine, histidine, and sphingomyelin metabolic pathways represent promising candidates for identifying patients at high risk for poor outcomes and investigation into their roles as markers or mediators of disease progression and RV adaptation.

NETWORK DIFFERENTIAL CONNECTIVITY ANALYSIS

Identifying differences in networks has become a canonical problem in many biological applications. Existing methods try to accomplish this goal by either directly comparing the estimated structures of two networks, or testing the null hypothesis that the covariance or inverse covariance matrices in two populations are identical. However, estimation approaches do not provide measures of uncertainty, e.g., p-values, whereas existing testing approaches could lead to misleading results, as we illustrate in this paper. To address these shortcomings, we propose a qualitative hypothesis testing framework, which tests whether the connectivity structures in the two networks are the same. our framework is especially appropriate if the goal is to identify nodes or edges that are differentially connected. No existing approach could test such hypotheses and provide corresponding measures of uncertainty. Theoretically, we show that under appropriate conditions, our proposal correctly controls the type-I error rate in testing the qualitative hypothesis. Empirically, we demonstrate the performance of our proposal using simulation studies and applications in cancer genomics.

Neural Granger Causality

While most classical approaches to Granger causality detection assume linear dynamics, many interactions in real-world applications, like neuroscience and genomics, are inherently nonlinear. In these cases, using linear models may lead to inconsistent estimation of Granger causal interactions. We propose a class of nonlinear methods by applying structured multilayer perceptrons (MLPs) or recurrent neural networks (RNNs) combined with sparsity-inducing penalties on the weights. By encouraging specific sets of weights to be zero-in particular, through the use of convex group-lasso penalties-we can extract the Granger causal structure. To further contrast with traditional approaches, our framework naturally enables us to efficiently capture long-range dependencies between series either via our RNNs or through an automatic lag selection in the MLP. We show that our neural Granger causality methods outperform state-of-the-art nonlinear Granger causality methods on the DREAM3 challenge data. This data consists of nonlinear gene expression and regulation time courses with only a limited number of time points. The successes we show in this challenging dataset provide a powerful example of how deep learning can be useful in cases that go beyond prediction on large datasets. We likewise illustrate our methods in detecting nonlinear interactions in a human motion capture dataset.

Proteomics and Population Biology in the Cardiovascular Health Study (CHS): design of a study with mentored access and active data sharing

BACKGROUND

In the last decade, genomic studies have identified and replicated thousands of genetic associations with measures of health and disease and contributed to the understanding of the etiology of a variety of health conditions. Proteins are key biomarkers in clinical medicine and often drug-therapy targets. Like genomics, proteomics can advance our understanding of biology.

METHODS AND RESULTS

In the setting of the Cardiovascular Health Study (CHS), a cohort study of older adults, an aptamer-based method that has high sensitivity for low-abundance proteins was used to assay 4979 proteins in frozen, stored plasma from 3188 participants (61% women, mean age 74 years). CHS provides active support, including central analysis, for seven phenotype-specific working groups (WGs). Each CHS WG is led by one or two senior investigators and includes 10 to 20 early or mid-career scientists. In this setting of mentored access, the proteomic data and analytic methods are widely shared with the WGs and investigators so that they may evaluate associations between baseline levels of circulating proteins and the incidence of a variety of health outcomes in prospective cohort analyses. We describe the design of CHS, the CHS Proteomics Study, characteristics of participants, quality control measures, and structural characteristics of the data provided to CHS WGs. We additionally highlight plans for validation and replication of novel proteomic associations.

CONCLUSION

The CHS Proteomics Study offers an opportunity for collaborative data sharing to improve our understanding of the etiology of a variety of health conditions in older adults.

Generalized score matching for general domains

Estimation of density functions supported on general domains arises when the data are naturally restricted to a proper subset of the real space. This problem is complicated by typically intractable normalizing constants. Score matching provides a powerful tool for estimating densities with such intractable normalizing constants but as originally proposed is limited to densities on [Formula: see text] and [Formula: see text]. In this paper, we offer a natural generalization of score matching that accommodates densities supported on a very general class of domains. We apply the framework to truncated graphical and pairwise interaction models and provide theoretical guarantees for the resulting estimators. We also generalize a recently proposed method from bounded to unbounded domains and empirically demonstrate the advantages of our method.

Network structure mediates functional reorganization induced by optogenetic stimulation of non-human primate sensorimotor cortex

Because aberrant network-level functional connectivity underlies a variety of neural disorders, the ability to induce targeted functional reorganization would be a profound development toward therapies for neural disorders. Brain stimulation has been shown to induce large-scale network-wide functional connectivity changes (FCC), but the mapping from stimulation to the induced changes is unclear. Here, we develop a model which jointly considers the stimulation protocol and the cortical network structure to accurately predict network-wide FCC in response to optogenetic stimulation of non-human primate primary sensorimotor cortex. We observe that the network structure has a much stronger effect than the stimulation protocol on the resulting FCC. We also observe that the mappings from these input features to the FCC diverge over frequency bands and successive stimulations. Our framework represents a paradigm shift for targeted neural stimulation and can be used to interrogate, improve, and develop stimulation-based interventions for neural disorders.

•

Optogenetic stimulation drives connectivity changes over sensorimotor cortex

•

Nonparametric model informed by protocol and network features predicts the changes

•

The underlying network is the primary mediator of the induced changes

•

The mappings governing the changes diverge over time and frequency bands

Granger Causality: A Review and Recent Advances

Introduced more than a half-century ago, Granger causality has become a popular tool for analyzing time series data in many application domains, from economics and finance to genomics and neuroscience. Despite this popularity, the validity of this framework for inferring causal relationships among time series has remained the topic of continuous debate. Moreover, while the original definition was general, limitations in computational tools have constrained the applications of Granger causality to primarily simple bivariate vector autoregressive processes. Starting with a review of early developments and debates, this article discusses recent advances that address various shortcomings of the earlier approaches, from models for high-dimensional time series to more recent developments that account for nonlinear and non-Gaussian observations and allow for subsampled and mixed-frequency time series.

Nonparametric Causal Structure Learning in High Dimensions

The PC and FCI algorithms are popular constraint-based methods for learning the structure of directed acyclic graphs (DAGs) in the absence and presence of latent and selection variables, respectively. These algorithms (and their order-independent variants, PC-stable and FCI-stable) have been shown to be consistent for learning sparse high-dimensional DAGs based on partial correlations. However, inferring conditional independences from partial correlations is valid if the data are jointly Gaussian or generated from a linear structural equation model—an assumption that may be violated in many applications. To broaden the scope of high-dimensional causal structure learning, we propose nonparametric variants of the PC-stable and FCI-stable algorithms that employ the conditional distance covariance (CdCov) to test for conditional independence relationships. As the key theoretical contribution, we prove that the high-dimensional consistency of the PC-stable and FCI-stable algorithms carry over to general distributions over DAGs when we implement CdCov-based nonparametric tests for conditional independence. Numerical studies demonstrate that our proposed algorithms perform nearly as good as the PC-stable and FCI-stable for Gaussian distributions, and offer advantages in non-Gaussian graphical models.

Islet Autoimmunity is Highly Prevalent and Associated With Diminished β-Cell Function in Patients With Type 2 Diabetes in the Grade Study

Islet autoimmunity may contribute to β-cell dysfunction in type 2 diabetes (T2D). Its prevalence and clinical significance have not been rigorously determined. In this ancillary study to the Glycemia Reduction Approaches in Diabetes-A Comparative Effectiveness (GRADE) Study, we investigated the prevalence of cellular and humoral islet autoimmunity in patients with T2D duration 4·0±3·0 y, HbA1c 7·5±0·5% on metformin alone. We measured T cell autoreactivity against islet proteins, islet autoantibodies against GAD65, IA2, ZnT8, and β-cell function. Cellular islet autoimmunity was present in 41·3%, humoral islet autoimmunity in 13·5%, and both in 5·3%. β-cell function calculated as iAUC-CG and ΔC-peptide(0- 30)/Δglucose(0-30) from an oral glucose tolerance test was lower among T cell-positives (T+) than T cell-negatives (T-) using two different adjustments for insulin sensitivity (iAUC-CG: 13·2% [95% CI 0·3, 24·4%] or 11·4% [95% CI 0·4, 21·2%] lower; ΔC-peptide(0-30)/Δglucose(0-30)) 19% [95% CI 3·1, 32·3%] or 17·7% [95% CI 2·6, 30·5%] lower). T+ patients had 17% higher HbA1c (95% CI 0·07, 0·28) and 7·7 mg/dL higher fasting plasma glucose levels (95% CI 0·2,15·3) than T- patients. We conclude that islet autoimmunity is much more prevalent in T2D patients than previously reported. T cell-mediated autoimmunity is associated with diminished β-cell function and worse glycemic control.

Penalized estimation of threshold auto-regressive models with many components and thresholds

Thanks to their simplicity and interpretable structure, autoregressive processes are widely used to model time series data. However, many real time series data sets exhibit non-linear patterns, requiring nonlinear modeling. The threshold Auto-Regressive (TAR) process provides a family of non-linear auto-regressive time series models in which the process dynamics are specific step functions of a thresholding variable. While estimation and inference for low-dimensional TAR models have been investigated, high-dimensional TAR models have received less attention. In this article, we develop a new framework for estimating high-dimensional TAR models, and propose two different sparsity-inducing penalties. The first penalty corresponds to a natural extension of classical TAR model to high-dimensional settings, where the same threshold is enforced for all model parameters. Our second penalty develops a more flexible TAR model, where different thresholds are allowed for different auto-regressive coefficients. We show that both penalized estimation strategies can be utilized in a three-step procedure that consistently learns both the thresholds and the corresponding auto-regressive coefficients. However, our theoretical and empirical investigations show that the direct extension of the TAR model is not appropriate for high-dimensional settings and is better suited for moderate dimensions. In contrast, the more flexible extension of the TAR model leads to consistent estimation and superior empirical performance in high dimensions.

Generalized Sparse Additive Models

We present a unified framework for estimation and analysis of generalized additive models in high dimensions. The framework defines a large class of penalized regression estimators, encompassing many existing methods. An efficient computational algorithm for this class is presented that easily scales to thousands of observations and features. We prove minimax optimal convergence bounds for this class under a weak compatibility condition. In addition, we characterize the rate of convergence when this compatibility condition is not met. Finally, we also show that the optimal penalty parameters for structure and sparsity penalties in our framework are linked, allowing cross-validation to be conducted over only a single tuning parameter. We complement our theoretical results with empirical studies comparing some existing methods within this framework.

Microglia subtype transcriptomes differ between Alzheimer Disease and control human postmortem brain samples

BACKGROUND

Microglia-mediated neuroinflammation is hypothesized to contribute to disease progression in neurodegenerative diseases such as Alzheimer's Disease (AD). Microglia subtypes are complex, with beneficial and harmful phenotypes. Understanding the gene expression networks which define the spectrum of microglia phenotypes is critical to identifying specific targets for neuroinflammation modulating therapies.

METHOD

Our study utilized post-mortem brain tissue from 22 total (7 male) participants; 12 (3 male) had significant AD neuropathic change. Nuclei isolated from prefrontal cortex were sorted for the myeloid marker PU.1 using fluorescence activated nucleus sorting (FANS). The FANS approach yields larger numbers of nuclei annotated as microglia with high quality sequence from each individual. We performed single-nucleus RNA-seq using the 10X Genomics Chromium platform.

RESULTS

We isolated more than 120,000 microglia nuclei, facilitating group comparisons based on disease state. Unbiased clustering revealed 10 microglia clusters and improved resolution of microglia heterogeneity compared to standard single-cell approaches. We identify clusters of microglia enriched for biological pathways implicating defined myeloid roles including interferon-stimulated, endo/lysosomal, neurodegenerative with a "disease-associated microglia" (DAM) signature, as well as a metabolically active and autophagic cluster. Interestingly, the cluster proportionately enriched for AD individuals' nuclei is not the DAM cluster but instead one of the clusters in which endo/lysosomal genes are highly upregulated. Furthermore, many of the genes in known AD risk loci are strongly differentially regulated in this AD associated cluster. We also identify a cluster of microglia that is proportionately enriched for control samples with upregulated cell cycle and proliferation genes. Trajectory analysis suggests that the paths AD and control nuclei take from unactivated "homeostatic" to various phenotypic states are also distinct.

CONCLUSION

Using human AD tissue collected with uniform protocols we characterize the transcriptomic profiles of microglia subtypes in human brain. By enriching for myeloid cells prior to analysis we can resolve microglia subtypes revealing the diversity of microglia which are "inflammatory" as well as other microglia subtypes responding with induction of metabolic and lysosomal pathways. Our data identifies subtypes of microglia that are unique to AD and control individuals. These results support the possibility of pharmacological targeting of specific subtypes of microglia to alter AD progression.

Causal Discovery in High-Dimensional Point Process Networks with Hidden Nodes

Thanks to technological advances leading to near-continuous time observations, emerging multivariate point process data offer new opportunities for causal discovery. However, a key obstacle in achieving this goal is that many relevant processes may not be observed in practice. Naïve estimation approaches that ignore these hidden variables can generate misleading results because of the unadjusted confounding. To plug this gap, we propose a deconfounding procedure to estimate high-dimensional point process networks with only a subset of the nodes being observed. Our method allows flexible connections between the observed and unobserved processes. It also allows the number of unobserved processes to be unknown and potentially larger than the number of observed nodes. Theoretical analyses and numerical studies highlight the advantages of the proposed method in identifying causal interactions among the observed processes.

In Defense of the Indefensible: A Very Naïve Approach to High-Dimensional Inference

A great deal of interest has recently focused on conducting inference on the parameters in a high-dimensional linear model. In this paper, we consider a simple and very naïve two-step procedure for this task, in which we (i) fit a lasso model in order to obtain a subset of the variables, and (ii) fit a least squares model on the lasso-selected set. Conventional statistical wisdom tells us that we cannot make use of the standard statistical inference tools for the resulting least squares model (such as confidence intervals and p-values), since we peeked at the data twice: once in running the lasso, and again in fitting the least squares model. However, in this paper, we show that under a certain set of assumptions, with high probability, the set of variables selected by the lasso is identical to the one selected by the noiseless lasso and is hence deterministic. Consequently, the naïve two-step approach can yield asymptotically valid inference. We utilize this finding to develop the naïve confidence interval, which can be used to draw inference on the regression coefficients of the model selected by the lasso, as well as the naïve score test, which can be used to test the hypotheses regarding the full-model regression coefficients.

REHE: Fast variance components estimation for linear mixed models

Linear mixed models are widely used in ecological and biological applications, especially in genetic studies. Reliable estimation of variance components is crucial for using linear mixed models. However, standard methods, such as the restricted maximum likelihood (REML), are computationally inefficient in large samples and may be unstable with small samples. Other commonly used methods, such as the Haseman-Elston (HE) regression, may yield negative estimates of variances. Utilizing regularized estimation strategies, we propose the restricted Haseman-Elston (REHE) regression and REHE with resampling (reREHE) estimators, along with an inference framework for REHE, as fast and robust alternatives that provide nonnegative estimates with comparable accuracy to REML. The merits of REHE are illustrated using real data and benchmark simulation studies.

CorDiffViz: an R package for visualizing multi-omics differential correlation networks

BACKGROUND

Differential correlation networks are increasingly used to delineate changes in interactions among biomolecules. They characterize differences between omics networks under two different conditions, and can be used to delineate mechanisms of disease initiation and progression.

RESULTS

We present a new R package, CorDiffViz, that facilitates the estimation and visualization of differential correlation networks using multiple correlation measures and inference methods. The software is implemented in R, HTML and Javascript, and is available at https://github.com/sqyu/CorDiffViz . Visualization has been tested for the Chrome and Firefox web browsers. A demo is available at https://diffcornet.github.io/CorDiffViz/demo.html .

CONCLUSIONS

Our software offers considerable flexibility by allowing the user to interact with the visualization and choose from different estimation methods and visualizations. It also allows the user to easily toggle between correlation networks for samples under one condition and differential correlations between samples under two conditions. Moreover, the software facilitates integrative analysis of cross-correlation networks between two omics data sets.

netgsa: Fast computation and interactive visualization for topology-based pathway enrichment analysis

Existing software tools for topology-based pathway enrichment analysis are either computationally inefficient, have undesirable statistical power, or require expert knowledge to leverage the methods' capabilities. To address these limitations, we have overhauled NetGSA, an existing topology-based method, to provide a computationally-efficient user-friendly tool that offers interactive visualization. Pathway enrichment analysis for thousands of genes can be performed in minutes on a personal computer without sacrificing statistical power. The new software also removes the need for expert knowledge by directly curating gene-gene interaction information from multiple external databases. Lastly, by utilizing the capabilities of Cytoscape, the new software also offers interactive and intuitive network visualization.

Differential Network Analysis: A Statistical Perspective

Networks effectively capture interactions among components of complex systems, and have thus become a mainstay in many scientific disciplines. Growing evidence, especially from biology, suggest that networks undergo changes over time, and in response to external stimuli. In biology and medicine, these changes have been found to be predictive of complex diseases. They have also been used to gain insight into mechanisms of disease initiation and progression. Primarily motivated by biological applications, this article provides a review of recent statistical machine learning methods for inferring networks and identifying changes in their structures.

Integer Programming for Learning Directed Acyclic Graphs from Continuous Data

Learning directed acyclic graphs (DAGs) from data is a challenging task both in theory and in practice, because the number of possible DAGs scales superexponentially with the number of nodes. In this paper, we study the problem of learning an optimal DAG from continuous observational data. We cast this problem in the form of a mathematical programming model that can naturally incorporate a superstructure to reduce the set of possible candidate DAGs. We use a negative log-likelihood score function with both [Formula: see text] and [Formula: see text] penalties and propose a new mixed-integer quadratic program, referred to as a layered network (LN) formulation. The LN formulation is a compact model that enjoys as tight an optimal continuous relaxation value as the stronger but larger formulations under a mild condition. Computational results indicate that the proposed formulation outperforms existing mathematical formulations and scales better than available algorithms that can solve the same problem with only [Formula: see text] regularization. In particular, the LN formulation clearly outperforms existing methods in terms of computational time needed to find an optimal DAG in the presence of a sparse superstructure.

The Convex Mixture Distribution: Granger Causality for Categorical Time Series

We present a framework for learning Granger causality networks for multivariate categorical time series based on the mixture transition distribution (MTD) model. Traditionally, MTD is plagued by a nonconvex objective, non-identifiability, and presence of local optima. To circumvent these problems, we recast inference in the MTD as a convex problem. The new formulation facilitates the application of MTD to high-dimensional multivariate time series. As a baseline, we also formulate a multi-output logistic autoregressive model (mLTD), which while a straightforward extension of autoregressive Bernoulli generalized linear models, has not been previously applied to the analysis of multivariate categorial time series. We establish identifiability conditions of the MTD model and compare them to those for mLTD. We further devise novel and efficient optimization algorithms for MTD based on our proposed convex formulation, and compare the MTD and mLTD in both simulated and real data experiments. Finally, we establish consistency of the convex MTD in high dimensions. Our approach simultaneously provides a comparison of methods for network inference in categorical time series and opens the door to modern, regularized inference with the MTD model.

Plasma lipidomic profiles after a low and high glycemic load dietary pattern in a randomized controlled crossover feeding study

BACKGROUND

Dietary patterns low in glycemic load are associated with reduced risk of cardiometabolic diseases. Improvements in serum lipid concentrations may play a role in these observed associations.

OBJECTIVE

We investigated how dietary patterns differing in glycemic load affect clinical lipid panel measures and plasma lipidomics profiles.

METHODS

In a crossover, controlled feeding study, 80 healthy participants (n = 40 men, n = 40 women), 18-45 y were randomized to receive low-glycemic load (LGL) or high glycemic load (HGL) diets for 28 days each with at least a 28-day washout period between controlled diets. Fasting plasma samples were collected at baseline and end of each diet period. Lipids on a clinical panel including total-, VLDL-, LDL-, and HDL-cholesterol and triglycerides were measured using an auto-analyzer. Lipidomics analysis using mass-spectrometry provided the concentrations of 863 species. Linear mixed models and lipid ontology enrichment analysis were implemented.

RESULTS

Lipids from the clinical panel were not significantly different between diets. Univariate analysis showed that 67 species on the lipidomics panel, predominantly in the triacylglycerol class, were higher after the LGL diet compared to the HGL (FDR < 0.05). Three species with FA 17:0 were lower after LGL diet with enrichment analysis (FDR < 0.05).

CONCLUSION

In the context of controlled eucaloric diets with similar macronutrient distribution, these results suggest that there are relative shifts in lipid species, but the overall pool does not change. Further studies are needed to better understand in which compartment the different lipid species are transported in blood, and how these shifts are related to health outcomes. This trial was registered at clinicaltrials.gov as NCT00622661.

Statistical significance in high-dimensional linear mixed models

This paper concerns the development of an inferential framework for high-dimensional linear mixed effect models. These are suitable models, for instance, when we have n repeated measurements for M subjects. We consider a scenario where the number of fixed effects p is large (and may be larger than M), but the number of random effects q is small. Our framework is inspired by a recent line of work that proposes de-biasing penalized estimators to perform inference for high-dimensional linear models with fixed effects only. In particular, we demonstrate how to correct a 'naive' ridge estimator in extension of work by Bühlmann (2013) to build asymptotically valid confidence intervals for mixed effect models. We validate our theoretical results with numerical experiments, in which we show our method outperforms those that fail to account for correlation induced by the random effects. For a practical demonstration we consider a riboflavin production dataset that exhibits group structure, and show that conclusions drawn using our method are consistent with those obtained on a similar dataset without group structure.

Genetic and metabolomic architecture of variation in diet restriction-mediated lifespan extension in Drosophila

In most organisms, dietary restriction (DR) increases lifespan. However, several studies have found that genotypes within the same species vary widely in how they respond to DR. To explore the mechanisms underlying this variation, we exposed 178 inbred Drosophila melanogaster lines to a DR or ad libitum (AL) diet, and measured a panel of 105 metabolites under both diets. Twenty four out of 105 metabolites were associated with the magnitude of the lifespan response. These included proteinogenic amino acids and metabolites involved in α-ketoglutarate (α-KG)/glutamine metabolism. We confirm the role of α-KG/glutamine synthesis pathways in the DR response through genetic manipulations. We used covariance network analysis to investigate diet-dependent interactions between metabolites, identifying the essential amino acids threonine and arginine as “hub” metabolites in the DR response. Finally, we employ a novel metabolic and genetic bipartite network analysis to reveal multiple genes that influence DR lifespan response, some of which have not previously been implicated in DR regulation. One of these is CCHa2R, a gene that encodes a neuropeptide receptor that influences satiety response and insulin signaling. Across the lines, variation in an intronic single nucleotide variant of CCHa2R correlated with variation in levels of five metabolites, all of which in turn were correlated with DR lifespan response. Inhibition of adult CCHa2R expression extended DR lifespan of flies, confirming the role of CCHa2R in lifespan response. These results provide support for the power of combined genomic and metabolomic analysis to identify key pathways underlying variation in this complex quantitative trait.

Dietary restriction extends lifespan across most organisms in which it has been tested. However, several studies have now demonstrated that this effect can vary dramatically across different genotypes within a population. Within a population, dietary restriction might be beneficial for some, yet detrimental for others. Here, we measure the metabolome of 178 genetically characterized fly strains on fully fed and restricted diets. The fly strains vary widely in their lifespan response to dietary restriction. We then use information about each strain’s genome and metabolome (a measure of small molecules circulating in flies) to pinpoint cellular pathways that govern this variation in response. We identify a novel pathway involving the gene CCHa2R, which encodes a neuropeptide receptor that has not previously been implicated in dietary restriction or age-related signaling pathways. This study demonstrates the power of leveraging systems biology and network biology methods to understand how and why different individuals vary in their response to health and lifespan-extending interventions.

Joint Structural Break Detection and Parameter Estimation in High-Dimensional Non-Stationary VAR Models

Assuming stationarity is unrealistic in many time series applications. A more realistic alternative is to assume piecewise stationarity, where the model can change at potentially many change points. We propose a three-stage procedure for simultaneous estimation of change points and parameters of high-dimensional piecewise vector autoregressive (VAR) models. In the first step, we reformulate the change point detection problem as a high-dimensional variable selection one, and solve it using a penalized least square estimator with a total variation penalty. We show that the penalized estimation method over-estimates the number of change points, and propose a selection criterion to identify the change points. In the last step of our procedure, we estimate the VAR parameters in each of the segments. We prove that the proposed procedure consistently detects the number and location of change points, and provides consistent estimates of VAR parameters. The performance of the method is illustrated through several simulated and real data examples.

The Generalized Matrix Decomposition Biplot and Its Application to Microbiome Data

Exploratory analysis of human microbiome data is often based on dimension-reduced graphical displays derived from similarities based on non-Euclidean distances, such as UniFrac or Bray-Curtis. However, a display of this type, often referred to as the principal-coordinate analysis (PCoA) plot, does not reveal which taxa are related to the observed clustering because the configuration of samples is not based on a coordinate system in which both the samples and variables can be represented. The reason is that the PCoA plot is based on the eigen-decomposition of a similarity matrix and not the singular value decomposition (SVD) of the sample-by-abundance matrix. We propose a novel biplot that is based on an extension of the SVD, called the generalized matrix decomposition biplot (GMD-biplot), which involves an arbitrary matrix of similarities and the original matrix of variable measures, such as taxon abundances. As in a traditional biplot, points represent the samples, and arrows represent the variables. The proposed GMD-biplot is illustrated by analyzing multiple real and simulated data sets which demonstrate that the GMD-biplot provides improved clustering capability and a more meaningful relationship between the arrows and points.IMPORTANCE Biplots that simultaneously display the sample clustering and the important taxa have gained popularity in the exploratory analysis of human microbiome data. Traditional biplots, assuming Euclidean distances between samples, are not appropriate for microbiome data, when non-Euclidean distances are used to characterize dissimilarities among microbial communities. Thus, incorporating information from non-Euclidean distances into a biplot becomes useful for graphical displays of microbiome data. The proposed GMD-biplot accounts for any arbitrary non-Euclidean distances and provides a robust and computationally efficient approach for graphical visualization of microbiome data. In addition, the proposed GMD-biplot displays both the samples and taxa with respect to the same coordinate system, which further allows the configuration of future samples.

A comparative study of topology-based pathway enrichment analysis methods

BACKGROUND

Pathway enrichment extensively used in the analysis of Omics data for gaining biological insights into the functional roles of pre-defined subsets of genes, proteins and metabolites. A large number of methods have been proposed in the literature for this task. The vast majority of these methods use as input expression levels of the biomolecules under study together with their membership in pathways of interest. The latest generation of pathway enrichment methods also leverages information on the topology of the underlying pathways, which as evidence from their evaluation reveals, lead to improved sensitivity and specificity. Nevertheless, a systematic empirical comparison of such methods is still lacking, making selection of the most suitable method for a specific experimental setting challenging. This comparative study of nine network-based methods for pathway enrichment analysis aims to provide a systematic evaluation of their performance based on three real data sets with different number of features (genes/metabolites) and number of samples.

RESULTS

The findings highlight both methodological and empirical differences across the nine methods. In particular, certain methods assess pathway enrichment due to differences both across expression levels and in the strength of the interconnectedness of the members of the pathway, while others only leverage differential expression levels. In the more challenging setting involving a metabolomics data set, the results show that methods that utilize both pieces of information (with NetGSA being a prototypical one) exhibit superior statistical power in detecting pathway enrichment.

CONCLUSION

The analysis reveals that a number of methods perform equally well when testing large size pathways, which is the case with genomic data. On the other hand, NetGSA that takes into consideration both differential expression of the biomolecules in the pathway, as well as changes in the topology exhibits a superior performance when testing small size pathways, which is usually the case for metabolomics data.

Bayesian Hidden Markov Models for Dependent Large-Scale Multiple Testing

An optimal and flexible multiple hypotheses testing procedure is constructed for dependent data based on Bayesian techniques, aiming at handling two challenges, namely dependence structure and non-null distribution specification. Ignoring dependence among hypotheses tests may lead to loss of efficiency and bias in decision. Misspecification in the non-null distribution, on the other hand, can result in both false positive and false negative errors. Hidden Markov models are used to accommodate the dependence structure among the tests. Dirichlet mixture process prior is applied on the non-null distribution to overcome the potential pitfalls in distribution misspecification. The testing algorithm based on Bayesian techniques optimizes the false negative rate (FNR) while controlling the false discovery rate (FDR). The procedure is applied to pointwise and clusterwise analysis. Its performance is compared with existing approaches using both simulated and real data examples.

Plasma metabolomics profiles suggest beneficial effects of a low-glycemic load dietary pattern on inflammation and energy metabolism

BACKGROUND

Low-glycemic load dietary patterns, characterized by consumption of whole grains, legumes, fruits, and vegetables, are associated with reduced risk of several chronic diseases.

METHODS

Using samples from a randomized, controlled, crossover feeding trial, we evaluated the effects on metabolic profiles of a low-glycemic whole-grain dietary pattern (WG) compared with a dietary pattern high in refined grains and added sugars (RG) for 28 d. LC-MS-based targeted metabolomics analysis was performed on fasting plasma samples from 80 healthy participants (n = 40 men, n = 40 women) aged 18-45 y. Linear mixed models were used to evaluate differences in response between diets for individual metabolites. Kyoto Encyclopedia of Genes and Genomes (KEGG)-defined pathways and 2 novel data-driven analyses were conducted to consider differences at the pathway level.

RESULTS

There were 121 metabolites with detectable signal in >98% of all plasma samples. Eighteen metabolites were significantly different between diets at day 28 [false discovery rate (FDR) < 0.05]. Inositol, hydroxyphenylpyruvate, citrulline, ornithine, 13-hydroxyoctadecadienoic acid, glutamine, and oxaloacetate were higher after the WG diet than after the RG diet, whereas melatonin, betaine, creatine, acetylcholine, aspartate, hydroxyproline, methylhistidine, tryptophan, cystamine, carnitine, and trimethylamine were lower. Analyses using KEGG-defined pathways revealed statistically significant differences in tryptophan metabolism between diets, with kynurenine and melatonin positively associated with serum C-reactive protein concentrations. Novel data-driven methods at the metabolite and network levels found correlations among metabolites involved in branched-chain amino acid (BCAA) degradation, trimethylamine-N-oxide production, and β oxidation of fatty acids (FDR < 0.1) that differed between diets, with more favorable metabolic profiles detected after the WG diet. Higher BCAAs and trimethylamine were positively associated with homeostasis model assessment-insulin resistance.

CONCLUSIONS

These exploratory metabolomics results support beneficial effects of a low-glycemic load dietary pattern characterized by whole grains, legumes, fruits, and vegetables, compared with a diet high in refined grains and added sugars on inflammation and energy metabolism pathways. This trial was registered at clinicaltrials.gov as NCT00622661.

Nonparametric regression with adaptive truncation via a convex hierarchical penalty

We consider the problem of nonparametric regression with a potentially large number of covariates. We propose a convex, penalized estimation framework that is particularly well suited to high-dimensional sparse additive models and combines the appealing features of finite basis representation and smoothing penalties. In the case of additive models, a finite basis representation provides a parsimonious representation for fitted functions but is not adaptive when component functions possess different levels of complexity. In contrast, a smoothing spline-type penalty on the component functions is adaptive but does not provide a parsimonious representation. Our proposal simultaneously achieves parsimony and adaptivity in a computationally efficient way. We demonstrate these properties through empirical studies and show that our estimator converges at the minimax rate for functions within a hierarchical class. We further establish minimax rates for a large class of sparse additive models. We also develop an efficient algorithm that scales similarly to the lasso with the number of covariates and sample size.

Plasma Lipidomics Profiles After a Diet Characterized by Whole Grains Compared to a Diet High in Refined Grains and Added Sugars (FS03-07-19)

Objectives

Dietary patterns high in fiber from sources including whole grains, legumes, fruits, vegetables, nuts and seeds, are associated with lower risk of chronic disease, such as cardiovascular disease and cancer. We investigated how plasma lipidomics profiles differed between a diet high in whole grains (WG) versus a diet high in refined grains and added sugars (RG).

Methods

Using a randomized, crossover, controlled feeding study, 80 healthy participants (n = 40 men, n = 40 women, 40 normal weight, 40 overweight/obese), 18–45 y, were randomized to receive either a WG or RG diet for 28 days. After a 28-day washout period where participants resumed their habitual diet, they crossed over to the other diet. Targeted, differential mobility mass spectrometry was performed on fasting plasma samples collected at the baseline and end of each diet period and quantified the concentrations of 863 lipids from 13 classes. Paired t-tests and pairwise partial least squares-discriminant analysis (PLS-DA) were used to evaluate differences in lipid profiles between the two diets.

Results

At a class level, only ceramides were significantly different when comparing the two diets. After removing lipid species with > 20% missing values or CVs < 25%, 606 were retained for species analysis. Sixty-seven lipid species were significantly different between diets at day 28 (FDR < 0.05): 38 of 414 detected triglycerides, 9 of 59 phosphatidylethanolamines, 9 of 63 phosphatidylcholines, 4 of 22 cholesterol esters, 3 of 11 sphingomyelins, 2 of 13 lysophosphatidylcholines, and 1 of 5 ceramides. The majority of significant lipids were higher in plasma after the WG diet. PLSDA analysis showed the first and second components explaining 49% and 8.4%, respectively. Based on the selected components, lipidomic profiles showed fair separation for the two groups of diet. R2 values were 0.07 and 0.43, and Q2 values were -0.03 and 0.04 for components 1 and 2, respectively.

Conclusions

Higher concentrations of some lipid species such as cholesterol ester 12:0, a carrier of high-density lipoprotein, could indicate a favorable shift in lipid profiles. Further investigation using more complex models are being conducted.

Funding Sources

National Cancer Institute - National Institutes of Health.

Generalized Score Matching for Non-Negative Data

A common challenge in estimating parameters of probability density functions is the intractability of the normalizing constant. While in such cases maximum likelihood estimation may be implemented using numerical integration, the approach becomes computationally intensive. The score matching method of Hyvärinen (2005) avoids direct calculation of the normalizing constant and yields closed-form estimates for exponential families of continuous distributions over R m . Hyvärinen (2007) extended the approach to distributions supported on the non-negative orthant, R + m . In this paper, we give a generalized form of score matching for non-negative data that improves estimation efficiency. As an example, we consider a general class of pairwise interaction models. Addressing an overlooked inexistence problem, we generalize the regularized score matching method of Lin et al. (2016) and improve its theoretical guarantees for non-negative Gaussian graphical models.

The Reduced PC-Algorithm: Improved Causal Structure Learning in Large Random Networks

We consider the task of estimating a high-dimensional directed acyclic graph, given observations from a linear structural equation model with arbitrary noise distribution. By exploiting properties of common random graphs, we develop a new algorithm that requires conditioning only on small sets of variables. The proposed algorithm, which is essentially a modified version of the PC-Algorithm, offers significant gains in both computational complexity and estimation accuracy. In particular, it results in more efficient and accurate estimation in large networks containing hub nodes, which are common in biological systems. We prove the consistency of the proposed algorithm, and show that it also requires a less stringent faithfulness assumption than the PC-Algorithm. Simulations in low and high-dimensional settings are used to illustrate these findings. An application to gene expression data suggests that the proposed algorithm can identify a greater number of clinically relevant genes than current methods.

Combining Supervised and Unsupervised Learning for Improved miRNA Target Prediction

MicroRNAs (miRNAs) are short non-coding RNAs which bind to mRNAs and regulate their expression. MiRNAs have been found to be associated with initiation and progression of many complex diseases. Investigating miRNAs and their targets can thus help develop new therapies by designing anti-miRNA oligonucleotides. While existing computational approaches can predict miRNA targets, these predictions have low accuracy. In this paper, we propose a two-step approach to refine the results of sequence-based prediction algorithms. The first step, which is based on our previous work, uses an ensemble learning approach that combines multiple existing methods. The second step utilizes support vector machine (SVM) classifiers in one- and two-class modes to infer miRNA-mRNA interactions based on both binding features, as well as network features extracted from gene regulatory network. Experimental results using two real data sets from TCGA indicate that the use of two-class SVM classification significantly improves the precision of miRNA-mRNA prediction.