A compendium of human genes regulating feeding behavior and body weight, its functional characterization and identification of GWAS genes involved in brain-specific PPI network

BACKGROUND: Obesity is heritable. It predisposes to many diseases. The objectives of this study were to create a compendium of genes relevant to feeding behavior (FB) and/or body weight (BW) regulation; to construct and to analyze networks formed by associations between genes/proteins; and to identify the most significant genes, biological processes/pathways, and tissues/organs involved in BW regulation. RESULTS: The compendium of genes controlling FB or BW includes 578 human genes. Candidate genes were identified from various sources, including previously published original research and review articles, GWAS meta-analyses, and OMIM (Online Mendelian Inheritance in Man). All genes were ranked according to knowledge about their biological role in body weight regulation and classified according to expression patterns or functional characteristics. Substantial and overrepresented numbers of genes from the compendium encoded cell surface receptors, signaling molecules (hormones, neuropeptides, cytokines), transcription factors, signal transduction proteins, cilium and BBSome components, and lipid binding proteins or were present in the brain-specific list of tissue-enriched genes identified with TSEA tool. We identified 27 pathways from KEGG, REACTOME and BIOCARTA whose genes were overrepresented in the compendium. Networks formed by physical interactions or homological relationships between proteins or interactions between proteins involved in biochemical/signaling pathways were reconstructed and analyzed. Subnetworks and clusters identified by the MCODE tool included genes/proteins associated with cilium morphogenesis, signal transduction proteins (particularly, G protein-coupled receptors, kinases or proteins involved in response to insulin stimulus) and transcription regulation (particularly nuclear receptors). We ranked GWAS genes according to the number of neighbors in three networks and revealed 22 GWAS genes involved in the brain-specific PPI network. On the base of the most reliable PPIs functioning in the brain tissue, new regulatory schemes interpreting relevance to BW regulation are proposed for three GWAS genes (ETV5, LRP1B, and NDUFS3). CONCLUSIONS: A compendium comprising 578 human genes controlling FB or BW was designed, and the most significant functional groups of genes, biological processes/pathways, and tissues/organs involved in BW regulation were revealed. We ranked genes from the GWAS meta-analysis set according to the number and quality of associations in the networks and then according to their involvement in the brain-specific PPI network and proposed new regulatory schemes involving three GWAS genes (ETV5, LRP1B, and NDUFS3) in BW regulation. The compendium is expected to be useful for pathology risk estimation and for design of new pharmacological approaches in the treatment of human obesity

Comparing Three Approaches

Hybridization-based target enrichment protocols require relatively large starting amounts of genomic DNA, which is not always available. Here, we tested three approaches to pre-capture library preparation starting from 10 ng of genomic DNA: (i and ii) whole-genome amplification of DNA samples with REPLI-g (Qiagen) and GenomePlex (Sigma) kits followed by standard library preparation, and (iii) library construction with a low input oriented ThruPLEX kit (Rubicon Genomics). Exome capture with Agilent SureSelectXT2 Human AllExon v4+UTRs capture probes, and HiSeq2000 sequencing were performed for test libraries along with the control library prepared from 1 µg of starting DNA. Tested protocols were characterized in terms of mapping efficiency, enrichment ratio, coverage of the target region, and reliability of SNP genotyping. REPLI-g- and ThruPLEX-FD-based protocols seem to be adequate solutions for exome sequencing of low input sample

Haplotype analysis of APOE intragenic SNPs

BACKGROUND: APOE epsilon4 allele is most common genetic risk factor for Alzheimer\u27s disease (AD) and cognitive decline. However, it remains poorly understood why only some carriers of APOE epsilon4 develop AD and how ethnic variabilities in APOE locus contribute to AD risk. Here, to address the role of APOE haplotypes, we reassessed the diversity of APOE locus in major ethnic groups and in Alzheimer\u27s Disease Neuroimaging Initiative (ADNI) dataset on patients with AD, and subjects with mild cognitive impairment (MCI), and control non-demented individuals. RESULTS: We performed APOE gene haplotype analysis for a short block of five SNPs across the gene using the ADNI whole genome sequencing dataset. The compilation of ADNI data with 1000 Genomes identified the APOE epsilon4 linked haplotypes, which appeared to be distant for the Asian, African and European populations. The common European epsilon4-bearing haplotype is associated with AD but not with MCI, and the Africans lack this haplotype. Haplotypic inference revealed alleles that may confer protection against AD. By assessing the DNA methylation profile of the APOE haplotypes, we found that the AD-associated haplotype features elevated APOE CpG content, implying that this locus can also be regulated by genetic-epigenetic interactions. CONCLUSIONS: We showed that SNP frequency profiles within APOE locus are highly skewed to population-specific haplotypes, suggesting that the ancestral background within different sites at APOE gene may shape the disease phenotype. We propose that our results can be utilized for more specific risk assessment based on population descent of the individuals and on higher specificity of five site haplotypes associated with AD

Genome analysis of American minks reveals link of mutations in Ras-related protein-38 gene to Moyle brown coat phenotype

Over 35 fur colours have been described in American mink (Neovison vison), only six of which have been previously linked to specific genes. Moyle fur colour belongs to a wide group of brownish colours that are highly similar to each other, which complicates selection and breeding procedures. We performed whole genome sequencing for two American minks with Moyle (m/m) and Violet (a/a m/m /p/p) phenotypes. We identified two frame-shift mutations in the gene encoding Ras-related protein-38 (RAB38), which regulates the trafficking of tyrosinase-containing vesicles to maturing melanosomes. The results highlight the role of RAB38 in the biogenesis of melanosomes and melanin and the genetic mechanism contributing to hair colour variety and intensity. These data are also useful for tracking economically valuable fur traits in mink breeding programmes

Mutational re-modeling of di-aspartyl intramembrane proteases: uncoupling physiologically-relevant activities from those associated with Alzheimer\u27s disease

The intramembrane proteolytic activities of presenilins (PSEN1/PS1 and PSEN2/PS2) underlie production of beta-amyloid, the key process in Alzheimer\u27s disease (AD). Dysregulation of presenilin-mediated signaling is linked to cancers. Inhibition of the gamma-cleavage activities of PSENs that produce Abeta, but not the epsilon-like cleavage activity that release physiologically essential transcription activators, is a potential approach for the development of rational therapies for AD. In order to identify whether different activities of PSEN1 can be dissociated, we designed multiple mutations in the evolutionary conserved sites of PSEN1. We tested them in vitro and in vivo assays and compared their activities with mutant isoforms of presenilin-related intramembrane di-aspartyl protease (IMPAS1 (IMP1)/signal peptide peptidase (SPP)). PSEN1 auto-cleavage was more resistant to the mutation remodeling than the epsilon-like proteolysis. PSEN1 with a G382A or a P433A mutation in evolutionary invariant sites retains functionally important APP epsilon- and Notch S3- cleavage activities, but G382A inhibits APP gamma-cleavage and Abeta production and a P433A elevates Abeta. The G382A variant cannot restore the normal cellular ER Ca(2+) leak in PSEN1/PSEN2 double knockout cells, but efficiently rescues the loss-of-function (Egl) phenotype of presenilin in C. elegans. We found that, unlike in PSEN1 knockout cells, endoplasmic reticulum (ER) Ca(2+) leak is not changed in the absence of IMP1/SPP. IMP1/SPP with the analogous mutations retained efficiency in cleavage of transmembrane substrates and rescued the lethality of Ce-imp-2 knockouts. In summary, our data show that mutations near the active catalytic sites of intramembrane di-aspartyl proteases have different consequences on proteolytic and signaling functions

Neurodevelopmental Syndrome with Intellectual Disability, Speech Impairment, and Quadrupedia Is Associated with Glutamate Receptor Delta 2 Gene Defect

Bipedalism, speech, and intellect are the most prominent traits that emerged in the evolution of; Homo sapiens; . Here, we describe a novel genetic cause of an "involution" phenotype in four patients, who are characterized by quadrupedal locomotion, intellectual impairment, the absence of speech, small stature, and hirsutism, observed in a consanguineous Brazilian family. Using whole-genome sequencing analysis and homozygous genetic mapping, we identified genes bearing homozygous genetic variants and found a homozygous 36.2 kb deletion in the gene of glutamate receptor delta 2 (; GRID2; ) in the patients, resulting in the lack of a coding region from the fifth to the seventh exons. The; GRID2; gene is highly expressed in the cerebellum cortex from prenatal development to adulthood, specifically in Purkinje neurons. Deletion in this gene leads to the loss of the alpha chain in the extracellular amino-terminal protein domain (ATD), essential in protein folding and transport from the endoplasmic reticulum (ER) to the cell surface. Then, we studied the evolutionary trajectories of the; GRID2; gene. There was no sign of strong selection of the highly conservative; GRID2; gene in ancient hominids (Neanderthals and Denisovans) or modern humans; however, according to in silico tests using the Mfold tool, the; GRID; 2 gene possibly gained human-specific mutations that increased the stability of; GRID2; mRNA

Chromatin profiling of cortical neurons identifies individual epigenetic signatures in schizophrenia

Both heritability and environment contribute to risk for schizophrenia. However, the molecular mechanisms of interactions between genetic and non-genetic factors remain unclear. Epigenetic regulation of neuronal genome may be a presumable mechanism in pathogenesis of schizophrenia. Here, we performed analysis of open chromatin landscape of gene promoters in prefrontal cortical (PFC) neurons from schizophrenic patients. We cataloged cell-type-based epigenetic signals of transcriptional start sites (TSS) marked by histone H3-K4 trimethylation (H3K4me3) across the genome in PFC from multiple schizophrenia subjects and age-matched control individuals. One of the top-ranked chromatin alterations was found in the major histocompatibility (MHC) locus on chromosome 6 highlighting the overlap between genetic and epigenetic risk factors in schizophrenia. The chromosome conformation capture (3C) analysis in human brain cells revealed the architecture of multipoint chromatin interactions between the schizophrenia-associated genetic and epigenetic polymorphic sites and distantly located HLA-DRB5 and BTNL2 genes. In addition, schizophrenia-specific chromatin modifications in neurons were particularly prominent for non-coding RNA genes, including an uncharacterized LINC01115 gene and recently identified BNRNA_052780. Notably, protein-coding genes with altered epigenetic state in schizophrenia are enriched for oxidative stress and cell motility pathways. Our results imply the rare individual epigenetic alterations in brain neurons are involved in the pathogenesis of schizophrenia

Whole exome sequencing links dental tumor to an autosomal-dominant mutation in ANO5 gene associated with gnathodiaphyseal dysplasia and muscle dystrophies

Tumors of the jaws may represent different human disorders and frequently associate with pathologic bone fractures. In this report, we analyzed two affected siblings from a family of Russian origin, with a history of dental tumors of the jaws, in correspondence to original clinical diagnosis of cementoma consistent with gigantiform cementoma (GC, OMIM: 137575). Whole exome sequencing revealed the heterozygous missense mutation c.1067G \u3e A (p.Cys356Tyr) in ANO5 gene in these patients. To date, autosomal-dominant mutations have been described in the ANO5 gene for gnathodiaphyseal dysplasia (GDD, OMIM: 166260), and multiple recessive mutations have been described in the gene for muscle dystrophies (OMIM: 613319, 611307); the same amino acid (Cys) at the position 356 is mutated in GDD. These genetic data and similar clinical phenotypes demonstrate that the GC and GDD likely represent the same type of bone pathology. Our data illustrate the significance of mutations in single amino-acid position for particular bone tissue pathology. Modifying role of genetic variations in another gene on the severity of the monogenic trait pathology is also suggested. Finally, we propose the model explaining the tissue-specific manifestation of clinically distant bone and muscle diseases linked to mutations in one gene

Human-specific histone methylation signatures at transcription start sites in prefrontal neurons

Cognitive abilities and disorders unique to humans are thought to result from adaptively driven changes in brain transcriptomes, but little is known about the role of cis-regulatory changes affecting transcription start sites (TSS). Here, we mapped in human, chimpanzee, and macaque prefrontal cortex the genome-wide distribution of histone H3 trimethylated at lysine 4 (H3K4me3), an epigenetic mark sharply regulated at TSS, and identified 471 sequences with human-specific enrichment or depletion. Among these were 33 loci selectively methylated in neuronal but not non-neuronal chromatin from children and adults, including TSS at DPP10 (2q14.1), CNTN4 and CHL1 (3p26.3), and other neuropsychiatric susceptibility genes. Regulatory sequences at DPP10 and additional loci carried a strong footprint of hominid adaptation, including elevated nucleotide substitution rates and regulatory motifs absent in other primates (including archaic hominins), with evidence for selective pressures during more recent evolution and adaptive fixations in modern populations. Chromosome conformation capture at two neurodevelopmental disease loci, 2q14.1 and 16p11.2, revealed higher order chromatin structures resulting in physical contact of multiple human-specific H3K4me3 peaks spaced 0.5-1 Mb apart, in conjunction with a novel cis-bound antisense RNA linked to Polycomb repressor proteins and downregulated DPP10 expression. Therefore, coordinated epigenetic regulation via newly derived TSS chromatin could play an important role in the emergence of human-specific gene expression networks in brain that contribute to cognitive functions and neurological disease susceptibility in modern day humans