A Principal Components Analysis and Functional Annotation of Differentially Expressed Genes in Brain Regions of Gray Rats Selected for Tame or Aggressive Behavior.

The process of domestication, despite its short duration as it compared with the time scale of the natural evolutionary process, has caused rapid and substantial changes in the phenotype of domestic animal species. Nonetheless, the genetic mechanisms underlying these changes remain poorly understood. The present study deals with an analysis of the transcriptomes from four brain regions of gray rats (Rattus norvegicus), serving as an experimental model object of domestication. We compared gene expression profiles in the hypothalamus, hippocampus, periaqueductal gray matter, and the midbrain tegmental region between tame domesticated and aggressive gray rats and revealed subdivisions of differentially expressed genes by principal components analysis that explain the main part of differentially gene expression variance. Functional analysis (in the DAVID (Database for Annotation, Visualization and Integrated Discovery) Bioinformatics Resources database) of the differentially expressed genes allowed us to identify and describe the key biological processes that can participate in the formation of the different behavioral patterns seen in the two groups of gray rats. Using the STRING- DB (search tool for recurring instances of neighboring genes) web service, we built a gene association network. The genes engaged in broad network interactions have been identified. Our study offers data on the genes whose expression levels change in response to artificial selection for behavior during animal domestication.

AtSNP_TATAdb: Candidate Molecular Markers of Plant Advantages Related to Single Nucleotide Polymorphisms within Proximal Promoters of Arabidopsis thaliana L.

The mainstream of the post-genome target-assisted breeding in crop plant species includes biofortification such as high-throughput phenotyping along with genome-based selection. Therefore, in this work, we used the Web-service Plant_SNP_TATA_Z-tester, which we have previously developed, to run a uniform in silico analysis of the transcriptional alterations of 54,013 protein-coding transcripts from 32,833 Arabidopsis thaliana L. genes caused by 871,707 SNPs located in the proximal promoter region. The analysis identified 54,993 SNPs as significantly decreasing or increasing gene expression through changes in TATA-binding protein affinity to the promoters. The existence of these SNPs in highly conserved proximal promoters may be explained as intraspecific diversity kept by the stabilizing natural selection. To support this, we hand-annotated papers on some of the Arabidopsis genes possessing these SNPs or on their orthologs in other plant species and demonstrated the effects of changes in these gene expressions on plant vital traits. We integrated in silico estimates of the TBP-promoter affinity in the AtSNP_TATAdb knowledge base and showed their significant correlations with independent in vivo experimental data. These correlations appeared to be robust to variations in statistical criteria, genomic environment of TATA box regions, plants species and growing conditions.

Differentially Expressed Genes and Molecular Susceptibility to Human Age-Related Diseases.

Mainstream transcriptome profiling of susceptibility versus resistance to age-related diseases (ARDs) is focused on differentially expressed genes (DEGs) specific to gender, age, and pathogeneses. This approach fits in well with predictive, preventive, personalized, participatory medicine and helps understand how, why, when, and what ARDs one can develop depending on their genetic background. Within this mainstream paradigm, we wanted to find out whether the known ARD-linked DEGs available in PubMed can reveal a molecular marker that will serve the purpose in anyone's any tissue at any time. We sequenced the periaqueductal gray (PAG) transcriptome of tame versus aggressive rats, identified rat-behavior-related DEGs, and compared them with their known homologous animal ARD-linked DEGs. This analysis yielded statistically significant correlations between behavior-related and ARD-susceptibility-related fold changes (log2 values) in the expression of these DEG homologs. We found principal components, PC1 and PC2, corresponding to the half-sum and the half-difference of these log2 values, respectively. With the DEGs linked to ARD susceptibility and ARD resistance in humans used as controls, we verified these principal components. This yielded only one statistically significant common molecular marker for ARDs: an excess of Fcγ receptor IIb suppressing immune cell hyperactivation.

Candidate SNP Markers Significantly Altering the Affinity of TATA-Binding Protein for the Promoters of Human Hub Genes for Atherogenesis, Atherosclerosis and Atheroprotection.

Atherosclerosis is a systemic disease in which focal lesions in arteries promote the build-up of lipoproteins and cholesterol they are transporting. The development of atheroma (atherogenesis) narrows blood vessels, reduces the blood supply and leads to cardiovascular diseases. According to the World Health Organization (WHO), cardiovascular diseases are the leading cause of death, which has been especially boosted since the COVID-19 pandemic. There is a variety of contributors to atherosclerosis, including lifestyle factors and genetic predisposition. Antioxidant diets and recreational exercises act as atheroprotectors and can retard atherogenesis. The search for molecular markers of atherogenesis and atheroprotection for predictive, preventive and personalized medicine appears to be the most promising direction for the study of atherosclerosis. In this work, we have analyzed 1068 human genes associated with atherogenesis, atherosclerosis and atheroprotection. The hub genes regulating these processes have been found to be the most ancient. In silico analysis of all 5112 SNPs in their promoters has revealed 330 candidate SNP markers, which statistically significantly change the affinity of the TATA-binding protein (TBP) for these promoters. These molecular markers have made us confident that natural selection acts against underexpression of the hub genes for atherogenesis, atherosclerosis and atheroprotection. At the same time, upregulation of the one for atheroprotection promotes human health.

Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing.

Studies on hereditary fixation of the tame-behavior phenotype during animal domestication remain relevant and important because they are of both basic research and applied significance. In model animals, gray rats Rattus norvegicus bred for either an enhancement or reduction in defensive response to humans, for the first time, we used high-throughput RNA sequencing to investigate differential expression of genes in tissue samples from the tegmental region of the midbrain in 2-month-old rats showing either tame or aggressive behavior. A total of 42 differentially expressed genes (DEGs; adjusted p-value < 0.01 and fold-change > 2) were identified, with 20 upregulated and 22 downregulated genes in the tissue samples from tame rats compared with aggressive rats. Among them, three genes encoding transcription factors (TFs) were detected: Ascl3 was upregulated, whereas Fos and Fosb were downregulated in tissue samples from the brains of tame rats brain. Other DEGs were annotated as associated with extracellular matrix components, transporter proteins, the neurotransmitter system, signaling molecules, and immune system proteins. We believe that these DEGs encode proteins that constitute a multifactorial system determining the behavior for which the rats have been artificially selected. We demonstrated that several structural subtypes of E-box motifs­known as binding sites for many developmental TFs of the bHLH class, including the ASCL subfamily of TFs­are enriched in the set of promoters of the DEGs downregulated in the tissue samples of tame rats'. Because ASCL3 may act as a repressor on target genes of other developmental TFs of the bHLH class, we hypothesize that the expression of TF gene Ascl3 in tame rats indicates longer neurogenesis (as compared to aggressive rats), which is a sign of neoteny and domestication. Thus, our domestication model shows a new function of TF ASCL3: it may play the most important role in behavioral changes in animals.

Plant_SNP_TATA_Z-Tester: A Web Service That Unequivocally Estimates the Impact of Proximal Promoter Mutations on Plant Gene Expression.

Synthetic targeted optimization of plant promoters is becoming a part of progress in mainstream postgenomic agriculture along with hybridization of cultivated plants with wild congeners, as well as marker-assisted breeding. Therefore, here, for the first time, we compiled all the experimental data-on mutational effects in plant proximal promoters on gene expression-that we could find in PubMed. Some of these datasets cast doubt on both the existence and the uniqueness of the sought solution, which could unequivocally estimate effects of proximal promoter mutation on gene expression when plants are grown under various environmental conditions during their development. This means that the inverse problem under study is ill-posed. Furthermore, we found experimental data on in vitro interchangeability of plant and human TATA-binding proteins allowing the application of Tikhonov's regularization, making this problem well-posed. Within these frameworks, we created our Web service Plant_SNP_TATA_Z-tester and then determined the limits of its applicability using those data that cast doubt on both the existence and the uniqueness of the sought solution. We confirmed that the effects (of proximal promoter mutations on gene expression) predicted by Plant_SNP_TATA_Z-tester correlate statistically significantly with all the experimental data under study. Lastly, we exemplified an application of Plant_SNP_TATA_Z-tester to agriculturally valuable mutations in plant promoters.

Stress Reactivity, Susceptibility to Hypertension, and Differential Expression of Genes in Hypertensive Compared to Normotensive Patients.

Although half of hypertensive patients have hypertensive parents, known hypertension-related human loci identified by genome-wide analysis explain only 3% of hypertension heredity. Therefore, mainstream transcriptome profiling of hypertensive subjects addresses differentially expressed genes (DEGs) specific to gender, age, and comorbidities in accordance with predictive preventive personalized participatory medicine treating patients according to their symptoms, individual lifestyle, and genetic background. Within this mainstream paradigm, here, we determined whether, among the known hypertension-related DEGs that we could find, there is any genome-wide hypertension theranostic molecular marker applicable to everyone, everywhere, anytime. Therefore, we sequenced the hippocampal transcriptome of tame and aggressive rats, corresponding to low and high stress reactivity, an increase of which raises hypertensive risk; we identified stress-reactivity-related rat DEGs and compared them with their known homologous hypertension-related animal DEGs. This yielded significant correlations between stress reactivity-related and hypertension-related fold changes (log2 values) of these DEG homologs. We found principal components, PC1 and PC2, corresponding to a half-difference and half-sum of these log2 values. Using the DEGs of hypertensive versus normotensive patients (as the control), we verified the correlations and principal components. This analysis highlighted downregulation of ß-protocadherins and hemoglobin as whole-genome hypertension theranostic molecular markers associated with a wide vascular inner diameter and low blood viscosity, respectively.

A Bioinformatics Model of Human Diseases on the Basis of Differentially Expressed Genes (of Domestic Versus Wild Animals) That Are Orthologs of Human Genes Associated with Reproductive-Potential Changes.

Earlier, after our bioinformatic analysis of single-nucleotide polymorphisms of TATA-binding protein-binding sites within gene promoters on the human Y chromosome, we suggested that human reproductive potential diminishes during self-domestication. Here, we implemented bioinformatics models of human diseases using animal in vivo genome-wide RNA-Seq data to compare the effect of co-directed changes in the expression of orthologous genes on human reproductive potential and during the divergence of domestic and wild animals from their nearest common ancestor (NCA). For example, serotonin receptor 3A (HTR3A) deficiency contributes to sudden death in pregnancy, consistently with Htr3a underexpression in guinea pigs (Cavia porcellus) during their divergence from their NCA with cavy (C. aperea). Overall, 25 and three differentially expressed genes (hereinafter, DEGs) in domestic animals versus 11 and 17 DEGs in wild animals show the direction consistent with human orthologous gene-markers of reduced and increased reproductive potential. This indicates a reliable association between DEGs in domestic animals and human orthologous genes reducing reproductive potential (Pearson's χ2 test p < 0.001, Fisher's exact test p < 0.05, binomial distribution p < 0.0001), whereas DEGs in wild animals uniformly match human orthologous genes decreasing and increasing human reproductive potential (p > 0.1; binomial distribution), thus enforcing the norm (wild type).

Unannotated single nucleotide polymorphisms in the TATA box of erythropoiesis genes show in vitro positive involvements in cognitive and mental disorders.

BACKGROUND: Hemoglobin is a tetramer consisting of two α-chains and two ß-chains of globin. Hereditary aberrations in the synthesis of one of the globin chains are at the root of thalassemia, one of the most prevalent monogenic diseases worldwide. In humans, in addition to α- and ß-globins, embryonic zeta-globin and fetal γ-globin are expressed. Immediately after birth, the expression of fetal Aγ- and Gγ-globin ceases, and then adult ß-globin is mostly expressed. It has been shown that in addition to erythroid cells, hemoglobin is widely expressed in nonerythroid cells including neurons of the cortex, hippocampus, and cerebellum in rodents; embryonic and adult brain neurons in mice; and mesencephalic dopaminergic brain cells in humans, mice, and rats. Lately, there is growing evidence that different forms of anemia (changes in the number and quality of blood cells) may be involved in (or may accompany) the pathogenesis of various cognitive and mental disorders, such as Alzheimer's and Parkinson's diseases, depression of various severity levels, bipolar disorders, and schizophrenia. Higher hemoglobin concentrations in the blood may lead to hyperviscosity, hypovolemia, and lung diseases, which may cause brain hypoxia and anomalies of brain function, which may also result in cognitive deficits. METHODS: In this study, a search for unannotated single-nucleotide polymorphisms (SNPs) of erythroid genes was initially performed using our previously created and published SNP-TATA_Z-tester, which is a Web service for computational analysis of a given SNP for in silico estimation of its influence on the affinity of TATA-binding protein (TBP) for TATA and TATA-like sequences. The obtained predictions were finally verified in vitro by an electrophoretic mobility shift assay (EMSA). RESULTS: On the basis of these experimental in vitro results and literature data, we studied TATA box SNPs influencing both human erythropoiesis and cognitive abilities. For instance, TBP-TATA affinity in the HbZ promoter decreases 6.6-fold as a result of a substitution in the TATA box (rs113180943), thereby possibly disrupting stage-dependent events of "switching" of hemoglobin genes and thus causing erythroblastosis. Therefore, rs113180943 may be a candidate marker of severe hemoglobinopathies with comorbid cognitive and mental disorders associated with cerebral blood flow disturbances. CONCLUSIONS: The literature data and experimental and computations results suggest that the uncovered candidate SNP markers of erythropoiesis anomalies may also be studied in cohorts of patients with cognitive and/or mental disorders with comorbid erythropoiesis diseases in comparison to conventionally healthy volunteers. Research into the regulatory mechanisms by which the identified SNP markers contribute to the development of hemoglobinopathies and of the associated cognitive deficits will allow physicians not only to take timely and adequate measures against hemoglobinopathies but also to implement strategies preventing cognitive and mental disorders.

Disruptive natural selection by male reproductive potential prevents underexpression of protein-coding genes on the human Y chromosome as a self-domestication syndrome.

BACKGROUND: In population ecology, the concept of reproductive potential denotes the most vital indicator of chances to produce and sustain a healthy descendant until his/her reproductive maturity under the best conditions. This concept links quality of life and longevity of an individual with disease susceptibilities encoded by his/her genome. Female reproductive potential has been investigated deeply, widely, and comprehensively in the past, but the male one has not received an equal amount of attention. Therefore, here we focused on the human Y chromosome and found candidate single-nucleotide polymorphism (SNP) markers of male reproductive potential. RESULTS: Examining in silico (i.e., using our earlier created Web-service SNP_TATA_Z-tester) all 1206 unannotated SNPs within 70 bp proximal promoters of all 63 Y-linked genes, we found 261 possible male-reproductive-potential SNP markers that can significantly alter the binding affinity of TATA-binding protein (TBP) for these promoters. Among them, there are candidate SNP markers of spermatogenesis disorders (e.g., rs1402972626), pediatric cancer (e.g., rs1483581212) as well as male anxiety damaging family relationships and mother's and children's health (e.g., rs187456378). First of all, we selectively verified in vitro both absolute and relative values of the analyzed TBP-promoter affinity, whose Pearson's coefficients of correlation between predicted and measured values were r = 0.84 (significance p <  0.025) and r = 0.98 (p <  0.025), respectively. Next, we found that there are twofold fewer candidate SNP markers decreasing TBP-promoter affinity relative to those increasing it, whereas in the genome-wide norm, SNP-induced damage to TBP-promoter complexes is fourfold more frequent than SNP-induced improvement (p <  0.05, binomial distribution). This means natural selection against underexpression of these genes. Meanwhile, the numbers of candidate SNP markers of an increase and decrease in male reproductive potential were indistinguishably equal to each other (p <  0.05) as if male self-domestication could have happened, with its experimentally known disruptive natural selection. Because there is still not enough scientific evidence that this could have happened, we discuss the human diseases associated with candidate SNP markers of male reproductive potential that may correspond to domestication-related disorders in pets. CONCLUSIONS: Overall, our findings seem to support a self-domestication syndrome with disruptive natural selection by male reproductive potential preventing Y-linked underexpression of a protein.

Candidate SNP Markers of Atherogenesis Significantly Shifting the Affinity of TATA-Binding Protein for Human Gene Promoters show stabilizing Natural Selection as a Sum of Neutral Drift Accelerating Atherogenesis and Directional Natural Selection Slowing It.

(1) Background: The World Health Organization (WHO) regards atherosclerosis-related myocardial infarction and stroke as the main causes of death in humans. Susceptibility to atherogenesis-associated diseases is caused by single-nucleotide polymorphisms (SNPs). (2) Methods: Using our previously developed public web-service SNP_TATA_Comparator, we estimated statistical significance of the SNP-caused alterations in TATA-binding protein (TBP) binding affinity for 70 bp proximal promoter regions of the human genes clinically associated with diseases syntonic or dystonic with atherogenesis. Additionally, we did the same for several genes related to the maintenance of mitochondrial genome integrity, according to present-day active research aimed at retarding atherogenesis. (3) Results: In dbSNP, we found 1186 SNPs altering such affinity to the same extent as clinical SNP markers do (as estimated). Particularly, clinical SNP marker rs2276109 can prevent autoimmune diseases via reduced TBP affinity for the human MMP12 gene promoter and therefore macrophage elastase deficiency, which is a well-known physiological marker of accelerated atherogenesis that could be retarded nutritionally using dairy fermented by lactobacilli. (4) Conclusions: Our results uncovered SNPs near clinical SNP markers as the basis of neutral drift accelerating atherogenesis and SNPs of genes encoding proteins related to mitochondrial genome integrity and microRNA genes associated with instability of the atherosclerotic plaque as a basis of directional natural selection slowing atherogenesis. Their sum may be stabilizing the natural selection that sets the normal level of atherogenesis.

Candidate SNP markers of reproductive potential are predicted by a significant change in the affinity of TATA-binding protein for human gene promoters.

BACKGROUND: The progress of medicine, science, technology, education, and culture improves, year by year, quality of life and life expectancy of the populace. The modern human has a chance to further improve the quality and duration of his/her life and the lives of his/her loved ones by bringing their lifestyle in line with their sequenced individual genomes. With this in mind, one of genome-based developments at the junction of personalized medicine and bioinformatics will be considered in this work, where we used two Web services: (i) SNP_TATA_Comparator to search for alleles with a single nucleotide polymorphism (SNP) that alters the affinity of TATA-binding protein (TBP) for the TATA boxes of human gene promoters and (ii) PubMed to look for retrospective clinical reviews on changes in physiological indicators of reproductive potential in carriers of these alleles. RESULTS: A total of 126 SNP markers of female reproductive potential, capable of altering the affinity of TBP for gene promoters, were found using the two above-mentioned Web services. For example, 10 candidate SNP markers of thrombosis (e.g., rs563763767) can cause overproduction of coagulation inducers. In pregnant women, Hughes syndrome provokes thrombosis with a fatal outcome although this syndrome can be diagnosed and eliminated even at the earliest stages of its development. Thus, in women carrying any of the above SNPs, preventive treatment of this syndrome before a planned pregnancy can reduce the risk of death. Similarly, seven SNP markers predicted here (e.g., rs774688955) can elevate the risk of myocardial infarction. In line with Bowles' lifespan theory, women carrying any of these SNPs may modify their lifestyle to improve their longevity if they can take under advisement that risks of myocardial infarction increase with age of the mother, total number of pregnancies, in multiple pregnancies, pregnancies under the age of 20, hypertension, preeclampsia, menstrual cycle irregularity, and in women smokers. CONCLUSIONS: According to Bowles' lifespan theory-which links reproductive potential, quality of life, and life expectancy-the above information was compiled for those who would like to reduce risks of diseases corresponding to alleles in own sequenced genomes. Candidate SNP markers can focus the clinical analysis of unannotated SNPs, after which they may become useful for people who would like to bring their lifestyle in line with their sequenced individual genomes.

Genomic landscape of CpG rich elements in human.

BACKGROUND: The studies on CpG islands (CGI) and Alu elements functions, evolution, and distribution in the genome started since the discovery in nineteen eighties (1981, 1986, correspondingly). Their highly skewed genome wide distribution implies the non-random retrotransposition pattern. Besides CGIs in gene promoters, CGIs clusters were observed in the homeobox gene regions and in the macrosatellites, but the whole picture of their distribution specifics was not grasped. Attempts to identify any causative features upon their (genome wide) distribution, such as the DNA context mediated preferred insertion sites of Alu repeats, have been made to ascribe their clusters location. METHODS: Recent emergence of high resolution 3D map of human genome allowed segregating the genome into the large scale chromatin domains of naturally observable nuclear subcompartments, or Topologically Associated Domains (TADs), designated by spatial chromatin distribution. We utilized the chromatin map to elucidate relations between large scale chromatin state and CpG rich elements landscape. In the course of analysis it was confirmed that genes, Alu and CGI clusters maintain obvious, albeit different in strength, preference for open chromatin. For the first time it was clearly shown that the clusters density of the Alu and CGIs monotonically depend on the chromatin accessibility rate. In particular, the highest density of these elements is found in A1 euchromatin regions characterized by a high density of small length genes replicating in the early S-phase. It implies that these elements mediate (CGIs) or are a side element (Alus) of chromatin accessibility. RESULTS: We elucidated that both methylated and non-methylated CGIs display the affinity to chromatin accessibility. As a part of comparative genomics section, we elucidated that the dog's genome non-canonical structure, outstanding in mammals for its high CGIs abundance compared to gene number, is explained by the presence of dense tandem CGI extended hotspots (500 kb on average) in subtelomeric and pericentromeric regions with highly skewed CG content, and not by CGIs global distribution pattern shift. CONCLUSIONS: The study underlines the close association of CG-rich elements distribution with the newly introduced large scale chromatin state map, proposing a refined standpoint on interrelation of aforementioned genome elements and the chromatin state. To our expertise, the TAD-associated partition model employed in the study is likely the most substantial one regarding CpG rich clusters distribution among the whole genome chromatin/isochores maps available.

Candidate SNP markers of aggressiveness-related complications and comorbidities of genetic diseases are predicted by a significant change in the affinity of TATA-binding protein for human gene promoters.

BACKGROUND: Aggressiveness in humans is a hereditary behavioral trait that mobilizes all systems of the body-first of all, the nervous and endocrine systems, and then the respiratory, vascular, muscular, and others-e.g., for the defense of oneself, children, family, shelter, territory, and other possessions as well as personal interests. The level of aggressiveness of a person determines many other characteristics of quality of life and lifespan, acting as a stress factor. Aggressive behavior depends on many parameters such as age, gender, diseases and treatment, diet, and environmental conditions. Among them, genetic factors are believed to be the main parameters that are well-studied at the factual level, but in actuality, genome-wide studies of aggressive behavior appeared relatively recently. One of the biggest projects of the modern science-1000 Genomes-involves identification of single nucleotide polymorphisms (SNPs), i.e., differences of individual genomes from the reference genome. SNPs can be associated with hereditary diseases, their complications, comorbidities, and responses to stress or a drug. Clinical comparisons between cohorts of patients and healthy volunteers (as a control) allow for identifying SNPs whose allele frequencies significantly separate them from one another as markers of the above conditions. Computer-based preliminary analysis of millions of SNPs detected by the 1000 Genomes project can accelerate clinical search for SNP markers due to preliminary whole-genome search for the most meaningful candidate SNP markers and discarding of neutral and poorly substantiated SNPs. RESULTS: Here, we combine two computer-based search methods for SNPs (that alter gene expression) {i} Web service SNP_TATA_Comparator (DNA sequence analysis) and {ii} PubMed-based manual search for articles on aggressiveness using heuristic keywords. Near the known binding sites for TATA-binding protein (TBP) in human gene promoters, we found aggressiveness-related candidate SNP markers, including rs1143627 (associated with higher aggressiveness in patients undergoing cytokine immunotherapy), rs544850971 (higher aggressiveness in old women taking lipid-lowering medication), and rs10895068 (childhood aggressiveness-related obesity in adolescence with cardiovascular complications in adulthood). CONCLUSIONS: After validation of these candidate markers by clinical protocols, these SNPs may become useful for physicians (may help to improve treatment of patients) and for the general population (a lifestyle choice preventing aggressiveness-related complications).

Disruptive Selection of Human Immunostimulatory and Immunosuppressive Genes Both Provokes and Prevents Rheumatoid Arthritis, Respectively, as a Self-Domestication Syndrome.

Using our previously published Web service SNP_TATA_Comparator, we conducted a genome-wide study of single-nucleotide polymorphisms (SNPs) within core promoters of 68 human rheumatoid arthritis (RA)-related genes. Using 603 SNPs within 25 genes clinically associated with RA-comorbid disorders, we predicted 84 and 70 candidate SNP markers for overexpression and underexpression of these genes, respectively, among which 58 and 96 candidate SNP markers, respectively, can relieve and worsen RA as if there is a neutral drift toward susceptibility to RA. Similarly, we predicted natural selection toward susceptibility to RA for 8 immunostimulatory genes (e.g., IL9R) and 10 genes most often associated with RA (e.g., NPY). On the contrary, using 25 immunosuppressive genes, we predicted 70 and 109 candidate SNP markers aggravating and relieving RA, respectively (e.g., IL1R2 and TGFB2), suggesting that natural selection can simultaneously additionally yield resistance to RA. We concluded that disruptive natural selection of human immunostimulatory and immunosuppressive genes is concurrently elevating and reducing the risk of RA, respectively. So, we hypothesize that RA in human could be a self-domestication syndrome referring to evolution patterns in domestic animals. We tested this hypothesis by means of public RNA-Seq data on 1740 differentially expressed genes (DEGs) of pets vs. wild animals (e.g., dogs vs. wolves). The number of DEGs in the domestic animals corresponding to worsened RA condition in humans was significantly larger than that in the related wild animals (10 vs. 3). Moreover, much less DEGs in the domestic animals were accordant to relieved RA condition in humans than those in the wild animals (1 vs. 8 genes). This indicates that the anthropogenic environment, in contrast to a natural one, affects gene expression across the whole genome (e.g., immunostimulatory and immunosuppressive genes) in a manner that likely contributes to RA. The difference in gene numbers is statistically significant as confirmed by binomial distribution (p < 0.01), Pearson's χ2 (p < 0.01), and Fisher's exact test (p < 0.05). This allows us to propose RA as a candidate symptom within a self-domestication syndrome. Such syndrome might be considered as a human's payment with health for the benefits received during evolution.

Domestication Explains Two-Thirds of Differential-Gene-Expression Variance between Domestic and Wild Animals; The Remaining One-Third Reflects Intraspecific and Interspecific Variation.

Belyaev's concept of destabilizing selection during domestication was a major achievement in the XX century. Its practical value has been realized in commercial colors of the domesticated fox that never occur in the wild and has been confirmed in a wide variety of pet breeds. Many human disease models involving animals allow to test drugs before human testing. Perhaps this is why investigators doing transcriptomic profiling of domestic versus wild animals have searched for breed-specific patterns. Here we sequenced hypothalamic transcriptomes of tame and aggressive rats, identified their differentially expressed genes (DEGs), and, for the first time, applied principal component analysis to compare them with all the known DEGs of domestic versus wild animals that we could find. Two principal components, PC1 and PC2, respectively explained 67% and 33% of differential-gene-expression variance (hereinafter: log2 value) between domestic and wild animals. PC1 corresponded to multiple orthologous DEGs supported by homologs; these DEGs kept the log2 value sign from species to species and from tissue to tissue (i.e., a common domestication pattern). PC2 represented stand-alone homologous DEG pairs reversing the log2 value sign from one species to another and from tissue to tissue (i.e., representing intraspecific and interspecific variation).

A Rat Model of Human Behavior Provides Evidence of Natural Selection Against Underexpression of Aggressiveness-Related Genes in Humans.

Aggressiveness is a hereditary behavioral pattern that forms a social hierarchy and affects the individual social rank and accordingly quality and duration of life. Thus, genome-wide studies of human aggressiveness are important. Nonetheless, the aggressiveness-related genome-wide studies have been conducted on animals rather than humans. Recently, in our genome-wide study, we uncovered natural selection against underexpression of human aggressiveness-related genes and proved it using F1 hybrid mice. Simultaneously, this natural selection equally supports two opposing traits in humans (dominance and subordination) as if self-domestication could have happened with its disruptive natural selection. Because there is still not enough scientific evidence that this could happen, here, we verified this natural selection pattern using quantitative PCR and two outbred rat lines (70 generations of artificial selection for aggressiveness or tameness, hereinafter: domestication). We chose seven genes-Cacna2d3, Gad2, Gria2, Mapk1, Nos1, Pomc, and Syn1-over- or underexpression of which corresponds to aggressive or domesticated behavior (in humans or mice) that has the same direction as natural selection. Comparing aggressive male rats with domesticated ones, we found that these genes are overexpressed statistically significantly in the hypothalamus (as a universal behavior regulator), not in the periaqueductal gray, where there was no aggressiveness-related expression of the genes in males. Database STRING showed statistically significant associations of the human genes homologous to these rat genes with long-term depression, circadian entrainment, Alzheimer's disease, and the central nervous system disorders during chronic IL-6 overexpression. This finding more likely supports positive perspectives of further studies on self-domestication syndromes.

Natural Selection Equally Supports the Human Tendencies in Subordination and Domination: A Genome-Wide Study With in silico Confirmation and in vivo Validation in Mice.

We proposed the following heuristic decision-making rule: "IF {an excess of a protein relating to the nervous system is an experimentally known physiological marker of low pain sensitivity, fast postinjury recovery, or aggressive, risk/novelty-seeking, anesthetic-like, or similar agonistic-intolerant behavior} AND IF {a single nucleotide polymorphism (SNP) causes overexpression of the gene encoding this protein} THEN {this SNP can be a SNP marker of the tendency in dominance} WHILE {underexpression corresponds to subordination} AND vice versa." Using this decision-making rule, we analyzed 231 human genes of neuropeptidergic, non-neuropeptidergic, and neurotrophinergic systems that encode neurotrophic and growth factors, interleukins, neurotransmitters, receptors, transporters, and enzymes. These proteins are known as key factors of human social behavior. We analyzed all the 5,052 SNPs within the 70 bp promoter region upstream of the position where the protein-coding transcript starts, which were retrieved from databases Ensembl and dbSNP using our previously created public Web service SNP_TATA_Comparator (http://beehive.bionet.nsc.ru/cgi-bin/mgs/tatascan/start.pl). This definition of the promoter region includes all TATA-binding protein (TBP)-binding sites. A total of 556 and 552 candidate SNP markers contributing to the dominance and the subordination, respectively, were uncovered. On this basis, we determined that 231 human genes under study are subject to natural selection against underexpression (significance p < 0.0005), which equally supports the human tendencies in domination and subordination such as the norm of a reaction (plasticity) of the human social hierarchy. These findings explain vertical transmission of domination and subordination traits previously observed in rodent models. Thus, the results of this study equally support both sides of the century-old unsettled scientific debate on whether both aggressiveness and the social hierarchy among humans are inherited (as suggested by Freud and Lorenz) or are due to non-genetic social education, when the children are influenced by older individuals across generations (as proposed by Berkowitz and Fromm).

Candidate SNP Markers of Familial and Sporadic Alzheimer's Diseases Are Predicted by a Significant Change in the Affinity of TATA-Binding Protein for Human Gene Promoters.

While year after year, conditions, quality, and duration of human lives have been improving due to the progress in science, technology, education, and medicine, only eight diseases have been increasing in prevalence and shortening human lives because of premature deaths according to the retrospective official review on the state of US health, 1990-2010. These diseases are kidney cancer, chronic kidney diseases, liver cancer, diabetes, drug addiction, poisoning cases, consequences of falls, and Alzheimer's disease (AD) as one of the leading pathologies. There are familial AD of hereditary nature (~4% of cases) and sporadic AD of unclear etiology (remaining ~96% of cases; i.e., non-familial AD). Therefore, sporadic AD is no longer a purely medical problem, but rather a social challenge when someone asks oneself: "What can I do in my own adulthood to reduce the risk of sporadic AD at my old age to save the years of my lifespan from the destruction caused by it?" Here, we combine two computational approaches for regulatory SNPs: Web service SNP_TATA_Comparator for sequence analysis and a PubMed-based keyword search for articles on the biochemical markers of diseases. Our purpose was to try to find answers to the question: "What can be done in adulthood to reduce the risk of sporadic AD in old age to prevent the lifespan reduction caused by it?" As a result, we found 89 candidate SNP markers of familial and sporadic AD (e.g., rs562962093 is associated with sporadic AD in the elderly as a complication of stroke in adulthood, where natural marine diets can reduce risks of both diseases in case of the minor allele of this SNP). In addition, rs768454929, and rs761695685 correlate with sporadic AD as a comorbidity of short stature, where maximizing stature in childhood and adolescence as an integral indicator of health can minimize (or even eliminate) the risk of sporadic AD in the elderly. After validation by clinical protocols, these candidate SNP markers may become interesting to the general population [may help to choose a lifestyle (in childhood, adolescence, and adulthood) that can reduce the risks of sporadic AD, its comorbidities, and complications in the elderly].

SNP_TATA_Comparator: genomewide landmarks for preventive personalized medicine.

Year after year, conditions, quality, and duration of human lives have been improving due to the progress of science, technology, education, and medicine, which however has a downside. Owing to improvement in children's nutrition, developmental acceleration occurs that imbalances a child's system. Because of virtual worlds of the Internet, social experience of teenagers expands and clashes with puberty of adolescents. Due to the comfort of cities, urbanization emerges and causes stress to adults because of artificial light, noise, pollution, violations of personal space, and family disruption. At old age, all these factors taken together contribute to loneliness, cancer, diabetes, drug addiction, and sporadic Alzheimer's disease, which shorten the lifespan, as reviewed in the US, 1990-2010. That is why, a person may ask oneself: "What can I do now to keep my health in my old age?" To help them, we provide this comprehensive review on predictive preventive personalized medicine. This branch of molecular medicine uses single nucleotide polymorphisms to prevent diseases on the basis of the difference between the individual and reference human genomes.