Genomic Signatures of Positive Selection in Human Populations of the OXT, OXTR, AVP, AVPR1A and AVR1B Gene Variants Related to the Regulation of Psychoemotional Response.

The neurobiological systems of maintenance and control of behavioral responses result from natural selection. We have analyzed the selection signatures for single nucleotide variants (SNV) of the genes of oxytocin (OXT, OXTR) and vasopressin (AVP, AVPR1A, AVPR1B) systems, which are associated with the regulation of social and emotional behavior in distinct populations. The analysis was performed using original WGS (whole genome sequencing) data on Eastern Slavs (SlEast), as well as publicly available data from the 1000 Genomes Project on GBR, FIN, IBR, PUR, BEB, CHB, and ACB populations (the latter were taken as reference). To identify selection signatures, we rated the integrated haplotype scores (iHS), the numbers of segregating sites by length (nSl), and the integrated haplotype homozygosity pooled (iHH12) measures; the fixation index Fst was implemented to assess genetic differentiation between populations. We revealed that the strongest genetic differentiation of populations was found with respect to the AVPR1B gene, with the greatest differentiation observed in GRB (Fst = 0.316) and CHB (Fst = 0.325) in comparison to ACB. Also, high Fst values were found for SNVs of the AVPR1B gene rs28499431, rs33940624, rs28477649, rs3883899, and rs28452187 in most of the populations. Selection signatures have also been identified in the AVP, AVPR1A, OXT, and OXTR genes. Our analysis shows that the OXT, OXTR, AVP, AVPR1A, and AVPR1B genes were subject to positive selection in a population-specific process, which was likely contributing to the diversity of adaptive emotional response types and social function realizations.

GWAS reveals genetic basis of a predisposition to severe COVID-19 through in silico modeling of the FYCO1 protein.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is heavily reliant on its natural ability to "hack" the host's genetic and biological pathways. The genetic susceptibility of the host is a key factor underlying the severity of the disease. Polygenic risk scores are essential for risk assessment, risk stratification, and the prevention of adverse outcomes. In this study, we aimed to assess and analyze the genetic predisposition to severe COVID-19 in a large representative sample of the Russian population as well as to build a reliable but simple polygenic risk score model with a lower margin of error. Another important goal was to learn more about the pathogenesis of severe COVID-19. We examined the tertiary structure of the FYCO1 protein, the only gene with mutations in its coding region and discovered changes in the coiled-coil domain. Our findings suggest that FYCO1 may accelerate viral intracellular replication and excessive exocytosis and may contribute to an increased risk of severe COVID-19. We found significant associations between COVID-19 and LZTFL1, FYCO1, XCR1, CCR9, TMLHE-AS1, and SCYL2 at 3p21.31. Our findings further demonstrate the polymorphic nature of the severe COVID-19 phenotype.

Genetics of psycho-emotional well-being: genome-wide association study and polygenic risk score analysis.

Background: Psycho-emotional well-being is essential for living a life of satisfaction and fulfillment. However, depression and anxiety have become the leading mental health issues worldwide, according to the World Health Organization. Both disorders have been linked to stress and other psychological factors. Their genetic basis remains understudied. Methods: In 2020-2021, the psycho-emotional well-being of 30,063 Russians with no known psychiatric history was assessed using the Hospital Anxiety and Depression Scale (HADS) for general mental health and the HADS subscale A (anxiety) for anxiety. Following the original instructions, an anxiety score of ≥11 points was used as the anxiety threshold. A genome-wide association study was performed to find associations between anxiety and HADS/HADS-A scores using linear and logistic regressions based on HADS/HADS-A scores as binary and continuous variables, respectively. In addition, the links between anxiety, sociodemographic factors (such as age, sex, and employment), lifestyle (such as physical activity, sleep duration, and smoking), and markers of caffeine and alcohol metabolism were analyzed. To assess the risk of anxiety, polygenic risk score modeling was carried out using open-access software and principal component analysis (PCA) to simplify the calculations (ROC AUC = 89.4 ± 2.2% on the test set). Results: There was a strong positive association between HADS/HADS-A scores and sociodemographic factors and lifestyle. New single-nucleotide polymorphisms (SNPs) with genome-wide significance were discovered, which had not been associated with anxiety or other stress-related conditions but were located in genes previously associated with bipolar disorder, schizophrenia, or emotional instability. The CACNA1C variant rs1205787230 was associated with clinical anxiety (a HADS-A score of ≥11 points). There was an association between anxiety levels (HADS-A scores) and genes involved in the activity of excitatory neurotransmitters: PTPRN2 (rs3857647), DLGAP4 (rs8114927), and STK24 (rs9517326). Conclusion: Our results suggest that calcium channels and monoamine neurotransmitters, as well as SNPs in genes directly or indirectly affecting neurogenesis and synaptic functions, may be involved in the development of increased anxiety. The role of some non-genetic factors and the clinical significance of physiological markers such as lifestyle were also demonstrated.

Single-Cell Gene Expression Analysis Revealed Immune Cell Signatures of Delta COVID-19.

The coronavirus disease 2019 (COVID-19) is accompanied by a cytokine storm with the release of many proinflammatory factors and development of respiratory syndrome. Several SARS-CoV-2 lineages have been identified, and the Delta variant (B.1.617), linked with high mortality risk, has become dominant in many countries. Understanding the immune responses associated with COVID-19 lineages may therefore aid the development of therapeutic and diagnostic strategies. Multiple single-cell gene expression studies revealed innate and adaptive immunological factors and pathways correlated with COVID-19 severity. Additional investigations covering host-pathogen response characteristics for infection caused by different lineages are required. Here, we performed single-cell transcriptome profiling of blood mononuclear cells from the individuals with different severity of the COVID-19 and virus lineages to uncover variant specific molecular factors associated with immunity. We identified significant changes in lymphoid and myeloid cells. Our study highlights that an abundant population of monocytes with specific gene expression signatures accompanies Delta lineage of SARS-CoV-2 and contributes to COVID-19 pathogenesis inferring immune components for targeted therapy.

A Splice Variant of the MYH7 Gene Is Causative in a Family with Isolated Left Ventricular Noncompaction Cardiomyopathy.

Variants of the MYH7 gene have been associated with a number of primary cardiac conditions, including left ventricular noncompaction cardiomyopathy (LVNC). Most cases of MYH7-related diseases are associated with such variant types as missense substitutions and in-frame indels. Thus, truncating variants in MYH7 (MYH7tv) and associated mechanism of haploinsufficiency are usually considered not pathogenic in these disorders. However, recent large-scale studies demonstrated evidence of the significance of MYH7tv for LVNC and gave rise to an assumption that haploinsufficiency may be the causal mechanism for LVNC. In this article, we present a family with isolated LVNC and a heterozygous splice variant of the MYH7 gene, analyze possible consequences of this variant and conclude that not all variants that are predicted truncating really act through haploinsufficiency. This study can highlight the importance of a precise assessment of MYH7 splicing variants and their participation in the development of LVNC.

Prostaglandin H synthase kinetics in the two-phase aqueous-micellar system.

Reaction mixture for PGHS (prostaglandin-H-synthase) is a two-phase system including micellar hydrophobic phase and hydrophilic aqueous phase. Reagents added to the mixture are distributed between phases, thus concentrations of reagents dissolved in phases can differ significantly from their overall contents. Using dynamic light scattering we found that the hydrophobic phase produced by tween-20 consists of micelles, which radius (4-5nm) does not depend on either tween-20 overall content (0.1%-1% v/v) or arachidonic acid (AA) addition (10-1000µM) or PGHS addition (1µM). Tween-20 overall content changing from 0.1% to 2% v/v dramatically affected COX kinetic, but accounting AA distribution between phases allowed us to estimate "true" parameters, independent of the tween-20 overall content and the concentration of another substrate: KM(Ox) equals 9.8µM O2 in the aqueous phase or 0.0074bar in the gaseous phase, KM(AA) equals 5400µM AA in the phase of tween-20 micelles and 5400/PµM AA in the aqueous phase (P is the distribution ratio for the AA between the aqueous phase and the hydrophobic phase (P≫1000)). This approach allowed to evaluate PS, the distribution ratio for the AA between the hydrophobic phase and the PGHS active center (PS ~310). This coefficient indicates the AA selectivity toward the cyclooxygenase active center.

New words in human mutagenesis.

BACKGROUND: The substitution rates within different nucleotide contexts are subject to varying levels of bias. The most well known example of such bias is the excess of C to T (C > T) mutations in CpG (CG) dinucleotides. The molecular mechanisms underlying this bias are important factors in human genome evolution and cancer development. The discovery of other nucleotide contexts that have profound effects on substitution rates can improve our understanding of how mutations are acquired, and why mutation hotspots exist. RESULTS: We compared rates of inherited mutations in 1-4 bp nucleotide contexts using reconstructed ancestral states of human single nucleotide polymorphisms (SNPs) from intergenic regions. Chimp and orangutan genomic sequences were used as outgroups. We uncovered 3.5 and 3.3-fold excesses of T > C mutations in the second position of ATTG and ATAG words, respectively, and a 3.4-fold excess of A > C mutations in the first position of the ACAA word. CONCLUSIONS: Although all the observed biases are less pronounced than the 5.1-fold excess of C > T mutations in CG dinucleotides, the three 4 bp mutation contexts mentioned above (and their complementary contexts) are well distinguished from all other mutation contexts. This provides a challenge to discover the underlying mechanisms responsible for the observed excesses of mutations.

Exclusive sequences of different genomes.

We studied the distribution of 1-7 bp words in a dataset that includes 139 complete eukaryotic genomes, 33 masked eukaryotic genomes and coding regions from 35 genomes. We tested different statistical models to determine over- and under-represented words. The method described by Karlin et al. has the strongest predictive power compared to other methods. Using this method we identified over- and under-represented words consistent within a large array of taxonomic groups. Some of those words have not yet been described as exclusive. For example, CGCG is over-represented in CG-deficient organisms. We also describe exceptions for widely known exclusive words, such as CG and TA.