A Human-Specific Schizophrenia Risk Tandem Repeat Affects Alternative Splicing of a Human-Unique Isoform AS3MTd2d3 and Mushroom Dendritic Spine Density

Multi-omic quantitative trait loci link tandem repeat size variation to gene regulation in human brain

Tandem repeat (TR) size variation is implicated in ~50 neurological disorders, yet its impact on gene regulation in the human brain remains largely unknown. In the present study, we quantified the impact of TR size variation on brain gene regulation across distinct molecular phenotypes, based on 4,412 multi-omics samples from 1,597 donors, including 1,586 newly sequenced ones. We identified ~2.2 million TR molecular quantitative trait loci (TR-xQTLs), linking ~139,000 unique TRs to nearby molecular phenotypes, including many known disease-risk TRs, such as the G2C4 expansion in C9orf72 associated with amyotrophic lateral sclerosis. Fine-mapping revealed ~18,700 TRs as potential causal variants. Our in vitro experiments further confirmed the causal and independent regulatory effects of three TRs. Additional colocalization analysis indicated the potential causal role of TR variation in brain-related phenotypes, highlighted by a 3'-UTR TR in NUDT14 linked to cortical surface area and a TG repeat in PLEKHA1, associated with Alzheimer's disease.

Decreased CNNM2 expression in prefrontal cortex affects sensorimotor gating function, cognition, dendritic spine morphogenesis and risk of schizophrenia

Genome-wide association studies (GWASs) have reported multiple single nucleotide polymorphisms (SNPs) associated with schizophrenia, yet the underlying molecular mechanisms are largely unknown. In this study, we aimed to identify schizophrenia relevant genes showing alterations in mRNA and protein expression associated with risk SNPs at the 10q24.32-33 GWAS locus. We carried out the quantitative trait loci (QTL) and summary data-based Mendelian randomization (SMR) analyses, using the PsychENCODE dorsolateral prefrontal cortex (DLPFC) expression QTL (eQTL) database, as well as the ROSMAP and Banner DLPFC protein QTL (pQTL) datasets. The gene CNNM2 (encoding a magnesium transporter) at 10q24.32-33 was identified to be a robust schizophrenia risk gene, and was highly expressed in human neurons according to single cell RNA-seq (scRNA-seq) data. We further revealed that reduced Cnnm2 in the mPFC of mice led to impaired cognition and compromised sensorimotor gating function, and decreased Cnnm2 in primary cortical neurons altered dendritic spine morphogenesis, confirming the link between CNNM2 and endophenotypes of schizophrenia. Proteomics analyses showed that reduced Cnnm2 level changed expression of proteins associated with neuronal structure and function. Together, these results identify a robust gene in the pathogenesis of schizophrenia.

Identification of a psychiatric risk gene NISCH at 3p21.1 GWAS locus mediating dendritic spine morphogenesis and cognitive function

BACKGROUND

Schizophrenia and bipolar disorder (BD) are believed to share clinical symptoms, genetic risk, etiological factors, and pathogenic mechanisms. We previously reported that single nucleotide polymorphisms spanning chromosome 3p21.1 showed significant associations with both schizophrenia and BD, and a risk SNP rs2251219 was in linkage disequilibrium with a human specific Alu polymorphism rs71052682, which showed enhancer effects on transcriptional activities using luciferase reporter assays in U251 and U87MG cells.

METHODS

CRISPR/Cas9-directed genome editing, real-time quantitative PCR, and public Hi-C data were utilized to investigate the correlation between the Alu polymorphism rs71052682 and NISCH. Primary neuronal culture, immunofluorescence staining, co-immunoprecipitation, lentiviral vector production, intracranial stereotaxic injection, behavioral assessment, and drug treatment were used to examine the physiological impacts of Nischarin (encoded by NISCH).

RESULTS

Deleting the Alu sequence in U251 and U87MG cells reduced mRNA expression of NISCH, the gene locates 180 kb from rs71052682, and Hi-C data in brain tissues confirmed the extensive chromatin contacts. These data suggested that the genetic risk of schizophrenia and BD predicted elevated NISCH expression, which was also consistent with the observed higher NISCH mRNA levels in the brain tissues from psychiatric patients compared with controls. We then found that overexpression of NISCH resulted in a significantly decreased density of mushroom dendritic spines with a simultaneously increased density of thin dendritic spines in primary cultured neurons. Intriguingly, elevated expression of this gene in mice also led to impaired spatial working memory in the Y-maze. Given that Nischarin is the target of anti-hypertensive agents clonidine and tizanidine, which have shown therapeutic effects in patients with schizophrenia and patients with BD in preliminary clinical trials, we demonstrated that treatment with those antihypertensive drugs could reduce NISCH mRNA expression and rescue the impaired working memory in mice.

CONCLUSIONS

We identify a psychiatric risk gene NISCH at 3p21.1 GWAS locus influencing dendritic spine morphogenesis and cognitive function, and Nischarin may have potentials for future therapeutic development.

De novo truncating NOVA2 variants affect alternative splicing and lead to heterogeneous neurodevelopmental phenotypes

Alternative splicing (AS) is crucial for cell-type-specific gene transcription and plays a critical role in neuronal differentiation and synaptic plasticity. De novo frameshift variants in NOVA2, encoding a neuron-specific key splicing factor, have been recently associated with a new neurodevelopmental disorder (NDD) with hypotonia, neurological features, and brain abnormalities. We investigated eight unrelated individuals by exome sequencing (ES) and identified seven novel pathogenic NOVA2 variants, including two with a novel localization at the KH1 and KH3 domains. In addition to a severe NDD phenotype, novel clinical features included psychomotor regression, attention deficit-hyperactivity disorder (ADHD), dyspraxia, and urogenital and endocrinological manifestations. To test the effect of the variants on splicing regulation, we transfected HeLa cells with wildtype and mutant NOVA2 complementary DNA (cDNA). The novel variants NM_002516.4:c.754_756delCTGinsTT p.(Leu252Phefs*144) and c.1329dup p.(Lys444Glnfs*82) all negatively affected AS events. The distal p.(Lys444Glnfs*82) variant, causing a partial removal of the KH3 domain, had a milder functional effect leading to an intermediate phenotype. Our findings expand the molecular and phenotypic spectrum of NOVA2-related NDD, supporting the pathogenic role of AS disruption by truncating variants and suggesting that this is a heterogeneous condition with variable clinical course.

Functional variant rs2270363 on 16p13.3 confers schizophrenia risk by regulating NMRAL1

Recent genome-wide association studies have reported multiple schizophrenia risk loci, yet the functional variants and their roles in schizophrenia remain to be characterized. Here we identify a functional single nucleotide polymorphism (rs2270363: G>A) at the schizophrenia risk locus 16p13.3. rs2270363 lies in the E-box element of the promoter of NMRAL1 and disrupts binding of the basic helix-loop-helix leucine zipper family proteins, including USF1, MAX and MXI1. We validated the regulatory effects of rs2270363 using reporter gene assays and electrophoretic mobility shift assay. Besides, expression quantitative trait loci analysis showed that the risk allele (A) of rs2270363 was significantly associated with elevated NMRAL1 expression in the human brain. Transcription factors knockdown and CRISPR-Cas9-mediated editing further confirmed the regulatory effects of the genomic region containing rs2270363 on NMRAL1. Intriguingly, NMRAL1 was significantly downregulated in the brain of schizophrenia patients compared with healthy subjects, and knockdown of Nmral1 expression affected proliferation and differentiation of mouse neural stem cells, as well as genes and pathways associated with brain development and synaptic transmission. Of note, Nmral1 knockdown resulted in significant decrease of dendritic spine density, revealing the potential pathophysiological mechanisms of NMRAL1 in schizophrenia. Finally, we independently confirmed the association between rs2270363 and schizophrenia in the Chinese population and found that the risk allele of rs2270363 was the same in European and Chinese populations. These lines of evidence suggest that rs2270363 may confer schizophrenia risk by regulating NMRAL1, a gene whose expression dysregulation might be involved in the pathogenesis of schizophrenia by affecting neurodevelopment and synaptic plasticity.

Integration of multidimensional splicing data and GWAS summary statistics for risk gene discovery

A common strategy for the functional interpretation of genome-wide association study (GWAS) findings has been the integrative analysis of GWAS and expression data. Using this strategy, many association methods (e.g., PrediXcan and FUSION) have been successful in identifying trait-associated genes via mediating effects on RNA expression. However, these approaches often ignore the effects of splicing, which can carry as much disease risk as expression. Compared to expression data, one challenge to detect associations using splicing data is the large multiple testing burden due to multidimensional splicing events within genes. Here, we introduce a multidimensional splicing gene (MSG) approach, which consists of two stages: 1) we use sparse canonical correlation analysis (sCCA) to construct latent canonical vectors (CVs) by identifying sparse linear combinations of genetic variants and splicing events that are maximally correlated with each other; and 2) we test for the association between the genetically regulated splicing CVs and the trait of interest using GWAS summary statistics. Simulations show that MSG has proper type I error control and substantial power gains over existing multidimensional expression analysis methods (i.e., S-MultiXcan, UTMOST, and sCCA+ACAT) under diverse scenarios. When applied to the Genotype-Tissue Expression Project data and GWAS summary statistics of 14 complex human traits, MSG identified on average 83%, 115%, and 223% more significant genes than sCCA+ACAT, S-MultiXcan, and UTMOST, respectively. We highlight MSG's applications to Alzheimer's disease, low-density lipoprotein cholesterol, and schizophrenia, and found that the majority of MSG-identified genes would have been missed from expression-based analyses. Our results demonstrate that aggregating splicing data through MSG can improve power in identifying gene-trait associations and help better understand the genetic risk of complex traits.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants

BACKGROUND

Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging.

METHODS

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding-disrupting SNPs (hereafter referred to as "functional SNPs") through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density.

RESULTS

We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding-disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD.

CONCLUSIONS

Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

An alternative splicing hypothesis for neuropathology of schizophrenia: evidence from studies on historical candidate genes and multi-omics data

Alternative splicing of schizophrenia risk genes, such as DRD2, GRM3, and DISC1, has been extensively described. Nevertheless, the alternative splicing characteristics of the growing number of schizophrenia risk genes identified through genetic analyses remain relatively opaque. Recently, transcriptomic analyses in human brains based on short-read RNA-sequencing have discovered many "local splicing" events (e.g., exon skipping junctions) associated with genetic risk of schizophrenia, and further molecular characterizations have identified novel spliced isoforms, such as AS3MT^d2d3 and ZNF804A^E3E4. In addition, long-read sequencing analyses of schizophrenia risk genes (e.g., CACNA1C and NRXN1) have revealed multiple previously unannotated brain-abundant isoforms with therapeutic potentials, and functional analyses of KCNH2-3.1 and Ube3a1 have provided examples for investigating such spliced isoforms in vitro and in vivo. These findings suggest that alternative splicing may be an essential molecular mechanism underlying genetic risk of schizophrenia, however, the incomplete annotations of human brain transcriptomes might have limited our understanding of schizophrenia pathogenesis, and further efforts to elucidate these transcriptional characteristics are urgently needed to gain insights into the illness-correlated brain physiology and pathology as well as to translate genetic discoveries into novel therapeutic targets.

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.

Proteome-wide Association Study Provides Insights Into the Genetic Component of Protein Abundance in Psychiatric Disorders

BACKGROUND

Genome-wide association studies have identified multiple risk variants for psychiatric disorders. Nevertheless, how the risk variants confer risk of psychiatric disorders remains largely unknown.

METHODS

We performed proteome-wide association studies to identify genes whose cis-regulated protein abundance change in the human brain were associated with psychiatric disorders.

RESULTS

By integrating genome-wide associations of four common psychiatric disorders and two independent brain proteomes (n = 376 and n = 152, respectively) from the dorsolateral prefrontal cortex, we identified 61 genes (including 48 genes for schizophrenia, 12 genes for bipolar disorder, 5 genes for depression, and 2 genes for attention-deficit/hyperactivity disorder) whose genetically regulated protein abundance levels were associated with risk of psychiatric disorders. Comparison with transcriptome-wide association studies identified 18 overlapping genes that showed significant associations with psychiatric disorders at both proteome-wide and transcriptome-wide levels, suggesting that genetic risk variants likely confer risk of psychiatric disorders by regulating messenger RNA expression and protein abundance of these genes.

CONCLUSIONS

Our study not only provides new insights into the genetic component of protein abundance in psychiatric disorders but also highlights several high-confidence risk proteins (including CNNM2 and CTNND1) for schizophrenia and depression. These high-confidence risk proteins represent promising therapeutic targets for future drug development.

AS3MT Polymorphism: A Risk Factor for Epilepsy Susceptibility and Adverse Drug Reactions to Valproic Acid and Oxcarbazepine Treatment in Children From South China

Epilepsy is a common neurologic disorder characterized by intractable seizures, involving genetic factors. There is a need to develop reliable genetic markers to predict the risk of epilepsy and design effective therapies. Arsenite methyltransferase (AS3MT) catalyzes the biomethylation of arsenic and hence regulates arsenic metabolism. AS3MT variation has been linked to the progression of various diseases including schizophrenia and attention deficit or hyperactivity disorder. Whether genetic polymorphism of AS3MT contributes to epilepsy remains unclear. In this study, we investigated the association of AS3MT gene polymorphism with susceptibility to epilepsy in children from south China. We also explored the effect of AS3MT variation on the safety of antiepileptic drugs. Genotypic analysis for AS3MT rs7085104 was performed using samples from a Chinese cohort of 200 epileptic children and 244 healthy individuals. The results revealed a genetic association of AS3MT rs7085104 with susceptibility to pediatric epilepsy. Mutant homozygous GG genotype exhibited a lower susceptibility to childhood epilepsy than AA genotype. Carriers of AS3MT rs7085104 AA genotype exhibited a higher risk of digestive adverse drug reactions (dADRs) in children when treated with valproic acid (VPA) or oxcarbazepine (OXC). Additionally, bioinformatics analysis identified eight AS3MT target genes related to epilepsy and three AS3MT-associated genes in VPA-related dADRs. The effects of AS3MT on epilepsy might involve multiple targets including CNNM2, CACNB2, TRIM26, MTHFR, GSTM1, CYP17A1, NT5C2, and YBX3. This study reveals that AS3MT may be a new gene contributing to epileptogenesis. Hence, analysis of AS3MT polymorphisms will help to evaluate susceptibility to pediatric epilepsy and drug safety.

Joint-Tissue Integrative Analysis Identified Hundreds of Schizophrenia Risk Genes

Genome-wide association studies (GWAS) have identified a large number of schizophrenia risk variants, and most of them are mapped to noncoding regions. By leveraging multiple joint-tissue gene expression data and GWAS data, we herein performed a transcriptome-wide association study (TWAS) and Mendelian randomization (MR) analysis and identified 144 genes whose mRNA levels were related to genetic risk of schizophrenia. Most of these genes exhibited diametrically opposite trends of expression in prenatal and postnatal brain tissues, despite that their expression levels in dorsolateral prefrontal cortex (DLPFC) tissues did not significantly differ between schizophrenics and healthy controls. We then found significant enrichment of these genes in dopamine-related pathways that were repeatedly implicated in schizophrenia pathogenesis and in the action of antipsychotic drugs. Gene expression analysis using single cell RNA-sequencing (scRNA-seq) data of mid-gestation fetal brains further revealed enrichment of these genes in glutamatergic excitatory neurons and cycling progenitors. These lines of evidence, in consistency with previous findings, confirmed the polygenic nature of schizophrenia and highlighted involvement of early neurodevelopment aberrations in this disorder. Further investigations using advanced algorithms in both bulk brain tissues and in single cells and at different developmental stages are necessary to characterize transcriptomic features of schizophrenia pathogenesis along brain development.

Downregulation by CNNM2 of ATP5MD expression in the 10q24.32 schizophrenia-associated locus involved in impaired ATP production and neurodevelopment

Genome-wide association studies (GWAS) have accelerated the discovery of numerous genetic variants associated with schizophrenia. However, most risk variants show a small effect size (odds ratio (OR) <1.2), suggesting that more functional risk variants remain to be identified. Here, we employed region-based multi-marker analysis of genomic annotation (MAGMA) to identify additional risk loci containing variants with large OR value from Psychiatry Genomics Consortium (PGC2) schizophrenia GWAS data and then employed summary-data-based mendelian randomization (SMR) to prioritize schizophrenia susceptibility genes. The top-ranked susceptibility gene ATP5MD, encoding an ATP synthase membrane subunit, is observed to be downregulated in schizophrenia by the risk allele of CNNM2-rs1926032 in the schizophrenia-associated 10q24.32 locus. The Atp5md knockout (KO) in mice was associated with abnormal startle reflex and gait, and ATP5MD knockdown (KD) in human induced pluripotent stem cell-derived neurons disrupted the neural development and mitochondrial respiration and ATP production. Moreover, CNNM2-rs1926032 KO could induce downregulation of ATP5MD expression and disruptions of mitochondrial respiration and ATP production. This study constitutes an important mechanistic component that links schizophrenia-associated CNNM2 regions to disruption in energy adenosine system modulation and neuronal function by long-distance chromatin domain downregulation of ATP5MD. This pathogenic mechanism provides therapeutic implications for schizophrenia.

Transcriptomic analyses of humans and mice provide insights into depression

Accumulating studies have been conducted to identify risk genes and relevant biological mechanisms underlying major depressive disorder (MDD). In particular, transcriptomic analyses in brain regions engaged in cognitive and emotional processes, e.g., the dorsolateral prefrontal cortex (DLPFC), have provided essential insights. Based on three independent DLPFC RNA-seq datasets of 79 MDD patients and 75 healthy controls, we performed differential expression analyses using two alternative approaches for cross-validation. We also conducted transcriptomic analyses in mice undergoing chronic variable stress (CVS) and chronic social defeat stress (CSDS). We identified 12 differentially expressed genes (DEGs) through both analytical methods in MDD patients, the majority of which were also dysregulated in stressed mice. Notably, the mRNA level of the immediate early gene FOS ( Fos proto-oncogene) was significantly decreased in both MDD patients and CVS-exposed mice, and CSDS-susceptible mice exhibited a greater reduction in Fos expression compared to resilient mice. These findings suggest the potential key roles of this gene in the pathogenesis of MDD related to stress exposure. Altered transcriptomes in the DLPFC of MDD patients might be, at least partially, the result of stress exposure, supporting that stress is a primary risk factor for MDD.

Functional Genomics Identify a Regulatory Risk Variation rs4420550 in the 16p11.2 Schizophrenia-Associated Locus

BACKGROUND

Genome-wide association studies (GWASs) have reported hundreds of genomic loci associated with schizophrenia, yet identifying the functional risk variations is a key step in elucidating the underlying mechanisms.

METHODS

We applied multiple bioinformatics and molecular approaches, including expression quantitative trait loci analyses, epigenome signature identification, luciferase reporter assay, chromatin conformation capture, homology-directed genome editing by CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/Cas9), RNA sequencing, and ATAC-Seq (assay for transposase-accessible chromatin using sequencing).

RESULTS

We found that the schizophrenia GWAS risk variations at 16p11.2 were significantly associated with messenger RNA levels of multiple genes in human brain, and one of the leading expression quantitative trait loci genes, MAPK3, is located ∼200 kb away from these risk variations in the genome. Further analyses based on the epigenome marks in human brain and cell lines suggested that a noncoding single nucleotide polymorphism, rs4420550 (p = 2.36 × 10^-9 in schizophrenia GWAS), was within a DNA enhancer region, which was validated via in vitro luciferase reporter assays. The chromatin conformation capture experiment showed that the rs4420550 region physically interacted with the MAPK3 promoter and TAOK2 promoter. Precise CRISPR/Cas9 editing of a single base pair in cells followed by RNA sequencing further confirmed the regulatory effects of rs4420550 on the transcription of 16p11.2 genes, and ATAC-Seq demonstrated that rs4420550 affected chromatin accessibility at the 16p11.2 region. The rs4420550-[A/A] cells showed significantly higher proliferation rates compared with rs4420550-[G/G] cells.

CONCLUSIONS

These results together suggest that rs4420550 is a functional risk variation, and this study illustrates an example of comprehensive functional characterization of schizophrenia GWAS risk loci.

Functional Analysis of Haplotypes in Bovine PSAP Gene and Their Relationship with Beef Cattle Production Traits

Simple Summary

With the rapid development of information technology and molecular biotechnology, animal molecular breeding technology is playing an increasingly important role in beef cattle breeding. Prosaposin (PSAP) is involved in regulating the growth and development of animals, and it is reported that PSAP is an important marker-assisted selection (MAS) in cattle herd. The purpose of this study was to explore the novel variants in 3’ UTR of cattle PSAP and evaluate their effects on the morphological traits of four Chinese cattle breeds. In this study, 13 variants were identified in the PSAP 3’ UTR from 501 individuals belonging to four cattle breeds. In Nanyang cattle, the distribution of haplotypes was different from the other three breeds. Two groups of haplotypes had association with morphological traits by changing the secondary structures of PSAP 3’ UTR rather than the miR-184 target sites. This study not only expands the genetic variation spectrum of cattle PSAP but also contributes to MAS genetics and breeding of Chinese cattle breeds.

Abstract

The purpose of this study was to explore functional variants in the prosaposin (PSAP) three prime untranslated region (3’ UTR) and clarify the relationship between the variants and morphological traits. Through Sanger sequencing, 13 variations were identified in bovine PSAP in four Chinese cattle breeds, with six of them being loci in 3’ UTR. In particular, Nanyang (NY) cattle had a special genotype and haplotype distribution compared to the other three breeds. NY cattle with ACATG and GCGTG haplotypes had higher morphological traits than GTACA and GTACG haplotypes. The results of dual-luciferase reporter assay showed that ACATG and GCGTG haplotypes affected the morphological traits of NY cattle by altering the secondary structure of PSAP 3’ UTR rather than the miR-184 target sites. The findings of this study could be an evidence of a complex and varying mechanism between variants and animal morphological traits and could be used to complement candidate genes for molecular breeding.