Immunoglobulin genes expressed in lymphoblastoid cell lines discern and predict lithium response in bipolar disorder patients.

Bipolar disorder (BD) is a neuropsychiatric mood disorder manifested by recurrent episodes of mania and depression. More than half of BD patients are non-responsive to lithium, the first-line treatment drug, complicating BD clinical management. Given its unknown etiology, it is pertinent to understand the genetic signatures that lead to variability in lithium response. We discovered a set of differentially expressed genes (DEGs) from the lymphoblastoid cell lines (LCLs) of 10 controls and 19 BD patients belonging mainly to the immunoglobulin gene family that can be used as potential biomarkers to diagnose and treat BD. Importantly, we trained machine learning algorithms on our datasets that predicted the lithium response of BD subtypes with minimal errors, even when used on a different cohort of 24 BD patients acquired by a different laboratory. This proves the scalability of our methodology for predicting lithium response in BD and for a prompt and suitable decision on therapeutic interventions.

Enriching iPSC research diversity: Harnessing human biobank collections for improved ethnic representation.

Biobanks of human biosamples and cell lines are indispensable for biomedical research on human health and disease and for drug development projects. Many human cell line biobanks worldwide hold collections of lymphoblastoid cell lines (LCLs), representing thousands of affected and control donors from diverse ethnic/ancestry groups. In recent years, induced human pluripotent stem cells (iPSCs) and differentiated human cells derived from these iPSCs have become indispensable for applied biomedical research. Establishing iPSCs remains a laborious and costly step towards generating differentiated human cells. To address this research need, several non-profit and commercial biobanks have established iPSC collections for distribution to researchers, thereby serving as a resource for generating differentiated human cells. The most common starting materials for generation of iPSCs are a skin biopsy for harvesting fibroblasts, or a blood sample for collection of peripheral blood mononuclear cells. However untapped resources include the large established collections of biobanked human LCLs which can be reprogrammed to iPSCs using a variety of published protocols including the use of non-integrating episomal vectors. Many biobanks curate LCLs from diverse ethnic/ancestry populations, an aspect largely absent in most established iPSC biobanks which tend to primarily reflect populations from developed countries. Here, we call upon researchers across the breadth of iPSC research to tap the unique resource of existing and diverse human LCL collections for establishing biobanked iPSC panels that better represent the varied human ethnic (and hence genomic) diversity, thereby benefiting precision medicine and drug development research on a global scale.

Cannabinoid receptor gene CNR1 is downregulated in subcortical brain samples and upregulated in blood samples of individuals with schizophrenia: A participant data systematic meta-analysis.

Cannabis use leads to symptom exacerbation in schizophrenia patients, and endocannabinoid ligands have been studied as tentative schizophrenia therapeutics. Here, we aimed to characterise the connection between schizophrenia and the cannabinoid receptor 1 gene (CNR1) and explore possible mechanisms affecting its expression in schizophrenia. We performed a participant data systematic meta-analysis of CNR1 gene expression and additional endocannabinoid system genes in both brain (subcortical areas) and blood samples. We integrated eight brain sample datasets (overall 316 samples; 149 schizophrenia and 167 controls) and two blood sample datasets (overall 90 samples; 53 schizophrenia and 37 controls) while following the PRISMA meta-analysis guidelines. CNR1 was downregulated in subcortical regions and upregulated in blood samples of patients with schizophrenia. CNR2 and genes encoding endocannabinoids synthesis and degradation did not show differential expression in the brain or blood, except fatty acid amide hydrolase (FAAH), which showed a downregulation trend in blood. In addition, the brain expression levels of CNR1 and three GABA receptor genes, GABRA1, GABRA6 and GABRG2, were positively correlated (R = .57, .36, .54; p = 2.7 × 10-14 , 6.9 × 10-6 and 1.1 × 10-12 , respectively). Brain CNR1 downregulation and the positive correlation with three GABA receptor genes suggest an association with GABA neurotransmission and possible effects on negative schizophrenia symptoms. Further studies are required for clarifying the opposite CNR1 dysregulation in the brain and blood of schizophrenia patients and the potential of endocannabinoid ligands as schizophrenia therapeutics.

High purity and integrity RNA from human cell lines stored in liquid nitrogen for over 20 years.

Biobanks are a key resource for obtaining human cell lines for biomedical research, including for drug development projects. Such projects often include comparative RNA-sequencing of large panels of human cell lines from individuals affected by certain disorders and healthy controls, or from individuals with different drug response phenotypes. RNA extractions are typically done from growing cell cultures, a process that may take several weeks. However, maintaining large numbers of cell lines in parallel increases the project workload. Here, we show that extracting RNAs directly from frozen vials of human cell lines stored for over 20 years in a liquid nitrogen freezer yields RNAs with the high purity and integrity parameters that conform to those required for optimal RNA-sequencing and are closely similar to those obtained for RNAs extracted from growing human cell lines.

Machine Learning-Based Blood RNA Signature for Diagnosis of Autism Spectrum Disorder.

Early diagnosis of autism spectrum disorder (ASD) is crucial for providing appropriate treatments and parental guidance from an early age. Yet, ASD diagnosis is a lengthy process, in part due to the lack of reliable biomarkers. We recently applied RNA-sequencing of peripheral blood samples from 73 American and Israeli children with ASD and 26 neurotypically developing (NT) children to identify 10 genes with dysregulated blood expression levels in children with ASD. Machine learning (ML) analyzes data by computerized analytical model building and may be applied to building diagnostic tools based on the optimization of large datasets. Here, we present several ML-generated models, based on RNA expression datasets collected during our recently published RNA-seq study, as tentative tools for ASD diagnosis. Using the random forest classifier, two of our proposed models yield an accuracy of 82% in distinguishing children with ASD and NT children. Our proof-of-concept study requires refinement and independent validation by studies with far larger cohorts of children with ASD and NT children and should thus be perceived as starting point for building more accurate ML-based tools. Eventually, such tools may potentially provide an unbiased means to support the early diagnosis of ASD.

Nuclear αvß3 integrin expression, post translational modifications and regulation in hematological malignancies.

αvß3 integrin, a plasma membrane protein, is amply expressed on an array of tumors. We identified nuclear αvß3 pool in ovarian cancer cells and were interested to explore this phenomenon in two rare and aggressive types of leukemia, T-cell acute lymphoblastic leukemia (T-ALL) and Mast cell leukemia (MCL) using Jurkat and HMC-1 cell lines, respectively. Moreover, we collected primary cells from patients with chronic lymphocytic leukemia (CLL, n = 11), the most common chronic adult leukemia and used human lymphoblastoid cell lines (LCL) generated from normal B cells. Nuclear αvß3 integrin was assessed by Western blots, confocal microscopy, and the ImageStream technology which combines flow-cytometry with microscopy. We further examined post translational modifications (phosphorylation/glycosylation), nuclear trafficking regulation using inhibitors for MAPK (U0126) and PI3K (LY294002), as well as nuclear interactions by performing Co-immunoprecipitation (Co-IP). αvß3 integrin was identified in all cell models within the nucleus and is N-glycosylated. In primary CLL cells the ß3 integrin monomer is tyrosine Y759 phosphorylated, suggesting an active receptor conformation. MAPK and PI3K inhibition in Jurkat and CLL cells led to αvß3 enhancement in the nucleus and a reduction in the membrane. The nuclear αvß3 integrin interacts with ERK, Histone H3 and Lamin B1 in Jurkat, Histone H3 in CLL cells, but not in control LCL cells. To conclude, this observational study provides the identification of nuclear αvß3 in hematological malignancies and lays the basis for novel cancer-relevant actions, which may be independent from the membrane functions.

Discovery of autism/intellectual disability somatic mutations in Alzheimer's brains: mutated ADNP cytoskeletal impairments and repair as a case study.

With Alzheimer's disease (AD) exhibiting reduced ability of neural stem cell renewal, we hypothesized that de novo mutations controlling embryonic development, in the form of brain somatic mutations instigate the disease. A leading gene presenting heterozygous dominant de novo autism-intellectual disabilities (ID) causing mutations is activity-dependent neuroprotective protein (ADNP), with intact ADNP protecting against AD-tauopathy. We discovered a genomic autism ADNP mutation (c.2188C>T) in postmortem AD olfactory bulbs and hippocampi. RNA-Seq of olfactory bulbs also identified a novel ADNP hotspot mutation, c.2187_2188insA. Altogether, 665 mutations in 596 genes with 441 mutations in AD patients (389 genes, 38% AD-exclusive mutations) and 104 genes presenting disease-causing mutations (OMIM) were discovered. OMIM AD mutated genes converged on cytoskeletal mechanisms, autism and ID causing mutations (about 40% each). The number and average frequencies of AD-related mutations per subject were higher in AD subjects compared to controls. RNA-seq datamining (hippocampus, dorsolateral prefrontal cortex, fusiform gyrus and superior frontal gyrus-583 subjects) yielded similar results. Overlapping all tested brain areas identified unique and shared mutations, with ADNP singled out as a gene associated with autism/ID/AD and presenting several unique aging/AD mutations. The large fusiform gyrus library (117 subjects) with high sequencing coverage correlated the c.2187_2188insA ADNP mutation frequency to Braak stage (tauopathy) and showed more ADNP mutations in AD specimens. In cell cultures, the ADNP-derived snippet NAP inhibited mutated-ADNP-microtubule (MT) toxicity and enhanced Tau-MT association. We propose a paradigm-shifting concept in the perception of AD whereby accumulating mosaic somatic mutations promote brain pathology.

Introducing ADNP and SIRT1 as new partners regulating microtubules and histone methylation.

Activity-dependent neuroprotective protein (ADNP) is essential for brain formation and function. As such, de novo mutations in ADNP lead to the autistic ADNP syndrome and somatic ADNP mutations may drive Alzheimer's disease (AD) tauopathy. Sirtuin 1 (SIRT1) is positively associated with aging, the major risk for AD. Here, we revealed two key interaction sites for ADNP and SIRT1. One, at the microtubule end-binding protein (EB1 and EB3) Tau level, with EB1/EB3 serving as amplifiers for microtubule dynamics, synapse formation, axonal transport, and protection against tauopathy. Two, on the DNA/chromatin site, with yin yang 1, histone deacetylase 2, and ADNP, sharing a DNA binding motif and regulating SIRT1, ADNP, and EB1 (MAPRE1). This interaction was linked to sex- and age-dependent altered histone modification, associated with ADNP/SIRT1/WD repeat-containing protein 5, which mediates the assembly of histone modification complexes. Single-cell RNA and protein expression analyses as well as gene expression correlations placed SIRT1-ADNP and either MAPRE1 (EB1), MAPRE3 (EB3), or both in the same mouse and human cell; however, while MAPRE1 seemed to be similarly regulated to ADNP and SIRT1, MAPRE3 seemed to deviate. Finally, we demonstrated an extremely tight correlation for the gene transcripts described above, including related gene products. This correlation was specifically abolished in affected postmortem AD and Parkinson's disease brain select areas compared to matched controls, while being maintained in blood samples. Thus, we identified an ADNP-SIRT1 complex that may serve as a new target for the understanding of brain degeneration.

d-Galactose treatment increases ACE2, TMPRSS2, and FURIN and reduces SERPINA1 mRNA expression in A549 human lung epithelial cells.

Several comorbidities including diabetes, immune deficiency, and chronic respiratory disorders increase the risk of severe Covid-19 and fatalities among SARS-CoV-2 infected individuals. Severe Covid-19 risk among diabetes patients may reflect reduced immune response to viral infections. SARS-CoV-2 initially infects respiratory tract epithelial cells by binding to the host cell membrane ACE2, followed by proteolytic priming for cell entry by the host cell membrane serine protease TMPRSS2. Additionally, the protease FURIN facilitates cell exit of mature SARS-CoV-2 virions. Alpha-1 antitrypsin (AAT), the major plasma serine protease inhibitor, encoded by SERPINA1, is known to promote immune response to viral infections. AAT inhibits neutrophil elastase, a key inflammatory serine protease implicated in alveolar cell damage during respiratory infections, and AAT deficiency is associated with susceptibility to lung infections. AAT is implicated in Covid-19 as it inhibits TMPRSS2, a protease essential for SARS-CoV-2 cell entry. Here we show that treatment of A549 human lung epithelial cells for 7 days with 25 mM d-galactose, an inducer of diabetic-like and oxidative stress cellular phenotypes, leads to increased mRNA levels of ACE2, TMPRSS2, and FURIN, along with reduced SERPINA1 mRNA. Together, the dysregulated transcription of these genes following d-galactose treatment suggests that chronic diabetic-like conditions may facilitate SARS-CoV-2 infection of lung epithelial cells. Our findings may in part explain the higher severe Covid-19 risk in diabetes, and highlight the need to develop special treatment protocols for diabetic patients.

Estradiol reduces ACE2 and TMPRSS2 mRNA levels in A549 human lung epithelial cells.

Epidemiologic studies suggest slightly higher risk of severe Covid-19 symptoms and fatalities following SARS-CoV-2 infection in men compared with women from similar age groups. This bias was suggested to reflect differences in the male and female immune system regulation, driven by different sex hormone levels in men and women, in particular, higher plasma estradiol in women. SARS-CoV-2 infects respiratory tract epithelial cells by binding to their cell membrane ACE2, followed by priming for cell entry by the host cell membrane serine protease TMPRSS2. The cell protease FURIN facilitates cell exit of mature SARS-CoV-2 virions. Our study examined the effects of in vitro treatment of A549 human lung epithelial cells with 17-ß-estradiol on mRNA expression of genes coding for these proteins. Treatment of A549 human lung epithelial cells with 17-ß-estradiol reduced the cellular mRNA levels of ACE2 and TMPRSS2 mRNA, while not affecting FURIN expression. Our findings suggest that 17-ß-estradiol may reduce SARS-CoV-2 infection of lung epithelial cells, which may in part explain the reduced incidence of severe Covid-19 and fatalities among women compared with men of similar age. Studies into the molecular pathways by which 17-ß-estradiol reduces ACE2 and TMPRSS2 mRNA expression in lung epithelial cells are needed for assessing its potential protective value against severe Covid-19.

Higher ATM expression in lymphoblastoid cell lines from centenarian compared with younger women.

With increased life expectancies in developed countries, cancer rates are becoming more common among the elderly. Cancer is typically driven by a combination of germline and somatic mutations accumulating during an individual's lifetime. Yet, many centenarians reach exceptionally old age without experiencing cancer. It was suggested that centenarians have more robust DNA repair and mitochondrial function, allowing improved maintenance of DNA stability. In this study, we applied real-time quantitative PCR to examine the expression of ATM in lymphoblastoid cell lines (LCLs) from 15 healthy female centenarians and 24 younger female donors aged 21-88 years. We observed higher ATM mRNA expression of in LCLs from female centenarians compared with both women aged 21-48 years (FD = 2.0, p = .0016) and women aged 56-88 years (FD = 1.8, p = .0094. Positive correlation was found between ATM mRNA expression and donors age (p = .0028). Levels of hsa-miR-181a-5p, which targets ATM, were lower in LCLs from centenarians compared with younger women. Our findings suggest a role for ATM in protection from age-related diseases, possibly reflecting more effective DNA repair, thereby reducing somatic mutation accumulation during aging. Further studies are required for analyzing additional DNA repair pathways in biosamples from centenarians and younger age men and women.

Blood RNA Sequencing Indicates Upregulated BATF2 and LY6E and Downregulated ISG15 and MT2A Expression in Children with Autism Spectrum Disorder.

Mutations in over 100 genes are implicated in autism spectrum disorder (ASD). DNA SNPs, CNVs, and epigenomic modifications also contribute to ASD. Transcriptomics analysis of blood samples may offer clues for pathways dysregulated in ASD. To expand and validate published findings of RNA-sequencing (RNA-seq) studies, we performed RNA-seq of whole blood samples from an Israeli discovery cohort of eight children with ASD compared with nine age- and sex-matched neurotypical children. This revealed 10 genes with differential expression. Using quantitative real-time PCR, we compared RNAs from whole blood samples of 73 Israeli and American children with ASD and 26 matched neurotypical children for the 10 dysregulated genes detected by RNA-seq. This revealed higher expression levels of the pro-inflammatory transcripts BATF2 and LY6E and lower expression levels of the anti-inflammatory transcripts ISG15 and MT2A in the ASD compared to neurotypical children. BATF2 was recently reported as upregulated in blood samples of Japanese adults with ASD. Our findings support an involvement of these genes in ASD phenotypes, independent of age and ethnicity. Upregulation of BATF2 and downregulation of ISG15 and MT2A were reported to reduce cancer risk. Implications of the dysregulated genes for pro-inflammatory phenotypes, immunity, and cancer risk in ASD are discussed.

Ethnic differences in alpha-1 antitrypsin deficiency allele frequencies may partially explain national differences in COVID-19 fatality rates.

Infection rates, severity, and fatalities due to COVID-19, the pandemic mediated by SARS-CoV-2, vary greatly between countries. With few exceptions, these are lower in East and Southeast Asian and Sub-Saharan African countries compared with other regions. Epidemiological differences may reflect differences in border closures, lockdowns, and social distancing measures taken by each county, and by cultural differences, such as common use of face masks in East and Southeast Asian countries. The plasma serine protease inhibitor alpha-1 antitrypsin was suggested to protect from COVID-19 by inhibiting TMPRSS2, a cell surface serine protease essential for the SARS-CoV-2 cell entry. Here, we present evidence that population differences in alpha-1 antitrypsin deficiency allele frequencies may partially explain national differences in the COVID-19 epidemiology. Our study compared reported national estimates for the major alpha-1 antitrypsin deficiency alleles PiZ and PiS (SERPINA1 rs28929474 and rs17580, respectively) with the Johns Hopkins University Coronavirus Resource Center dataset. We found a significant positive correlation (R = .54, P = 1.98e-6) between the combined frequencies of the alpha-1 antitrypsin PiZ and PiS deficiency alleles in 67 countries and their reported COVID-19 mortality rates. Our observations suggest that alpha-1 antitrypsin deficiency alleles may contribute to national differences in COVID-19 infection, severity, and mortality rates. Population-wide screening for carriers of alpha-1 antitrypsin deficiency alleles should be considered for prioritizing individuals for stricter social distancing measures and for receiving a SARS-CoV-2 vaccine once it becomes available.

Identification of thioredoxin-interacting protein (TXNIP) as a downstream target for IGF1 action.

Laron syndrome (LS), or primary growth hormone (GH) insensitivity, is the best-characterized entity among the congenital insulin-like growth factor 1 (IGF1) deficiencies. Life-long exposure to minute endogenous IGF1 levels is linked to low stature as well as a number of endocrine and metabolic abnormalities. While elevated IGF1 is correlated with increased cancer incidence, epidemiological studies revealed that patients with LS do not develop tumors. The mechanisms associated with cancer protection in LS are yet to be discovered. Recent genomic analyses identified a series of metabolic genes that are overrepresented in patients with LS. Given the augmented expression of these genes in a low IGF1 milieu, we hypothesized that they may constitute targets for IGF1 action. Thioredoxin-interacting protein (TXNIP) plays a critical role in cellular redox control by thioredoxin. TXNIP serves as a glucose and oxidative stress sensor, being commonly silenced by genetic or epigenetic events in cancer cells. Consistent with its enhanced expression in LS, we provide evidence that TXNIP gene expression is negatively regulated by IGF1. These results were corroborated in animal studies. In addition, we show that oxidative and glucose stresses led to marked increases in TXNIP expression. Supplementation of IGF1 attenuated TXNIP levels, suggesting that IGF1 exerts its antiapoptotic effect via inhibition of TXNIP Augmented TXNIP expression in LS may account for cancer protection in this condition. Finally, TXNIP levels could be potentially useful in the clinic as a predictive or diagnostic biomarker for IGF1R-targeted therapies.

Angiotensin receptor blockers as tentative SARS-CoV-2 therapeutics.

At the time of writing this commentary (February 2020), the coronavirus COVID-19 epidemic has already resulted in more fatalities compared with the SARS and MERS coronavirus epidemics combined. Therapeutics that may assist to contain its rapid spread and reduce its high mortality rates are urgently needed. Developing vaccines against the SARS-CoV-2 virus may take many months. Moreover, vaccines based on viral-encoded peptides may not be effective against future coronavirus epidemics, as virus mutations could make them futile. Indeed, new Influenza virus strains emerge every year, requiring new immunizations. A tentative suggestion based on existing therapeutics, which would likely be resistant to new coronavirus mutations, is to use available angiotensin receptor 1 (AT1R) blockers, such as losartan, as therapeutics for reducing the aggressiveness and mortality from SARS-CoV-2 virus infections. This idea is based on observations that the angiotensin-converting enzyme 2 (ACE2) very likely serves as the binding site for SARS-CoV-2, the strain implicated in the current COVID-19 epidemic, similarly to strain SARS-CoV implicated in the 2002-2003 SARS epidemic. This commentary elaborates on the idea of considering AT1R blockers as tentative treatment for SARS-CoV-2 infections, and proposes a research direction based on datamining of clinical patient records for assessing its feasibility.

Repurposing current therapeutics for treating COVID-19: A vital role of prescription records data mining.

Since its outbreak in late 2019, the SARS-Cov-2 pandemic already infected over 3.7 million people and claimed more than 250,000 lives globally. At least 1 year may take for an approved vaccine to be in place, and meanwhile millions more could be infected, some with fatal outcome. Over thousand clinical trials with COVID-19 patients are already listed in ClinicalTrials.com, some of them for assessing the utility of therapeutics approved for other conditions. However, clinical trials take many months, and are typically done with small cohorts. A much faster and by far more efficient method for rapidly identifying approved therapeutics that can be repurposed for treating COVID-19 patients is data mining their past and current electronic health and prescription records for identifying drugs that may protect infected individuals from severe COVID-19 symptoms. Examples are discussed for applying health and prescription records for assessing the potential repurposing (repositioning) of angiotensin receptor blockers, estradiol, or antiandrogens for reducing COVID-19 morbidity and fatalities. Data mining of prescription records of COVID-19 patients will not cancel the need for conducting controlled clinical trials, but could substantially assist in trial design, drug choice, inclusion and exclusion criteria, and prioritization. This approach requires a strong commitment of health provides for open collaboration with the biomedical research community, as health provides are typically the sole owners of retrospective drug prescription records.

Smoking and COVID-19: Similar bronchial ACE2 and TMPRSS2 expression and higher TMPRSS4 expression in current versus never smokers.

Uncertainties remain concerning the pathophysiology, epidemiology, and potential therapeutics for COVID-19. Among unsettled controversies is whether tobacco smoking increases or protects from severe COVID-19. Several epidemiological studies reported reduced COVID-19 hospitalizations among smokers, while other studies reported the opposite trend. Some authors assumed that smokers have elevated airway expression of ACE2, the cell recognition site of the SARS-Cov-2 spike protein, but this suggestion remains unverified. We therefore performed data mining of two independent NCBI GEO genome-wide RNA expression files (GSE7894 and GSE994) and report that in both data sets, current smokers and never smokers have, on average, closely similar bronchial epithelial cell mRNA levels of ACE2, as well as TMPRSS2, coding for a serine protease priming SARS-Cov-2 for cell entry, and ADAM17, coding for a protease implicated in ACE2 membrane shedding. In contrast, the expression levels of TMPRSS4, coding for a protease that primes SARS-CoV-2 for cell entry similarly to TMPRSS2, were elevated in bronchial epithelial cells from current smokers compared with never smokers, suggesting that higher bronchial TMPRSS4 levels in smokers might put them at higher SARS-Cov-2 infection risk. The effects of smoking on COVID-19 severity need clarification with larger studies. Additionally, the postulated protective effects of nicotine and nitric oxide, which may presumably reduce the risk of a "cytokine storm" in infected individuals, deserve assessment by controlled clinical trials.

SCN11A mRNA levels in female bipolar disorder PBMCs as tentative biomarker for distinct patient sub-phenotypes.

Bipolar disorder (BD) is a complex neuropsychiatric disorder characterized by recurrent mania and depression episodes and requiring lifelong treatment with mood stabilizing drugs. Several lines of evidence, including with BD patient iPSC-derived neurons, suggest that neuronal hyperexcitability may underlie the key clinical symptoms of BD. Indeed, higher mRNA levels of SCN11A, coding for the voltage-gated sodium channel NaV 1.9 implicated in nociception, were detected in iPSC-derived neurons from BD patients, and were normalized by in vitro lithium. Here we studied SCN11A expression in peripheral blood mononuclear cells (PBMCs) from well-phenotyped female BD patients and controls and evaluated their association with several clinical sub-phenotypes. We observed higher mRNA levels of SCN11A in PBMCs from female BD patients with no records of alcohol dependence (p = .0050), no records of psychosis (p = .0097), or no records of suicide attempts (p = .0409). A trend was observed for higher SCN11A expression (FD = 1.91; p = .052) in BD PBMCs compared with controls. Datamining of published postmortem gene expression datasets indicated higher SCN11A expression in dorsolateral prefrontal cortex and orbitofrontal cortex tissues from BD patients compared with controls. Higher phenotype-associated expression levels in PBMC from BD patients were also observed for ID2 (alcohol dependence, suicide attempts) and HDGFRP3 (seasonal BD pattern). Our findings suggest that higher PBMC SCN11A expression levels may be associated with certain behavioral BD sub-phenotypes, including lack of alcohol dependence and psychosis, among BD patients. The NaV 1.9 voltage-gated sodium channel thus deserves consideration as a tentative phenotype modifier in BD.