Temporal associations of plasma levels of the secreted phospholipase A2 family and mortality in severe COVID-19.

Previous research suggests that group IIA-secreted phospholipase A2 (sPLA2-IIA) plays a role in and predicts lethal COVID-19 disease. The current study reanalyzed a longitudinal proteomic data set to determine the temporal relationship between levels of several members of a family of sPLA2 isoforms and the severity of COVID-19 in 214 ICU patients. The levels of six secreted PLA2 isoforms, sPLA2-IIA, sPLA2-V, sPLA2-X, sPLA2-IB, sPLA2-IIC, and sPLA2-XVI, increased over the first 7 ICU days in those who succumbed to the disease but attenuated over the same time period in survivors. In contrast, a reversed pattern in sPLA2-IID and sPLA2-XIIB levels over 7 days suggests a protective role of these two isoforms. Furthermore, decision tree models demonstrated that sPLA2-IIA outperformed top-ranked cytokines and chemokines as a predictor of patient outcome. Taken together, proteomic analysis revealed temporal sPLA2 patterns that reflect the critical roles of sPLA2 isoforms in severe COVID-19 disease.

Candesartan restores blood-brain barrier dysfunction, mitigates aberrant gene expression, and extends lifespan in a knockin mouse model of epileptogenesis.

Blockade of Angiotensin type 1 receptor (AT1R) has potential therapeutic utility in the treatment of numerous detrimental consequences of epileptogenesis, including oxidative stress, neuroinflammation, and blood-brain barrier (BBB) dysfunction. We have recently shown that many of these pathological processes play a critical role in seizure onset and propagation in the Scn8a-N1768D mouse model. Here we investigate the efficacy and potential mechanism(s) of action of candesartan (CND), an FDA-approved angiotensin receptor blocker (ARB) indicated for hypertension, in improving outcomes in this model of pediatric epilepsy. We compared length of lifespan, seizure frequency, and BBB permeability in juvenile (D/D) and adult (D/+) mice treated with CND at times after seizure onset. We performed RNAseq on hippocampal tissue to quantify differences in genome-wide patterns of transcript abundance and inferred beneficial and detrimental effects of canonical pathways identified by enrichment methods in untreated and treated mice. Our results demonstrate that treatment with CND gives rise to increased survival, longer periods of seizure freedom, and diminished BBB permeability. CND treatment also partially reversed or 'normalized' disease-induced genome-wide gene expression profiles associated with inhibition of NF-κB, TNFα, IL-6, and TGF-ß signaling in juvenile and adult mice. Pathway analyses reveal that efficacy of CND is due to its known dual mechanism of action as both an AT1R antagonist and a PPARγ agonist. The robust effectiveness of CND across ages, sexes and mouse strains is a positive indication for its translation to humans and its suitability of use for clinical trials in children with SCN8A epilepsy.

Sex differences in physiological response to increased neuronal excitability in a knockin mouse model of pediatric epilepsy.

BACKGROUND: Epilepsy is a common neurological disease; however, few if any of the currently marketed antiseizure medications prevent or cure epilepsy. Discovery of pathological processes in the early stages of epileptogenesis has been challenging given the common use of preclinical models that induce seizures in physiologically normal animals. Moreover, despite known sex dimorphism in neurological diseases, females are rarely included in preclinical epilepsy models. METHODS: We characterized sex differences in mice carrying a pathogenic knockin variant (p.N1768D) in the Scn8a gene that causes spontaneous tonic-clonic seizures (TCs) at ∼3 months of age and found that heterozygous females are more resilient than males in mortality and morbidity. To investigate the cellular mechanisms that underlie female resilience, we utilized blood-brain barrier (BBB) and hippocampal transcriptomic analyses in heterozygous mice before seizure onset (pre-TC) and in mice that experienced ∼20 TCs (post-TC). RESULTS: In the pre-TC latent phase, both sexes exhibited leaky BBB; however, patterns of gene expression were sexually dimorphic. Females exhibited enhanced oxidative phosphorylation and protein biogenesis, while males activated gliosis and CREB signaling. After seizure onset (chronic phase), females exhibited a metabolic switch to lipid metabolism, while males exhibited increased gliosis and BBB dysfunction and a strong activation of neuroinflammatory pathways. CONCLUSION: The results underscore the central role of oxidative stress and BBB permeability in the early stages of epileptogenesis, as well as sex dimorphism in response to increasing neuronal hyperexcitability. Our results also highlight the need to include both sexes in preclinical studies to effectively translate results of drug efficacy studies.

Clinical severity is correlated with age at seizure onset and biophysical properties of recurrent gain of function variants associated with SCN8A-related epilepsy.

OBJECTIVE: Genetic variants in the SCN8A gene underlie a wide spectrum of neurodevelopmental phenotypes including several distinct seizure types and a host of comorbidities. One of the major challenges facing clinicians and researchers alike is to identify genotype-phenotype (G-P) correlations that may improve prognosis, guide treatment decisions, and lead to precision medicine approaches. METHODS: We investigated G-P correlations among 270 participants harboring gain-of-function (GOF) variants enrolled in the International SCN8A Registry, a patient-driven online database. We performed correlation analyses stratifying the cohort by clinical phenotypes to identify diagnostic features that differ among patients with varying levels of clinical severity, and that differ among patients with distinct GOF variants. RESULTS: Our analyses confirm positive correlations between age at seizure onset and developmental skills acquisition (developmental quotient), rate of seizure freedom, and percentage of cohort with developmental delays, and identify negative correlations with number of current and weaned antiseizure medications. This set of features is more detrimentally affected in individuals with a priori expectations of more severe clinical phenotypes. Our analyses also reveal a significant correlation between a severity index combining clinical features of individuals with a particular highly recurrent variant and an independent electrophysiological score assigned to each variant based on in vitro testing. SIGNIFICANCE: This is one of the first studies to identify statistically significant G-P correlations for individual SCN8A variants with GOF properties. The results suggest that individual GOF variants (1) are predictive of clinical severity for individuals carrying those variants and (2) may underlie distinct clinical phenotypes of SCN8A disease, thus helping to explain the wide SCN8A-related epilepsy disease spectrum. These results also suggest that certain features present at initial diagnosis are predictive of clinical severity, and with more informed treatment plans, may serve to improve prognosis for patients with SCN8A GOF variants.

Greater female than male resilience to mortality and morbidity in the Scn8a mouse model of pediatric epilepsy.

AIMS: Females and males of all ages are affected by epilepsy; however, unlike many clinical studies, most preclinical research has focused on males. Genetic variants in the voltage-gated sodium channel gene, SCN8A, are associated with a broad spectrum of neurological and epileptic syndromes. Here we investigate sex differences in the natural history of the Scn8a-N1768D knockin mouse model of pediatric epilepsy. METHODS: We utilize 24/7 video to monitor juveniles and adults of both sexes to investigate variability in seizure activity (e.g. onset and frequency), mortality and morbidity, response to cannabinoids, and mode of death. We also monitor sleep architecture using a noninvasive piezoelectric method in order to identify factors that influence seizure severity and outcome. RESULTS: Both sexes had nearly 100% penetrance in seizure onset and early mortality. However, adult heterozygous (D/+) females were more resilient as exhibited by the ability to tolerate more seizures over a longer lifespan. Homozygous (D/D) juveniles did not exhibit a sex difference in overall survival. Female estrus cycle was disrupted before seizure onset, while sleep was disrupted in both sexes in association with seizure onset. Females typically died while in convulsive status epilepticus; however, a high proportion of males died while not experiencing behavioral seizures. Only juvenile and adult males benefited from cannabinoid administration. CONCLUSIONS: These results support the hypothesis that factors associated with sexual differentiation play a role in the neurobiology of epilepsy and point to the importance of including both sexes in the design of studies to identify new epilepsy therapies.

A novel nonlinear dimension reduction approach to infer population structure for low-coverage sequencing data.

BACKGROUND: Low-depth sequencing allows researchers to increase sample size at the expense of lower accuracy. To incorporate uncertainties while maintaining statistical power, we introduce MCPCA_PopGen to analyze population structure of low-depth sequencing data. RESULTS: The method optimizes the choice of nonlinear transformations of dosages to maximize the Ky Fan norm of the covariance matrix. The transformation incorporates the uncertainty in calling between heterozygotes and the common homozygotes for loci having a rare allele and is more linear when both variants are common. CONCLUSIONS: We apply MCPCA_PopGen to samples from two indigenous Siberian populations and reveal hidden population structure accurately using only a single chromosome. The MCPCA_PopGen package is available on https://github.com/yiwenstat/MCPCA_PopGen .

Genomic Evidence of Local Adaptation to Climate and Diet in Indigenous Siberians.

The indigenous inhabitants of Siberia live in some of the harshest environments on earth, experiencing extended periods of severe cold temperatures, dramatic variation in photoperiod, and limited and highly variable food resources. While the successful long-term settlement of this area by humans required multiple behavioral and cultural innovations, the nature of the underlying genetic changes has generally remained elusive. In this study, we used a three-part approach to identify putative targets of positive natural selection in Siberians. We first performed selection scans on whole exome and genome-wide single nucleotide polymorphism array data from multiple Siberian populations. We then annotated candidates in the tails of the empirical distributions, focusing on candidates with evidence linking them to biological processes and phenotypes previously identified as relevant to adaptation in circumpolar groups. The top candidates were then genotyped in additional populations to determine their spatial allele frequency distributions and associations with climate variables. Our analysis reveals missense mutations in three genes involved in lipid metabolism (PLA2G2A, PLIN1, and ANGPTL8) that exhibit genomic and spatial patterns consistent with selection for cold climate and/or diet. These variants are unified by their connection to brown adipose tissue and may help to explain previously observed physiological differences in Siberians such as low serum lipid levels and increased basal metabolic rate. These results support the hypothesis that indigenous Siberians have genetically adapted to their local environment by selection on multiple genes.

The Role of Phylogenetically Conserved Elements in Shaping Patterns of Human Genomic Diversity.

Evolutionary genetic studies have shown a positive correlation between levels of nucleotide diversity and either rates of recombination or genetic distance to genes. Both positive-directional and purifying selection have been offered as the source of these correlations via genetic hitchhiking and background selection, respectively. Phylogenetically conserved elements (CEs) are short (∼100 bp), widely distributed (comprising ∼5% of genome), sequences that are often found far from genes. While the function of many CEs is unknown, CEs also are associated with reduced diversity at linked sites. Using high coverage (>80×) whole genome data from two human populations, the Yoruba and the CEU, we perform fine scale evaluations of diversity, rates of recombination, and linkage to genes. We find that the local rate of recombination has a stronger effect on levels of diversity than linkage to genes, and that these effects of recombination persist even in regions far from genes. Our whole genome modeling demonstrates that, rather than recombination or GC-biased gene conversion, selection on sites within or linked to CEs better explains the observed genomic diversity patterns. A major implication is that very few sites in the human genome are predicted to be free of the effects of selection. These sites, which we refer to as the human "neutralome," comprise only 1.2% of the autosomes and 5.1% of the X chromosome. Demographic analysis of the neutralome reveals larger population sizes and lower rates of growth for ancestral human populations than inferred by previous analyses.

Influence of age at seizure onset on the acquisition of neurodevelopmental skills in an SCN8A cohort.

OBJECTIVE: To characterize a cohort of patients with SCN8A-related epilepsy and to perform analyses to identify correlations involving the acquisition of neurodevelopmental skills. METHODS: We analyzed patient data (n = 91) submitted to an online registry tailored to characteristics of children with SCN8A variants. Participants provided information on the history of their child's seizures, medications, comorbidities, and developmental skills based on the Denver II items. Spearman rank tests were utilized to test for correlations among a variety of aspects of seizures, medications, and neurodevelopmental progression. RESULTS: The 91 participants carried 71 missense variants (41 newly reported) and three truncating variants. Ages at seizure onset ranged from birth to >12 months of age (mean ± SD = 5 months 21 days ± 7 months 14 days). Multiple seizure types with multimodal onset times and developmental delay were observed as general features of this cohort. We found a positive correlation between a developmental score based upon percentage of acquired skills and the age at seizure onset, current seizure freedom, and initial febrile seizures. Analyses of cohort subgroups revealed clear distinctions between patients who had a single reported variant in SCN8A and those with an additional variant reported in a gene other than SCN8A, as well as between patients with different patterns of regression before and at seizure onset. SIGNIFICANCE: This is the first study of an SCN8A patient cohort of this size and for which correlations between age at seizure onset and neurodevelopment were investigated. Our correlation studies suggest that variants of uncertain significance should be considered in assessing children with SCN8A-related disorders. This study substantially improves the characterization of this patient population and our understanding of the neurodevelopmental effects associated with seizures for SCN8A patients, and provides a clinical context at initial presentation that may be prognostic for developmental outcome.

Exome Sequencing Provides Evidence of Polygenic Adaptation to a Fat-Rich Animal Diet in Indigenous Siberian Populations.

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.

A de novo missense mutation of GABRB2 causes early myoclonic encephalopathy.

BACKGROUND: Early myoclonic encephalopathy (EME), a disease with a devastating prognosis, is characterised by neonatal onset of seizures and massive myoclonus accompanied by a continuous suppression-burst EEG pattern. Three genes are associated with EMEs that have metabolic features. Here, we report a pathogenic mutation of an ion channel as a cause of EME for the first time. METHODS: Sequencing was performed for 214 patients with epileptic seizures using a gene panel with 109 genes that are known or suspected to cause epileptic seizures. Functional assessments were demonstrated by using electrophysiological experiments and immunostaining for mutant γ-aminobutyric acid-A (GABAA) receptor subunits in HEK293T cells. RESULTS: We discovered a de novo heterozygous missense mutation (c.859A>C [p.Thr287Pro]) in the GABRB2-encoded ß2 subunit of the GABAA receptor in an infant with EME. No GABRB2 mutations were found in three other EME cases or in 166 patients with infantile spasms. GABAA receptors bearing the mutant ß2 subunit were poorly trafficked to the cell membrane and prevented γ2 subunits from trafficking to the cell surface. The peak amplitudes of currents from GABAA receptors containing only mutant ß2 subunits were smaller than that of those from receptors containing only wild-type ß2 subunits. The decrease in peak current amplitude (96.4% reduction) associated with the mutant GABAA receptor was greater than expected, based on the degree to which cell surface expression was reduced (66% reduction). CONCLUSION: This mutation has complex functional effects on GABAA receptors, including reduction of cell surface expression and attenuation of channel function, which would significantly perturb GABAergic inhibition in the brain.

Clinical implications of SCN1A missense and truncation variants in a large Japanese cohort with Dravet syndrome.

OBJECTIVE: Two major classes of SCN1A variants are associated with Dravet syndrome (DS): those that result in haploinsufficiency (truncating) and those that result in an amino acid substitution (missense). The aim of this retrospective study was to describe the first large cohort of Japanese patients with SCN1A mutation-positive DS (n = 285), and investigate the relationship between variant (type and position) and clinical expression and response to treatment. METHODS: We sequenced all exons and intron-exon boundaries of SCN1A in our cohort, investigated differences in the distribution of truncating and missense variants, tested for associations between variant type and phenotype, and compared these patterns with those of cohorts with milder epilepsy and healthy individuals. RESULTS: Unlike truncation variants, missense variants are found at higher density in the S4 voltage sensor and pore loops and at lower density in the domain I-II and II-III linkers and the first three segments of domain II. Relative to healthy individuals, there is an increased frequency of truncating (but not missense) variants in the noncoding C-terminus. The rate of cognitive decline is more rapid for patients with truncation variants regardless of age at seizure onset, whereas age at onset is a predictor of the rate of cognitive decline for patients with missense variants. SIGNIFICANCE: We found significant differences in the distribution of truncating and missense variants across the SCN1A sequence among healthy individuals, patients with DS, and those with milder forms of SCN1A-variant positive epilepsy. Testing for associations with phenotype revealed that variant type can be predictive of rate of cognitive decline. Analysis of descriptive medication data suggests that in addition to conventional drug therapy in DS, bromide, clonazepam and topiramate may reduce seizure frequency.

Evidence for increased levels of positive and negative selection on the X chromosome versus autosomes in humans.

Partially recessive variants under positive selection are expected to go to fixation more quickly on the X chromosome as a result of hemizygosity, an effect known as faster-X. Conversely, purifying selection is expected to reduce substitution rates more effectively on the X chromosome. Previous work in humans contrasted divergence on the autosomes and X chromosome, with results tending to support the faster-X effect. However, no study has yet incorporated both divergence and polymorphism to quantify the effects of both purifying and positive selection, which are opposing forces with respect to divergence. In this study, we develop a framework that integrates previously developed theory addressing differential rates of X and autosomal evolution with methods that jointly estimate the level of purifying and positive selection via modeling of the distribution of fitness effects (DFE). We then utilize this framework to estimate the proportion of nonsynonymous substitutions fixed by positive selection (α) using exome sequence data from a West African population. We find that varying the female to male breeding ratio (ß) has minimal impact on the DFE for the X chromosome, especially when compared with the effect of varying the dominance coefficient of deleterious alleles (h). Estimates of α range from 46% to 51% and from 4% to 24% for the X chromosome and autosomes, respectively. While dependent on h, the magnitude of the difference between α values estimated for these two systems is highly statistically significant over a range of biologically realistic parameter values, suggesting faster-X has been operating in humans.

A haplotype at STAT2 Introgressed from neanderthals and serves as a candidate of positive selection in Papua New Guinea.

Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (~80 thousand years ago) and is found throughout Eurasia at an average frequency of ~5%. Interestingly, N is found in Melanesian populations at ~10-fold higher frequency (~54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2-all of which are part of the same 250 kb introgressive haplotype-as good candidates.

Does Management of Diabetic Ketoacidosis with Subcutaneous Rapid-acting Insulin Reduce the Need for Intensive Care Unit Admission?

BACKGROUND: In the last 20 years, rapid-acting insulin analogs have emerged on the market, including aspart and lispro, which may be efficacious in the management of diabetic ketoacidosis (DKA) when administered by non-intravenous (i.v.) routes. CLINICAL QUESTION: In patients with mild-to-moderate DKA without another reason for intensive care unit (ICU) admission, is the administration of a subcutaneous (s.c.) rapid-acting insulin analog a safe and effective alternative to a continuous infusion of i.v. regular insulin, and would such a strategy eliminate the need for ICU admission? EVIDENCE REVIEW: Five randomized controlled trials were identified and critically appraised. RESULTS: The outcomes suggest that there is no difference in the duration of therapy required to resolve DKA with either strategy. CONCLUSION: Current evidence supports DKA management with s.c. rapid-acting insulin analogs in a non-ICU setting in carefully selected patients.

Neandertal origin of genetic variation at the cluster of OAS immunity genes.

Analyses of ancient DNA from extinct humans reveal signals of at least two independent hybridization events in the history of non-African populations. To date, there are very few examples of specific genetic variants that have been rigorously identified as introgressive. Here, we survey DNA sequence variation in the OAS gene cluster on chromosome 12 and provide strong evidence that a haplotype extending for ~185 kb introgressed from Neandertals. This haplotype is nearly restricted to Eurasians and is estimated to have diverged from the Neandertal sequence ~125 kya. Despite the potential for novel functional variation, the observed frequency of this haplotype is consistent with neutral introgression. This is the second locus in the human genome, after STAT2, carrying distinct haplotypes that appear to have introgressed separately from both Neandertals and Denisova.

A stochastic model of translation with -1 programmed ribosomal frameshifting.

Many viruses produce multiple proteins from a single mRNA sequence by encoding overlapping genes. One mechanism to decode both genes, which reside in alternate reading frames, is -1 programmed ribosomal frameshifting. Although recognized for over 25 years, the molecular and physical mechanism of -1 frameshifting remains poorly understood. We have developed a mathematical model that treats mRNA translation and associated -1 frameshifting as a stochastic process in which the transition probabilities are based on the energetics of local molecular interactions. The model predicts both the location and efficiency of -1 frameshift events in HIV-1. Moreover, we compute -1 frameshift efficiencies upon mutations in the viral mRNA sequence and variations in relative tRNA abundances, predictions that are directly testable in experiment.

Genetic evidence for archaic admixture in Africa.

A long-debated question concerns the fate of archaic forms of the genus Homo: did they go extinct without interbreeding with anatomically modern humans, or are their genes present in contemporary populations? This question is typically focused on the genetic contribution of archaic forms outside of Africa. Here we use DNA sequence data gathered from 61 noncoding autosomal regions in a sample of three sub-Saharan African populations (Mandenka, Biaka, and San) to test models of African archaic admixture. We use two complementary approximate-likelihood approaches and a model of human evolution that involves recent population structure, with and without gene flow from an archaic population. Extensive simulation results reject the null model of no admixture and allow us to infer that contemporary African populations contain a small proportion of genetic material (≈ 2%) that introgressed ≈ 35 kya from an archaic population that split from the ancestors of anatomically modern humans ≈ 700 kya. Three candidate regions showing deep haplotype divergence, unusual patterns of linkage disequilibrium, and small basal clade size are identified and the distributions of introgressive haplotypes surveyed in a sample of populations from across sub-Saharan Africa. One candidate locus with an unusual segment of DNA that extends for >31 kb on chromosome 4 seems to have introgressed into modern Africans from a now-extinct taxon that may have lived in central Africa. Taken together our results suggest that polymorphisms present in extant populations introgressed via relatively recent interbreeding with hominin forms that diverged from the ancestors of modern humans in the Lower-Middle Pleistocene.

Machine learning models reveal distinct disease subgroups and improve diagnostic and prognostic accuracy for individuals with pathogenic SCN8A gain-of-function variants.

Distinguishing clinical subgroups for patients suffering with diseases characterized by a wide phenotypic spectrum is essential for developing precision therapies. Patients with gain-of-function (GOF) variants in the SCN8A gene exhibit substantial clinical heterogeneity, viewed historically as a linear spectrum ranging from mild to severe. To test for hidden clinical subgroups, we applied two machine-learning algorithms to analyze a dataset of patient features collected by the International SCN8A Patient Registry. We used two research methodologies: a supervised approach that incorporated feature severity cutoffs based on clinical conventions, and an unsupervised approach employing an entirely data-driven strategy. Both approaches found statistical support for three distinct subgroups and were validated by correlation analyses using external variables. However, distinguishing features of the three subgroups within each approach were not concordant, suggesting a more complex phenotypic landscape. The unsupervised approach yielded strong support for a model involving three partially ordered subgroups rather than a linear spectrum. Application of these machine-learning approaches may lead to improved prognosis and clinical management of individuals with SCN8A GOF variants and provide insights into the underlying mechanisms of the disease.

Identification of an Allele-Specific Transcription Factor Binding Interaction that May Regulate PLA2G2A Gene Expression.

The secreted phospholipase A2 (sPLA2) isoform, sPLA2-IIA, has been implicated in a variety of diseases and conditions, including bacteremia, cardiovascular disease, COVID-19, sepsis, adult respiratory distress syndrome, and certain cancers. Given its significant role in these conditions, understanding the regulatory mechanisms impacting its levels is crucial. Genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs), including rs11573156, that are associated with circulating levels of sPLA2-IIA. The work in the manuscript leveraged 4 publicly available datasets to investigate the mechanism by which rs11573156 influences sPLA2-IIA levels via bioinformatics and modeling analysis. Through genotype-tissue expression (GTEx), 234 expression quantitative trait loci (eQTLs) were identified for the gene that encodes for sPLA2-IIA, PLA2G2A. SNP2TFBS was used to ascertain the binding affinities between transcription factors (TFs) to both the reference and alternative alleles of identified eQTL SNPs. Subsequently, candidate TF-SNP interactions were cross-referenced with the ChIP-seq results in matched tissues from ENCODE. SP1-rs11573156 emerged as the significant TF-SNP pair in the liver. Further analysis revealed that the upregulation of PLA2G2A transcript levels through the rs11573156 variant was likely affected by tissue SP1 protein levels. Using an ordinary differential equation based on Michaelis-Menten kinetic assumptions, we modeled the dependence of PLA2G2A transcription on SP1 protein levels, incorporating the SNP influence. Collectively, our analysis strongly suggests that the difference in the binding dynamics of SP1 to different rs11573156 alleles may underlie the allele-specific PLA2G2A expression in different tissues, a mechanistic model that awaits future direct experimental validation. This mechanism likely contributes to the variation in circulating sPLA2-IIA protein levels in the human population, with implications for a wide range of human diseases.