The anatomical distribution of genetic associations

The Relevance of Reperfusion Stroke Therapy for miR-9-3p and miR-9-5p Expression in Acute Stroke-A Preliminary Study

Reperfusion stroke therapy is a modern treatment that involves thrombolysis and the mechanical removal of thrombus from the extracranial and/or cerebral arteries, thereby increasing penumbra reperfusion. After reperfusion therapy, 46% of patients are able to live independently 3 months after stroke onset. MicroRNAs (miRNAs) are essential regulators in the development of cerebral ischemia/reperfusion injury and the efficacy of the applied treatment. The first aim of this study was to examine the change in serum miRNA levels via next-generation sequencing (NGS) 10 days after the onset of acute stroke and reperfusion treatment. Next, the predictive values of the bioinformatics analysis of miRNA gene targets for the assessment of brain ischemic response to reperfusion treatment were explored. Human serum samples were collected from patients on days 1 and 10 after stroke onset and reperfusion treatment. The samples were subjected to NGS and then validated using qRT-PCR. Differentially expressed miRNAs (DEmiRNAs) were used for enrichment analysis. Hsa-miR-9-3p and hsa-miR-9-5p expression were downregulated on day 10 compared to reperfusion treatment on day 1 after stroke. The functional analysis of miRNA target genes revealed a strong association between the identified miRNA and stroke-related biological processes related to neuroregeneration signaling pathways. Hsa-miR-9-3p and hsa-miR-9-5p are potential candidates for the further exploration of reperfusion treatment efficacy in stroke patients.

Increased burden of rare protein-truncating variants in constrained, brain-specific and synaptic genes in extremely impulsively violent males with antisocial personality disorder

The genetic correlates of extreme impulsive violence are poorly understood, and there have been few studies that have characterized a large group of affected individuals both clinically and genetically. We performed whole exome sequencing (WES) in 290 males with the life-course-persistent, extremely impulsively violent form of antisocial personality disorder (APD) and analyzed the spectrum of rare protein-truncating variants (rPTVs). Comparisons were made with 314 male controls and publicly available genotype data. Functional annotation tools were used for biological interpretation. Participants were significantly more likely to harbor rPTVs in genes that are intolerant to loss-of-function variants (odds ratio [OR] 2.06; p < 0.001), specifically expressed in brain (OR 2.80; p = 0.036) and enriched for those involved in neurotransmitter transport and synaptic processes. In 60 individuals (20%), we identified rPTVs that we classified as clinically relevant based on their clinical associations, biological function and gene expression patterns. Of these, 37 individuals harbored rPTVs in 23 genes that are associated with a monogenic neurological disorder, and 23 individuals harbored rPTVs in 20 genes reportedly intolerant to loss-of-function variants. The analysis presents evidence in support of a model where presence of either one or several private, functionally relevant mutations contribute significantly to individual risk of life-course-persistent APD and reveals multiple individuals who could be affected by clinically unrecognized neuropsychiatric Mendelian disease. Thus, Mendelian diseases and increased rPTV burden may represent important factors for the development of extremely impulsive violent life-course-persistent forms of APD irrespective of their clinical presentation.

Sbp2l contributes to oligodendrocyte maturation through translational control in Tcf7l2 signaling

Oligodendrocytes (OLs) are the myelin-forming cells in the CNS that support neurons through the insulating sheath of axons. This unique feature and developmental processes are achieved by extrinsic and intrinsic gene expression programs, where RNA-binding proteins can contribute to dynamic and fine-tuned post-transcriptional regulation. Here, we identified SECIS-binding protein 2-like (Sbp2l), which is specifically expressed in OLs by integrated transcriptomics. Histological analysis revealed that Sbp2l is a molecular marker of OL maturation. Sbp2l knockdown (KD) led to suppression of matured OL markers, but not a typical selenoprotein, Gpx4. Transcriptome analysis demonstrated that Sbp2l KD decreased cholesterol-biosynthesis-related genes regulated by Tcf7l2 transcription factor. Indeed, we confirmed the downregulation of Tcf7l2 protein without changing its mRNA in Sbp2l KD OPCs. Furthermore, Sbp2l KO mice showed the decrease of Tcf7l2 protein and deficiency of OL maturation. These results suggest that Sbp2l contributes to OL maturation by translational control of Tcf7l2.

Microarray meta-analysis reveals comprehensive effects of 3,4,5-tricaffeolyquinic acid in cell differentiation and signaling

Caffeoylquinic acids (CQA) are polyphenolic compounds found in fruits, vegetables, coffee, and spices that have exhibited several beneficial activities, including antioxidant, antibacterial, neuroprotective, anti-inflammatory, anticancer, antiviral, antidiabetic, and cardiovascular effects. A derivative, TCQA (3,4,5-Tri-O-caffeoylquinic acid), has also shown both neurogenic and pigment differentiation potential. A transcriptomic-based meta-analysis was conducted to explore potential biochemical processes and molecular targets of TCQA. This approach involved integrating data from various cell and tissue types, including human amniotic stem cells, human neural stem cells, human dermal papilla cells, and the brain cortex of aging model mice. It offered a comprehensive perspective on the significant gene regulations in response to TCQA treatment. The objective was to uncover the mechanism and novel targets of TCQA, facilitating a further understanding of its functions. New areas of interest found were TCQA's effect on adipogenesis, heart, and muscle tissue development. In addition, significantly enhanced biological activities found through meta-analysis included cell cycle, VEGFA-VEGFR2 pathway, and BMP signaling. Overall, a comprehensive functional and visual analysis using available biological databases uncovered the multi-target potential of this natural compound.

Evolution of human genes encoding cell surface receptors involved in the regulation of appetite: an analysis based on the phylostratigraphic age and divergence indexes

Genes encoding cell surface receptors make up a significant portion of the human genome (more than a thousand genes) and play an important role in gene networks. Cell surface receptors are transmembrane proteins that interact with molecules (ligands) located outside the cell. This interaction activates signal transduction pathways in the cell. A large number of exogenous ligands of various origins, including drugs, are known for cell surface receptors, which accounts for interest in them from biomedical researchers. Appetite (the desire of the animal organism to consume food) is one of the most primitive instincts that contribute to survival. However, when the supply of nutrients is stable, the mechanism of adaptation to adverse factors acquired in the course of evolution turned out to be excessive, and therefore obesity has become one of the most serious public health problems of the twenty-first century. Pathological human conditions characterized by appetite violations include both hyperphagia, which inevitably leads to obesity, and anorexia nervosa induced by psychosocial stimuli, as well as decreased appetite caused by neurodegeneration, inflammation or cancer. Understanding the evolutionary mechanisms of human diseases, especially those related to lifestyle changes that have occurred over the past 100-200 years, is of fundamental and applied importance. It is also very important to identify relationships between the evolutionary characteristics of genes in gene networks and the resistance of these networks to changes caused by mutations. The aim of the current study is to identify the distinctive features of human genes encoding cell surface receptors involved in appetite regulation using the phylostratigraphic age index (PAI) and divergence index (DI). The values of PAI and DI were analyzed for 64 human genes encoding cell surface receptors, the orthologs of which were involved in the regulation of appetite in model animal species. It turned out that the set of genes under consideration contains an increased number of genes with the same phylostratigraphic age (PAI = 5, the stage of vertebrate divergence), and almost all of these genes (28 out of 31) belong to the superfamily of G-protein coupled receptors. Apparently, the synchronized evolution of such a large group of genes (31 genes out of 64) is associated with the development of the brain as a separate organ in the first vertebrates. When studying the distribution of genes from the same set by DI values, a significant enrichment with genes having a low DIs was revealed: eight genes (GPR26, NPY1R, GHSR, ADIPOR1, DRD1, NPY2R, GPR171, NPBWR1) had extremely low DIs (less than 0.05). Such low DI values indicate that most likely these genes are subjected to stabilizing selection. It was also found that the group of genes with low DIs was enriched with genes that had brain-specific patterns of expression. In particular, GPR26, which had the lowest DI, is in the group of brain-specific genes. Because the endogenous ligand for the GPR26 receptor has not yet been identified, this gene seems to be an extremely interesting object for further theoretical and experimental research. We believe that the features of the genes encoding cell surface receptors we have identified using the evolutionary metrics PAI and DI can be a starting point for further evolutionary analysis of the gene network regulating appetite.

Challenges in screening for de novo noncoding variants contributing to genetically complex phenotypes

Understanding the genetic basis for complex, heterogeneous disorders, such as autism spectrum disorder (ASD), is a persistent challenge in human medicine. Owing to their phenotypic complexity, the genetic mechanisms underlying these disorders may be highly variable across individual patients. Furthermore, much of their heritability is unexplained by known regulatory or coding variants. Indeed, there is evidence that much of the causal genetic variation stems from rare and de novo variants arising from ongoing mutation. These variants occur mostly in noncoding regions, likely affecting regulatory processes for genes linked to the phenotype of interest. However, because there is no uniform code for assessing regulatory function, it is difficult to separate these mutations into likely functional and nonfunctional subsets. This makes finding associations between complex diseases and potentially causal de novo single-nucleotide variants (dnSNVs) a difficult task. To date, most published studies have struggled to find any significant associations between dnSNVs from ASD patients and any class of known regulatory elements. We sought to identify the underlying reasons for this and present strategies for overcoming these challenges. We show that, contrary to previous claims, the main reason for failure to find robust statistical enrichments is not only the number of families sampled, but also the quality and relevance to ASD of the annotations used to prioritize dnSNVs, and the reliability of the set of dnSNVs itself. We present a list of recommendations for designing future studies of this sort that will help researchers avoid common pitfalls.

deCS: A Tool for Systematic Cell Type Annotations of Single-cell RNA Sequencing Data among Human Tissues

Single-cell RNA sequencing (scRNA-seq) is revolutionizing the study of complex and dynamic cellular mechanisms. However, cell type annotation remains a main challenge as it largely relies on a priori knowledge and manual curation, which is cumbersome and subjective. The increasing number of scRNA-seq datasets, as well as numerous published genetic studies, has motivated us to build a comprehensive human cell type reference atlas.Here, we present decoding Cell type Specificity (deCS), an automatic cell type annotation method augmented by a comprehensive collection of human cell type expression profiles and marker genes. We used deCS to annotate scRNA-seq data from various tissue types and systematically evaluated the annotation accuracy under different conditions, including reference panels, sequencing depth, and feature selection strategies. Our results demonstrate that expanding the references is critical for improving annotation accuracy. Compared to many existing state-of-the-art annotation tools, deCS significantly reduced computation time and increased accuracy. deCS can be integrated into the standard scRNA-seq analytical pipeline to enhance cell type annotation. Finally, we demonstrated the broad utility of deCS to identify trait-cell type associations in 51 human complex traits, providing deep insights into the cellular mechanisms underlying disease pathogenesis. All documents for deCS, including source code, user manual, demo data, and tutorials, are freely available at https://github.com/bsml320/deCS.

Molecular Landscape of Tourette's Disorder

Tourette's disorder (TD) is a highly heritable childhood-onset neurodevelopmental disorder and is caused by a complex interplay of multiple genetic and environmental factors. Yet, the molecular mechanisms underlying the disorder remain largely elusive. In this study, we used the available omics data to compile a list of TD candidate genes, and we subsequently conducted tissue/cell type specificity and functional enrichment analyses of this list. Using genomic data, we also investigated genetic sharing between TD and blood and cerebrospinal fluid (CSF) metabolite levels. Lastly, we built a molecular landscape of TD through integrating the results from these analyses with an extensive literature search to identify the interactions between the TD candidate genes/proteins and metabolites. We found evidence for an enriched expression of the TD candidate genes in four brain regions and the pituitary. The functional enrichment analyses implicated two pathways ('cAMP-mediated signaling' and 'Endocannabinoid Neuronal Synapse Pathway') and multiple biological functions related to brain development and synaptic transmission in TD etiology. Furthermore, we found genetic sharing between TD and the blood and CSF levels of 39 metabolites. The landscape of TD not only provides insights into the (altered) molecular processes that underlie the disease but, through the identification of potential drug targets (such as FLT3, NAALAD2, CX3CL1-CX3CR1, OPRM1, and HRH2), it also yields clues for developing novel TD treatments.

TP53 promotes lineage commitment of human embryonic stem cells through ciliogenesis and sonic hedgehog signaling

Aneuploidy, defective differentiation, and inactivation of the tumor suppressor TP53 all occur frequently during tumorigenesis. Here, we probe the potential links among these cancer traits by inactivating TP53 in human embryonic stem cells (hESCs). TP53-/- hESCs exhibit increased proliferation rates, mitotic errors, and low-grade structural aneuploidy; produce poorly differentiated immature teratomas in mice; and fail to differentiate into neural progenitor cells (NPCs) in vitro. Genome-wide CRISPR screen reveals requirements of ciliogenesis and sonic hedgehog (Shh) pathways for hESC differentiation into NPCs. TP53 deletion causes abnormal ciliogenesis in neural rosettes. In addition to restraining cell proliferation through CDKN1A, TP53 activates the transcription of BBS9, which encodes a ciliogenesis regulator required for proper Shh signaling and NPC formation. This developmentally regulated transcriptional program of TP53 promotes ciliogenesis, restrains Shh signaling, and commits hESCs to neural lineages.

A Catalog of Human Genes Associated With Pathozoospermia and Functional Characteristics of These Genes

Genetic causes of the global decline in male fertility are among the hot spots of scientific research in reproductive genetics. The most common way to evaluate male fertility in clinical trials is to determine semen quality. Lower semen quality is very often accompanied by subfertility or infertility, occurs in many diseases and can be caused by many factors, including genetic ones. The following forms of lowered semen quality (pathozoospermia) are known: azoospermia, oligozoospermia, asthenozoospermia, teratozoospermia, and some combined forms. To systematize information about the genetic basis of impaired spermatogenesis, we created a catalog of human genes associated with lowered semen quality (HGAPat) and analyzed their functional characteristics. The catalog comprises data on 126 human genes. Each entry of the catalog describes an association between an allelic variant of the gene and a particular form of lowered semen quality, extracted from the experimental study. Most genes included into the catalog are located on autosomes and are associated with such pathologies as non-obstructive azoospermia, oligozoospermia or asthenozoospermia. Slightly less than half of the included genes (43%) are expressed in the testes in a tissue-specific manner. Functional annotation of genes from the catalog showed that spermatogenic failure can be associated with mutations in genes that control biological processes essential for spermiogenesis (regulating DNA metabolism, cell division, formation of cellular structures, which provide cell movement) as well as with mutations in genes that control cellular responses to unfavorable conditions (stress factors, including oxidative stress and exposure to toxins).

Bromodomain protein BRDT directs ΔNp63 function and super-enhancer activity in a subset of esophageal squamous cell carcinomas

Esophageal squamous cell carcinoma (ESCC) is the predominant subtype of esophageal cancer with a particularly high prevalence in certain geographical regions and a poor prognosis with a 5-year survival rate of 15-25%. Despite numerous studies characterizing the genetic and transcriptomic landscape of ESCC, there are currently no effective targeted therapies. In this study, we used an unbiased screening approach to uncover novel molecular precision oncology targets for ESCC and identified the bromodomain and extraterminal (BET) family member bromodomain testis-specific protein (BRDT) to be uniquely expressed in a subgroup of ESCC. Experimental studies revealed that BRDT expression promotes migration but is dispensable for cell proliferation. Further mechanistic insight was gained through transcriptome analyses, which revealed that BRDT controls the expression of a subset of ΔNp63 target genes. Epigenome and genome-wide occupancy studies, combined with genome-wide chromatin interaction studies, revealed that BRDT colocalizes and interacts with ΔNp63 to drive a unique transcriptional program and modulate cell phenotype. Our data demonstrate that these genomic regions are enriched for super-enhancers that loop to critical ΔNp63 target genes related to the squamous phenotype such as KRT14, FAT2, and PTHLH. Interestingly, BET proteolysis-targeting chimera, MZ1, reversed the activation of these genes. Importantly, we observed a preferential degradation of BRDT by MZ1 compared with BRD2, BRD3, and BRD4. Taken together, these findings reveal a previously unknown function of BRDT in ESCC and provide a proof-of-concept that BRDT may represent a novel therapeutic target in cancer.

Multiomics dissection of molecular regulatory mechanisms underlying autoimmune-associated noncoding SNPs

More than 90% of autoimmune-associated variants are located in noncoding regions, leading to challenges in deciphering the underlying causal roles of functional variants and genes and biological mechanisms. Therefore, to reduce the gap between traditional genetic findings and mechanistic understanding of disease etiologies and clinical drug development, it is important to translate systematically the regulatory mechanisms underlying noncoding variants. Here, we prioritized functional noncoding SNPs with regulatory gene targets associated with 19 autoimmune diseases by incorporating hundreds of immune cell-specific multiomics data. The prioritized SNPs are associated with transcription factor (TF) binding, histone modification, or chromatin accessibility, indicating their allele-specific regulatory roles. Their target genes are significantly enriched in immunologically related pathways and other known immunologically related functions. We found that 90.1% of target genes are regulated by distal SNPs involving several TFs (e.g., the DNA-binding protein CCCTC-binding factor [CTCF]), suggesting the importance of long-range chromatin interaction in autoimmune diseases. Moreover, we predicted potential drug targets for autoimmune diseases, including 2 genes (NFKB1 and SH2B3) with known drug indications on other diseases, highlighting their potential drug repurposing opportunities. Taken together, these findings may provide useful information for future experimental follow-up and drug applications on autoimmune diseases.

Rare Genetic Diseases: Nature's Experiments on Human Development

Rare genetic diseases are the result of a continuous forward genetic screen that nature is conducting on humans. Here, we present epistemological and systems biology arguments highlighting the importance of studying these rare genetic diseases. We contend that the expanding catalog of mutations in ∼4,000 genes, which cause ∼6,500 diseases and their annotated phenotypes, offer a wide landscape for discovering fundamental mechanisms required for human development and involved in common diseases. Rare afflictions disproportionately affect the nervous system in children, but paradoxically, the majority of these disease-causing genes are evolutionarily ancient and ubiquitously expressed in human tissues. We propose that the biased prevalence of childhood rare diseases affecting nervous tissue results from the topological complexity of the protein interaction networks formed by ubiquitous and ancient proteins encoded by childhood disease genes. Finally, we illustrate these principles discussing Menkes disease, an example of the discovery power afforded by rare diseases.

The Integrated Landscape of Biological Candidate Causal Genes in Coronary Artery Disease

Background

Genome-wide association studies (GWASs) have identified more than 150 genetic loci that demonstrate robust association with coronary artery disease (CAD). In contrast to the success of GWAS, the translation from statistical signals to biological mechanism and exploration of causal genes for drug development remain difficult, owing to the complexity of gene regulatory and linkage disequilibrium patterns. We aim to prioritize the plausible causal genes for CAD at a genome-wide level.

Methods

We integrated the latest GWAS summary statistics with other omics data from different layers and utilized eight different computational methods to predict CAD potential causal genes. The prioritized candidate genes were further characterized by pathway enrichment analysis, tissue-specific expression analysis, and pathway crosstalk analysis.

Results

Our analysis identified 55 high-confidence causal genes for CAD, among which 15 genes (LPL, COL4A2, PLG, CDKN2B, COL4A1, FES, FLT1, FN1, IL6R, LPA, PCSK9, PSRC1, SMAD3, SWAP70, and VAMP8) ranked the highest priority because of consistent evidence from different data-driven approaches. GO analysis showed that these plausible causal genes were enriched in lipid metabolic and extracellular regions. Tissue-specific enrichment analysis revealed that these genes were significantly overexpressed in adipose and liver tissues. Further, KEGG and crosstalk analysis also revealed several key pathways involved in the pathogenesis of CAD.

Conclusion

Our study delineated the landscape of CAD potential causal genes and highlighted several biological processes involved in CAD pathogenesis. Further studies and experimental validations of these genes may shed light on mechanistic insights into CAD development and provide potential drug targets for future therapeutics.

A viral toolkit for recording transcription factor-DNA interactions in live mouse tissues

Transcription factors (TFs) enact precise regulation of gene expression through site-specific, genome-wide binding. Common methods for TF-occupancy profiling, such as chromatin immunoprecipitation, are limited by requirement of TF-specific antibodies and provide only end-point snapshots of TF binding. Alternatively, TF-tagging techniques, in which a TF is fused to a DNA-modifying enzyme that marks TF-binding events across the genome as they occur, do not require TF-specific antibodies and offer the potential for unique applications, such as recording of TF occupancy over time and cell type specificity through conditional expression of the TF-enzyme fusion. Here, we create a viral toolkit for one such method, calling cards, and demonstrate that these reagents can be delivered to the live mouse brain and used to report TF occupancy. Further, we establish a Cre-dependent calling cards system and, in proof-of-principle experiments, show utility in defining cell type-specific TF profiles and recording and integrating TF-binding events across time. This versatile approach will enable unique studies of TF-mediated gene regulation in live animal models.

Integrative Omics Approach to Identifying Genes Associated With Atrial Fibrillation

Rationale: GWAS (Genome-Wide Association Studies) have identified hundreds of genetic loci associated with atrial fibrillation (AF). However, these loci explain only a small proportion of AF heritability. Objective: To develop an approach to identify additional AF-related genes by integrating multiple omics data. Methods and Results: Three types of omics data were integrated: (1) summary statistics from the AFGen 2017 GWAS; (2) a whole blood EWAS (Epigenome-Wide Association Study) of AF; and (3) a whole blood TWAS (Transcriptome-Wide Association Study) of AF. The variant-level GWAS results were collapsed into gene-level associations using fast set-based association analysis. The CpG-level EWAS results were also collapsed into gene-level associations by an adapted SNP-set Kernel Association Test approach. Both GWAS and EWAS gene-based associations were then meta-analyzed with TWAS using a fixed-effects model weighted by the sample size of each data set. A tissue-specific network was subsequently constructed using the NetWAS (Network-Wide Association Study). The identified genes were then compared with the AFGen 2018 GWAS that contained more than triple the number of AF cases compared with AFGen 2017 GWAS. We observed that the multiomics approach identified many more relevant AF-related genes than using AFGen 2018 GWAS alone (1931 versus 206 genes). Many of these genes are involved in the development and regulation of heart- and muscle-related biological processes. Moreover, the gene set identified by multiomics approach explained much more AF variance than those identified by GWAS alone (10.4% versus 3.5%). Conclusions: We developed a strategy to integrate multiple omics data to identify AF-related genes. Our integrative approach may be useful to improve the power of traditional GWAS, which might be particularly useful for rare traits and diseases with limited sample size.

Genetics of schizophrenia in the South African Xhosa

Africa, the ancestral home of all modern humans, is the most informative continent for understanding the human genome and its contribution to complex disease. To better understand the genetics of schizophrenia, we studied the illness in the Xhosa population of South Africa, recruiting 909 cases and 917 age-, gender-, and residence-matched controls. Individuals with schizophrenia were significantly more likely than controls to harbor private, severely damaging mutations in genes that are critical to synaptic function, including neural circuitry mediated by the neurotransmitters glutamine, γ-aminobutyric acid, and dopamine. Schizophrenia is genetically highly heterogeneous, involving severe ultrarare mutations in genes that are critical to synaptic plasticity. The depth of genetic variation in Africa revealed this relationship with a moderate sample size and informed our understanding of the genetics of schizophrenia worldwide.

Biobank-driven genomic discovery yields new insight into atrial fibrillation biology

To identify genetic variation underlying atrial fibrillation, the most common cardiac arrhythmia, we performed a genome-wide association study of >1,000,000 people, including 60,620 atrial fibrillation cases and 970,216 controls. We identified 142 independent risk variants at 111 loci and prioritized 151 functional candidate genes likely to be involved in atrial fibrillation. Many of the identified risk variants fall near genes where more deleterious mutations have been reported to cause serious heart defects in humans (GATA4, MYH6, NKX2-5, PITX2, TBX5)¹, or near genes important for striated muscle function and integrity (for example, CFL2, MYH7, PKP2, RBM20, SGCG, SSPN). Pathway and functional enrichment analyses also suggested that many of the putative atrial fibrillation genes act via cardiac structural remodeling, potentially in the form of an 'atrial cardiomyopathy'², either during fetal heart development or as a response to stress in the adult heart.

The Expanding Toolkit of Translating Ribosome Affinity Purification

Translating ribosome affinity purification is a method initially developed for profiling mRNA from genetically defined cell types in complex tissues. It has been applied both to identify target molecules in cell types that are important for controlling a variety of behaviors in the brain, and to understand the molecular consequences on those cells due to experimental manipulations, ranging from drugs of abuse to disease-causing mutations. Since its inception, a variety of methodological advances are opening new avenues of investigation. These advances include a variety of new methods for targeting cells for translating ribosome affinity purification by features such as their projections or activity, additional tags and mouse reagents increasing the flexibility of the system, and new modifications of the method specifically focused on studying the regulation of translation. The latter includes methods to assess cell type-specific regulation of translation in specific subcellular compartments. Here, I provide a summary of these recent advances and resources, highlighting both new experimental opportunities and areas for future technical development.

A compendium and functional characterization of mammalian genes involved in adaptation to Arctic or Antarctic environments

Background

Many mammals are well adapted to surviving in extremely cold environments. These species have likely accumulated genetic changes that help them efficiently cope with low temperatures. It is not known whether the same genes related to cold adaptation in one species would be under selection in another species. The aims of this study therefore were: to create a compendium of mammalian genes related to adaptations to a low temperature environment; to identify genes related to cold tolerance that have been subjected to independent positive selection in several species; to determine promising candidate genes/pathways/organs for further empirical research on cold adaptation in mammals.

Results

After a search for publications containing keywords: “whole genome”, “transcriptome or exome sequencing data”, and “genome-wide genotyping array data” authors looked for information related to genetic signatures ascribable to positive selection in Arctic or Antarctic mammalian species. Publications related to Human, Arctic fox, Yakut horse, Mammoth, Polar bear, and Minke whale were chosen. The compendium of genes that potentially underwent positive selection in >1 of these six species consisted of 416 genes. Twelve of them showed traces of positive selection in three species. Gene ontology term enrichment analysis of 416 genes from the compendium has revealed 13 terms relevant to the scope of this study. We found that enriched terms were relevant to three major groups: terms associated with collagen proteins and the extracellular matrix; terms associated with the anatomy and physiology of cilium; terms associated with docking. We further revealed that genes from compendium were over-represented in the lists of genes expressed in the lung and liver.

Conclusions

A compendium combining mammalian genes involved in adaptation to cold environment was designed, based on the intersection of positively selected genes from six Arctic and Antarctic species. The compendium contained 416 genes that have been positively selected in at least two species. However, we did not reveal any positively selected genes that would be related to cold adaptation in all species from our list. But, our work points to several strong candidate genes involved in mechanisms and biochemical pathways related to cold adaptation response in different species.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0580-9) contains supplementary material, which is available to authorized users.

Lessons from ten years of genome-wide association studies of asthma

Twenty-five genome-wide association studies (GWAS) of asthma were published between 2007 and 2016, the largest with a sample size of 157242 individuals. Across these studies, 39 genetic variants in low linkage disequilibrium (LD) with each other were reported to associate with disease risk at a significance threshold of P<5 × 10⁻⁸, including 31 in populations of European ancestry. Results from analyses of the UK Biobank data (n=380 503) indicate that at least 28 of the 31 associations reported in Europeans represent true-positive findings, collectively explaining 2.5% of the variation in disease liability (median of 0.06% per variant). We identified 49 transcripts as likely target genes of the published asthma risk variants, mostly based on LD with expression quantitative trait loci (eQTL). Of these genes, 16 were previously implicated in disease pathophysiology by functional studies, including TSLP, TNFSF4, ADORA1, CHIT1 and USF1. In contrast, at present, there is limited or no functional evidence directly implicating the remaining 33 likely target genes in asthma pathophysiology. Some of these genes have a known function that is relevant to allergic disease, including F11R, CD247, PGAP3, AAGAB, CAMK4 and PEX14, and so could be prioritized for functional follow-up. We conclude by highlighting three areas of research that are essential to help translate GWAS findings into clinical research or practice, namely validation of target gene predictions, understanding target gene function and their role in disease pathophysiology and genomics-guided prioritization of targets for drug development.

Varying Intolerance of Gene Pathways to Mutational Classes Explain Genetic Convergence across Neuropsychiatric Disorders

Genetic susceptibility to intellectual disability (ID), autism spectrum disorder (ASD), and schizophrenia (SCZ) often arises from mutations in the same genes, suggesting that they share common mechanisms. We studied genes with de novo mutations in the three disorders and genes implicated in SCZ by genome-wide association study (GWAS). Using biological annotations and brain gene expression, we show that mutation class explains enrichment patterns more than specific disorder. Genes with loss-of-function mutations and genes with missense mutations were associated with different pathways across disorders. Conversely, gene expression patterns were specific for each disorder. ID genes were preferentially expressed in the cortex; ASD genes were expressed in the fetal cortex, cerebellum, and striatum; and genes associated with SCZ were expressed in the adolescent cortex. Our study suggests that convergence across neuropsychiatric disorders stems from common pathways that are consistently vulnerable to genetic variations but that spatiotemporal activity of genes contributes to specific phenotypes.

Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology

Asthma, hay fever (or allergic rhinitis) and eczema (or atopic dermatitis) often coexist in the same individuals1, partly because of a shared genetic origin2–4. To identify shared risk variants, we performed a genome-wide association study (GWAS, n=360,838) of a broad allergic disease phenotype that considers the presence of any one of these three diseases. We identified 136 independent risk variants (P<3x10^-8), including 73 not previously reported, which implicate 132 nearby genes in allergic disease pathophysiology. Disease-specific effects were detected for only six variants, confirming that most represent shared risk factors. Tissue-specific heritability and biological process enrichment analyses suggest that shared risk variants influence lymphocyte-mediated immunity. Six target genes provide an opportunity for drug repositioning, while for 36 genes CpG methylation was found to influence transcription independently of genetic effects. Asthma, hay fever and eczema partly coexist because they share many genetic risk variants that dysregulate the expression of immune-related genes.

Analysis of chromatin accessibility in human epidermis identifies putative barrier dysfunction-sensing enhancers

To identify putative gene regulatory regions that respond to epidermal injury, we mapped chromatin dynamics in a stratified human epidermis during barrier maturation and disruption. Engineered skin substitutes (ESS) cultured at the air-liquid interface were used as a model of developing human epidermis with incomplete barrier formation. The epidermal barrier stabilized following engraftment onto immunocompromised mice, and was compromised again upon injury. Modified formaldehyde-assisted isolation of regulatory elements (FAIRE) was used to identify accessible genomic regions characteristic of monolayer keratinocytes, ESS in vitro, grafted ESS, and tape-stripped ESS graft. We mapped differentiation- and maturation-associated changes in transcription factor binding sites enriched at each stage and observed overrepresentation of AP-1 gene family motifs in barrier-deficient samples. Transcription of TSLP, an important effector of immunological memory in response to allergen exposure, was dramatically elevated in our barrier-deficient samples. We identified dynamic DNA elements that correlated with TSLP induction and may contain enhancers that regulate TSLP. Two dynamic regions were located near the TSLP promoter and overlapped with allergy-associated SNPs rs17551370 and rs2289877, strongly implicating these loci in the regulation of TSLP expression in allergic disease. Additional dynamic chromatin regions ~250kb upstream of the TSLP promoter were found to be in high linkage disequilibrium with allergic disease SNPs. Taken together, these results define dynamic chromatin accessibility changes during epidermal development and dysfunction.

Transcriptomic Analysis of Ribosome-Bound mRNA in Cortical Neurites In Vivo

Localized translation in neurites helps regulate synaptic strength and development. Dysregulation of local translation is associated with many neurological disorders. However, due to technical limitations, study of this phenomenon has largely been limited to brain regions with laminar organization of dendrites such as the hippocampus or cerebellum. It has not been examined in the cortex, a region of importance for most neurological disorders, where dendrites of each neuronal population are densely intermingled with cell bodies of others. Therefore, we have developed a novel method, SynapTRAP, which combines synaptoneurosomal fractionation with translating ribosome affinity purification to identify ribosome-bound mRNA in processes of genetically defined cell types. We demonstrate SynapTRAP's efficacy and report local translation in the cortex of mice, where we identify a subset of mRNAs that are translated in dendrites by neuronal ribosomes. These mRNAs have disproportionately longer lengths, enrichment for FMRP binding and G-quartets, and their genes are under greater evolutionary constraint in humans. In addition, we show that alternative splicing likely regulates this phenomenon. Overall, SynapTRAP allows for rapid isolation of cell-type-specific localized translation and is applicable to classes of previously inaccessible neuronal and non-neuronal cells in vivoSIGNIFICANCE STATEMENT Instructions for making proteins are found in the genome, housed within the nucleus of each cell. These are then copied as RNA and exported to manufacture new proteins. However, in the brain, memory is thought to be encoded by strengthening individual connections (synapses) between neurons far from the nucleus. Thus, to efficiently make new proteins specifically where they are needed, neurons can transport RNAs to sites near synapses to locally produce proteins. Importantly, several mutations that cause autism disrupt this process. It has been assumed this process occurs in all brain regions, but has never been measured in the cortex. We applied a newly developed method measure to study, for the first time, local translation in cortical neurons.

Systems biology in the central nervous system: a brief perspective on essential recent advancements

As recent advances in human genetics have begun to more rapidly identify the individual genes contributing to risk of psychiatric disease, the spotlight now turns to understanding how disruption of these genes alters the brain, and thus behavior. Compared to other tissues, cellular complexity in the brain provides both a substantial challenge and a significant opportunity for systems biology approaches. Current methods are maturing that will allow for finally defining the 'parts list' for the functioning mouse and human brains, enabling new approaches to defining how the system goes awry in disorders of the CNS. However, the availability of tissue is certainly a challenge for systems biology of neuroscience, compared to systems biology of other tissues, where biopsy is feasible. This challenge is particularly notable for disorders caused by extremely rare genetic variants. Thus computational and systems biology approaches, as well as precise experimental models by way of genome editing, will play key roles in defining mechanisms for disorders, and their individual symptoms, across varied genetic etiologies. Here, we highlight recent progress in neurogenetics, postmortem genomics, cell-type specific profiling, and precision modeling toward defining mechanisms in disease.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0466-2) contains supplementary material, which is available to authorized users.

Phenotypically anchored transcriptome profiling of developmental exposure to the antimicrobial agent, triclosan, reveals hepatotoxicity in embryonic zebrafish

Triclosan (TCS) is an antimicrobial agent commonly found in a variety of personal care products and cosmetics. TCS readily enters the environment through wastewater and is detected in human plasma, urine, and breast milk due to its widespread use. Studies have implicated TCS as a disruptor of thyroid and estrogen signaling; therefore, research examining the developmental effects of TCS is warranted. In this study, we used embryonic zebrafish to investigate the developmental toxicity and potential mechanism of action of TCS. Embryos were exposed to graded concentrations of TCS from 6 to 120hours post-fertilization (hpf) and the concentration where 80% of the animals had mortality or morbidity at 120hpf (EC80) was calculated. Transcriptomic profiling was conducted on embryos exposed to the EC80 (7.37μM). We identified a total of 922 significant differentially expressed transcripts (FDR adjusted P-value≤0.05; fold change ≥2). Pathway and gene ontology enrichment analyses identified biological networks and transcriptional hubs involving normal liver functioning, suggesting TCS may be hepatotoxic in zebrafish. Tissue-specific gene enrichment analysis further supported the role of the liver as a target organ for TCS toxicity. We also examined the in vitro bioactivity profile of TCS reported by the ToxCast screening program. TCS had a diverse bioactivity profile and was a hit in 217 of the 385 assay endpoints we identified. We observed similarities in gene expression and hepatic steatosis assays; however, hit data for TCS were more concordant with the hypothesized CAR/PXR activity of TCS from rodent and human in vitro studies.