Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Characterization of human embryonic stem cell lines by the International Stem Cell Initiative

The International Stem Cell Initiative characterized 59 human embryonic stem cell lines from 17 laboratories worldwide. Despite diverse genotypes and different techniques used for derivation and maintenance, all lines exhibited similar expression patterns for several markers of human embryonic stem cells. They expressed the glycolipid antigens SSEA3 and SSEA4, the keratan sulfate antigens TRA-1-60, TRA-1-81, GCTM2 and GCT343, and the protein antigens CD9, Thy1 (also known as CD90), tissue-nonspecific alkaline phosphatase and class 1 HLA, as well as the strongly developmentally regulated genes NANOG, POU5F1 (formerly known as OCT4), TDGF1, DNMT3B, GABRB3 and GDF3. Nevertheless, the lines were not identical: differences in expression of several lineage markers were evident, and several imprinted genes showed generally similar allele-specific expression patterns, but some gene-dependent variation was observed. Also, some female lines expressed readily detectable levels of XIST whereas others did not. No significant contamination of the lines with mycoplasma, bacteria or cytopathic viruses was detected.

Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects

Copy number variants (CNVs) have been strongly implicated in the genetic etiology of schizophrenia (SCZ). However, genome-wide investigation of the contribution of CNV to risk has been hampered by limited sample sizes. We sought to address this obstacle by applying a centralized analysis pipeline to a SCZ cohort of 21,094 cases and 20,227 controls. A global enrichment of CNV burden was observed in cases (odds ratio (OR) = 1.11, P = 5.7 × 10-15), which persisted after excluding loci implicated in previous studies (OR = 1.07, P = 1.7 × 10-6). CNV burden was enriched for genes associated with synaptic function (OR = 1.68, P = 2.8 × 10-11) and neurobehavioral phenotypes in mouse (OR = 1.18, P = 7.3 × 10-5). Genome-wide significant evidence was obtained for eight loci, including 1q21.1, 2p16.3 (NRXN1), 3q29, 7q11.2, 15q13.3, distal 16p11.2, proximal 16p11.2 and 22q11.2. Suggestive support was found for eight additional candidate susceptibility and protective loci, which consisted predominantly of CNVs mediated by nonallelic homologous recombination.

Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa

Characterized primarily by a low body-mass index, anorexia nervosa is a complex and serious illness1, affecting 0.9-4% of women and 0.3% of men2-4, with twin-based heritability estimates of 50-60%5. Mortality rates are higher than those in other psychiatric disorders6, and outcomes are unacceptably poor7. Here we combine data from the Anorexia Nervosa Genetics Initiative (ANGI)8,9 and the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED) and conduct a genome-wide association study of 16,992 cases of anorexia nervosa and 55,525 controls, identifying eight significant loci. The genetic architecture of anorexia nervosa mirrors its clinical presentation, showing significant genetic correlations with psychiatric disorders, physical activity, and metabolic (including glycemic), lipid and anthropometric traits, independent of the effects of common variants associated with body-mass index. These results further encourage a reconceptualization of anorexia nervosa as a metabo-psychiatric disorder. Elucidating the metabolic component is a critical direction for future research, and paying attention to both psychiatric and metabolic components may be key to improving outcomes.

Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa

OBJECTIVE

The authors conducted a genome-wide association study of anorexia nervosa and calculated genetic correlations with a series of psychiatric, educational, and metabolic phenotypes.

METHOD

Following uniform quality control and imputation procedures using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, the authors performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression was used to calculate genome-wide common variant heritability (single-nucleotide polymorphism [SNP]-based heritability [h2SNP]), partitioned heritability, and genetic correlations (rg) between anorexia nervosa and 159 other phenotypes.

RESULTS

Results were obtained for 10,641,224 SNPs and insertion-deletion variants with minor allele frequencies >1% and imputation quality scores >0.6. The h2SNP of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. The authors identified one genome-wide significant locus on chromosome 12 (rs4622308) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high-density lipoprotein cholesterol, and significant negative genetic correlations were observed between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes.

CONCLUSIONS

Anorexia nervosa is a complex heritable phenotype for which this study has uncovered the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. The study results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Copy number variation in schizophrenia in Sweden

Schizophrenia (SCZ) is a highly heritable neuropsychiatric disorder of complex genetic etiology. Previous genome-wide surveys have revealed a greater burden of large, rare copy number variations (CNVs) in SCZ cases and identified multiple rare recurrent CNVs that increase risk of SCZ although with incomplete penetrance and pleiotropic effects. Identification of additional recurrent CNVs and biological pathways enriched for SCZ CNVs requires greater sample sizes. We conducted a genome-wide survey for CNVs associated with SCZ using a Swedish national sample (4719 cases and 5917 controls). High-confidence CNV calls were generated using genotyping array intensity data, and their effect on risk of SCZ was measured. Our data confirm increased burden of large, rare CNVs in SCZ cases as well as significant associations for recurrent 16p11.2 duplications, 22q11.2 deletions and 3q29 deletions. We report a novel association for 17q12 duplications (odds ratio=4.16, P=0.018), previously associated with autism and mental retardation but not SCZ. Intriguingly, gene set association analyses implicate biological pathways previously associated with SCZ through common variation and exome sequencing (calcium channel signaling and binding partners of the fragile X mental retardation protein). We found significantly increased burden of the largest CNVs (>500 kb) in genes present in the postsynaptic density, in genomic regions implicated via SCZ genome-wide association studies and in gene products localized to mitochondria and cytoplasm. Our findings suggest that multiple lines of genomic inquiry--genome-wide screens for CNVs, common variation and exonic variation--are converging on similar sets of pathways and/or genes.

The recombinational anatomy of a mouse chromosome

Among mammals, genetic recombination occurs at highly delimited sites known as recombination hotspots. They are typically 1–2 kb long and vary as much as a 1,000-fold or more in recombination activity. Although much is known about the molecular details of the recombination process itself, the factors determining the location and relative activity of hotspots are poorly understood. To further our understanding, we have collected and mapped the locations of 5,472 crossover events along mouse Chromosome 1 arising in 6,028 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. Crossovers were mapped to a minimum resolution of 225 kb, and those in the telomere-proximal 24.7 Mb were further mapped to resolve individual hotspots. Recombination rates were evolutionarily conserved on a regional scale, but not at the local level. There was a clear negative-exponential relationship between the relative activity and abundance of hotspot activity classes, such that a small number of the most active hotspots account for the majority of recombination. Females had 1.2× higher overall recombination than males did, although the sex ratio showed considerable regional variation. Locally, entirely sex-specific hotspots were rare. The initiation of recombination at the most active hotspot was regulated independently on the two parental chromatids, and analysis of reciprocal crosses indicated that parental imprinting has subtle effects on recombination rates. It appears that the regulation of mammalian recombination is a complex, dynamic process involving multiple factors reflecting species, sex, individual variation within species, and the properties of individual hotspots.

In most eukaryotic organisms, recombination—the exchange of genetic information between homologous chromosomes—ensures the proper recognition and segregation of chromosomes during meiosis. Recombination events in mammals are not randomly positioned along the chromosomes but occur in preferential 1–2-kilobase sequences termed hotspots. Different species such as humans and mice do not share hotspots, although the same principles almost certainly regulate their placement in the genome. Hotspot positions and activities depend on genetic background and show sex-specific differences. In this study, we present a detailed analysis of recombination activity along the largest mouse chromosome, finding that recombination is regulated on multiple levels, including regional positioning relative to the chromosomal ends, local gene content, sex-specific mechanisms of hotspot recognition, and parental origin. Our results will contribute to further understanding of one of the most fundamental biological processes and are likely to cast light on several aspects of population genetics and evolutionary biology, as well as enhance our practical ability to define the genetic components of human disease.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

Evidence that duplications of 22q11.2 protect against schizophrenia

A number of large, rare copy number variants (CNVs) are deleterious for neurodevelopmental disorders, but large, rare, protective CNVs have not been reported for such phenotypes. Here we show in a CNV analysis of 47 005 individuals, the largest CNV analysis of schizophrenia to date, that large duplications (1.5-3.0 Mb) at 22q11.2--the reciprocal of the well-known, risk-inducing deletion of this locus--are substantially less common in schizophrenia cases than in the general population (0.014% vs 0.085%, OR=0.17, P=0.00086). 22q11.2 duplications represent the first putative protective mutation for schizophrenia.

Crossover interference underlies sex differences in recombination rates

In many organisms, recombination rates differ between the two sexes. Here we show that in mice, this is because of a shorter genomic interference distance in females than in males, measured in Mb. However, the interference distance is the same in terms of bivalent length. We propose a model in which the interference distance in the two sexes reflects the compaction of chromosomes at the pachytene stage of meiosis.

CNV analysis in a large schizophrenia sample implicates deletions at 16p12.1 and SLC1A1 and duplications at 1p36.33 and CGNL1

Large and rare copy number variants (CNVs) at several loci have been shown to increase risk for schizophrenia. Aiming to discover novel susceptibility CNV loci, we analyzed 6882 cases and 11 255 controls genotyped on Illumina arrays, most of which have not been used for this purpose before. We identified genes enriched for rare exonic CNVs among cases, and then attempted to replicate the findings in additional 14 568 cases and 15 274 controls. In a combined analysis of all samples, 12 distinct loci were enriched among cases with nominal levels of significance (P < 0.05); however, none would survive correction for multiple testing. These loci include recurrent deletions at 16p12.1, a locus previously associated with neurodevelopmental disorders (P = 0.0084 in the discovery sample and P = 0.023 in the replication sample). Other plausible candidates include non-recurrent deletions at the glutamate transporter gene SLC1A1, a CNV locus recently suggested to be involved in schizophrenia through linkage analysis, and duplications at 1p36.33 and CGNL1. A burden analysis of large (>500 kb), rare CNVs showed a 1.2% excess in cases after excluding known schizophrenia-associated loci, suggesting that additional susceptibility loci exist. However, even larger samples are required for their discovery.

An imputed genotype resource for the laboratory mouse

We have created a high-density SNP resource encompassing 7.87 million polymorphic loci across 49 inbred mouse strains of the laboratory mouse by combining data available from public databases and training a hidden Markov model to impute missing genotypes in the combined data. The strong linkage disequilibrium found in dense sets of SNP markers in the laboratory mouse provides the basis for accurate imputation. Using genotypes from eight independent SNP resources, we empirically validated the quality of the imputed genotypes and demonstrated that they are highly reliable for most inbred strains. The imputed SNP resource will be useful for studies of natural variation and complex traits. It will facilitate association study designs by providing high-density SNP genotypes for large numbers of mouse strains. We anticipate that this resource will continue to evolve as new genotype data become available for laboratory mouse strains. The data are available for bulk download or query at http://cgd.jax.org /.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Endogenous retinoids in the pathogenesis of alopecia areata

Alopecia areata (AA) is an autoimmune disease that attacks anagen hair follicles. Gene array in graft-induced C3H/HeJ mice revealed that genes involved in retinoic acid (RA) synthesis were increased, whereas RA degradation genes were decreased in AA compared with sham controls. This was confirmed by immunohistochemistry in biopsies from patients with AA and both mouse and rat AA models. RA levels were also increased in C3H/HeJ mice with AA. C3H/HeJ mice were fed a purified diet containing one of the four levels of dietary vitamin A or an unpurified diet 2 weeks before grafting and disease progression followed. High vitamin A accelerated AA, whereas mice that were not fed vitamin A had more severe disease by the end of the study. More hair follicles were in anagen in mice fed high vitamin A. Both the number and localization of granzyme B-positive cells were altered by vitamin A. IFNγ was also the lowest and IL13 highest in mice fed high vitamin A. Other cytokines were reduced and chemokines increased as the disease progressed, but no additional effects of vitamin A were seen. Combined, these results suggest that vitamin A regulates both the hair cycle and immune response to alter the progression of AA.

An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100 kb), rare copy-number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1564 cases and 1748 controls all from Ireland, and further extended the analysis to include an additional 5196 UK controls. We found association with duplications at chr20p12.2 (P = 0.007) and evidence of replication in large independent European schizophrenia (P = 0.052) and UK bipolar disorder case-control cohorts (P = 0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11 707 cases and 10 carriers in 21 204 controls [meta-analysis Cochran-Mantel-Haenszel P-value = 2 × 10(-4); odds ratio (OR) = 11.3, 95% CI = 3.7, ∞]. Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68 Mb with similar breakpoints across samples. By haplotype analysis and sequencing, we identified a tandem ~149 kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P = 2.5 × 10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions

Background

Otitis media (OM) is a common worldwide pediatric health care problem that is known to be influenced by genetics. The objective of our study was to use linkage analysis to map possible OM susceptibility genes.

Methods

Using a stringent diagnostic model in which only those who underwent tympanostomy tube insertion at least once for recurrent/persistent OM are considered affected, we have carried out a genome-wide linkage scan using the 10K Affymetrix SNP panel. We genotyped 403 Caucasian families containing 1,431 genotyped individuals and 377 genotyped affected sib pairs, and 26 African American families containing 75 genotyped individuals and 27 genotyped affected sib pairs. After careful quality control, non-parametric linkage analysis was carried out using 8,802 SNPs.

Results

In the Caucasian-only data set, our most significant linkage peak is on chromosome 17q12 at rs226088 with a p-value of 0.00007. Other peaks of potential interest are on 10q22.3 (0.00181 at rs1878001), 7q33 (0.00105 at rs958408), 6p25.1 (0.00261 at rs554653), and 4p15.2 (0.00301 at rs2133507). In the combined Caucasian and African American dataset, the 10q22.3 peak becomes more significant, with a minimal p-value of 0.00026 at rs719871. Family-based association testing reveals signals near previously implicated genes: 513 kb from SFTPA2 (10q22.3), 48 kb from IFNG (12q14), and 870 kb from TNF (6p21.3).

Conclusion

Our scan does not provide evidence for linkage in the previously reported regions of 10q26.3 and 19q13.43. Our best-supported linkage regions may contain susceptibility genes that influence the risk for recurrent/persistent OM. Plausible candidates in 17q12 include AP2B1, CCL5, and a cluster of other CCL genes, and in 10q22.3, SFTPA2.

Chromosome neighborhood composition determines translocation outcomes after exposure to high-dose radiation in primary cells

Radiation exposure is an occupational hazard for military personnel, some health care professionals, airport security screeners, and medical patients, with some individuals at risk for acute, high-dose exposures. Therefore, the biological effects of radiation, especially the potential for chromosome damage, are major occupational and health concerns. However, the biophysical mechanisms of chromosome instability subsequent to radiation-induced DNA damage are poorly understood. It is clear that interphase chromosomes occupy discrete structural and functional subnuclear domains, termed chromosome territories (CT), which may be organized into 'neighborhoods' comprising groups of specific CTs. We directly evaluated the relationship between chromosome positioning, neighborhood composition, and translocation partner choice in primary lymphocytes, using a cell-based system in which we could induce multiple, concentrated DNA breaks via high-dose irradiation. We critically evaluated mis-rejoining profiles and tested whether breaks occurring nearby were more likely to fuse than breaks occurring at a distance. We show that CT neighborhoods comprise heterologous chromosomes, within which inter-CT distances directly relate to translocation partner choice. These findings demonstrate that interphase chromosome arrangement is a principal factor in genomic instability outcomes in primary lymphocytes, providing a structural context for understanding the biological effects of radiation exposure, and the molecular etiology of tumor-specific translocation patterns.

Improving detection of copy-number variation by simultaneous bias correction and read-depth segmentation

Structural variation is an important class of genetic variation in mammals. High-throughput sequencing (HTS) technologies promise to revolutionize copy-number variation (CNV) detection but present substantial analytic challenges. Converging evidence suggests that multiple types of CNV-informative data (e.g. read-depth, read-pair, split-read) need be considered, and that sophisticated methods are needed for more accurate CNV detection. We observed that various sources of experimental biases in HTS confound read-depth estimation, and note that bias correction has not been adequately addressed by existing methods. We present a novel read-depth-based method, GENSENG, which uses a hidden Markov model and negative binomial regression framework to identify regions of discrete copy-number changes while simultaneously accounting for the effects of multiple confounders. Based on extensive calibration using multiple HTS data sets, we conclude that our method outperforms existing read-depth-based CNV detection algorithms. The concept of simultaneous bias correction and CNV detection can serve as a basis for combining read-depth with other types of information such as read-pair or split-read in a single analysis. A user-friendly and computationally efficient implementation of our method is freely available.

Implication of LRRC4C and DPP6 in neurodevelopmental disorders

We performed whole-genome sequencing on an individual from a family with variable psychiatric phenotypes that had a sensory processing disorder, apraxia, and autism. The proband harbored a maternally inherited balanced translocation (46,XY,t(11;14)(p12;p12)mat) that disrupted LRRC4C, a member of the highly specialized netrin G family of axon guidance molecules. The proband also inherited a paternally derived chromosomal inversion that disrupted DPP6, a potassium channel interacting protein. Copy Number (CN) analysis in 14,077 cases with neurodevelopmental disorders and 8,960 control subjects revealed that 60% of cases with exonic deletions in LRRC4C had a second clinically recognizable syndrome associated with variable clinical phenotypes, including 16p11.2, 1q44, and 2q33.1 CN syndromes, suggesting LRRC4C deletion variants may be modifiers of neurodevelopmental disorders. In vitro, functional assessments modeling patient deletions in LRRC4C suggest a negative regulatory role of these exons found in the untranslated region of LRRC4C, which has a single, terminal coding exon. These data suggest that the proband's autism may be due to the inheritance of disruptions in both DPP6 and LRRC4C, and may highlight the importance of the netrin G family and potassium channel interacting molecules in neurodevelopmental disorders. © 2016 Wiley Periodicals, Inc.

Defective carbohydrate metabolism in mice homozygous for the tubby mutation

Tub is a member of a small gene family, the tubby-like proteins (TULPs), with predominant expression in neurons. Mice carrying a mutation in Tub develop retinal and cochlear degeneration as well as late-onset obesity with insulin resistance. During behavioral and metabolic testing, we found that homozygous C57BL/6J-Tub(tub) mice have a lower respiratory quotient than C57BL/6J controls before the onset of obesity, indicating that tubby homozygotes fail to activate carbohydrate metabolism and instead rely on fat metabolism for energy needs. In concordance with this, tubby mice show higher excretion of ketone bodies and accumulation of glycogen in the liver. Quantitation of liver mRNA levels shows that, during the transition from light to dark period, tubby mice fail to induce glucose-6-phosphate dehydrogenase (G6pdh), the rate-limiting enzyme in the pentose phosphate pathway that normally supplies NADPH for de novo fatty acid synthesis and glutathione reduction. Reduced G6PDH protein levels and enzymatic activity in tubby mice lead accordingly to lower levels of NADPH and reduced glutathione (GSH), respectively. mRNA levels for the lipolytic enzymes acetyl-CoA synthetase and carnitine palmitoyltransferase are increased during the dark cycle and decreased during the light period, and several citric acid cycle genes are dysregulated in tubby mice. Examination of hypothalamic gene expression showed high levels of preproorexin mRNA leading to accumulation of orexin peptide in the lateral hypothalamus. We hypothesize that abnormal hypothalamic orexin expression leads to changes in liver carbohydrate metabolism and may contribute to the moderate obesity observed in tubby mice.

A survey of airway responsiveness in 36 inbred mouse strains facilitates gene mapping studies and identification of quantitative trait loci

Airway hyper-responsiveness (AHR) is a critical phenotype of human asthma and animal models of asthma. Other studies have measured AHR in nine mouse strains, but only six strains have been used to identify genetic loci underlying AHR. Our goals were to increase the genetic diversity of available strains by surveying 27 additional strains, to apply haplotype association mapping to the 36-strain survey, and to identify new genetic determinants for AHR. We derived AHR from the increase in airway resistance in females subjected to increasing levels of methacholine concentrations. We used haplotype association mapping to identify associations between AHR and haplotypes on chromosomes 3, 5, 8, 12, 13, and 14. And we used bioinformatics techniques to narrow the identified region on chromosome 13, reducing the region to 29 candidate genes, with 11 of considerable interest. Our combined use of haplotype association mapping with bioinformatics tools is the first study of its kind for AHR on these 36 strains of mice. Our analyses have narrowed the possible QTL genes and will facilitate the discovery of novel genes that regulate AHR in mice.

Association of candidate genes with phenotypic traits relevant to anorexia nervosa

This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein, we investigate alternative phenotypes associated with AN. In 1762 females, using regression analyses, we examined the following: (i) lowest illness-related attained body mass index; (ii) age at menarche; (iii) drive for thinness; (iv) body dissatisfaction; (v) trait anxiety; (vi) concern over mistakes; and (vii) the anticipatory worry and pessimism versus uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes.