Comprehensive functional annotation of susceptibility SNPs prioritized 10 genes for schizophrenia

Subjective sleep assessment in individuals with SYNGAP1-associated syndrome

STUDY OBJECTIVES

Sleep disturbances are common in neurodevelopmental disorders, affecting patients and caregivers' quality of life. SYNGAP1-associated syndrome, a rare neurodevelopmental disorder, is marked by intellectual disability, developmental delay, epilepsy, and sleep issues. However, research on sleep quality in these individuals is limited. This study aimed to evaluate genetic variants, epilepsy, and sleep patterns in SYNGAP1-associated syndrome patients and their caregivers.

METHODS

An online survey was applied to 11 caregivers of individuals diagnosed with SYNGAP1-associated syndrome. Specific clinical inquiries were included, addressing childbirth, previous surgeries, and medication use. Inquiries about epilepsy included type of epilepsy, type and frequency of seizures, antiseizure medications, and complementary nonpharmacological treatments. Children's Sleep Habits Questionnaire was applied to assess the patients' sleep profile. Pittsburgh Sleep Quality Index was used to evaluate the sleep quality of caregivers.

RESULTS

Genetic analysis showed heterozygous mutations in SYNGAP1, often leading to loss of function. Epilepsy was present in 82% of participants, with 77.8% having drug-resistant seizures. Using the Children's Sleep Habits Questionnaire, 81.8% of patients exhibited poor sleep habits, including bedtime resistance, anxiety, night awakenings, parasomnias, and daytime sleepiness. Caregivers also reported poor sleep quality according to the Pittsburgh Sleep Quality Index.

CONCLUSIONS

This study highlights the high prevalence of epilepsy and sleep problems in SYNGAP1-associated syndrome, impacting both patients and caregivers. Further research is crucial to understand the syndrome's effects on sleep disturbances, emphasizing the need for targeted interventions to improve sleep quality in individuals with rare genetic syndromes and their caregivers.

CITATION

Mosini A, Moysés-Oliveira M, Adami L, et al. Subjective sleep assessment in individuals with SYNGAP1-associated syndrome. J Clin Sleep Med. 2024;20(12):1879-1885.

Transcriptomic analysis of the human habenula in schizophrenia

Pathophysiology of many neuropsychiatric disorders, including schizophrenia (SCZD), is linked to habenula (Hb) function. While pharmacotherapies and deep brain stimulation targeting the Hb are emerging as promising therapeutic treatments, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in SCZD. Here we define the molecular neuroanatomy of the human Hb and identify transcriptomic changes in individuals with SCZD compared to neurotypical controls. Utilizing Hb-enriched postmortem human brain tissue, we performed single nucleus RNA-sequencing (snRNA-seq; n=7 neurotypical donors) and identified 17 molecularly defined Hb cell types across 16,437 nuclei, including 3 medial and 7 lateral Hb populations, several of which were conserved between rodents and humans. Single molecule fluorescent in situ hybridization (smFISH; n=3 neurotypical donors) validated snRNA-seq Hb cell types and mapped their spatial locations. Bulk RNA-sequencing and cell type deconvolution in Hb-enriched tissue from 35 individuals with SCZD and 33 neurotypical controls yielded 45 SCZD-associated differentially expressed genes (DEGs, FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. eQTL analysis identified 717 independent SNP-gene pairs (FDR < 0.05), where either the SNP is a SCZD risk variant (16 pairs) or the gene is a SCZD DEG (7 pairs). eQTL and SCZD risk colocalization analysis identified 16 colocalized genes. These results identify topographically organized cell types with distinct molecular signatures in the human Hb and demonstrate unique genetic changes associated with SCZD, thereby providing novel molecular insights into the role of Hb in neuropsychiatric disorders.

One Sentence Summary

Transcriptomic analysis of the human habenula and identification of molecular changes associated with schizophrenia risk and illness state.

Exploration on the potential efficacy and mechanism of methyl salicylate glycosides in the treatment of schizophrenia based on bioinformatics, molecular docking and dynamics simulation

The etiological and therapeutic complexities of schizophrenia (SCZ) persist, prompting exploration of anti-inflammatory therapy as a potential treatment approach. Methyl salicylate glycosides (MSGs), possessing a structural parent nucleus akin to aspirin, are being investigated for their therapeutic potential in schizophrenia. Utilizing bioinformation mining, network pharmacology, molecular docking and dynamics simulation, the potential value and mechanism of MSGs (including MSTG-A, MSTG-B, and Gaultherin) in the treatment of SCZ, as well as the underlying pathogenesis of the disorder, were examined. 581 differentially expressed genes related to SCZ were identified in patients and healthy individuals, with 349 up-regulated genes and 232 down-regulated genes. 29 core targets were characterized by protein-protein interaction (PPI) network, with the top 10 core targets being BDNF, VEGFA, PVALB, KCNA1, GRIN2A, ATP2B2, KCNA2, APOE, PPARGC1A and SCN1A. The pathogenesis of SCZ primarily involves cAMP signaling, neurodegenerative diseases and other pathways, as well as regulation of ion transmembrane transport. Molecular docking analysis revealed that the three candidates exhibited binding activity with certain targets with binding affinities ranging from -4.7 to -109.2 kcal/mol. MSTG-A, MSTG-B and Gaultherin show promise for use in the treatment of SCZ, potentially through their ability to modulate the expression of multiple genes involved in synaptic structure and function, ion transport, energy metabolism. Molecular dynamics simulation revealed good binding abilities between MSTG-A, MSTG-B, Gaultherin and ATP2B2. It suggests new avenues for further investigation in this area.

Identifying potential risk genes and pathways for neuropsychiatric and substance use disorders using intermediate molecular mediator information

Neuropsychiatric and substance use disorders (NPSUDs) have a complex etiology that includes environmental and polygenic risk factors with significant cross-trait genetic correlations. Genome-wide association studies (GWAS) of NPSUDs yield numerous association signals. However, for most of these regions, we do not yet have a firm understanding of either the specific risk variants or the effects of these variants. Post-GWAS methods allow researchers to use GWAS summary statistics and molecular mediators (transcript, protein, and methylation abundances) infer the effect of these mediators on risk for disorders. One group of post-GWAS approaches is commonly referred to as transcriptome/proteome/methylome-wide association studies, which are abbreviated as T/P/MWAS (or collectively as XWAS). Since these approaches use biological mediators, the multiple testing burden is reduced to the number of genes (∼20,000) instead of millions of GWAS SNPs, which leads to increased signal detection. In this work, our aim is to uncover likely risk genes for NPSUDs by performing XWAS analyses in two tissues-blood and brain. First, to identify putative causal risk genes, we performed an XWAS using the Summary-data-based Mendelian randomization, which uses GWAS summary statistics, reference xQTL data, and a reference LD panel. Second, given the large comorbidities among NPSUDs and the shared cis-xQTLs between blood and the brain, we improved XWAS signal detection for underpowered analyses by performing joint concordance analyses between XWAS results i) across the two tissues and ii) across NPSUDs. All XWAS signals i) were adjusted for heterogeneity in dependent instruments (HEIDI) (non-causality) p-values and ii) used to test for pathway enrichment. The results suggest that there were widely shared gene/protein signals within the major histocompatibility complex region on chromosome 6 (BTN3A2 and C4A) and elsewhere in the genome (FURIN, NEK4, RERE, and ZDHHC5). The identification of putative molecular genes and pathways underlying risk may offer new targets for therapeutic development. Our study revealed an enrichment of XWAS signals in vitamin D and omega-3 gene sets. So, including vitamin D and omega-3 in treatment plans may have a modest but beneficial effect on patients with bipolar disorder.

G protein-coupled receptors in neurodegenerative diseases and psychiatric disorders

Neuropsychiatric disorders are multifactorial disorders with diverse aetiological factors. Identifying treatment targets is challenging because the diseases are resulting from heterogeneous biological, genetic, and environmental factors. Nevertheless, the increasing understanding of G protein-coupled receptor (GPCR) opens a new possibility in drug discovery. Harnessing our knowledge of molecular mechanisms and structural information of GPCRs will be advantageous for developing effective drugs. This review provides an overview of the role of GPCRs in various neurodegenerative and psychiatric diseases. Besides, we highlight the emerging opportunities of novel GPCR targets and address recent progress in GPCR drug development.

Lifetime cannabis use and childhood trauma associated with CNR1 genetic variants increase the risk of psychosis: findings from the STREAM study

OBJECTIVES

Gene-environment interactions increase the risk of psychosis. The objective of this study was to investigate gene-gene and gene-environment interactions in psychosis, including single nucleotide variants (SNVs) of dopamine-2 receptor (D2R), N-methyl-d-aspartate receptor (NMDAR), and cannabinoid receptor type 1 (CB1R), lifetime cannabis use, and childhood trauma.

METHODS

Twenty-three SNVs of genes encoding D2R (DRD2: rs1799978, rs7131056, rs6275), NMDAR (GRIN1: rs4880213, rs11146020; GRIN2A: rs1420040, rs11866328; GRIN2B: rs890, rs2098469, rs7298664), and CB1R (CNR1: rs806380, rs806379, rs1049353, rs6454674, rs1535255, rs2023239, rs12720071, rs6928499, rs806374, rs7766029, rs806378, rs10485170, rs9450898) were genotyped in 143 first-episode psychosis patients (FEPp) and 286 community-based controls by Illumina HumanCoreExome-24 BeadChip. Gene-gene and gene-environment associations were assessed using nonparametric Multifactor Dimensionality Reduction software.

RESULTS

Single-locus analyses among the 23 SNVs for psychosis and gene-gene interactions were not significant (p > 0.05 for all comparisons); however, both environmental risk factors showed an association with psychosis (p < 0.001). Moreover, gene-environment interactions were significant for an SNV in CNR1 and cannabis use. The best-performing model was the combination of CNR1 rs12720071 and lifetime cannabis use (p < 0.001), suggesting an increased risk of psychosis.

CONCLUSION

Our study supports the hypothesis of gene-environment interactions for psychosis involving T-allele carriers of CNR1 SNVs, childhood trauma, and cannabis use.

Functional connectivity signatures of NMDAR dysfunction in schizophrenia-integrating findings from imaging genetics and pharmaco-fMRI

Both, pharmacological and genome-wide association studies suggest N-methyl-D-aspartate receptor (NMDAR) dysfunction and excitatory/inhibitory (E/I)-imbalance as a major pathophysiological mechanism of schizophrenia. The identification of shared fMRI brain signatures of genetically and pharmacologically induced NMDAR dysfunction may help to define biomarkers for patient stratification. NMDAR-related genetic and pharmacological effects on functional connectivity were investigated by integrating three different datasets: (A) resting state fMRI data from 146 patients with schizophrenia genotyped for the disease-associated genetic variant rs7191183 of GRIN2A (encoding the NMDAR 2 A subunit) as well as 142 healthy controls. (B) Pharmacological effects of the NMDAR antagonist ketamine and the GABA-A receptor agonist midazolam were obtained from a double-blind, crossover pharmaco-fMRI study in 28 healthy participants. (C) Regional gene expression profiles were estimated using a postmortem whole-brain microarray dataset from six healthy donors. A strong resemblance was observed between the effect of the genetic variant in schizophrenia and the ketamine versus midazolam contrast of connectivity suggestive for an associated E/I-imbalance. This similarity became more pronounced for regions with high density of NMDARs, glutamatergic neurons, and parvalbumin-positive interneurons. From a functional perspective, increased connectivity emerged between striato-pallido-thalamic regions and cortical regions of the auditory-sensory-motor network, while decreased connectivity was observed between auditory (superior temporal gyrus) and visual processing regions (lateral occipital cortex, fusiform gyrus, cuneus). Importantly, these imaging phenotypes were associated with the genetic variant, the differential effect of ketamine versus midazolam and schizophrenia (as compared to healthy controls). Moreover, the genetic variant was associated with language-related negative symptomatology which correlated with disturbed connectivity between the left posterior superior temporal gyrus and the superior lateral occipital cortex. Shared genetic and pharmacological functional connectivity profiles were suggestive of E/I-imbalance and associated with schizophrenia. The identified brain signatures may help to stratify patients with a common molecular disease pathway providing a basis for personalized psychiatry.

DRUG-seq Provides Unbiased Biological Activity Readouts for Neuroscience Drug Discovery

Unbiased transcriptomic RNA-seq data has provided deep insights into biological processes. However, its impact in drug discovery has been narrow given high costs and low throughput. Proof-of-concept studies with Digital RNA with pertUrbation of Genes (DRUG)-seq demonstrated the potential to address this gap. We extended the DRUG-seq platform by subjecting it to rigorous testing and by adding an open-source analysis pipeline. The results demonstrate high reproducibility and ability to resolve the mechanism(s) of action for a diverse set of compounds. Furthermore, we demonstrate how this data can be incorporated into a drug discovery project aiming to develop therapeutics for schizophrenia using human stem cell-derived neurons. We identified both an on-target activation signature, induced by a set of chemically distinct positive allosteric modulators of the N-methyl-d-aspartate (NMDA) receptor, and independent off-target effects. Overall, the protocol and open-source analysis pipeline are a step toward industrializing RNA-seq for high-complexity transcriptomics studies performed at a saturating scale.

Integrative omics of schizophrenia: from genetic determinants to clinical classification and risk prediction

Schizophrenia (SCZ) is a debilitating neuropsychiatric disorder with high heritability and complex inheritance. In the past decade, successful identification of numerous susceptibility loci has provided useful insights into the molecular etiology of SCZ. However, applications of these findings to clinical classification and diagnosis, risk prediction, or intervention for SCZ have been limited, and elucidating the underlying genomic and molecular mechanisms of SCZ is still challenging. More recently, multiple Omics technologies - genomics, transcriptomics, epigenomics, proteomics, metabolomics, connectomics, and gut microbiomics - have all been applied to examine different aspects of SCZ pathogenesis. Integration of multi-Omics data has thus emerged as an approach to provide a more comprehensive view of biological complexity, which is vital to enable translation into assessments and interventions of clinical benefit to individuals with SCZ. In this review, we provide a broad survey of the single-omics studies of SCZ, summarize the advantages and challenges of different Omics technologies, and then focus on studies in which multiple omics data are integrated to unravel the complex pathophysiology of SCZ. We believe that integration of multi-Omics technologies would provide a roadmap to create a more comprehensive picture of interactions involved in the complex pathogenesis of SCZ, constitute a rich resource for elucidating the potential molecular mechanisms of the illness, and eventually improve clinical assessments and interventions of SCZ to address clinical translational questions from bench to bedside.

A preliminary genetic association study of GAD1 and GABAB receptor genes in patients with treatment-resistant schizophrenia

BACKGROUND

GABAergic system dysfunction has been implicated in the etiology of schizophrenia and of cognitive impairments in particular. Patients with treatment-resistant schizophrenia (TRS) generally suffer from profound cognitive impairments in addition to severe positive symptoms, suggesting that GABA system dysfunction could be involved more closely in patients with TRS.

METHODS AND RESULTS

In the present study, exome sequencing was conducted on fourteen TRS patients, whereby four SNPs were identified on GAD1, GABBR1 and GABBR2 genes. An association study for five SNPs including these 4 SNPs and rs3749034 on GAD1 as then performed among 357 patients with TRS, 682 non-TRS patients and 508 healthy controls (HC). The results revealed no significant differences in allelic and/or genetic distributions for any of the five SNPs. However, several subanalyses in comparisons between schizophrenia and HC groups, as well as between the three groups, showed nominal-level significance for rs3749034 on GAD1 and rs10985765/rs3750344 on GABBR2. In particular, in comparisons of female subjects, rigorous analysis for rs3749034 showed a statistical difference between the schizophrenia and HC groups and between the TRS and HC groups.

CONCLUSIONS

Several positive results in subanalyses suggested that genetic vulnerability in the GABA system to schizophrenia or TRS could be affected by sex or sampling area, and overall, that rs3749034 on GAD1 and rs10985765 on GABBR2 could be related to TRS. In the present study, only a few SNPs were examined; it is possible that other important genetic variants in other regions of GABA-related genes were not captured in this preliminary study.

A Selective Review of the Excitatory-Inhibitory Imbalance in Schizophrenia: Underlying Biology, Genetics, Microcircuits, and Symptoms

Schizophrenia is a chronic disorder characterized by specific positive and negative primary symptoms, social behavior disturbances and cognitive deficits (e.g., impairment in working memory and cognitive flexibility). Mounting evidence suggests that altered excitability and inhibition at the molecular, cellular, circuit and network level might be the basis for the pathophysiology of neurodevelopmental and neuropsychiatric disorders such as schizophrenia. In the past decades, human and animal studies have identified that glutamate and gamma-aminobutyric acid (GABA) neurotransmissions are critically involved in several cognitive progresses, including learning and memory. The purpose of this review is, by analyzing emerging findings relating to the balance of excitatory and inhibitory, ranging from animal models of schizophrenia to clinical studies in patients with early onset, first-episode or chronic schizophrenia, to discuss how the excitatory-inhibitory imbalance may relate to the pathophysiology of disease phenotypes such as cognitive deficits and negative symptoms, and highlight directions for appropriate therapeutic strategies.

Alternative models for transgenerational epigenetic inheritance: Molecular psychiatry beyond mice and man

Mental illness remains the greatest chronic health burden globally with few in-roads having been made despite significant advances in genomic knowledge in recent decades. The field of psychiatry is constantly challenged to bring new approaches and tools to address and treat the needs of vulnerable individuals and subpopulations, and that has to be supported by a continuous growth in knowledge. The majority of neuropsychiatric symptoms reflect complex gene-environment interactions, with epigenetics bridging the gap between genetic susceptibility and environmental stressors that trigger disease onset and drive the advancement of symptoms. It has more recently been demonstrated in preclinical models that epigenetics underpins the transgenerational inheritance of stress-related behavioural phenotypes in both paternal and maternal lineages, providing further supporting evidence for heritability in humans. However, unbiased prospective studies of this nature are practically impossible to conduct in humans so preclinical models remain our best option for researching the molecular pathophysiologies underlying many neuropsychiatric conditions. While rodents will remain the dominant model system for preclinical studies (especially for addressing complex behavioural phenotypes), there is scope to expand current research of the molecular and epigenetic pathologies by using invertebrate models. Here, we will discuss the utility and advantages of two alternative model organisms-Caenorhabditis elegans and Drosophila melanogaster-and summarise the compelling insights of the epigenetic regulation of transgenerational inheritance that are potentially relevant to human psychiatry.

SYNGAP1 and Its Related Epileptic Syndromes

Synaptic Ras GTPase-activating protein 1 (SYNGAP1) is abundantly expressed in the postsynaptic space in brain tissue and has a crucial role in the regulation of the excitatory/inhibitory balance and in brain development. It is estimated that SYNGAP1 loss of function variants have an incidence of 1 to 4/10,000 individuals, mostly occurring de novo, even if few cases of vertical transmission of mosaic mutations have been reported. Loss-of-function mutations within this gene have been related with an epileptic encephalopathy characterized by eyelid myoclonia with absences (EMA) and myoclonic-atonic seizures (MAE) with early onset, commonly resistant to antiepileptic drugs (AED). Epilepsy is often associated with other clinical features, including truncal and/or facial hypotonia and/or ataxia with a wide-based and unsteady gate. Other clinical signs are intellectual disability, developmental delay, and behavioral and speech impairment, in a context of a normal neuroimaging study. In selected cases, dysmorphic features, skeletal abnormalities, and eye involvement are also described. The diagnosis of the disorder is usually established by multigene panel and, in unsolved cases, by exome sequencing. Management of the affected individuals involves different specialists and is mainly symptomatic. No clinical trials about the efficacy of AED in SYNGAP1 encephalopathy have been performed yet and Lamotrigine and valproate are commonly prescribed. In more than half of cases, however, epilepsy is refractory to AED.

Differential methylation of imprinting genes and MHC locus in 22q11.2 deletion syndrome-related schizophrenia spectrum disorders

OBJECTIVES

22q11.2 deletion syndrome (DS) is the strongest known genetic risk for schizophrenia. Methylome screening was conducted to elucidate possible involvement of epigenetic alterations in the emergence of schizophrenia spectrum disorders (SZ-SD) in 22q11.2DS.

METHODS

Sixteen adult men with/without SZ-SD were recruited from a 22q11.2DS cohort and underwent genome-wide DNA methylation profile analysis. Differentially methylated probes (DMPs) and regions (DMRs) were analysed using the ChAMP software.

RESULTS

The DMPs (p-value <10^-6) and DMRs (p-valueArea <0.01) were enriched in two gene sets, 'imprinting genes' and 'chr6p21', a region overlapping the MHC locus. Most of the identified imprinting genes are involved in neurodevelopment and located in clusters under imprinting control region (ICR) regulation, including PEG10, SGCE (7q21.3), GNAS, GNAS-AS1 (20q13.32) and SNHG14, SNURF-SNRPN, SNORD115 (15q11.2). The differentially methylated genes from the MHC locus included immune HLA-genes and non-immune genes, RNF39, PPP1R18 and NOTCH4, implicated in neurodevelopment and synaptic plasticity. The most significant DMR is located in MHC locus and covered the transcription regulator ZFP57 that is required for control and maintenance of gene imprinting at multiple ICRs.

CONCLUSIONS

The differential methylation in imprinting genes and in chr6p21-22 indicate the neurodevelopmental nature of 22q11.2DS-related SZ and the major role of MHC locus in the risk to develop SZ.

Twenty Years of SynGAP Research: From Synapses to Cognition

SynGAP is a potent regulator of biochemical signaling in neurons and plays critical roles in neuronal function. It was first identified in 1998, and has since been extensively characterized as a mediator of synaptic plasticity. Because of its involvement in synaptic plasticity, SynGAP has emerged as a critical protein for normal cognitive function. In recent years, mutations in the SYNGAP1 gene have been shown to cause intellectual disability in humans and have been linked to other neurodevelopmental disorders, such as autism spectrum disorders and schizophrenia. While the structure and biochemical function of SynGAP have been well characterized, a unified understanding of the various roles of SynGAP at the synapse and its contributions to neuronal function remains to be achieved. In this review, we summarize and discuss the current understanding of the multifactorial role of SynGAP in regulating neuronal function gathered over the last two decades.

A genome-wide multiphenotypic association analysis identified candidate genes and gene ontology shared by four common risky behaviors

BACKGROUND

Risky behaviors can lead to huge economic and health losses. However, limited efforts are paid to explore the genetic mechanisms of risky behaviors.

RESULT

MASH analysis identified a group of target genes for risky behaviors, such as APBB2, MAPT and DCC. For GO enrichment analysis, FUMA detected multiple risky behaviors related GO terms and brain related diseases, such as regulation of neuron differentiation (adjusted P value = 2.84×10^-5), autism spectrum disorder (adjusted P value =1.81×10^-27) and intelligence (adjusted P value =5.89×10^-15).

CONCLUSION

We reported multiple candidate genes and GO terms shared by the four risky behaviors, providing novel clues for understanding the genetic mechanism of risky behaviors.

METHODS

Multivariate Adaptive Shrinkage (MASH) analysis was first applied to the GWAS data of four specific risky behaviors (automobile speeding, drinks per week, ever-smoker, number of sexual partners) to detect the common genetic variants shared by the four risky behaviors. Utilizing genomic functional annotation data of SNPs, the SNPs detected by MASH were then mapped to target genes. Finally, gene set enrichment analysis of the identified candidate genes were conducted by the FUMA platform to obtain risky behaviors related gene ontology (GO) terms as well as diseases and traits, respectively.

Single nucleotide polymorphism rs110861313 in the intergenic region of chromosome 23 is associated with the development of leukosis in the Russian Black Pied cattle

In recent years, using a genome-wide association study (GWAS), a number of single nucleotide polymorphisms (SNPs) have been suggested to be associated with susceptibility to leukemia in cattle. However, all studies have been done with purebred Holstein cows and their hybrids. In this regard, it is important to confirm the functional role of polymorphisms previously identified in a GWAS study in Russian cattle breeds. The aim of this study was to verify the association between rs110861313 in the intergenic region of bovine chromosome 23 and leukemia in the Russian Black Pied cattle. Based on the levels of bovine leukemia virus (BLV)-specific antibodies detected in serum using serodiagnostic techniques, animals were divided into three groups: healthy animals (n = 115), asymptomatic virus carriers (n = 145) and animals with leukemia (n = 107). Genotyping of rs110861313 was carried out using polymerase chain reaction followed by analysis of restriction fragment length polymorphisms. A significant decrease in the frequency of the A/A genotype (11.2 %) was revealed in animals with persistent lymphocytosis compared to virus carriers (27.6 %) (p < 0.002). At the same time, the frequency of animals with the C/C genotype in animals with persistent lymphocytosis (41.1 %) was significantly higher than that of virus carriers (21.4 %) (p < 0.001). In this case, asymptomatic virus carriers can be considered a more suitable control than healthy animals that have not been in contact with the virus. According to bioinformatics analysis, resistance to BLV can be due to the presence of the transcription factor FOXM1 binding site in the region of rs110861313. FOXM1 is expressed in immune cells and can potentially affect the expression of the neighboring genes (LY6G5B, GPANK1, ABHD16A, LY6G6F, LY6G6E, CSNK2B, ApoM). Thus, we found that SNP rs110861313 in the intergenic region of bovine chromosome 23 is associated with the development of leukemia following BLV infection in the Russian Black Pied cattle.