Sleep Problems and Autism Impairments in a Large Community Sample of Children and Adolescents

Sleep problems are common in individuals with autism spectrum disorder (ASD). How sleep problems reflect specific ASD phenotypes is unclear. We studied whether sleep problems indexed functional impairment in a heterogeneous community sample of individuals with ASD. We analyzed 977 probands (233 females; age = 11.27 ± 4.13 years) from the Rhode Island Consortium for Autism Research and Treatment dataset, a unique public-private-academic collaboration involving all major points of service for families in Rhode Island. We found that individuals with a confirmed diagnosis of ASD were more likely to have sleep problems. However, across the whole sample and above and beyond a formal diagnosis, sleep problems were dimensionally associated with worse social impairment and poorer adaptive functioning. By using a large dataset reflective of the diversity of presentations in the community, this study underscores the importance of considering sleep problems in clinical practice to improve adaptive functioning in individuals with ASD.

Identification and functional modelling of plausibly causative cis-regulatory variants in a highly-selected cohort with X-linked intellectual disability

Identifying causative variants in cis-regulatory elements (CRE) in neurodevelopmental disorders has proven challenging. We have used in vivo functional analyses to categorize rigorously filtered CRE variants in a clinical cohort that is plausibly enriched for causative CRE mutations: 48 unrelated males with a family history consistent with X-linked intellectual disability (XLID) in whom no detectable cause could be identified in the coding regions of the X chromosome (chrX). Targeted sequencing of all chrX CRE identified six rare variants in five affected individuals that altered conserved bases in CRE targeting known XLID genes and segregated appropriately in families. Two of these variants, FMR1CRE and TENM1CRE, showed consistent site- and stage-specific differences of enhancer function in the developing zebrafish brain using dual-color fluorescent reporter assay. Mouse models were created for both variants. In male mice Fmr1CRE induced alterations in neurodevelopmental Fmr1 expression, olfactory behavior and neurophysiological indicators of FMRP function. The absence of another likely causative variant on whole genome sequencing further supported FMR1CRE as the likely basis of the XLID in this family. Tenm1CRE mice showed no phenotypic anomalies. Following the release of gnomAD 2.1, reanalysis showed that TENM1CRE exceeded the maximum plausible population frequency of a XLID causative allele. Assigning causative status to any ultra-rare CRE variant remains problematic and requires disease-relevant in vivo functional data from multiple sources. The sequential and bespoke nature of such analyses renders them time-consuming and challenging to scale for routine clinical use.

DNA methylation analysis reveals epimutation hotspots in patients with dilated cardiomyopathy-associated laminopathies

BACKGROUND

Mutations in LMNA, encoding lamin A/C, lead to a variety of diseases known as laminopathies including dilated cardiomyopathy (DCM) and skeletal abnormalities. Though previous studies have investigated the dysregulation of gene expression in cells from patients with DCM, the role of epigenetic (gene regulatory) mechanisms, such as DNA methylation, has not been thoroughly investigated. Furthermore, the impact of family-specific LMNA mutations on DNA methylation is unknown. Here, we performed reduced representation bisulfite sequencing on ten pairs of fibroblasts and their induced pluripotent stem cell (iPSC) derivatives from two families with DCM due to distinct LMNA mutations, one of which also induces brachydactyly.

RESULTS

Family-specific differentially methylated regions (DMRs) were identified by comparing the DNA methylation landscape of patient and control samples. Fibroblast DMRs were found to enrich for distal regulatory features and transcriptionally repressed chromatin and to associate with genes related to phenotypes found in tissues affected by laminopathies. These DMRs, in combination with transcriptome-wide expression data and lamina-associated domain (LAD) organization, revealed the presence of inter-family epimutation hotspots near differentially expressed genes, most of which were located outside LADs redistributed in LMNA-related DCM. Comparison of DMRs found in fibroblasts and iPSCs identified regions where epimutations were persistent across both cell types. Finally, a network of aberrantly methylated disease-associated genes revealed a potential molecular link between pathways involved in bone and heart development.

CONCLUSIONS

Our results identified both shared and mutation-specific laminopathy epimutation landscapes that were consistent with lamin A/C mutation-mediated epigenetic aberrancies that arose in somatic and early developmental cell stages.

Revisiting the concept of activation in supercooled liquids

In this work, we revisit the description of dynamics based on the concepts of metabasins and activation in mildly supercooled liquids via the analysis of the dynamics of a paradigmatic glass former between its onset temperature [Formula: see text] and mode-coupling temperature [Formula: see text]. First, we provide measures that demonstrate that the onset of glassiness is indeed connected to the landscape, and that metabasin waiting time distributions are so broad that the system can remain stuck in a metabasin for times that exceed [Formula: see text] by orders of magnitude. We then reanalyze the transitions between metabasins, providing several indications that the standard picture of activated dynamics in terms of traps does not hold in this regime. Instead, we propose that here activation is principally driven by entropic instead of energetic barriers. In particular, we illustrate that activation is not controlled by the hopping of high energetic barriers and should more properly be interpreted as the entropic selection of nearly barrierless but rare pathways connecting metabasins on the landscape.

Quantifying cell transitions in C. elegans with data-fitted landscape models

Increasing interest has emerged in new mathematical approaches that simplify the study of complex differentiation processes by formalizing Waddington's landscape metaphor. However, a rational method to build these landscape models remains an open problem. Here we study vulval development in C. elegans by developing a framework based on Catastrophe Theory (CT) and approximate Bayesian computation (ABC) to build data-fitted landscape models. We first identify the candidate qualitative landscapes, and then use CT to build the simplest model consistent with the data, which we quantitatively fit using ABC. The resulting model suggests that the underlying mechanism is a quantifiable two-step decision controlled by EGF and Notch-Delta signals, where a non-vulval/vulval decision is followed by a bistable transition to the two vulval states. This new model fits a broad set of data and makes several novel predictions.

A comparison of neuronal population dynamics measured with calcium imaging and electrophysiology

Calcium imaging with fluorescent protein sensors is widely used to record activity in neuronal populations. The transform between neural activity and calcium-related fluorescence involves nonlinearities and low-pass filtering, but the effects of the transformation on analyses of neural populations are not well understood. We compared neuronal spikes and fluorescence in matched neural populations in behaving mice. We report multiple discrepancies between analyses performed on the two types of data, including changes in single-neuron selectivity and population decoding. These were only partially resolved by spike inference algorithms applied to fluorescence. To model the relation between spiking and fluorescence we simultaneously recorded spikes and fluorescence from individual neurons. Using these recordings we developed a model transforming spike trains to synthetic-imaging data. The model recapitulated the differences in analyses. Our analysis highlights challenges in relating electrophysiology and imaging data, and suggests forward modeling as an effective way to understand differences between these data.

Efficient calculation of heterogeneous non-equilibrium statistics in coupled firing-rate models

Understanding nervous system function requires careful study of transient (non-equilibrium) neural response to rapidly changing, noisy input from the outside world. Such neural response results from dynamic interactions among multiple, heterogeneous brain regions. Realistic modeling of these large networks requires enormous computational resources, especially when high-dimensional parameter spaces are considered. By assuming quasi-steady-state activity, one can neglect the complex temporal dynamics; however, in many cases the quasi-steady-state assumption fails. Here, we develop a new reduction method for a general heterogeneous firing-rate model receiving background correlated noisy inputs that accurately handles highly non-equilibrium statistics and interactions of heterogeneous cells. Our method involves solving an efficient set of nonlinear ODEs, rather than time-consuming Monte Carlo simulations or high-dimensional PDEs, and it captures the entire set of first and second order statistics while allowing significant heterogeneity in all model parameters.

ChIPulate: A comprehensive ChIP-seq simulation pipeline

ChIP-seq (Chromatin Immunoprecipitation followed by sequencing) is a high-throughput technique to identify genomic regions that are bound in vivo by a particular protein, e.g., a transcription factor (TF). Biological factors, such as chromatin state, indirect and cooperative binding, as well as experimental factors, such as antibody quality, cross-linking, and PCR biases, are known to affect the outcome of ChIP-seq experiments. However, the relative impact of these factors on inferences made from ChIP-seq data is not entirely clear. Here, via a detailed ChIP-seq simulation pipeline, ChIPulate, we assess the impact of various biological and experimental sources of variation on several outcomes of a ChIP-seq experiment, viz., the recoverability of the TF binding motif, accuracy of TF-DNA binding detection, the sensitivity of inferred TF-DNA binding strength, and number of replicates needed to confidently infer binding strength. We find that the TF motif can be recovered despite poor and non-uniform extraction and PCR amplification efficiencies. The recovery of the motif is, however, affected to a larger extent by the fraction of sites that are either cooperatively or indirectly bound. Importantly, our simulations reveal that the number of ChIP-seq replicates needed to accurately measure in vivo occupancy at high-affinity sites is larger than the recommended community standards. Our results establish statistical limits on the accuracy of inferences of protein-DNA binding from ChIP-seq and suggest that increasing the mean extraction efficiency, rather than amplification efficiency, would better improve sensitivity. The source code and instructions for running ChIPulate can be found at https://github.com/vishakad/chipulate.

DNA-binding proteins perform many key roles in biology, such as transcriptional regulation of gene expression and chromatin modification. ChIP-seq (Chromatin immunoprecipitation followed by high-throughput sequencing) is a widely used experimental technique to identify DNA-binding sites of specific proteins of interest, within cells, genome-wide. DNA fragments from genomic regions that are bound by a protein of interest, often a transcription factor (TF), are selectively extracted using specific antibodies, amplified using PCR, and sequenced. The sequences are mapped to the reference genome. Regions where many sequences map, called “peaks”, are used to infer the location of TF-bound loci (peaks), in vivo occupancy at those loci, and the sequence pattern (motif) to which the TF shows a binding affinity. But measurements of TF occupancy and motif inference are vulnerable to several biological and experimental sources of variation that are poorly understood and difficult to assess directly. Here, we simulate key steps of the ChIP-seq protocol with the aim of estimating the relative effects of various sources of variations on motif inference and binding affinity estimations. Besides providing specific insights and recommendations, we provide a general framework to simulate sequence reads in a ChIP-seq experiment, which should considerably aid in the development of software aimed at analyzing ChIP-seq data.

Successful attentional set-shifting in 2-year-olds with and without Autism Spectrum Disorder

The development of executive function is necessary for flexible and voluntary control of behavior. Deficits in executive function are purported to be a primary cause of behavioral inflexibility—a core clinical symptom—in Autism Spectrum Disorder (ASD). Attentional set-shifting has traditionally been measured with the Dimensional Change Card Sort, however, this task requires following verbal instructions. Here, we used a novel visual search task that does not require verbal instructions in conjunction with eye-tracking to test attentional set-shifting in 2-year-old toddlers diagnosed with ASD (N = 29) and chronological age-matched typically developing controls (N = 30). On each trial, a relevant and an irrelevant target were embedded in a set of feature-conjunction distractors, and toddlers were tasked with searching for the relevant target. Critically, after a set of trials the targets switched roles (i.e., the previously relevant target became irrelevant, and the previously relevant target became irrelevant). We measured visual search performance prior to and following a target switch. We found that both groups of toddlers could readily switch targets, and found strikingly similar performance between typically developing toddlers and toddlers with ASD. Our results challenge the centrality of deficits in attentional set-shifting to early behavioral inflexibility in ASD.

Modeling craniofacial development reveals spatiotemporal constraints on robust patterning of the mandibular arch

How does pattern formation occur accurately when confronted with tissue growth and stochastic fluctuations (noise) in gene expression? Dorso-ventral (D-V) patterning of the mandibular arch specifies upper versus lower jaw skeletal elements through a combination of Bone morphogenetic protein (Bmp), Endothelin-1 (Edn1), and Notch signaling, and this system is highly robust. We combine NanoString experiments of early D-V gene expression with live imaging of arch development in zebrafish to construct a computational model of the D-V mandibular patterning network. The model recapitulates published genetic perturbations in arch development. Patterning is most sensitive to changes in Bmp signaling, and the temporal order of gene expression modulates the response of the patterning network to noise. Thus, our integrated systems biology approach reveals non-intuitive features of the complex signaling system crucial for craniofacial development, including novel insights into roles of gene expression timing and stochasticity in signaling and gene regulation.

A constructive approach to the epistemological problem of emergence in complex systems

Emergent patterns in complex systems are related with many intriguing phenomena in modern science. One question that has sparked vigorous debates is if difficulties in the modelization of emergent behaviours are a consequence of ontological or epistemological limitations. To elucidate this question, we propose a novel approximation through constructive logic. Under this framework, experimental measurements will be considered conceptual building blocks from which we aim to achieve a description of the microstates ensemble mapping the macroscopic emergent observation. This procedure allow us to have full control of any information loss, thus making the analysis of different systems fairly comparable. In particular, we aim to look for compact descriptions of the constraints underlying a dynamical system, as a necessary a priori step to develop explanatory (mechanistic) models. We apply our proposal to a synthetic system to show that the number and scope of the system’s constraints hinder our ability to build compact descriptions, being those systems under global constraints a limiting case in which such a description is unreachable. This result clearly links the epistemological limits of the framework selected with an ontological feature of the system, leading us to propose a definition of emergence strength which we make compatible with the scientific method through the active intervention of the observer on the system, following the spirit of Granger causality. We think that our approximation clarifies previous discrepancies found in the literature, reconciles distinct attempts to classify emergent processes, and paves the way to understand other challenging concepts such as downward causation.

Inferring decoding strategies for multiple correlated neural populations

Studies of neuron-behaviour correlation and causal manipulation have long been used separately to understand the neural basis of perception. Yet these approaches sometimes lead to drastically conflicting conclusions about the functional role of brain areas. Theories that focus only on choice-related neuronal activity cannot reconcile those findings without additional experiments involving large-scale recordings to measure interneuronal correlations. By expanding current theories of neural coding and incorporating results from inactivation experiments, we demonstrate here that it is possible to infer decoding weights of different brain areas at a coarse scale without precise knowledge of the correlation structure. We apply this technique to neural data collected from two different cortical areas in macaque monkeys trained to perform a heading discrimination task. We identify two opposing decoding schemes, each consistent with data depending on the nature of correlated noise. Our theory makes specific testable predictions to distinguish these scenarios experimentally without requiring measurement of the underlying noise correlations.

The neocortex is structurally organized into distinct brain areas. The role of specific brain areas in sensory perception is typically studied using two kinds of laboratory experiments: those that measure correlations between neural activity and reported percepts, and those that inactivate a brain region and measure the resulting changes in percepts. The two types of experiments have generally been interpreted in isolation, in part because no theory has been able combine their outcomes. Here, we describe a mathematical framework that synthesizes both kinds of results, giving us a new way to assess how different brain areas contribute to perception. When we apply our framework to experiments on behaving monkeys, we discover two models that can explain the perplexing finding that one brain area can predict an animal’s reported percepts, even though the percepts are not affected when that brain area is inactivated. The two models ascribe dramatically different efficiencies to brain computation. We show that these two models could be distinguished by a proposed experiment that measures correlations while inactivating different brain areas.

Complete functional mapping of infection- and vaccine-elicited antibodies against the fusion peptide of HIV

Eliciting broadly neutralizing antibodies (bnAbs) targeting envelope (Env) is a major goal of HIV vaccine development, but cross-clade breadth from immunization has only sporadically been observed. Recently, Xu et al (2018) elicited cross-reactive neutralizing antibody responses in a variety of animal models using immunogens based on the epitope of bnAb VRC34.01. The VRC34.01 antibody, which was elicited by natural human infection, targets the N terminus of the Env fusion peptide, a critical component of the virus entry machinery. Here we precisely characterize the functional epitopes of VRC34.01 and two vaccine-elicited murine antibodies by mapping all single amino-acid mutations to the BG505 Env that affect viral neutralization. While escape from VRC34.01 occurred via mutations in both fusion peptide and distal interacting sites of the Env trimer, escape from the vaccine-elicited antibodies was mediated predominantly by mutations in the fusion peptide. Cryo-electron microscopy of four vaccine-elicited antibodies in complex with Env trimer revealed focused recognition of the fusion peptide and provided a structural basis for development of neutralization breadth. Together, these functional and structural data suggest that the breadth of vaccine-elicited antibodies targeting the fusion peptide can be enhanced by specific interactions with additional portions of Env. Thus, our complete maps of viral escape both delineate pathways of resistance to these fusion peptide-directed antibodies and provide a strategy to improve the breadth or potency of future vaccine-induced antibodies against Env’s fusion peptide.

A major goal of HIV-1 vaccine design is to elicit antibodies that neutralize diverse strains of HIV-1. Recently, some of us elicited such antibodies in animal models using immunogens based on the epitope of a broad antibody (VRC34.01) isolated from an infected individual. Further improving these vaccine-elicited antibody responses will require a detailed understanding of how the resulting antibodies target HIV’s envelope protein (Env). Here, we used mutational antigenic profiling to precisely map the epitopes of two vaccine-elicited antibodies and the template VRC34.01 antibody. We did this by quantifying the effect of all possible amino acid mutations to Env on antibody neutralization. Although all antibodies target a similar region of Env, we found clear differences in the functional interaction of Env with these vaccine- and infection-elicited antibodies. We combined these functional data with structural analyses to identify antibody–Env interactions that may contribute to the relatively greater breadth of the infection-elicited antibody and could improve the breadth of vaccine-elicited antibodies. These data thereby help to refine vaccination schemes to achieve broader responses.

Computable performance guarantees for compressed sensing matrices

The null space condition for ℓ1 minimization in compressed sensing is a necessary and sufficient condition on the sensing matrices under which a sparse signal can be uniquely recovered from the observation data via ℓ1 minimization. However, verifying the null space condition is known to be computationally challenging. Most of the existing methods can provide only upper and lower bounds on the proportion parameter that characterizes the null space condition. In this paper, we propose new polynomial-time algorithms to establish upper bounds of the proportion parameter. We leverage on these techniques to find upper bounds and further develop a new procedure-tree search algorithm-that is able to precisely and quickly verify the null space condition. Numerical experiments show that the execution speed and accuracy of the results obtained from our methods far exceed those of the previous methods which rely on linear programming (LP) relaxation and semidefinite programming (SDP).

Interactions between species introduce spurious associations in microbiome studies

Microbiota contribute to many dimensions of host phenotype, including disease. To link specific microbes to specific phenotypes, microbiome-wide association studies compare microbial abundances between two groups of samples. Abundance differences, however, reflect not only direct associations with the phenotype, but also indirect effects due to microbial interactions. We found that microbial interactions could easily generate a large number of spurious associations that provide no mechanistic insight. Using techniques from statistical physics, we developed a method to remove indirect associations and applied it to the largest dataset on pediatric inflammatory bowel disease. Our method corrected the inflation of p-values in standard association tests and showed that only a small subset of associations is directly linked to the disease. Direct associations had a much higher accuracy in separating cases from controls and pointed to immunomodulation, butyrate production, and the brain-gut axis as important factors in the inflammatory bowel disease.

Temporal dynamics of neurogenomic plasticity in response to social interactions in male threespined sticklebacks

Animals exhibit dramatic immediate behavioral plasticity in response to social interactions, and brief social interactions can shape the future social landscape. However, the molecular mechanisms contributing to behavioral plasticity are unclear. Here, we show that the genome dynamically responds to social interactions with multiple waves of transcription associated with distinct molecular functions in the brain of male threespined sticklebacks, a species famous for its behavioral repertoire and evolution. Some biological functions (e.g., hormone activity) peaked soon after a brief territorial challenge and then declined, while others (e.g., immune response) peaked hours afterwards. We identify transcription factors that are predicted to coordinate waves of transcription associated with different components of behavioral plasticity. Next, using H3K27Ac as a marker of chromatin accessibility, we show that a brief territorial intrusion was sufficient to cause rapid and dramatic changes in the epigenome. Finally, we integrate the time course brain gene expression data with a transcriptional regulatory network, and link gene expression to changes in chromatin accessibility. This study reveals rapid and dramatic epigenomic plasticity in response to a brief, highly consequential social interaction.

Social interactions provoke changes in the brain and behavior but their underlying molecular mechanisms remain obscure. Male sticklebacks are small fish whose fitness depends on their ability to defend a territory. Here, by measuring the time course of gene expression in response to a territorial challenge in two brain regions, we show that a single brief territorial intrusion provoked waves of gene expression that persisted for hours afterwards, with waves of transcription associated with distinct biological processes. Moreover, a single territorial challenge caused dramatic changes to the epigenome. Changes in chromatin accessibility corresponded to changes in gene expression, and to the activity of transcription factors operating within gene regulatory networks. This study reveals rapid and dramatic epigenomic plasticity in response to a brief, highly consequential social interaction. These results suggest that meaningful social interactions (even brief ones) can provoke waves of transcription and changes to the epigenome which lead to changes in neural functioning, and those changes are a mechanism by which animals update their assessment of their social world.

Linking structure and activity in nonlinear spiking networks

Recent experimental advances are producing an avalanche of data on both neural connectivity and neural activity. To take full advantage of these two emerging datasets we need a framework that links them, revealing how collective neural activity arises from the structure of neural connectivity and intrinsic neural dynamics. This problem of structure-driven activity has drawn major interest in computational neuroscience. Existing methods for relating activity and architecture in spiking networks rely on linearizing activity around a central operating point and thus fail to capture the nonlinear responses of individual neurons that are the hallmark of neural information processing. Here, we overcome this limitation and present a new relationship between connectivity and activity in networks of nonlinear spiking neurons by developing a diagrammatic fluctuation expansion based on statistical field theory. We explicitly show how recurrent network structure produces pairwise and higher-order correlated activity, and how nonlinearities impact the networks’ spiking activity. Our findings open new avenues to investigating how single-neuron nonlinearities—including those of different cell types—combine with connectivity to shape population activity and function.

Neuronal networks, like many biological systems, exhibit variable activity. This activity is shaped by both the underlying biology of the component neurons and the structure of their interactions. How can we combine knowledge of these two things—that is, models of individual neurons and of their interactions—to predict the statistics of single- and multi-neuron activity? Current approaches rely on linearizing neural activity around a stationary state. In the face of neural nonlinearities, however, these linear methods can fail to predict spiking statistics and even fail to correctly predict whether activity is stable or pathological. Here, we show how to calculate any spike train cumulant in a broad class of models, while systematically accounting for nonlinear effects. We then study a fundamental effect of nonlinear input-rate transfer–coupling between different orders of spiking statistic–and how this depends on single-neuron and network properties.

Robust information propagation through noisy neural circuits

Sensory neurons give highly variable responses to stimulation, which can limit the amount of stimulus information available to downstream circuits. Much work has investigated the factors that affect the amount of information encoded in these population responses, leading to insights about the role of covariability among neurons, tuning curve shape, etc. However, the informativeness of neural responses is not the only relevant feature of population codes; of potentially equal importance is how robustly that information propagates to downstream structures. For instance, to quantify the retina’s performance, one must consider not only the informativeness of the optic nerve responses, but also the amount of information that survives the spike-generating nonlinearity and noise corruption in the next stage of processing, the lateral geniculate nucleus. Our study identifies the set of covariance structures for the upstream cells that optimize the ability of information to propagate through noisy, nonlinear circuits. Within this optimal family are covariances with “differential correlations”, which are known to reduce the information encoded in neural population activities. Thus, covariance structures that maximize information in neural population codes, and those that maximize the ability of this information to propagate, can be very different. Moreover, redundancy is neither necessary nor sufficient to make population codes robust against corruption by noise: redundant codes can be very fragile, and synergistic codes can—in some cases—optimize robustness against noise.

Information about the outside world, which originates in sensory neurons, propagates through multiple stages of processing before reaching the neural structures that control behavior. While much work in neuroscience has investigated the factors that affect the amount of information contained in peripheral sensory areas, very little work has asked how much of that information makes it through subsequent processing stages. That’s the focus of this paper, and it’s an important issue because information that fails to propagate cannot be used to affect decision-making. We find a tradeoff between information content and information transmission: neural codes which contain a large amount of information can transmit that information poorly to subsequent processing stages. Thus, the problem of robust information propagation—which has largely been overlooked in previous research—may be critical for determining how our sensory organs communicate with our brains. We identify the conditions under which information propagates well—or poorly—through multiple stages of neural processing.

Neurogenetic analysis of childhood disintegrative disorder

BACKGROUND

Childhood disintegrative disorder (CDD) is a rare form of autism spectrum disorder (ASD) of unknown etiology. It is characterized by late-onset regression leading to significant intellectual disability (ID) and severe autism. Although there are phenotypic differences between CDD and other forms of ASD, it is unclear if there are neurobiological differences.

METHODS

We pursued a multidisciplinary study of CDD (n = 17) and three comparison groups: low-functioning ASD (n = 12), high-functioning ASD (n = 50), and typically developing (n = 26) individuals. We performed whole-exome sequencing (WES), copy number variant (CNV), and gene expression analyses of CDD and, on subsets of each cohort, non-sedated functional magnetic resonance imaging (fMRI) while viewing socioemotional (faces) and non-socioemotional (houses) stimuli and eye tracking while viewing emotional faces.

RESULTS

We observed potential differences between CDD and other forms of ASD. WES and CNV analyses identified one or more rare de novo, homozygous, and/or hemizygous (mother-to-son transmission on chrX) variants for most probands that were not shared by unaffected sibling controls. There were no clearly deleterious variants or highly recurrent candidate genes. Candidate genes that were found to be most conserved at variant position and most intolerant of variation, such as TRRAP, ZNF236, and KIAA2018, play a role or may be involved in transcription. Using the human BrainSpan transcriptome dataset, CDD candidate genes were found to be more highly expressed in non-neocortical regions than neocortical regions. This expression profile was similar to that of an independent cohort of ASD probands with regression. The non-neocortical regions overlapped with those identified by fMRI as abnormally hyperactive in response to viewing faces, such as the thalamus, cerebellum, caudate, and hippocampus. Eye-tracking analysis showed that, among individuals with ASD, subjects with CDD focused on eyes the most when shown pictures of faces.

CONCLUSIONS

Given that cohort sizes were limited by the rarity of CDD, and the challenges of conducting non-sedated fMRI and eye tracking in subjects with ASD and significant ID, this is an exploratory study designed to investigate the neurobiological features of CDD. In addition to reporting the first multimodal analysis of CDD, a combination of fMRI and eye-tracking analyses are being presented for the first time for low-functioning individuals with ASD. Our results suggest differences between CDD and other forms of ASD on the neurobiological as well as clinical level.

Exaggerated CpH methylation in the autism-affected brain

BACKGROUND

The etiology of autism, a complex, heritable, neurodevelopmental disorder, remains largely unexplained. Given the unexplained risk and recent evidence supporting a role for epigenetic mechanisms in the development of autism, we explored the role of CpG and CpH (H = A, C, or T) methylation within the autism-affected cortical brain tissue.

METHODS

Reduced representation bisulfite sequencing (RRBS) was completed, and analysis was carried out in 63 post-mortem cortical brain samples (Brodmann area 19) from 29 autism-affected and 34 control individuals. Analyses to identify single sites that were differentially methylated and to identify any global methylation alterations at either CpG or CpH sites throughout the genome were carried out.

RESULTS

We report that while no individual site or region of methylation was significantly associated with autism after multi-test correction, methylated CpH dinucleotides were markedly enriched in autism-affected brains (~2-fold enrichment at p < 0.05 cutoff, p = 0.002).

CONCLUSIONS

These results further implicate epigenetic alterations in pathobiological mechanisms that underlie autism.

Quality filtering of Illumina index reads mitigates sample cross-talk

BACKGROUND

Multiplexing multiple samples during Illumina sequencing is a common practice and is rapidly growing in importance as the throughput of the platform increases. Misassignments during de-multiplexing, where sequences are associated with the wrong sample, are an overlooked error mode on the Illumina sequencing platform. This results in a low rate of cross-talk among multiplexed samples and can cause detrimental effects in studies requiring the detection of rare variants or when multiplexing a large number of samples.

RESULTS

We observed rates of cross-talk averaging 0.24 % when multiplexing 14 different samples with unique i5 and i7 index sequences. This cross-talk rate corresponded to 254,632 misassigned reads on a single lane of the Illumina HiSeq 2500. Notably, all types of misassignment occur at similar rates: incorrect i5, incorrect i7, and incorrect sequence reads. We demonstrate that misassignments can be nearly eliminated by quality filtering of index reads while preserving about 90 % of the original sequences.

CONCLUSIONS

Cross-talk among multiplexed samples is a significant error mode on the Illumina platform, especially if samples are only separated by a single unique index. Quality filtering of index sequences offers an effective solution to minimizing cross-talk among samples. Furthermore, we propose a straightforward method for verifying the extent of cross-talk between samples and optimizing quality score thresholds that does not require additional control samples and can even be performed post hoc on previous runs.

Unconventional height functions in simultaneous Diophantine approximation

Simultaneous Diophantine approximation is concerned with the approximation of a point x ∈ R d by points r ∈ Q d , with a view towards jointly minimizing the quantities ‖ x - r ‖ and H ( r ) . Here H ( r ) is the so-called "standard height" of the rational point r . In this paper the authors ask: What changes if we replace the standard height function by a different one? As it turns out, this change leads to dramatic differences from the classical theory and requires the development of new methods. We discuss three examples of nonstandard height functions, computing their exponents of irrationality as well as giving more precise results. A list of open questions is also given.

Multi-Layer Mechanical Model of Glagov Remodeling in Coronary Arteries: Differences between In-Vivo and Ex-Vivo Measurements

When blood vessels undergo remodeling because of the buildup of atherosclerotic plaque, it is thought that they first undergo compensatory or outward remodeling, followed by inward remodeling: the lumen area stays roughly constant or increases slightly and then decreases rapidly. The second phase of remodeling is supposed to start after the plaque burden exceeds about 40%. These changes in the vessel were first observed by S. Glagov who examined cross-sections of coronary arteries at different stages of the disease. In this paper, we use a mathematical model based on growth and elasticity theory to verify the main aspects of Glagov’s result. However, both our model and curve-fitting to the data suggest that the critical stenosis is around 20% rather than 40%. Our model and data from the PROSPECT trial also show that Glagov remodeling is qualitatively different depending on whether measurements are taken ex-vivo or in-vivo. Our results suggest that the first outward phase of “Glagov remodeling” is largely absent for in-vivo measurements: that is, the lumen area always decreases as plaque builds up. We advocate that care must be taken when infering how in-vivo vessels remodel from ex-vivo data.

Constrained Allocation Flux Balance Analysis

New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-wide constraint. Its roots lie in the experimentally observed pattern of proteome allocation for metabolic functions, allowing to bridge regulation and metabolism in a transparent way under the principle of growth-rate maximization. We provide a simple method to solve CAFBA efficiently and propose an “ensemble averaging” procedure to account for unknown protein costs. Applying this approach to modeling E. coli metabolism, we find that, as the growth rate increases, CAFBA solutions cross over from respiratory, growth-yield maximizing states (preferred at slow growth) to fermentative states with carbon overflow (preferred at fast growth). In addition, CAFBA allows for quantitatively accurate predictions on the rate of acetate excretion and growth yield based on only 3 parameters determined by empirical growth laws.

The intracellular protein levels of exponentially growing bacteria are known to vary strongly with growth conditions, as described by quantitative “growth laws”. This work introduces a computational genome-scale framework (Constrained Allocation Flux Balance Analysis, CAFBA) which incorporates growth laws into canonical Flux Balance Analysis. Upon introducing 3 parameters based on established growth laws for E. coli, CAFBA accurately reproduces empirical results on the growth-rate dependent rate of carbon overflow and growth yield, and generates testable predictions about cellular energetic strategies and protein expression levels. CAFBA therefore provides a simple, quantitative approach to balancing the trade-off between growth and its associated biosynthetic costs at genome-scale, without the burden of tuning many inaccessible parameters.

Neuroprosthetic Decoder Training as Imitation Learning

Neuroprosthetic brain-computer interfaces function via an algorithm which decodes neural activity of the user into movements of an end effector, such as a cursor or robotic arm. In practice, the decoder is often learned by updating its parameters while the user performs a task. When the user's intention is not directly observable, recent methods have demonstrated value in training the decoder against a surrogate for the user's intended movement. Here we show that training a decoder in this way is a novel variant of an imitation learning problem, where an oracle or expert is employed for supervised training in lieu of direct observations, which are not available. Specifically, we describe how a generic imitation learning meta-algorithm, dataset aggregation (DAgger), can be adapted to train a generic brain-computer interface. By deriving existing learning algorithms for brain-computer interfaces in this framework, we provide a novel analysis of regret (an important metric of learning efficacy) for brain-computer interfaces. This analysis allows us to characterize the space of algorithmic variants and bounds on their regret rates. Existing approaches for decoder learning have been performed in the cursor control setting, but the available design principles for these decoders are such that it has been impossible to scale them to naturalistic settings. Leveraging our findings, we then offer an algorithm that combines imitation learning with optimal control, which should allow for training of arbitrary effectors for which optimal control can generate goal-oriented control. We demonstrate this novel and general BCI algorithm with simulated neuroprosthetic control of a 26 degree-of-freedom model of an arm, a sophisticated and realistic end effector.