Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Genome-wide association study of population-standardised cognitive performance phenotypes in a rural South African community

Cognitive function is an indicator for global physical and mental health, and cognitive impairment has been associated with poorer life outcomes and earlier mortality. A standard cognition test, adapted to a rural-dwelling African community, and the Oxford Cognition Screen-Plus were used to capture cognitive performance as five continuous traits (total cognition score, verbal episodic memory, executive function, language, and visuospatial ability) for 2,246 adults in this population of South Africans. A novel common variant, rs73485231, reached genome-wide significance for association with episodic memory using data for ~14 million markers imputed from the H3Africa genotyping array data. Window-based replication of previously implicated variants and regions of interest support the discovery of African-specific associated variants despite the small population size and low allele frequency. This African genome-wide association study identifies suggestive associations with general cognition and domain-specific cognitive pathways and lays the groundwork for further genomic studies on cognition in Africa.

Multimodal bioinformatic analyses of the neurodegenerative disease-associated TECPR2 gene reveal its diverse roles

BACKGROUND

Loss of tectonin β-propeller repeat-containing 2 (TECPR2) function has been implicated in an array of neurodegenerative disorders, yet its physiological function remains largely unknown. Understanding TECPR2 function is essential for developing much needed precision therapeutics for TECPR2-related diseases.

METHODS

We leveraged considerable amounts of functional data to obtain a comprehensive perspective of the role of TECPR2 in health and disease. We integrated expression patterns, population variation, phylogenetic profiling, protein-protein interactions and regulatory network data for a minimally biased multimodal functional analysis. Genes and proteins linked to TECPR2 via multiple lines of evidence were subject to functional enrichment analyses to identify molecular mechanisms involving TECPR2.

RESULTS

TECPR2 was found to be part of a tight neurodevelopmental gene expression programme that includes KIF1A, ATXN1, TOM1L2 and FA2H, all implicated in neurological diseases. Functional enrichment analyses of TECPR2-related genes converged on a role in late autophagy and ribosomal processes. Large-scale population variation data demonstrated that this role is non-redundant.

CONCLUSIONS

TECPR2 might serve as an indicator for the energy balance between protein synthesis and autophagy, and a marker for diseases associated with their imbalance, such as Alzheimer's disease and Huntington's disease. Specifically, we speculate that TECPR2 plays an important role as a proteostasis regulator during synaptogenesis, highlighting its importance in developing neurons. By advancing our understanding of TECPR2 function, this work provides an essential stepping stone towards the development of precision diagnostics and targeted treatment options for TECPR2-related disorders.

Expression analyses of Rac3, a Rho family small GTPase, during mouse brain development

Rac3 is a member of Rho family small GTPases which regulates cellular signaling and cytoskeletal dynamics. The RAC3 gene abnormalities have been shown to cause neurodevelopmental disorders with structural brain anomalies, including polymicrogyria/dysgyria, callosal abnormalities, brainstem anomalies, and cerebellar dysplasia. Although this evidence indicates that Rac3 is essential in brain development, not only its molecular mechanism but also the expression profile is yet to be elucidated. In this study, we carried out expression analyses of Rac3 with mouse brain tissues. In immunoblotting, Rac3 exhibited a tissue-dependent expression profile in the young adult mouse and was expressed in a developmental stage-dependent manner in brain. In primary cultured hippocampal neurons, while Rac3 was distributed mainly in the cytoplasm, it was visualized in axon and dendrites with partial localization at synapses, in consistent with the observation in biochemical fractionation analyses. In immunofluorescence analyses with brain slices, Rac3 was distributed strongly and moderately in the axon and cytoplasm, respectively, of cerebral cortex at postnatal day (P) 2 and P18. Similar distribution profile was also observed in hippocampus. Taken together, the results obtained strongly suggest that Rac3 plays an important physiological role in neuronal tissues during corticogenesis, and defects in the Rac3 function induce structural brain anomalies leading to pathogenesis of neurodevelopmental disorders.

Physiological significance of WDR45, a responsible gene for β-propeller protein associated neurodegeneration (BPAN), in brain development

WDR45 plays an essential role in the early stage of autophagy. De novo heterozygous mutations in WDR45 have been known to cause β-propeller protein-associated neurodegeneration (BPAN), a subtype of neurodegeneration with brain iron accumulation (NBIA). Although BPAN patients display global developmental delay with intellectual disability, the neurodevelopmental pathophysiology of BPAN remains largely unknown. In the present study, we analyzed the physiological role of Wdr45 and pathophysiological significance of the gene abnormality during mouse brain development. Morphological and biochemical analyses revealed that Wdr45 is expressed in a developmental stage-dependent manner in mouse brain. Wdr45 was also found to be located in excitatory synapses by biochemical fractionation. Since WDR45 mutations are thought to cause protein degradation, we conducted acute knockdown experiments by in utero electroporation in mice to recapitulate the pathophysiological conditions of BPAN. Knockdown of Wdr45 caused abnormal dendritic development and synaptogenesis during corticogenesis, both of which were significantly rescued by co-expression with RNAi-resistant version of Wdr45. In addition, terminal arbors of callosal axons were less developed in Wdr45-deficient cortical neurons of adult mouse when compared to control cells. These results strongly suggest a pathophysiological significance of WDR45 gene abnormalities in neurodevelopmental aspects of BPAN.

The synaptic scaffolding protein CNKSR2 interacts with CYTH2 to mediate hippocampal granule cell development

CNKSR2 is a synaptic scaffolding molecule that is encoded by the CNKSR2 gene located on the X chromosome. Heterozygous mutations to CNKSR2 in humans are associated with intellectual disability and epileptic seizures, yet the cellular and molecular roles for CNKSR2 in nervous system development and disease remain poorly characterized. Here, we identify a molecular complex comprising CNKSR2 and the guanine nucleotide exchange factor (GEF) for ARF small GTPases, CYTH2, that is necessary for the proper development of granule neurons in the mouse hippocampus. Notably, we show that CYTH2 binding prevents proteasomal degradation of CNKSR2. Furthermore, to explore the functional significance of coexpression of CNKSR2 and CYTH2 in the soma of granule cells within the hippocampal dentate gyrus, we transduced mouse granule cell precursors in vivo with small hairpin RNAs (shRNAs) to silence CNKSR2 or CYTH2 expression. We found that such manipulations resulted in the abnormal localization of transduced cells at the boundary between the granule cell layer and the hilus. In both cases, CNKSR2-knockdown and CYTH2-knockdown cells exhibited characteristics of immature granule cells, consistent with their putative roles in neuron differentiation. Taken together, our results demonstrate that CNKSR2 and its molecular interaction partner CYTH2 are necessary for the proper development of dentate granule cells within the hippocampus through a mechanism that involves the stabilization of a complex comprising these proteins.

Transcription factors in neurodevelopmental and associated psychiatric disorders: A potential convergence for genetic and environmental risk factors

Neurodevelopmental disorders (NDDs) are a heterogeneous and highly prevalent group of psychiatric conditions marked by impairments in the nervous system. Their onset occurs during gestation, and the alterations are observed throughout the postnatal life. Although many genetic and environmental risk factors have been described in this context, the interactions between them challenge the understanding of the pathways associated with NDDs. Transcription factors (TFs)-a group of over 1,600 proteins that can interact with DNA, regulating gene expression through modulation of RNA synthesis-represent a point of convergence for different risk factors. In addition, TFs organize critical processes like angiogenesis, blood-brain barrier formation, myelination, neuronal migration, immune activation, and many others in a time and location-dependent way. In this review, we summarize important TF alterations in NDD and associated disorders, along with specific impairments observed in animal models, and, finally, establish hypotheses to explain how these proteins may be critical mediators in the context of genome-environment interactions.

Expression Analyses of Mediator Complex Subunit 13-Like: A Responsible Gene for Neurodevelopmental Disorders during Mouse Brain Development

MED13L (mediator complex subunit 13-like) is a component of the mediator complex, which functions as a regulator for gene transcription. Since gene abnormalities in MED13L are responsible for neurodevelopmental disorders, MED13L is presumed to play an essential role in brain development. In this study, we prepared a specific antibody against MED13L, anti-MED13L, and analyzed its expression profile in mouse tissues with focusing on the central nervous system. In Western blotting, MED13L exhibited a tissue-dependent expression profile in the adult mouse and was expressed in a developmental stage-dependent manner in brain. In immunofluorescence analyses, MED13L was at least partially colocalized with pre- and post-synaptic markers, synaptophysin, and PSD95, in primary cultured hippocampal neurons. Immunohistochemical analyses revealed that MED13L was relatively highly expressed in ventricular zone surface of cerebral cortex, and was also located both in the cytoplasm and nucleus of neurons in the cortical plate at embryonic day 14. Then, MED13L showed diffuse cytoplasmic distribution throughout the cerebral cortex at the postnatal day (P) 30. In addition, MED13L appeared to be localized in cell type- and developmental stage-specific manners in the hippocampus and cerebellum. These results suggest that MED13L is involved in the development of the central nervous system and synaptic function.

Behavioural Characterisation of Macrod1 and Macrod2 Knockout Mice

Adenosine diphosphate ribosylation (ADP-ribosylation; ADPr), the addition of ADP-ribose moieties onto proteins and nucleic acids, is a highly conserved modification involved in a wide range of cellular functions, from viral defence, DNA damage response (DDR), metabolism, carcinogenesis and neurobiology. Here we study MACROD1 and MACROD2 (mono-ADP-ribosylhydrolases 1 and 2), two of the least well-understood ADPr-mono-hydrolases. MACROD1 has been reported to be largely localized to the mitochondria, while the MACROD2 genomic locus has been associated with various neurological conditions such as autism, attention deficit hyperactivity disorder (ADHD) and schizophrenia; yet the potential significance of disrupting these proteins in the context of mammalian behaviour is unknown. Therefore, here we analysed both Macrod1 and Macrod2 gene knockout (KO) mouse models in a battery of well-defined, spontaneous behavioural testing paradigms. Loss of Macrod1 resulted in a female-specific motor-coordination defect, whereas Macrod2 disruption was associated with hyperactivity that became more pronounced with age, in combination with a bradykinesia-like gait. These data reveal new insights into the importance of ADPr-mono-hydrolases in aspects of behaviour associated with both mitochondrial and neuropsychiatric disorders.

Expression analyses of PLEKHG2, a Rho family-specific guanine nucleotide exchange factor, during mouse brain development

Abnormalities of PLEKHG2 gene, encoding a Rho family-specific guanine nucleotide exchange factor, are involved in microcephaly with intellectual disability. However, not only the role of PLEKHG2 in the developmental process but also its expression profile is unknown. In this study, we prepared a specific antibody against PLEKHG2 and carried out expression analyses with mouse tissues. In western blotting, PLEKHG2 exhibited a tissue-dependent expression profile in adult mouse and was expressed in a developmental stage-dependent manner in brain. Then, in immunohistochemical analyses, while PLEKHG2 was observed in the cortical plate and ventricular zone surface of the cerebral cortex at embryonic day 14, it came to be distributed throughout the cerebral cortex in layer II/III and V during corticogenesis. PLEKHG2 was also detected mainly in the nucleus of neurons in the hippocampal CA regions and dentate gyrus at P7. Notably, the nuclear accumulation disappeared at P30 and PLEKHG2 came to be located at the axons and/or dendrites at this time point. Moreover, in vitro immunofluorescence revealed that PLEKHG2 was at least partially localized at both excitatory and inhibitory synapses in primary cultured hippocampal neurons. These results suggest roles of PLEKHG2 in the development of the central nervous tissue and synaptic function.

A behavioral model for mapping the genetic architecture of gut-microbiota networks

The gut microbiota may play an important role in affecting human health. To explore the genetic mechanisms underlying microbiota-host relationships, many genome-wide association studies have begun to identify host genes that shape the microbial composition of the gut. It is becoming increasingly clear that the gut microbiota impacts host processes not only through the action of individual microbes but also their interaction networks. However, a systematic characterization of microbial interactions that occur in densely packed aggregates of the gut bacteria has proven to be extremely difficult. We develop a computational rule of thumb for addressing this issue by integrating ecological behavioral theory and genetic mapping theory. We introduce behavioral ecology theory to derive mathematical descriptors of how each microbe interacts with every other microbe through a web of cooperation and competition. We estimate the emergent properties of gut-microbiota networks reconstructed from these descriptors and map host-driven mutualism, antagonism, aggression, and altruism QTLs. We further integrate path analysis and mapping theory to detect and visualize how host genetic variants affect human diseases by perturbing the internal workings of the gut microbiota. As the proof of concept, we apply our model to analyze a published dataset of the gut microbiota, showing its usefulness and potential to gain new insight into how microbes are organized in human guts. The new model provides an analytical tool for revealing the "endophenotype" role of microbial networks in linking genotype to end-point phenotypes.

Brief Report: Genetic Links Between Autism and Suicidal Behavior-A Preliminary Investigation

Evidence suggests there may be increased risk for suicidal behavior among individuals with autism spectrum disorder (ASD). An emerging body of research explores social factors that may contribute to increased risk, however little is known about the potential role of biological factors. The current project addresses this knowledge gap through a preliminary study of genes associated with both ASD and suicidal behavior. Gene set enrichment tests of eight genes strongly associated with both ASD and suicidal behavior revealed overrepresentation of nine biological processes, including cognition and synapse function, and 14 cellular components, including the neuron, the synapse, and the synaptic and postsynaptic membrane. These results can be used to inform future investigations of the biological underpinnings of suicidal behavior and ASD.

Heterotrimeric G-protein, Gi1, is involved in the regulation of proliferation, neuronal migration, and dendrite morphology during cortical development in vivo

Heterotrimeric G-proteins are composed of α, β, and γ subunits, and function as signal transducers. Critical roles of the α-subunits of Gi/o family heterotrimeric G-proteins, Gαi2, and Gαo1, have so far been reported in brain development and neurodevelopmental disorders. In this study, we tried to clarify the role of Gαi1, α-subunit of another Gi/o family member Gi1, during corticogenesis, based on the recent identification of its gene abnormalities in neurodevelopmental disorders. In western blot analyses, Gαi1 was found to be expressed in mouse brain in a developmental stage-dependent manner. Morphological analyses revealed that Gαi1 was broadly distributed in cerebral cortex with relatively high expression in the ventricular zone (VZ) at embryonic day (E) 14. Meanwhile, Gαi1 was enriched in membrane area of yet unidentified early mitotic cells in the VZ and the marginal zone at E14. Acute knockdown of Gαi1 with in utero electroporation in cerebral cortex caused cell cycle elongation of the neural progenitor cells and promoted their cell cycle exit. Gαi1-deficient cortical neurons also exhibited delayed radial migration during corticogenesis, with abnormally elongated leading processes and hampered nucleokinesis. In addition, silencing of Gαi1 prevented basal dendrite development. The migration and dendritic phenotypes were at least partially rescued by an RNAi-resistant version of Gαi1. Collectively, these results strongly suggest a crucial role of Gi1 in cortical development, and disturbance of its function may cause deficits in synaptic network formation, leading to neurodevelopmental disorders.

Neuropathophysiological significance of the c.1449T>C/p.(Tyr64Cys) mutation in the CDC42 gene responsible for Takenouchi-Kosaki syndrome

Takenouchi-Kosaki syndrome (TKS) is an autosomal dominant congenital syndrome, of which pathogenesis is not well understood. Recently, a heterozygous mutation c.1449T > C/p.(Tyr64Cys) in the CDC42 gene, encoding a Rho family small GTPase, has been demonstrated to contribute to the TKS clinical features, including developmental delay with intellectual disability (ID). However, specific molecular mechanisms underlying the neuronal pathophysiology of TKS remain largely unknown. In this study, biochemical analyses revealed that the mutation moderately activates Cdc42. In utero electroporation-based acute expression of Cdc42-Y64C in ventricular zone progenitor cells in embryonic mice cerebral cortex resulted in migration defects and cluster formation of excitatory neurons. Expression the mutant in primary cultured hippocampal neurons caused impaired axon elongation. These data suggest that the c.1449T > C/p.(Tyr64Cys) mutation causes altered CDC42 function and results in defects in neuronal morphology and migration during brain development, which is likely to be responsible for pathophysiology of psychomotor delay and ID in TKS.

Comparative analysis of MACROD1, MACROD2 and TARG1 expression, localisation and interactome

The posttranslational modification ADP-ribosylation is involved in many cellular processes, with distinct roles for poly- and mono(ADP-ribosyl)ation (PAR- and MARylation, respectively). Reversibility of intracellular MARylation was demonstrated with the discovery of MACROD1, MACROD2 and TARG1, three macrodomain-containing enzymes capable of reversing MARylation of proteins and RNA. While the three enzymes have identical activities in vitro, their roles in cells are unclear and published data are partially contradictory, possibly due to a lack of validated reagents. We developed monoclonal antibodies to study these proteins and analysed their tissue distribution and intracellular localisation. MACROD1 is most prevalent in mitochondria of skeletal muscle, MACROD2 localises to nucleo- and cytoplasm and is found so far only in neuroblastoma cells, whereas the more ubiquitously expressed TARG1 is present in nucleoplasm, nucleolus and stress granules. Loss of MACROD1 or loss of TARG1 leads to disruption of mitochondrial or nucleolar morphology, respectively, hinting at their importance for these organelles. To start elucidating the underlying mechanisms, we have mapped their interactomes using BioID. The cellular localisation of interactors supports the mitochondrial, nucleolar and stress granule localisation of MACROD1 and TARG1, respectively. Gene ontology analysis suggests an involvement of MACROD1 and TARG1 in RNA metabolism in their respective compartments. The detailed description of the hydrolases’ expression, localisation and interactome presented here provides a solid basis for future work addressing their physiological function in more detail.

Biochemical and morphological characterization of SEPT1 in mouse brain

Septins are a highly conserved family of GTPases which are identified in diverse organisms ranging from yeast to humans. In mammals, nervous tissues abundantly contain septins and associations of septins with neurological disorders such as Alzheimer's disease and Parkinson's disease have been reported. However, roles of septins in the brain development have not been fully understood. In this study, we produced a specific antibody against mouse SEPT1 and carried out biochemical and morphological characterization of SEPT1. When the expression profile of SEPT1 during mouse brain development was analyzed by western blotting, we found that SEPT1 expression began to increase after birth and the increase continued until postnatal day 22. Subcellular fractionation of mouse brain and subsequent western blot analysis revealed the distribution of SEPT1 in synaptic fractions. Immunofluorescent analyses showed the localization of SEPT1 at synapses in primary cultured mouse hippocampal neurons. We also found the distribution of SEPT1 at synapses in mouse brain by immunohistochemistry. These results suggest that SEPT1 participates in various synaptic events such as the signaling, the neurotransmitter release, and the synapse formation/maintenance.

The Controversial Roles of ADP-Ribosyl Hydrolases MACROD1, MACROD2 and TARG1 in Carcinogenesis

Post-translational modifications (PTM) of proteins are crucial for fine-tuning a cell's response to both intracellular and extracellular cues. ADP-ribosylation is a PTM, which occurs in two flavours: modification of a target with multiple ADP-ribose moieties (poly(ADP-ribosyl)ation or PARylation) or with only one unit (MARylation), which are added by the different enzymes of the PARP family (also known as the ARTD family). PARylation has been relatively well-studied, particularly in the DNA damage response. This has resulted in the development of PARP inhibitors such as olaparib, which are increasingly employed in cancer chemotherapeutic approaches. Despite the fact that the majority of PARP enzymes catalyse MARylation, MARylation is not as well understood as PARylation. MARylation is a dynamic process: the enzymes reversing intracellular MARylation of acidic amino acids (MACROD1, MACROD2, and TARG1) were discovered in 2013. Since then, however, little information has been published about their physiological function. MACROD1, MACROD2, and TARG1 have a 'macrodomain' harbouring the catalytic site, but no other domains have been identified. Despite the lack of information regarding their cellular roles, there are a number of studies linking them to cancer. However, some of these publications oppose each other, some rely on poorly-characterised antibodies, or on aberrant localisation of overexpressed rather than native protein. In this review, we critically assess the available literature on a role for the hydrolases in cancer and find that, currently, there is limited evidence for a role for MACROD1, MACROD2, or TARG1 in tumorigenesis.

(ADP-ribosyl)hydrolases: structure, function, and biology

ADP-ribosylation is an intricate and versatile posttranslational modification involved in the regulation of a vast variety of cellular processes in all kingdoms of life. Its complexity derives from the varied range of different chemical linkages, including to several amino acid side chains as well as nucleic acids termini and bases, it can adopt. In this review, we provide an overview of the different families of (ADP-ribosyl)hydrolases. We discuss their molecular functions, physiological roles, and influence on human health and disease. Together, the accumulated data support the increasingly compelling view that (ADP-ribosyl)hydrolases are a vital element within ADP-ribosyl signaling pathways and they hold the potential for novel therapeutic approaches as well as a deeper understanding of ADP-ribosylation as a whole.

Expression Analyses of POGZ, A Responsible Gene for Neurodevelopmental Disorders, during Mouse Brain Development

POGZ is a heterochromatin protein 1 α-binding protein and regulates gene expression. On the other hand, accumulating pieces of evidence indicate that the POGZ gene abnormalities are involved in various neurodevelopmental disorders. In this study, we prepared a specific antibody against POGZ, anti-POGZ, and carried out biochemical and morphological characterization with mouse brain tissues. Western blotting analyses revealed that POGZ is expressed strongly at embryonic day 13 and then gradually decreased throughout the brain development process. In immunohistochemical analyses, POGZ was found to be enriched in cerebrocortical and hippocampal neurons in the early developmental stage. The nuclear expression was also detected in Purkinje cells in cerebellum at postnatal day (P)7 and P15 but disappeared at P30. In primary cultured hippocampal neurons, while POGZ was distributed mainly in the nucleus, it was also visualized in axon and dendrites with partial localization at synapses in consistency with the results obtained in biochemical fractionation analyses. The obtained results suggest that POGZ takes part in the regulation of synaptic function as well as gene expression during brain development.