Identification of complex molecular networks underlying common human phenotypes is a major challenge of modern genetics. In this study, we develop a method for network-based analysis of genetic associations (NETBAG). We use NETBAG to identify a large biological network of genes affected by rare de novo CNVs in autism. The genes forming the network are primarily related to synapse development, axon targeting, and neuron motility. The identified network is strongly related to genes previously implicated in autism and intellectual disability phenotypes. Our results are also consistent with the hypothesis that significantly stronger functional perturbations are required to trigger the autistic phenotype in females compared to males. Overall, the presented analysis of de novo variants supports the hypothesis that perturbed synaptogenesis is at the heart of autism. More generally, our study provides proof of the principle that networks underlying complex human phenotypes can be identified by a network-based functional analysis of rare genetic variants.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

Human exome resequencing using commercial target capture kits has been and is being used for sequencing large numbers of individuals to search for variants associated with various human diseases. We rigorously evaluated the capabilities of two solution exome capture kits. These analyses help clarify the strengths and limitations of those data as well as systematically identify variables that should be considered in the use of those data.

Ribosome profiling (Ribo-seq) is a powerful technique for measuring protein translation; however, sampling errors and biological biases are prevalent and poorly understood. Addressing these issues, we present Scikit-ribo (https://github.com/schatzlab/scikit-ribo), an open-source analysis package for accurate genome-wide A-site prediction and translation efficiency (TE) estimation from Ribo-seq and RNA sequencing data. Scikit-ribo accurately identifies A-site locations and reproduces codon elongation rates using several digestion protocols (r = 0.99). Next, we show that the commonly used reads per kilobase of transcript per million mapped reads-derived TE estimation is prone to biases, especially for low-abundance genes. Scikit-ribo introduces a codon-level generalized linear model with ridge penalty that correctly estimates TE, while accommodating variable codon elongation rates and mRNA secondary structure. This corrects the TE errors for over 2,000 genes in S. cerevisiae, which we validate using mass spectrometry of protein abundances (r = 0.81), and allows us to determine the Kozak-like sequence directly from Ribo-seq. We conclude with an analysis of coverage requirements needed for robust codon-level analysis and quantify the artifacts that can occur from cycloheximide treatment.

Human Phenotype Ontology (HPO) has risen as a useful tool for precision medicine by providing a standardized vocabulary of phenotypic abnormalities to describe presentations of human pathologies; however, there have been relatively few reports combining whole genome sequencing (WGS) and HPO, especially in the context of structural variants.

Nonalcoholic fatty liver disease (NAFLD) is a very common disorder affecting between 20 and 30% of adults in the United States. However, there is no effective pharmacotherapy for treating NAFLD. Niacin, a water-soluble vitamin (B3), at pharmacological doses, decreases hepatic triglyceride (TG) content in NAFLD through inhibition of diacylglycerol acyltransferase 2, a key enzyme that catalyzes the final step in TG synthesis. Alternatively, some studies indicate that niacin induces fatty liver in high-fat diet (HFD)-fed rats. Therefore, in this study we investigated whether niacin is beneficial in treating NAFLD in two strains of mice, C57BL/6J (B6) and B6129SF2/J (B6129) mice, with 20 weeks of HFD feeding. Niacin treatment was started from week 5 until the end of the study. Niacin treatment increased normalized liver weight, hepatic TG content and NAFLD score in HFD-fed B6129 mice but had no impact on B6 mice. Metabolomics analysis revealed that in B6129 mice, 4-hydroxyphenylpyruvic acid (4-HPP), which is associated with fatty acid oxidation, did not change with HFD feeding but significantly decreased with niacin treatment. Lipidomics analysis discovered that the abundance of phosphocholine (PC), which is critical for very low-density lipoprotein (VLDL)-TG production and secretion, was decreased in HFD-fed B6129 with niacin treatment. In conclusion, niacin had no impact on diet-induced NAFLD development in B6 mice but potentiated hepatic steatosis in HFD-fed B6129 mice due to impaired fatty acid oxidation and decreased VLDL-TG production and secretion.

Type 2 diabetes mellitus (T2DM) and anemia are related to some cardiovascular diseases and can predict poor outcomes. Both of them can damage the heart in their own ways, but their combined effects have not been well explored. This study aimed to explore the combined effects of T2DM and anemia and the interaction between left atrial (LA) and left ventricular (LV) function by cardiac magnetic resonance (CMR).

INDELs, especially those disrupting protein-coding regions of the genome, have been strongly associated with human diseases. However, there are still many errors with INDEL variant calling, driven by library preparation, sequencing biases, and algorithm artifacts.

We describe an X-linked genetic syndrome associated with mutations in TAF1 and manifesting with global developmental delay, intellectual disability (ID), characteristic facial dysmorphology, generalized hypotonia, and variable neurologic features, all in male individuals. Simultaneous studies using diverse strategies led to the identification of nine families with overlapping clinical presentations and affected by de novo or maternally inherited single-nucleotide changes. Two additional families harboring large duplications involving TAF1 were also found to share phenotypic overlap with the probands harboring single-nucleotide changes, but they also demonstrated a severe neurodegeneration phenotype. Functional analysis with RNA-seq for one of the families suggested that the phenotype is associated with downregulation of a set of genes notably enriched with genes regulated by E-box proteins. In addition, knockdown and mutant studies of this gene in zebrafish have shown a quantifiable, albeit small, effect on a neuronal phenotype. Our results suggest that mutations in TAF1 play a critical role in the development of this X-linked ID syndrome.

We present an open-source algorithm, Scalpel (http://scalpel.sourceforge.net/), which combines mapping and assembly for sensitive and specific discovery of insertions and deletions (indels) in exome-capture data. A detailed repeat analysis coupled with a self-tuning k-mer strategy allows Scalpel to outperform other state-of-the-art approaches for indel discovery, particularly in regions containing near-perfect repeats. We analyzed 593 families from the Simons Simplex Collection and demonstrated Scalpel's power to detect long (≥30 bp) transmitted events and enrichment for de novo likely gene-disrupting indels in autistic children.

Physicochemical properties of compounds have been instrumental in selecting lead compounds with increased drug-likeness. However, the relationship between physicochemical properties of constituent drugs and the tendency to exhibit drug interaction has not been systematically studied. We assembled physicochemical descriptors for a set of antifungal compounds ("drugs") previously examined for interaction. Analyzing the relationship between molecular weight, lipophilicity, H-bond donor, and H-bond acceptor values for drugs and their propensity to show pairwise antifungal drug synergy, we found that combinations of two lipophilic drugs had a greater tendency to show drug synergy. We developed a more refined decision tree model that successfully predicted drug synergy in stringent cross-validation tests based on only lipophilicity of drugs. Our predictions achieved a precision of 63% and allowed successful prediction for 58% of synergistic drug pairs, suggesting that this phenomenon can extend our understanding for a substantial fraction of synergistic drug interactions. We also generated and analyzed a large-scale synergistic human toxicity network, in which we observed that combinations of lipophilic compounds show a tendency for increased toxicity. Thus, lipophilicity, a simple and easily determined molecular descriptor, is a powerful predictor of drug synergy. It is well established that lipophilic compounds (i) are promiscuous, having many targets in the cell, and (ii) often penetrate into the cell via the cellular membrane by passive diffusion. We discuss the positive relationship between drug lipophilicity and drug synergy in the context of potential drug synergy mechanisms.

The ribosomal p70 S6 Kinase 1 (S6K1) has been implicated in the etiology of complex neurological diseases including autism, depression and dementia. Though no major gene disruption has been reported in humans in RPS6KB1, single nucleotide variants (SNVs) causing missense mutations have been identified, which have not been assessed for their impact on protein function. These S6K1 mutations have the potential to influence disease progression and treatment response. We mined the Simon Simplex Collection (SSC) and SPARK autism database to find inherited SNVs in S6K1 and characterized the effect of two missense SNVs, Asp14Asn (allele frequency = 0.03282%) and Glu44Gln (allele frequency = 0.0008244%), on S6K1 function in HEK293, human ES cells and primary neurons. Expressing Asp14Asn in HEK293 cells resulted in increased basal phosphorylation of downstream targets of S6K1 and increased de novo translation. This variant also showed blunted response to the specific S6K1 inhibitor, FS-115. In human embryonic cell line Shef4, Asp14Asn enhanced spontaneous neural fate specification in the absence of differentiating growth factors. In addition to enhanced translation, neurons expressing Asp14Asn exhibited impaired dendritic arborization and increased levels of phosphorylated ERK 1/2. Finally, in the SSC families tracked, Asp14Asn segregated with lower IQ scores when found in the autistic individual rather than the unaffected sibling. The Glu44Gln mutation showed a milder, but opposite phenotype in HEK cells as compared to Asp14Asn. Although the Glu44Gln mutation displayed increased neuronal translation, it had no impact on neuronal morphology. Our results provide the first characterization of naturally occurring human S6K1 variants on cognitive phenotype, neuronal morphology and maturation, underscoring again the importance of translation control in neural development and plasticity.

Obesity is an immunometabolic disease associated with chronic inflammation and the dysregulation of pro- and anti-inflammatory cytokines. One hallmark of obesity is reduced concentrations of the anti-inflammatory adipokine, adiponectin. Pharmacologic doses of niacin produce multiple metabolic benefits, including attenuating high-fat diet (HFD)-induced adipose tissue inflammation and increasing adiponectin concentrations. To determine if adiponectin mediates the anti-inflammatory effects of niacin, male C57BL/6J (WT) and adiponectin null (Adipoq-/-) mice were maintained on a low-fat diet (LFD) or HFD for 6 weeks, before being administered either vehicle or niacin (360 mg/kg/day) for 5 weeks. HFD-fed mice had increased expression of genes associated with macrophage recruitment (Ccl2) and number (Cd68), and increased crown-like structure (CLS) number in adipose tissue. While niacin attenuated Ccl2 expression, there were no effects on Cd68 or CLS number. The absence of adiponectin did not hinder the ability of niacin to reduce Ccl2 expression. HFD feeding increased gene expression of inflammatory markers in the adipose tissue of WT and Adipoq-/- mice. While niacin tended to decrease the expression of inflammatory markers in WT mice, niacin increased their expression in HFD-fed Adipoq-/- mice. Therefore, our results indicate that the absence of adiponectin alters the effects of niacin on markers of adipose tissue inflammation in HFD-fed mice, suggesting that the effects of niacin on tissue cytokines may involve adiponectin.

Background. In recent years, there has been an explosion in the number of technical and medical diagnostic platforms being developed. This has greatly improved our ability to more accurately, and more comprehensively, explore and characterize human biological systems on the individual level. Large quantities of biomedical data are now being generated and archived in many separate research and clinical activities, but there exists a paucity of studies that integrate the areas of clinical neuropsychiatry, personal genomics and brain-machine interfaces. Methods. A single person with severe mental illness was implanted with the Medtronic Reclaim(®) Deep Brain Stimulation (DBS) Therapy device for Obsessive Compulsive Disorder (OCD), targeting his nucleus accumbens/anterior limb of the internal capsule. Programming of the device and psychiatric assessments occurred in an outpatient setting for over two years. His genome was sequenced and variants were detected in the Illumina Whole Genome Sequencing Clinical Laboratory Improvement Amendments (CLIA)-certified laboratory. Results. We report here the detailed phenotypic characterization, clinical-grade whole genome sequencing (WGS), and two-year outcome of a man with severe OCD treated with DBS. Since implantation, this man has reported steady improvement, highlighted by a steady decline in his Yale-Brown Obsessive Compulsive Scale (YBOCS) score from ∼38 to a score of ∼25. A rechargeable Activa RC neurostimulator battery has been of major benefit in terms of facilitating a degree of stability and control over the stimulation. His psychiatric symptoms reliably worsen within hours of the battery becoming depleted, thus providing confirmatory evidence for the efficacy of DBS for OCD in this person. WGS revealed that he is a heterozygote for the p.Val66Met variant in BDNF, encoding a member of the nerve growth factor family, and which has been found to predispose carriers to various psychiatric illnesses. He carries the p.Glu429Ala allele in methylenetetrahydrofolate reductase (MTHFR) and the p.Asp7Asn allele in ChAT, encoding choline O-acetyltransferase, with both alleles having been shown to confer an elevated susceptibility to psychoses. We have found thousands of other variants in his genome, including pharmacogenetic and copy number variants. This information has been archived and offered to this person alongside the clinical sequencing data, so that he and others can re-analyze his genome for years to come. Conclusions. To our knowledge, this is the first study in the clinical neurosciences that integrates detailed neuropsychiatric phenotyping, deep brain stimulation for OCD and clinical-grade WGS with management of genetic results in the medical treatment of one person with severe mental illness. We offer this as an example of precision medicine in neuropsychiatry including brain-implantable devices and genomics-guided preventive health care.

Naa10 is an Nα -terminal acetyltransferase that, in a complex with its auxiliary subunit Naa15, co-translationally acetylates the α-amino group of newly synthetized proteins as they emerge from the ribosome. Roughly 40-50% of the human proteome is acetylated by Naa10, rendering this an enzyme one of the most broad substrate ranges known. Recently, we reported an X-linked disorder of infancy, Ogden syndrome, in two families harbouring a c.109 T > C (p.Ser37Pro) variant in NAA10. In the present study we performed in-depth characterization of a yeast model of Ogden syndrome. Stress tests and proteomic analyses suggest that the S37P mutation disrupts Naa10 function and reduces cellular fitness during heat shock, possibly owing to dysregulation of chaperone expression and accumulation. Microarray and RNA-seq revealed a pseudo-diploid gene expression profile in ΔNaa10 cells, probably responsible for a mating defect. In conclusion, the data presented here further support the disruptive nature of the S37P/Ogden mutation and identify affected cellular processes potentially contributing to the severe phenotype seen in Ogden syndrome. Data are available via GEO under identifier GSE86482 or with ProteomeXchange under identifier PXD004923. © 2016 The Authors. Yeast published by John Wiley & Sons, Ltd.

Autism arises in high and low-risk families. De novo mutation contributes to autism incidence in low-risk families as there is a higher incidence in the affected of the simplex families than in their unaffected siblings. But the extent of contribution in low-risk families cannot be determined solely from simplex families as they are a mixture of low and high-risk. The rate of de novo mutation in nearly pure populations of high-risk families, the multiplex families, has not previously been rigorously determined. Moreover, rates of de novo mutation have been underestimated from studies based on low resolution microarrays and whole exome sequencing. Here we report on findings from whole genome sequence (WGS) of both simplex families from the Simons Simplex Collection (SSC) and multiplex families from the Autism Genetic Resource Exchange (AGRE). After removing the multiplex samples with excessive cell-line genetic drift, we find that the contribution of de novo mutation in multiplex is significantly smaller than the contribution in simplex. We use WGS to provide high resolution CNV profiles and to analyze more than coding regions, and revise upward the rate in simplex autism due to an excess of de novo events targeting introns. Based on this study, we now estimate that de novo events contribute to 52-67% of cases of autism arising from low risk families, and 30-39% of cases of all autism.

Dietary protein restriction and dietary methionine restriction (MR) produce a comparable series of behavioral, physiological, biochemical, and transcriptional responses. Both dietary regimens produce a similar reduction in intake of sulfur amino acids (e.g., methionine and cystine), and both diets increase expression and release of hepatic FGF21. Given that FGF21 is an essential mediator of the metabolic phenotype produced by both diets, an important unresolved question is whether dietary protein restriction represents de facto methionine restriction. Using diets formulated from either casein or soy protein with matched reductions in sulfur amino acids, we compared the ability of the respective diets to recapitulate the metabolic phenotype produced by methionine restriction using elemental diets. Although the soy-based control diets supported faster growth compared to casein-based control diets, casein-based protein restriction and soy-based protein restriction produced comparable reductions in body weight and fat deposition, and similar increases in energy intake, energy expenditure, and water intake. In addition, the prototypical effects of dietary MR on hepatic and adipose tissue target genes were similarly regulated by casein- and soy-based protein restriction. The present findings support the feasibility of using restricted intake of diets from various protein sources to produce therapeutically effective implementation of dietary methionine restriction.

Lung adenocarcinoma (LUAD) is a common type of lung cancer with a high risk of metastasis, but the exact molecular mechanisms of metastasis are not yet understood.

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-seq data in GTEx v8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased WGS data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Autism spectrum disorders (ASDs) are a group of developmental disabilities that affect social interaction and communication and are characterized by repetitive behaviors. There is now a large body of evidence that suggests a complex role of genetics in ASDs, in which many different loci are involved. Although many current population-scale genomic studies have been demonstrably fruitful, these studies generally focus on analyzing a limited part of the genome or use a limited set of bioinformatics tools. These limitations preclude the analysis of genome-wide perturbations that may contribute to the development and severity of ASD-related phenotypes. To overcome these limitations, we have developed and utilized an integrative clinical and bioinformatics pipeline for generating a more complete and reliable set of genomic variants for downstream analyses. Our study focuses on the analysis of three simplex autism families consisting of one affected child, unaffected parents, and one unaffected sibling. All members were clinically evaluated and widely phenotyped. Genotyping arrays and whole-genome sequencing were performed on each member, and the resulting sequencing data were analyzed using a variety of available bioinformatics tools. We searched for rare variants of putative functional impact that were found to be segregating according to de novo, autosomal recessive, X-linked, mitochondrial, and compound heterozygote transmission models. The resulting candidate variants included three small heterozygous copy-number variations (CNVs), a rare heterozygous de novo nonsense mutation in MYBBP1A located within exon 1, and a novel de novo missense variant in LAMB3. Our work demonstrates how more comprehensive analyses that include rich clinical data and whole-genome sequencing data can generate reliable results for use in downstream investigations.