Genetic variants in miRNAs differentially expressed during brain development and their relevance to psychiatric disorders susceptibility

OBJECTIVES

MicroRNAs (miRNAs) play an important regulatory role in the expression of genes involved in brain functions during development. Genetic variants in miRNA genes may impact their regulatory function and lead to psychiatric disorders. To evaluate the role of genetic variants in genes of miRNAs differentially expressed during neurodevelopment on autism spectrum disorder (ASD), attention deficit/hyperactivity disorder (ADHD), schizophrenia (SCZ), and major depressive disorder (MDD).

METHODS

The miRNAs were identified in the literature. Summary statistics from the most recent genome-wide association studies to date were used to evaluate the association between the selected polymorphisms and each disorder in a look-up approach. In a global analysis, we compared the standardised risk effect of variants in neurodevelopment-related miRNAs with those in the remaining miRNAs from miRBase.

RESULTS

The global analysis showed that variants in neurodevelopment-related miRNAs had higher risk effects compared to the other miRNAs for SCZ (p = 0.010) and ADHD (p = 0.001). MIR33B, MIR29B2, MIR29C, MIR137, and MIR135A1 were significantly associated with SCZ, while 55.9% of the miRNAs were at least nominally associated with one or more psychiatric disorders (p < 0.05).

CONCLUSIONS

Genetic variants in neurodevelopment-related miRNAs play an important role in the genetic susceptibility of psychiatric disorders, mainly SCZ and ADHD.

Comparative Genomics Identifies Putative Interspecies Mechanisms Underlying Crbn-Sall4-Linked Thalidomide Embryopathy

The identification of thalidomide–Cereblon-induced SALL4 degradation has brought new understanding for thalidomide embryopathy (TE) differences across species. Some questions, however, regarding species variability, still remain. The aim of this study was to detect sequence divergences between species, affected or not by TE, and to evaluate the regulated gene co-expression in a murine model. Here, we performed a comparative analysis of proteins experimentally established as affected by thalidomide exposure, evaluating 14 species. The comparative analysis, regarding synteny, neighborhood, and protein conservation, was performed in 42 selected genes. Differential co-expression analysis was performed, using a publicly available assay, GSE61306, which evaluated mouse embryonic stem cells (mESC) exposed to thalidomide. The comparative analyses evidenced 20 genes in the upstream neighborhood of NOS3, which are different between the species who develop, or not, the classic TE phenotype. Considering protein sequence alignments, RECQL4, SALL4, CDH5, KDR, and NOS2 proteins had the biggest number of variants reported in unaffected species. In co-expression analysis, Crbn was a gene identified as a driver of the co-expression of other genes implicated in genetic, non-teratogenic, limb reduction defects (LRD), such as Tbx5, Esco2, Recql4, and Sall4; Crbn and Sall4 were shown to have a moderate co-expression correlation, which is affected after thalidomide exposure. Hence, even though the classic TE phenotype is not identified in mice, a deregulatory Crbn-induced mechanism is suggested in this animal. Functional studies are necessary, especially evaluating the genes responsible for LRD syndromes and their interaction with thalidomide–Cereblon.

Methylphenidate effects on mice odontogenesis and connections with human odontogenesis

The purpose of this study is to evaluate the effects of Methylphenidate exposure on mice odontogenesis and connect them by bioinformatics with human odontogenesis. Thirty-two pregnant Swiss mice were divided into treated group and control group, which received, respectively, 5 mg/kg of Methylphenidate and saline solution from the 5th to the 17th day of pregnancy. The mouse embryos tooth germs were analyzed through optical microscopy, and the data collected were analyzed statistically by Fisher's exact test. The presence and similarity of Methylphenidate-associated genes (Pharmgkb database) in both organisms and their interaction with dental development genes (AmiGO2 database) were verified on STRING database. Rates of tooth germ malformations were higher in treated than in control group (Control: 18; Treated: 27; p = 0.035). Mouse embryo malformations were connected with 238 interactions between 69 dental development genes with 35 Methylphenidate genes. Fourteen interactions for four Methylphenidate genes with four dental development genes, with human experimental data, were connected with mouse phenotype data. By homology, the interactions and conservation of proteins/genes may indicate similar outcomes for both organisms. The exposure to Methylphenidate during pregnancy affected odontogenesis in mouse embryos and may affect human odontogenesis. The study of malformations in mice, with a bioinformatics approach, could contribute to understanding of the Methylphenidate effect on embryo development. These results may provide novel hypotheses for further testing and reinforce the FDA protocol: as Methylphenidate is included in category C, its use during pregnancy should be considered if the benefits outweigh the risks.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

FOXP in Tetrapoda: Intrinsically Disordered Regions, Short Linear Motifs and their evolutionary significance

The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

The role of protein intrinsic disorder in major psychiatric disorders

Although new candidate genes for Autism Spectrum Disorder (ASD), Schizophrenia (SCZ), Attention-Deficit/Hyperactivity Disorder (ADHD), and Bipolar Disorder (BD) emerged from genome-wide association studies (GWAS), their underlying molecular mechanisms remain poorly understood. Evidences of the involvement of intrinsically disordered proteins in diseases have grown in the last decade. These proteins lack tridimensional structure under physiological conditions and are involved in important cellular functions such as signaling, recognition and regulation. The aim of the present study was to identify the role and abundance of intrinsically disordered proteins in a set of psychiatric diseases and to test whether diseases are different regarding protein intrinsic disorder. Our hypothesis is that differences across psychiatric illnesses phenotypes and symptoms may arise from differences in intrinsic protein disorder content and properties of each group. A bioinformatics prediction of intrinsic disorder was performed in proteins retrieved based on top findings from GWAS, Copy Number Variation and candidate gene investigations for each disease. This approach revealed that about 80% of studied proteins presented long stretches of disorder. This amount was significantly higher than that observed in general eukaryotic proteins, and those involved in cardiovascular diseases. These results suggest that proteins with intrinsic disorder are a common feature of neurodevelopment and synaptic transmission processes which are potentially involved in the etiology of psychiatric diseases. Moreover, we identified differences between ADHD and ASD when the binary prediction of structure and putative binding sites were compared. These differences may be related to variation in symptom complexity between both diseases. © 2016 Wiley Periodicals, Inc.

The contribution of protein intrinsic disorder to understand the role of genetic variants uncovered by autism spectrum disorders exome studies

Several autism spectrum disorders (ASD) exome studies suggest that coding single nucleotide variants (SNVs) play an important role on ASD etiology. Usually, the pathogenic effect of missense mutations is estimated through predictors that lose accuracy for those SNVs placed in intrinsically disordered regions of protein. Here, we used bioinformatics tools to investigate the effect of mutations described in ASD published exome studies (549 mutations) in protein disorder, considering post-translational modification, PEST and Molecular Recognition Features (MoRFs) motifs. Schizophrenia and type 2 diabetes (T2D) datasets were created for comparison purposes. The frequency of mutations predicted as disordered was comparable among the three datasets (38.1% in ASD, 35.7% in schizophrenia, 46.4% in T2D). However, the frequency of SNVs predicted to lead a gain or loss of functional sites or change intrinsic disorder tendencies was higher in ASD and schizophrenia than T2D (46.9%, 36.4%, and 23.1%, respectively). The results obtained by SIFT and PolyPhen-2 indicated that 38.9% and 34.4% of the mutations predicted, respectively, as tolerated and benign showed functional alterations in disorder properties. Given the frequency of mutations placed in IDRs and their functional impact, this study suggests that alterations in intrinsic disorder properties might play a role in ASD and schizophrenia etiologies. They should be taken into consideration when researching the pathogenicity of mutations in neurodevelopmental and psychiatric diseases. Finally, mutations with functional alterations in disorder properties must be potential targets for in vitro and in vivo functional studies.

Ocular and craniofacial phenotypes in a large Brazilian family with congenital aniridia

Congenital aniridia is a rare genetic disorder characterized by varying degrees of iris hypoplasia that are associated with additional ocular abnormalities. More than 90% of the causal mutations identified are found in the PAX6 gene, a transcription factor of critical importance in the process of neurogenesis and ocular development. Here, we investigate clinical, molecular, and craniofacial features of a large Brazilian family with congenital aniridia. Among the 56 eyes evaluated, phenotype variation encompassed bilateral total aniridia to mild iris defects with extensive variation between eyes of the same individual. PAX6 molecular screening indicated a heterozygous splice mutation (c.141 + 1G>A). Thus, we hypothesize that this splicing event may cause variation in the expression of the wild-type transcript, which may lead to the observed variation in phenotype. Affected individuals were more brachycephalic, even though their face height and cephalic circumference were not significantly different when compared to those of non-affected relatives. From this, we infer that the head shape of affected subjects may also be a result of the PAX6 splice-site mutation. Our data summarize the clinical variability associated with the ocular phenotype in a large family with aniridia, and help shed light on the role of PAX6 in neurocranial development.

Evolutionary history of chordate PAX genes: dynamics of change in a complex gene family

Paired box (PAX) genes are transcription factors that play important roles in embryonic development. Although the PAX gene family occurs in animals only, it is widely distributed. Among the vertebrates, its 9 genes appear to be the product of complete duplication of an original set of 4 genes, followed by an additional partial duplication. Although some studies of PAX genes have been conducted, no comprehensive survey of these genes across the entire taxonomic unit has yet been attempted. In this study, we conducted a detailed comparison of PAX sequences from 188 chordates, which revealed restricted variation. The absence of PAX4 and PAX8 among some species of reptiles and birds was notable; however, all 9 genes were present in all 74 mammalian genomes investigated. A search for signatures of selection indicated that all genes are subject to purifying selection, with a possible constraint relaxation in PAX4, PAX7, and PAX8. This result indicates asymmetric evolution of PAX family genes, which can be associated with the emergence of adaptive novelties in the chordate evolutionary trajectory.

The role of a mineralocorticoid receptor gene functional polymorphism in the symptom dimensions of persistent ADHD

Attention-deficit/hyperactivity disorder (ADHD) affects approximately 5 % of school-aged children and 2.5 % of adults. Genetic studies in ADHD have pointed to genes in different neurobiological systems, with relatively small individual effects. The mineralocorticoid receptor is the main receptor involved in the initial triggering of stress response. Therefore, its encoding gene (NR3C2) is a candidate for psychiatric disorder studies, including ADHD, and behavioral phenotypes. There is evidence that the Val allele of the MRI180V polymorphism (rs5522) increases the risk of depression, attention and cognitive deficits. We investigated the possible role of the mineralocorticoid receptor gene in the symptom dimensions and susceptibility to persistent ADHD. We compared genotype and allele frequencies in 478 adult patients with ADHD and 597 controls and symptom dimensions in 449 patients and 132 controls. Diagnoses were based on the DSM-IV criteria. ADHD symptom dimensions were investigated with SNAP-IV for ADHD severity and Barkley scales for severity and impairment. Carriers of the Val allele presented higher inattention, hyperactivity/impulsivity and impairment scores, while genotype and allele frequencies did not differ between patients and controls. These results are consistent with a possible link between genetic variations in the HPA axis and inattention and hyperactivity measures.

Molecular evolution of the melanocortin 1-receptor pigmentation gene in rodents

Adaptive variation in the melanocortin 1-receptor gene (MC1R), a key locus in melanogenesis, has been identified in some species of rodents. However, in others, MC1R has no causative role in pigmentation phenotypes despite their coat color variation. In this study, we characterized the rates and patterns of MC1R nucleotide and amino acid sequence evolution and, particularly, selective pressures in the separated domains of the protein using a comparative analysis of 43 species representing three major lineages of rodents with variable coat colors. We found high amino acid variation (44% of sites) throughout the protein. Most substitutions were observed in extracellular and transmembrane domains; the intracellular segment was conserved across species. Pairwise non-synonymous substitutions did not vary significantly in different domains among the rodent lineages - i.e., variation was not associated with phylogenetic distance. Phylogeny-based likelihood analysis suggested that purifying selection has mostly shaped the evolutionary course of MC1R. However, a high proportion of sites (27%) were under relaxation of functional constraints (ω = 0.38), and four sites (3, 14, 26, and 251) clearly evolved under positive selection (ω ≅ 2.9). Thus, our data indicate a high proportion of sites evolving under relaxed evolutionary constraints, which might indicate the evolvability of the system in the generation of adaptive changes in specific taxa in rodent lineages.

The TP53 fertility network

The TP53 gene, first described in 1979, was identified as a tumor suppressor gene in 1989, when it became clear that its product, the p53 nuclear phosphoprotein, was frequently inactivated in many different forms of cancers. Nicknamed "guardian of the genome", TP53 occupies a central node in stress response networks. The p53 protein has a key role as transcription factor in limiting oncogenesis through several growth suppressive functions, such as initiating apoptosis, senescence, or cell cycle arrest. The p53 protein is directly inactivated in about 50% of all tumors as a result of somatic gene mutations or deletions, and over 80% of tumors demonstrate dysfunctional p53 signaling. Beyond the undeniable importance of p53 as a tumor suppressor, an increasing number of new functions for p53 have been reported, including its ability to regulate energy metabolism, to control autophagy, and to participate in various aspects of differentiation and development. Recently, studies on genetic variations in TP53 among different populations have led to the notion that the p53 protein might play an important role in regulating fertility. This review summarizes current knowledge on the basic functions of different genes of the TP53 family and TP53 pathway with respect to fertility. We also provide original analyses based on genomic and genotype databases, providing further insights into the possible roles of the TP53 pathway in human reproduction.

Linking dopamine neurotransmission and neurogenesis: The evolutionary history of the NTAD (NCAM1-TTC12-ANKK1-DRD2) gene cluster

Genetic studies have long suggested the important role of the DRD2 gene in psychiatric disorders and behavior. Further research has shown a conjoined effect of genes in the Chr11q22–23 region, which includes the NCAM1, TTC12, ANKK1 and DRD2 genes, or NTAD cluster. Despite a growing need to unravel the role of this cluster, few studies have taken into account interspecies and evolutionary approaches. This study shows that behaviorally relevant SNPs from the NTAD cluster, such as rs1800497 (Taq1A) and rs6277, are ancient polymorphisms that date back to the common ancestor between modern humans and Neanderthals/Denisovans. Conserved synteny and neighborhood indicate the NTAD cluster seems to have been established at least 400 million years ago, when the first Sarcopterygians emerged. The NTAD genes are apparently co-regulated and this could be attributed to adaptive functional properties, including those that emerged when the central nervous system became more complex. Finally, our findings indicate that NTAD genes, which are related to neurogenesis and dopaminergic neurotransmission, should be approached as a unit in behavioral and psychiatric genetic studies.

Homo sapiens, Homo neanderthalensis and the Denisova specimen: New insights on their evolutionary histories using whole-genome comparisons

After a brief review of the most recent findings in the study of human evolution, an extensive comparison of the complete genomes of our nearest relative, the chimpanzee (Pan troglodytes), of extant Homo sapiens, archaic Homo neanderthalensis and the Denisova specimen were made. The focus was on non-synonymous mutations, which consequently had an impact on protein levels and these changes were classified according to degree of effect. A total of 10,447 non-synonymous substitutions were found in which the derived allele is fixed or nearly fixed in humans as compared to chimpanzee. Their most frequent location was on chromosome 21. Their presence was then searched in the two archaic genomes. Mutations in 381 genes would imply radical amino acid changes, with a fraction of these related to olfaction and other important physiological processes. Eight new alleles were identified in the Neanderthal and/or Denisova genetic pools. Four others, possibly affecting cognition, occured both in the sapiens and two other archaic genomes. The selective sweep that gave rise to Homo sapiens could, therefore, have initiated before the modern/archaic human divergence.

Evolutionary responses to a constructed niche: ancient Mesoamericans as a model of gene-culture coevolution

Culture and genetics rely on two distinct but not isolated transmission systems. Cultural processes may change the human selective environment and thereby affect which individuals survive and reproduce. Here, we evaluated whether the modes of subsistence in Native American populations and the frequencies of the ABCA1*Arg230Cys polymorphism were correlated. Further, we examined whether the evolutionary consequences of the agriculturally constructed niche in Mesoamerica could be considered as a gene-culture coevolution model. For this purpose, we genotyped 229 individuals affiliated with 19 Native American populations and added data for 41 other Native American groups (n = 1905) to the analysis. In combination with the SNP cluster of a neutral region, this dataset was then used to unravel the scenario involved in 230Cys evolutionary history. The estimated age of 230Cys is compatible with its origin occurring in the American continent. The correlation of its frequencies with the archeological data on Zea pollen in Mesoamerica/Central America, the neutral coalescent simulations, and the F_ST-based natural selection analysis suggest that maize domestication was the driving force in the increase in the frequencies of 230Cys in this region. These results may represent the first example of a gene-culture coevolution involving an autochthonous American allele.

Glutamate and synaptic plasticity systems and smoking behavior: results from a genetic association study

Smoking behavior is a multifactorial phenotype with significant heritability. Identifying the specific loci that influence smoking behavior could provide important etiological insights and facilitate the development of treatments to further reduce smoking related mortality. Although several studies pointed to different candidate genes for smoking, there is still a need for replication especially in samples from different countries. In the present study, we investigated whether 21 positive signals for smoking behavior from these studies are replicated in a sample of 531 blood donors from the Brazilian population. The polymorphisms were chosen based on their representativeness of different candidate biologic systems, strength of previous evidence, location and allele frequencies. By genotyping with the Sequenom MassARRAY iPLEX platform and subsequent statistical analysis using Plink software, we show that two of the SNPs studied, in the SLC1A2 (rs1083658) and ACTN1 (rs2268983) genes, were associated with smoking behavior in our study population. These genes are involved in crucial aspects of nicotine dependence, glutamate system and synaptic plasticity, and as such, are biologically plausible candidates that merit further molecular analyses so as to clarify their potential role in smoking behavior.

PAX9 and MSX1 transcription factor genes in non-syndromic dental agenesis

OBJECTIVE

The molecular variation of paired domain box gene 9 (PAX9) was previously investigated by our research group and a high degree of evolutionary conservation in coding and non-coding regions was observed except in exon 3. PAX9 is a transcription factor important in tooth development, and we wanted to verify its role in dental agenesis in detail. Since dental development is a complex trait we also decided to examine the influence of another transcription factor, muscle segment homeodomain-homeobox 1 (MSX1) on it.

DESIGN

A total of 360 consecutively ascertained patients seeking orthodontic treatment were screened for tooth agenesis and 33% of them were found to have it. Thirty-five of those with agenesis and 15 controls had their DNA studied for PAX9 exons 2, 3, 4 and adjacent regions (total of 1476 base pairs, bp) as well as MSX1 exon 2 (698bp). A trio (a proband and her parents) was also studied.

RESULTS

Six polymorphic sites were found, three in PAX9 exon 3 and three in MSX1 exon2. MSX1 rs1095 derived allele occurred in individuals with agenesis only, and two other mutations in this gene had been earlier associated with tooth agenesis. Homozygosity for the PAX9 Ala240Pro mutation was studied in a family (proband and her parents), suggesting recessive inheritance with variable expressivity for the dental agenesis found.

CONCLUSION

Common variants located out of the DNA binding domain of the two PAX9 and MSX1 genes can also be related to tooth agenesis.