A VNTR Regulates miR-137 Expression Through Novel Alternative Splicing and Contributes to Risk for Schizophrenia

Schizophrenia risk-associated SNPs affect expression of microRNA 137 host gene: a postmortem study

Common variants in the MicroRNA 137 host gene MIR137HG and its adjacent gene DPYD have been associated with schizophrenia risk and the latest Psychiatric Genomics Consortium (PGC). Genome-Wide Association Study on schizophrenia has confirmed and extended these findings. To elucidate the association of schizophrenia risk-associated SNPs in this genomic region, we examined the expression of both mature and immature transcripts of the miR-137 host gene (MIR137HG) in the dorsolateral prefrontal cortex (DLPFC) and subgenual anterior cingulate cortex (sgACC) of postmortem brain samples of donors with schizophrenia and psychiatrically-unaffected controls using qPCR and RNA-Seq approaches. No differential expression of miR-137, MIR137HG, or its transcripts was observed. Two schizophrenia risk-associated SNPs identified in the PGC study, rs11165917 (DLPFC: P = 2.0e-16; sgACC: P = 6.4e-10) and rs4274102 (DLPFC: P = 0.036; sgACC: P = 0.002), were associated with expression of the MIR137HG long non-coding RNA transcript MIR137HG-203 (ENST00000602672.2) in individuals of European ancestry. Carriers of the minor (risk) allele of rs11165917 had significantly lower expression of MIR137HG-203 compared with those carrying the major allele. However, we were unable to validate this result by short-read sequencing of RNA extracted from DLPFC or sgACC tissue. This finding suggests that immature transcripts of MIR137HG may contribute to genetic risk for schizophrenia.

Rediscovering tandem repeat variation in schizophrenia: challenges and opportunities

Tandem repeats (TRs) are prevalent throughout the genome, constituting at least 3% of the genome, and often highly polymorphic. The high mutation rate of TRs, which can be orders of magnitude higher than single-nucleotide polymorphisms and indels, indicates that they are likely to make significant contributions to phenotypic variation, yet their contribution to schizophrenia has been largely ignored by recent genome-wide association studies (GWAS). Tandem repeat expansions are already known causative factors for over 50 disorders, while common tandem repeat variation is increasingly being identified as significantly associated with complex disease and gene regulation. The current review summarizes key background concepts of tandem repeat variation as pertains to disease risk, elucidating their potential for schizophrenia association. An overview of next-generation sequencing-based methods that may be applied for TR genome-wide identification is provided, and some key methodological challenges in TR analyses are delineated.

Genome-wide identification of tandem repeats associated with splicing variation across 49 tissues in humans

Tandem repeats (TRs) are one of the largest sources of polymorphism, and their length is associated with gene regulation. Although previous studies reported several tandem repeats regulating gene splicing in cis (spl-TRs), no large-scale study has been conducted. In this study, we established a genome-wide catalog of 9537 spl-TRs with a total of 58,290 significant TR-splicing associations across 49 tissues (false discovery rate 5%) by using Genotype-Tissue expression (GTex) Project data. Regression models explaining splicing variation by using spl-TRs and other flanking variants suggest that at least some of the spl-TRs directly modulate splicing. In our catalog, two spl-TRs are known loci for repeat expansion diseases, spinocerebellar ataxia 6 (SCA6) and 12 (SCA12). Splicing alterations by these spl-TRs were compatible with those observed in SCA6 and SCA12. Thus, our comprehensive spl-TR catalog may help elucidate the pathomechanism of genetic diseases.

RNA epitranscriptomics dysregulation: A major determinant for significantly increased risk of ASD pathogenesis

Autism spectrum disorders (ASDs) are perhaps the most severe, intractable and challenging child psychiatric disorders. They are complex, pervasive and highly heterogeneous and depend on multifactorial neurodevelopmental conditions. Although the pathogenesis of autism remains unclear, it revolves around altered neurodevelopmental patterns and their implications for brain function, although these cannot be specifically linked to symptoms. While these affect neuronal migration and connectivity, little is known about the processes that lead to the disruption of specific laminar excitatory and inhibitory cortical circuits, a key feature of ASD. It is evident that ASD has multiple underlying causes and this multigenic condition has been considered to also dependent on epigenetic effects, although the exact nature of the factors that could be involved remains unclear. However, besides the possibility for differential epigenetic markings directly affecting the relative expression levels of individual genes or groups of genes, there are at least three mRNA epitranscriptomic mechanisms, which function cooperatively and could, in association with both genotypes and environmental conditions, alter spatiotemporal proteins expression patterns during brain development, at both quantitative and qualitative levels, in a tissue-specific, and context-dependent manner. As we have already postulated, sudden changes in environmental conditions, such as those conferred by maternal inflammation/immune activation, influence RNA epitranscriptomic mechanisms, with the combination of these processes altering fetal brain development. Herein, we explore the postulate whereby, in ASD pathogenesis, RNA epitranscriptomics might take precedence over epigenetic modifications. RNA epitranscriptomics affects real-time differential expression of receptor and channel proteins isoforms, playing a prominent role in central nervous system (CNS) development and functions, but also RNAi which, in turn, impact the spatiotemporal expression of receptors, channels and regulatory proteins irrespective of isoforms. Slight dysregulations in few early components of brain development, could, depending upon their extent, snowball into a huge variety of pathological cerebral alterations a few years after birth. This may very well explain the enormous genetic, neuropathological and symptomatic heterogeneities that are systematically associated with ASD and psychiatric disorders at large.

MicroRNA binding site variation is enriched in psychiatric disorders

Psychiatric disorders have a polygenic architecture, often associated with dozens or hundreds of independent genomic loci. Most associated loci impact noncoding regions of the genome, suggesting that the majority of disease heritability originates from the disruption of regulatory sequences. While most research has focused on variants that modify regulatory DNA elements, those affecting cis-acting RNA sequences, such as miRNA binding sites, are also likely to have a significant impact. We intersected genome-wide association study (GWAS) summary statistics with the dbMTS database of predictions for miRNA binding site variants (MBSVs). We compared the distributions of MBSV association statistics to non-MBSVs within brain-expressed 3'UTR regions. We aggregated GWAS p values at the gene, pathway, and miRNA family levels to investigate cellular functions and miRNA families strongly associated with each trait. We performed these analyses in several psychiatric disorders as well as nonpsychiatric traits for comparison. We observed significant enrichment of MBSVs in schizophrenia, depression, bipolar disorder, and anorexia nervosa, particularly in genes targeted by several miRNA families, including miR-335-5p, miR-21-5p/590-5p, miR-361-5p, and miR-557, and a nominally significant association between miR-323b-3p MBSVs and schizophrenia risk. We identified evidence for the association between MBSVs in synaptic gene sets in schizophrenia and bipolar disorder. We also observed a significant association of MBSVs in other complex traits including type 2 diabetes. These observations support the role of miRNA in the pathophysiology of psychiatric disorders and suggest that MBSVs are an important class of regulatory variants that have functional implications for many disorders, as well as other complex human traits.

A polymorphic transcriptional regulatory domain in the amyotrophic lateral sclerosis risk gene CFAP410 correlates with differential isoform expression

We describe the characterisation of a variable number tandem repeat (VNTR) domain within intron 1 of the amyotrophic lateral sclerosis (ALS) risk gene CFAP410 (Cilia and flagella associated protein 410) (previously known as C21orf2), providing insight into how this domain could support differential gene expression and thus be a modulator of ALS progression or risk. We demonstrated the VNTR was functional in a reporter gene assay in the HEK293 cell line, exhibiting both the properties of an activator domain and a transcriptional start site, and that the differential expression was directed by distinct repeat number in the VNTR. These properties embedded in the VNTR demonstrated the potential for this VNTR to modulate CFAP410 expression. We extrapolated these findings in silico by utilisation of tagging SNPs for the two most common VNTR alleles to establish a correlation with endogenous gene expression. Consistent with in vitro data, CFAP410 isoform expression was found to be variable in the brain. Furthermore, although the number of matched controls was low, there was evidence for one specific isoform being correlated with lower expression in those with ALS. To address if the genotype of the VNTR was associated with ALS risk, we characterised the variation of the CFAP410 VNTR in ALS cases and matched controls by PCR analysis of the VNTR length, defining eight alleles of the VNTR. No significant difference was observed between cases and controls, we noted, however, the cohort was unlikely to contain sufficient power to enable any firm conclusion to be drawn from this analysis. This data demonstrated that the VNTR domain has the potential to modulate CFAP410 expression as a regulatory element that could play a role in its tissue-specific and stimulus-inducible regulation that could impact the mechanism by which CFAP410 is involved in ALS.

Downregulating EVA1C exerts the potential to promote neuron growth after neonatal hypoxic-ischemic encephalopathy injury associated with alternative splicing

Neonatal hypoxic-ischemic encephalopathy (NHIE) is one of the major diseases in newborns during the perinatal stage, which globally is the main reason for children's morbidity and mortality. However, the mechanism of NHIE still remains poorly clear. In this study, the 7-day-old rats were subjected to hypoxic-ischemia (HI), then brain damage was detected. Afterward, the expression of eva-1 homolog C (EVA1C) was measured in vitro by establishing the oxygen-glucose deprivation (OGD) model in SHSY5Y cells and human fetal neurons. Subsequently, the potential function and mechanism of EVA1C were explored by silencing EVA1C and alternative splicing prediction. As a result, obvious neurobehavioral impairment and brain infarction were detected through Zea-Longa score and TTC staining; meanwhile, neuron injury was tested by HE and Nissl staining post HI. Moreover, it was found that the expression of EVA1C was notably upregulated in SHSY5Y cells and human fetal neurons after OGD. In addition, cell survival and growth were increased after silencing EVA1C, which might be associated with alternative splicing. In conclusion, EVA1C interference exhibited potential in promoting neuron survival and growth, associated with exon skipping with the alternative splicing site in 34613318:34687258, which may provide the basis for the therapeutic target and mechanism research of NHIE.

Significant Association of Variable Number Tandem Repeat Polymorphism rs58335419 in the MIR137 Gene With the Risk of Gastric and Colon Cancers

The purpose of the article: The MIR137 gene acts as a tumor-suppressor gene in colon and gastric cancers. The aim of this study was to investigate the association of functional variable number tandem repeat (VNTR) polymorphism rs58335419 locating in the upstream of the MIR137 gene with the risk of colon and gastric cancers.

Materials and methods: Totally, 429 individuals were contributed in the study, including 154 colon and 120 gastric cancer patients and 155 healthy controls. The target VNTR was genotyped using PCR and electrophoresis for all samples. Statistical analysis was performed using SPSS 21.0 software and by T, χ2 and logistic regression tests.

Results: Excluding the rare genotypes, our results showed that genotype 3/5 (95% CI = 1.08–3.73, OR = 2.01, p = 0.026) significantly increased the risk of colon cancer but not gastric cancer (95% CI = 0.88–3.30, OR = 1.70, p = 0.114). Also, in the stratification analysis for VNTRs and sex, genotypes 3/4 (95% CI = 1.00–6.07, OR = 2.46, p = 0.049) and 3/5 (95% CI = 1.25–7.18, OR = 2.99, p = 0.014) significantly increased the risk of colon cancer in men but not in women. In addition, all genotypes including the rare genotypes as a group, significantly increase the risk of gastric (95% CI = 1.14–3.00, OR = 1.85, p = 0.012) and colon (95% CI = 1.38–3.43, OR = 2.17, p = 0.001) cancers compared to the genotype 3/3 as a reference.

Conclusion: The results show that increasing the copy of VNTR in the MIR137 gene, increases the risk of colon and gastric cancers and can serve as a marker for susceptibility to colon and gastric cancers.

Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

Multifaceted Regulation of MicroRNA Biogenesis: Essential Roles and Functional Integration in Neuronal and Glial Development

MicroRNAs (miRNAs) are small, non-coding RNAs that function as endogenous gene silencers. Soon after the discovery of miRNAs, a subset of brain-enriched and brain-specific miRNAs were identified and significant advancements were made in delineating miRNA function in brain development. However, understanding the molecular mechanisms that regulate miRNA biogenesis in normal and diseased brains has become a prevailing challenge. Besides transcriptional regulation of miRNA host genes, miRNA processing intermediates are subjected to multifaceted regulation by canonical miRNA processing enzymes, RNA binding proteins (RBPs) and epitranscriptomic modifications. Further still, miRNA activity can be regulated by the sponging activity of other non-coding RNA classes, namely circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs). Differential abundance of these factors in neuronal and glial lineages partly underlies the spatiotemporal expression and function of lineage-specific miRNAs. Here, we review the continuously evolving understanding of the regulation of neuronal and glial miRNA biogenesis at the transcriptional and posttranscriptional levels and the cooperativity of miRNA species in targeting key mRNAs to drive lineage-specific development. In addition, we review dysregulation of neuronal and glial miRNAs and the detrimental impacts which contribute to developmental brain disorders.

An Insight into the microRNAs Associated with Arteriovenous and Cavernous Malformations of the Brain

Background: Brain arteriovenous malformations (BAVMs) and cerebral cavernous malformations (CCMs) are rare developmental anomalies of the intracranial vasculature, with an irregular tendency to rupture, and as of yet incompletely deciphered pathophysiology. Because of their variety in location, morphology, and size, as well as unpredictable natural history, they represent a management challenge. MicroRNAs (miRNAs) are strands of non-coding RNA of around 20 nucleotides that are able to modulate the expression of target genes by binding completely or partially to their respective complementary sequences. Recent breakthroughs have been made on elucidating their contribution to BAVM and CCM occurrence, growth, and evolution; however, there are still countless gaps in our understanding of the mechanisms involved. Methods: We have searched the Medline (PubMed; PubMed Central) database for pertinent articles on miRNAs and their putative implications in BAVMs and CCMs. To this purpose, we employed various permutations of the terms and idioms: ‘arteriovenous malformation’, ‘AVM’, and ‘BAVM’, or ‘cavernous malformation’, ‘cavernoma’, and ‘cavernous angioma’ on the one hand; and ‘microRNA’, ‘miRNA’, and ‘miR’ on the other. Using cross-reference search; we then investigated additional articles concerning the individual miRNAs identified in other cerebral diseases. Results: Seven miRNAs were discovered to play a role in BAVMs, three of which were downregulated (miR-18a, miR-137, and miR-195*) and four upregulated (miR-7-5p, miR-199a-5p, miR-200b-3p, and let-7b-3p). Similarly, eight miRNAs were identified in CCM in humans and experimental animal models, two being upregulated (miR-27a and mmu-miR-3472a), and six downregulated (miR-125a, miR-361-5p, miR-370-3p, miR-181a-2-3p, miR-95-3p, and let-7b-3p). Conclusions: The following literature review endeavored to address the recent discoveries related to the various implications of miRNAs in the formation and growth of BAVMs and CCMs. Additionally, by presenting other cerebral pathologies correlated with these miRNAs, it aimed to emphasize the potential directions of upcoming research and biological therapies.

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

Variable number tandem repeats - Their emerging role in sickness and health

Understanding the mechanisms regulating tissue specific and stimulus inducible regulation is at the heart of understanding human biology and how this translates to wellbeing, the ageing process, and disease progression. Polymorphic DNA variation is superimposed as an extra layer of complexity in such processes which underpin our individuality and are the focus of personalized medicine. This review focuses on the role and action of repetitive DNA, specifically variable number tandem repeats and SINE-VNTR-Alu domains, highlighting their role in modification of gene structure and gene expression in addition to their polymorphic nature being a genetic modifier of disease risk and progression. Although the literature focuses on their role in disease, it illustrates their potential to be major contributors to normal physiological function. To date, these elements have been under-reported in genomic analysis due to the difficulties in their characterization with short read DNA sequencing methods. However, recent advances in long read sequencing methods should resolve these problems allowing for a greater understanding of their contribution to a host of genomic and functional mechanisms underlying physiology and disease.

Translational genomics and beyond in bipolar disorder

Genome-wide association studies (GWAS) have revealed multiple genomic loci conferring risk of bipolar disorder (BD), providing hints for its underlying pathobiology. However, there are still remaining questions to answer. For example, discordance exists between BD heritability estimated with earlier epidemiological evidence and that calculated based on common GWAS variations. Where is the "missing heritability"? How can we explain the biology of the disease based on genetic findings? In this review, we summarize the accomplishments and limitations of current BD GWAS, and discuss potential reasons for the "missing heritability." In addition, progresses of research for the biological mechanisms underlying BD genetic risk using brain tissues, reprogrammed cells, and model animals are reviewed. While our knowledge of BD genetic basis is significantly promoted by these efforts, the complexities of gene regulation in the genome, the spatial-temporal heterogeneity during brain development, and the limitations of different experimental models should always be considered. Notably, several genes have been widely studied given their relatively well-characterized involvement in BD (e.g., CACAN1C and ANK3), and findings of these genes are summarized to both outline possible biological mechanisms of BD and describe examples of translating GWAS discoveries into the pathophysiology.

The MIR137 VNTR rs58335419 Is Associated With Cognitive Impairment in Schizophrenia and Altered Cortical Morphology

Genome-wide association studies (GWAS) of schizophrenia have strongly implicated a risk locus in close proximity to the gene for miR-137. While there are candidate single-nucleotide polymorphisms (SNPs) with functional implications for the microRNA's expression encompassed by the common haplotype tagged by rs1625579, there are likely to be others, such as the variable number tandem repeat (VNTR) variant rs58335419, that have no proxy on the SNP genotyping platforms used in GWAS to date. Using whole-genome sequencing data from schizophrenia patients (n = 299) and healthy controls (n = 131), we observed that the MIR137 4-repeats VNTR (VNTR4) variant was enriched in a cognitive deficit subtype of schizophrenia and associated with altered brain morphology, including thicker left inferior temporal gyrus and deeper right postcentral sulcus. These findings suggest that the MIR137 VNTR4 may impact neuroanatomical development that may, in turn, influence the expression of more severe cognitive symptoms in patients with schizophrenia.

Enhancer Function of MicroRNA-3681 Derived from Long Terminal Repeats Represses the Activity of Variable Number Tandem Repeats in the 3' UTR of SHISA7

microRNAs (miRNAs) are non-coding RNA molecules involved in the regulation of gene expression. miRNAs inhibit gene expression by binding to the 3' untranslated region (UTR) of their target gene. miRNAs can originate from transposable elements (TEs), which comprise approximately half of the eukaryotic genome and one type of TE, called the long terminal repeat (LTR) is found in class of retrotransposons. Amongst the miRNAs derived from LTR, hsa-miR-3681 was chosen and analyzed using bioinformatics tools and experimental analysis. Studies on hsa-miR-3681 have been scarce and this study provides the relative expression analysis of hsa-miR-3681-5p from humans, chimpanzees, crab-eating monkeys, and mice. Luciferase assay for hsa-miR-3681-5p and its target gene SHISA7 supports our hypothesis that the number of miRNA binding sites affects target gene expression. Especially, the variable number tandem repeat (VNTR) and hsa-miR-3681-5p share the binding sites in the 3' UTR of SHISA7, which leads the enhancer function of hsa-miR-3681-5p to inhibit the activity of VNTR. In conclusion, hsa-miR-3681-5p acts as a super-enhancer and the enhancer function of hsa-miR-3681-5p acts as a repressor of VNTR activity in the 3' UTR of SHISA7.

miR-137: A Novel Therapeutic Target for Human Glioma

MicroRNA (miR)-137 is highly expressed in the brain and plays a crucial role in the development and prognosis of glioma. In this review, we aim to summarize the latest findings regarding miR-137 in glioma cell apoptosis, proliferation, migration, invasion, angiogenesis, drug resistance, and cancer treatment. In addition, we focus on the identified miR-137 targets and pathways in the occurrence and development of glioma. Finally, future implications for the diagnostic and therapeutic potential of miR-137 in glioma were discussed.

Neurodevelopmental concepts of schizophrenia in the genome-wide association era: AKT/mTOR signaling as a pathological mediator of genetic and environmental programming during development

Normative brain development is contingent on the complex interplay between genes and environment. Schizophrenia (SCZ) is considered a highly polygenic, neurodevelopmental disorder associated with impaired neural circuit development, neurocognitive function and variations in neurotransmitter signaling systems, including dopamine. Significant evidence, accumulated over the last 30 years indicates a role for the in utero environment in SCZ pathophysiology. Emerging data suggests that changes in placental programming and function may mediate the link between genetic risk, early life complications (ELC) and adverse neurodevelopmental outcomes, with risk highlighted in key developmental drivers that converge on AKT/mTOR signaling. In this article we overview select risk genes identified through recent genome-wide association studies of SCZ including AKT3, miR-137, DRD2, and AKT1 itself. We propose that through convergence on AKT/mTOR signaling, these genes are critical factors directing both placentation and neurodevelopment, influencing risk for SCZ through dysregulation of placental function, metabolism and early brain development. We discuss association of risk genes in the context of their known roles in neurodevelopment, placental expression and their possible mechanistic links to SCZ in the broad context of the 'developmental origins of adult disease' construct. Understanding how common genetic variation impacts early fetal programming may advance our knowledge of disease etiology and identify early critical developmental windows for prevention and intervention.

A single-nucleotide polymorphism induced alternative splicing in Tacr3 involves in hypoxic-ischemic brain damage

Single-nucleotide polymorphism (SNP) and Alternative splicing (AS) were found to be implicated in certain diseases, nevertheless, the contributions of mRNA SNPs and AS to pathogenesis in developing rat brains with hypoxic-ischemic encephalopathy (HIE) remained largely vague. Additionally, the disease associated with Tacr3 was normosmic congenital hypogonadotropic hypogonadism, while the relationship between HIE and Tacr3 remained largely elusive. The current study was designed to investigate the differentially expressed mRNAs and related SNPs as well as AS in neonatal rats subjected to HIE to identify if the exhibition of AS was associated with SNPs under pathological condition. Firstly, we used postnatal day 7 Sprague-Dawley rats to construct neonatal HIE model, and analyzed the expression profiles of SNP mRNA in hypoxic-ischemic (HI) and sham brains by using RNA sequencing. Then four genes, including Mdfic, Lpp, Bag3 and Tacr3, connecting with HIE and exhibiting SNPs and AS were identified by bioinformatics analysis. Moreover, combined with exonic splicing enhancer (ESE) and alternative splice site predictor (ASSP) analysis, we found that Tacr3 is associated specifically with HIE through 258547789 G > A SNP in inside the Alt First Exon and 258548573 G > A SNP in outside the Alt First Exon. Taken together, our study provides new evidence to understand the role of Tacr3 in HIE and it is possibly a potential target for the treatment of HIE in future clinic trial.