Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

A sterol-binding protein integrates endosomal lipid metabolism with TOR signaling and nitrogen sensing.

Kes1, and other oxysterol-binding protein superfamily members, are involved in membrane and lipid trafficking through trans-Golgi network (TGN) and endosomal systems. We demonstrate that Kes1 represents a sterol-regulated antagonist of TGN/endosomal phosphatidylinositol-4-phosphate signaling. This regulation modulates TOR activation by amino acids and dampens gene expression driven by Gcn4, the primary transcriptional activator of the general amino acid control regulon. Kes1-mediated repression of Gcn4 transcription factor activity is characterized by nonproductive Gcn4 binding to its target sequences, involves TGN/endosome-derived sphingolipid signaling, and requires activity of the cyclin-dependent kinase 8 (CDK8) module of the enigmatic "large Mediator" complex. These data describe a pathway by which Kes1 integrates lipid metabolism with TORC1 signaling and nitrogen sensing.

Gene copy number variation and pediatric mental health/neurodevelopment in a general population.

We assessed the relationship of gene copy number variation (CNV) in mental health/neurodevelopmental traits and diagnoses, physical health and cognition in a community sample of 7100 unrelated children and youth of European or East Asian ancestry (Spit for Science). Clinically significant or susceptibility CNVs were present in 3.9% of participants and were associated with elevated scores on a continuous measure of attention-deficit/hyperactivity disorder (ADHD) traits (P = 5.0 × 10-3), longer response inhibition (a cognitive deficit found in several mental health and neurodevelopmental disorders; P = 1.0 × 10-2) and increased prevalence of mental health diagnoses (P = 1.9 × 10-6, odds ratio: 3.09), specifically ADHD, autism spectrum disorder anxiety and learning problems/learning disorder (P's < 0.01). There was an increased burden of rare deletions in gene-sets related to brain function or expression in brain associated with more ADHD traits. With the current mental health crisis, our data established a baseline for delineating genetic contributors in pediatric-onset conditions.

Oxidized phospholipids are proinflammatory and proatherogenic in hypercholesterolaemic mice.

Oxidized phospholipids (OxPL) are ubiquitous, are formed in many inflammatory tissues, including atherosclerotic lesions, and frequently mediate proinflammatory changes 1 . Because OxPL are mostly the products of non-enzymatic lipid peroxidation, mechanisms to specifically neutralize them are unavailable and their roles in vivo are largely unknown. We previously cloned the IgM natural antibody E06, which binds to the phosphocholine headgroup of OxPL, and blocks the uptake of oxidized low-density lipoprotein (OxLDL) by macrophages and inhibits the proinflammatory properties of OxPL2-4. Here, to determine the role of OxPL in vivo in the context of atherogenesis, we generated transgenic mice in the Ldlr-/- background that expressed a single-chain variable fragment of E06 (E06-scFv) using the Apoe promoter. E06-scFv was secreted into the plasma from the liver and macrophages, and achieved sufficient plasma levels to inhibit in vivo macrophage uptake of OxLDL and to prevent OxPL-induced inflammatory signalling. Compared to Ldlr-/- mice, Ldlr -/- E06-scFv mice had 57-28% less atherosclerosis after 4, 7 and even 12 months of 1% high-cholesterol diet. Echocardiographic and histologic evaluation of the aortic valves demonstrated that E06-scFv ameliorated the development of aortic valve gradients and decreased aortic valve calcification. Both cholesterol accumulation and in vivo uptake of OxLDL were decreased in peritoneal macrophages, and both peritoneal and aortic macrophages had a decreased inflammatory phenotype. Serum amyloid A was decreased by 32%, indicating decreased systemic inflammation, and hepatic steatosis and inflammation were also decreased. Finally, the E06-scFv prolonged life as measured over 15 months. Because the E06-scFv lacks the functional effects of an intact antibody other than the ability to bind OxPL and inhibit OxLDL uptake in macrophages, these data support a major proatherogenic role of OxLDL and demonstrate that OxPL are proinflammatory and proatherogenic, which E06 counteracts in vivo. These studies suggest that therapies inactivating OxPL may be beneficial for reducing generalized inflammation, including the progression of atherosclerosis, aortic stenosis and hepatic steatosis.

Publisher Correction: Oxidized phospholipids are proinflammatory and proatherogenic in hypercholesterolaemic mice.

In this Letter, affiliation number 1 was originally missing from the HTML; the affiliations were missing for author Ming-Yow Hung in the HTML; and the Fig. 4 legend erroneously referred to panels a-h, instead of a-g. These errors have been corrected online.

Rare copy number variation in posttraumatic stress disorder.

Posttraumatic stress disorder (PTSD) is a heritable (h2 = 24-71%) psychiatric illness. Copy number variation (CNV) is a form of rare genetic variation that has been implicated in the etiology of psychiatric disorders, but no large-scale investigation of CNV in PTSD has been performed. We present an association study of CNV burden and PTSD symptoms in a sample of 114,383 participants (13,036 cases and 101,347 controls) of European ancestry. CNVs were called using two calling algorithms and intersected to a consensus set. Quality control was performed to remove strong outlier samples. CNVs were examined for association with PTSD within each cohort using linear or logistic regression analysis adjusted for population structure and CNV quality metrics, then inverse variance weighted meta-analyzed across cohorts. We examined the genome-wide total span of CNVs, enrichment of CNVs within specified gene-sets, and CNVs overlapping individual genes and implicated neurodevelopmental regions. The total distance covered by deletions crossing over known neurodevelopmental CNV regions was significant (beta = 0.029, SE = 0.005, P = 6.3 × 10-8). The genome-wide neurodevelopmental CNV burden identified explains 0.034% of the variation in PTSD symptoms. The 15q11.2 BP1-BP2 microdeletion region was significantly associated with PTSD (beta = 0.0206, SE = 0.0056, P = 0.0002). No individual significant genes interrupted by CNV were identified. 22 gene pathways related to the function of the nervous system and brain were significant in pathway analysis (FDR q < 0.05), but these associations were not significant once NDD regions were removed. A larger sample size, better detection methods, and annotated resources of CNV are needed to explore this relationship further.

Direct Detection of Glutathione Biosynthesis, Conjugation, Depletion and Recovery in Intact Hepatoma Cells.

Nuclear magnetic resonance (NMR) spectroscopy was used to monitor glutathione metabolism in alginate-encapsulated JM-1 hepatoma cells perfused with growth media containing [3,3'-13C2]-cystine. After 20 h of perfusion with labeled medium, the 13C NMR spectrum is dominated by the signal from the 13C-labeled glutathione. Once 13C-labeled, the high intensity of the glutathione resonance allows the acquisition of subsequent spectra in 1.2 min intervals. At this temporal resolution, the detailed kinetics of glutathione metabolism can be monitored as the thiol alkylating agent monobromobimane (mBBr) is added to the perfusate. The addition of a bolus dose of mBBr results in rapid diminution of the resonance for 13C-labeled glutathione due to a loss of this metabolite through alkylation by mBBr. As the glutathione resonance decreases, a new resonance due to the production of intracellular glutathione-bimane conjugate is detectable. After clearance of the mBBr dose from the cells, intracellular glutathione repletion is then observed by a restoration of the 13C-glutathione signal along with wash-out of the conjugate. These data demonstrate that standard NMR techniques can directly monitor intracellular processes such as glutathione depletion with a time resolution of approximately < 2 min.

A Comprehensive Workflow for Read Depth-Based Identification of Copy-Number Variation from Whole-Genome Sequence Data.

A remaining hurdle to whole-genome sequencing (WGS) becoming a first-tier genetic test has been accurate detection of copy-number variations (CNVs). Here, we used several datasets to empirically develop a detailed workflow for identifying germline CNVs >1 kb from short-read WGS data using read depth-based algorithms. Our workflow is comprehensive in that it addresses all stages of the CNV-detection process, including DNA library preparation, sequencing, quality control, reference mapping, and computational CNV identification. We used our workflow to detect rare, genic CNVs in individuals with autism spectrum disorder (ASD), and 120/120 such CNVs tested using orthogonal methods were successfully confirmed. We also identified 71 putative genic de novo CNVs in this cohort, which had a confirmation rate of 70%; the remainder were incorrectly identified as de novo due to false positives in the proband (7%) or parental false negatives (23%). In individuals with an ASD diagnosis in which both microarray and WGS experiments were performed, our workflow detected all clinically relevant CNVs identified by microarrays, as well as additional potentially pathogenic CNVs < 20 kb. Thus, CNVs of clinical relevance can be discovered from WGS with a detection rate exceeding microarrays, positioning WGS as a single assay for genetic variation detection.

A copy number variation map of the human genome.

A major contribution to the genome variability among individuals comes from deletions and duplications - collectively termed copy number variations (CNVs) - which alter the diploid status of DNA. These alterations may have no phenotypic effect, account for adaptive traits or can underlie disease. We have compiled published high-quality data on healthy individuals of various ethnicities to construct an updated CNV map of the human genome. Depending on the level of stringency of the map, we estimated that 4.8-9.5% of the genome contributes to CNV and found approximately 100 genes that can be completely deleted without producing apparent phenotypic consequences. This map will aid the interpretation of new CNV findings for both clinical and research applications.

De novo and rare inherited copy-number variations in the hemiplegic form of cerebral palsy.

PurposeHemiplegia is a subtype of cerebral palsy (CP) in which one side of the body is affected. Our earlier study of unselected children with CP demonstrated de novo and clinically relevant rare inherited genomic copy-number variations (CNVs) in 9.6% of participants. Here, we examined the prevalence and types of CNVs specifically in hemiplegic CP.MethodsWe genotyped 97 unrelated probands with hemiplegic CP and their parents. We compared their CNVs to those of 10,851 population controls, in order to identify rare CNVs (<0.1% frequency) that might be relevant to CP. We also sequenced exomes of "CNV-positive" trios.ResultsWe detected de novo CNVs and/or sex chromosome abnormalities in 7/97 (7.2%) of probands, impacting important developmental genes such as GRIK2, LAMA1, DMD, PTPRM, and DIP2C. In 18/97 individuals (18.6%), rare inherited CNVs were found, affecting loci associated with known genomic disorders (17p12, 22q11.21) or involving genes linked to neurodevelopmental disorders.ConclusionWe found an increased rate of de novo CNVs in the hemiplegic CP subtype (7.2%) compared to controls (1%). This result is similar to that for an unselected CP group. Combined with rare inherited CNVs, the genomic data impacts the understanding of the potential etiology of hemiplegic CP in 23/97 (23.7%) of participants.

The Personal Genome Project Canada: findings from whole genome sequences of the inaugural 56 participants.

BACKGROUND: The Personal Genome Project Canada is a comprehensive public data resource that integrates whole genome sequencing data and health information. We describe genomic variation identified in the initial recruitment cohort of 56 volunteers. METHODS: Volunteers were screened for eligibility and provided informed consent for open data sharing. Using blood DNA, we performed whole genome sequencing and identified all possible classes of DNA variants. A genetic counsellor explained the implication of the results to each participant. RESULTS: Whole genome sequencing of the first 56 participants identified 207 662 805 sequence variants and 27 494 copy number variations. We analyzed a prioritized disease-associated data set (n = 1606 variants) according to standardized guidelines, and interpreted 19 variants in 14 participants (25%) as having obvious health implications. Six of these variants (e.g., in BRCA1 or mosaic loss of an X chromosome) were pathogenic or likely pathogenic. Seven were risk factors for cancer, cardiovascular or neurobehavioural conditions. Four other variants - associated with cancer, cardiac or neurodegenerative phenotypes - remained of uncertain significance because of discrepancies among databases. We also identified a large structural chromosome aberration and a likely pathogenic mitochondrial variant. There were 172 recessive disease alleles (e.g., 5 individuals carried mutations for cystic fibrosis). Pharmacogenomics analyses revealed another 3.9 potentially relevant genotypes per individual. INTERPRETATION: Our analyses identified a spectrum of genetic variants with potential health impact in 25% of participants. When also considering recessive alleles and variants with potential pharmacologic relevance, all 56 participants had medically relevant findings. Although access is mostly limited to research, whole genome sequencing can provide specific and novel information with the potential of major impact for health care.

Association of IMMP2L deletions with autism spectrum disorder: A trio family study and meta-analysis.

IMMP2L, the gene encoding the inner mitochondrial membrane peptidase subunit 2-like protein, has been reported as a candidate gene for Tourette syndrome, autism spectrum disorder (ASD) and additional neurodevelopmental disorders. Here we genotyped 100 trio families with an index proband with autism spectrum disorder in Han Chinese population and found three cases with rare exonic IMMP2L deletions. We have conducted a comprehensive meta-analysis to quantify the association of IMMP2L deletions with ASD using 5,568 cases and 10,279 controls. While the IMMP2L deletions carried non-recurrent breakpoints, in contrast to previous reports, our meta-analysis found no evidence of association (P > 0.05) between IMMP2L deletions and ASD. We also observed common exonic deletions impacting IMMP2L in a separate control (5,971 samples) cohort where subjects were screened for psychiatric conditions. This is the first systematic review and meta-analysis regarding the effect of IMMP2L deletions on ASD, but further investigations in different populations, especially Chinese population may be still needed to confirm our results.

Ribosomal protein-Mdm2-p53 pathway coordinates nutrient stress with lipid metabolism by regulating MCD and promoting fatty acid oxidation.

The tumor suppressor p53 has recently been shown to regulate energy metabolism through multiple mechanisms. However, the in vivo signaling pathways related to p53-mediated metabolic regulation remain largely uncharacterized. By using mice bearing a single amino acid substitution at cysteine residue 305 of mouse double minute 2 (Mdm2(C305F)), which renders Mdm2 deficient in binding ribosomal proteins (RPs) RPL11 and RPL5, we show that the RP-Mdm2-p53 signaling pathway is critical for sensing nutrient deprivation and maintaining liver lipid homeostasis. Although the Mdm2(C305F) mutation does not significantly affect growth and development in mice, this mutation promotes fat accumulation under normal feeding conditions and hepatosteatosis under acute fasting conditions. We show that nutrient deprivation inhibits rRNA biosynthesis, increases RP-Mdm2 interaction, and induces p53-mediated transactivation of malonyl-CoA decarboxylase (MCD), which catalyzes the degradation of malonyl-CoA to acetyl-CoA, thus modulating lipid partitioning. Fasted Mdm2(C305F) mice demonstrate attenuated MCD induction and enhanced malonyl-CoA accumulation in addition to decreased oxidative respiration and increased fatty acid accumulation in the liver. Thus, the RP-Mdm2-p53 pathway appears to function as an endogenous sensor responsible for stimulating fatty acid oxidation in response to nutrient depletion.

G-protein coupled receptor-mediated nutrient sensing and developmental control in Aspergillus nidulans.

Nutrient sensing and utilisation are fundamental for all life forms. As heterotrophs, fungi have evolved a diverse range of mechanisms for sensing and taking up various nutrients. Despite its importance, only a limited number of nutrient receptors and their corresponding ligands have been identified in fungi. G-protein coupled receptors (GPCRs) are the largest family of transmembrane receptors. The Aspergillus nidulans genome encodes 16 putative GPCRs, but only a few have been functionally characterised. Our previous study showed the increased expression of an uncharacterised putative GPCR, gprH, during carbon starvation. GprH appears conserved throughout numerous filamentous fungi. Here, we reveal that GprH is a putative receptor involved in glucose and tryptophan sensing. The absence of GprH results in a reduction in cAMP levels and PKA activity upon adding glucose or tryptophan to starved cells. GprH is pre-formed in conidia and is increasingly active during carbon starvation, where it plays a role in glucose uptake and the recovery of hyphal growth. GprH also represses sexual development under conditions favouring sexual fruiting and during carbon starvation in submerged cultures. In summary, the GprH nutrient-sensing system functions upstream of the cAMP-PKA pathway, influences primary metabolism and hyphal growth, while represses sexual development in A. nidulans.

Disruption of the ASTN2/TRIM32 locus at 9q33.1 is a risk factor in males for autism spectrum disorders, ADHD and other neurodevelopmental phenotypes.

Rare copy number variants (CNVs) disrupting ASTN2 or both ASTN2 and TRIM32 have been reported at 9q33.1 by genome-wide studies in a few individuals with neurodevelopmental disorders (NDDs). The vertebrate-specific astrotactins, ASTN2 and its paralog ASTN1, have key roles in glial-guided neuronal migration during brain development. To determine the prevalence of astrotactin mutations and delineate their associated phenotypic spectrum, we screened ASTN2/TRIM32 and ASTN1 (1q25.2) for exonic CNVs in clinical microarray data from 89 985 individuals across 10 sites, including 64 114 NDD subjects. In this clinical dataset, we identified 46 deletions and 12 duplications affecting ASTN2. Deletions of ASTN1 were much rarer. Deletions near the 3' terminus of ASTN2, which would disrupt all transcript isoforms (a subset of these deletions also included TRIM32), were significantly enriched in the NDD subjects (P = 0.002) compared with 44 085 population-based controls. Frequent phenotypes observed in individuals with such deletions include autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), speech delay, anxiety and obsessive compulsive disorder (OCD). The 3'-terminal ASTN2 deletions were significantly enriched compared with controls in males with NDDs, but not in females. Upon quantifying ASTN2 human brain RNA, we observed shorter isoforms expressed from an alternative transcription start site of recent evolutionary origin near the 3' end. Spatiotemporal expression profiling in the human brain revealed consistently high ASTN1 expression while ASTN2 expression peaked in the early embryonic neocortex and postnatal cerebellar cortex. Our findings shed new light on the role of the astrotactins in psychopathology and their interplay in human neurodevelopment.

Origins and functional impact of copy number variation in the human genome.

Structural variations of DNA greater than 1 kilobase in size account for most bases that vary among human genomes, but are still relatively under-ascertained. Here we use tiling oligonucleotide microarrays, comprising 42 million probes, to generate a comprehensive map of 11,700 copy number variations (CNVs) greater than 443 base pairs, of which most (8,599) have been validated independently. For 4,978 of these CNVs, we generated reference genotypes from 450 individuals of European, African or East Asian ancestry. The predominant mutational mechanisms differ among CNV size classes. Retrotransposition has duplicated and inserted some coding and non-coding DNA segments randomly around the genome. Furthermore, by correlation with known trait-associated single nucleotide polymorphisms (SNPs), we identified 30 loci with CNVs that are candidates for influencing disease susceptibility. Despite this, having assessed the completeness of our map and the patterns of linkage disequilibrium between CNVs and SNPs, we conclude that, for complex traits, the heritability void left by genome-wide association studies will not be accounted for by common CNVs.

The Database of Genomic Variants: a curated collection of structural variation in the human genome.

Over the past decade, the Database of Genomic Variants (DGV; http://dgv.tcag.ca/) has provided a publicly accessible, comprehensive curated catalogue of structural variation (SV) found in the genomes of control individuals from worldwide populations. Here, we describe updates and new features, which have expanded the utility of DGV for both the basic research and clinical diagnostic communities. The current version of DGV consists of 55 published studies, comprising >2.5 million entries identified in >22,300 genomes. Studies included in DGV are selected from the accessioned data sets in the archival SV databases dbVar (NCBI) and DGVa (EBI), and then further curated for accuracy and validity. The core visualization tool (gbrowse) has been upgraded with additional functions to facilitate data analysis and comparison, and a new query tool has been developed to provide flexible and interactive access to the data. The content from DGV is regularly incorporated into other large-scale genome reference databases and represents a standard data resource for new product and database development, in particular for copy number variation testing in clinical labs. The accurate cataloguing of variants in DGV will continue to enable medical genetics and genome sequencing research.

A high-resolution copy-number variation resource for clinical and population genetics.

PURPOSE: Chromosomal microarray analysis to assess copy-number variation has become a first-tier genetic diagnostic test for individuals with unexplained neurodevelopmental disorders or multiple congenital anomalies. More than 100 cytogenetic laboratories worldwide use the new ultra-high resolution Affymetrix CytoScan-HD array to genotype hundreds of thousands of samples per year. Our aim was to develop a copy-number variation resource from a new population sample that would enable more accurate interpretation of clinical genetics data on this microarray platform and others. METHODS: Genotyping of 1,000 adult volunteers who are broadly representative of the Ontario population (as obtained from the Ontario Population Genomics Platform) was performed with the CytoScan-HD microarray system, which has 2.7 million probes. Four independent algorithms were applied to detect copy-number variations. Reproducibility and validation metrics were quantified using sample replicates and quantitative-polymerase chain reaction, respectively. RESULTS: DNA from 873 individuals passed quality control and we identified 71,178 copy-number variations (81 copy-number variations/individual); 9.8% (6,984) of these copy-number variations were previously unreported. After applying three layers of filtering criteria, from our highest confidence copy-number variation data set we obtained >95% reproducibility and >90% validation rates (73% of these copy-number variations overlapped at least one gene). CONCLUSION: The genotype data and annotated copy-number variations for this largely Caucasian population will represent a valuable public resource enabling clinical genetics research and diagnostics.

Biocompatibility of Tygon® tubing in microfluidic cell culture.

Growth of the MDA-MB-231 breast cancer cell line in microfluidic channels was inhibited when culture media was delivered to the channels via microbore Tygon® tubing. Culture media incubated within this tubing also inhibited growth of these cells in conventional 96-well plates. These detrimental effects were not due to depletion of critical nutrients due to adsorption of media components onto the tubing surface. A pH change was also ruled out as a cause. Nuclear magnetic resonance spectroscopy of the cell growth media before and after incubation in the tubing confirmed no detectable loss of media components but did detect the presence of additional unidentified signals in the aliphatic region of the spectrum. These results indicate leaching of a chemical species from microbore Tygon® tubing that can affect cell growth in microfluidic devices.

Genome assembly comparison identifies structural variants in the human genome.

Numerous types of DNA variation exist, ranging from SNPs to larger structural alterations such as copy number variants (CNVs) and inversions. Alignment of DNA sequence from different sources has been used to identify SNPs and intermediate-sized variants (ISVs). However, only a small proportion of total heterogeneity is characterized, and little is known of the characteristics of most smaller-sized (<50 kb) variants. Here we show that genome assembly comparison is a robust approach for identification of all classes of genetic variation. Through comparison of two human assemblies (Celera's R27c compilation and the Build 35 reference sequence), we identified megabases of sequence (in the form of 13,534 putative non-SNP events) that were absent, inverted or polymorphic in one assembly. Database comparison and laboratory experimentation further demonstrated overlap or validation for 240 variable regions and confirmed >1.5 million SNPs. Some differences were simple insertions and deletions, but in regions containing CNVs, segmental duplication and repetitive DNA, they were more complex. Our results uncover substantial undescribed variation in humans, highlighting the need for comprehensive annotation strategies to fully interpret genome scanning and personalized sequencing projects.