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

An experimental study of the effects of SNPs in the TATA boxes of the <i>GRIN1, ASCL3</i> and <i>NOS1</i> genes on interactions with the TATA-binding protein

The GRIN1, ASCL3, and NOS1 genes are associated with various phenotypes of neuropsychiatric disorders. For instance, these genes contribute to the development of schizophrenia, Alzheimer’s and Parkinson’s diseases, and epilepsy. These genes are also associated with various cancers. For example, ASCL3 is overexpressed in breast cancer, and NOS1, in ovarian cancer cell lines. Based on our findings and literature data, we had previously obtained results suggesting that the single-nucleotide polymorphisms (SNPs) that disrupt erythropoiesis are highly likely to be associated with cognitive and neuropsychiatric disorders in humans. In the present work, using SNP_TATA_Z-tester, we investigated the influence of unannotated SNPs in the TATA boxes of the promoters of the GRIN1, ASCL3, and NOS1 genes (which are involved in neuropsychiatric disorders and cancers) on the interaction of the TATA boxes with the TATA-binding protein (TBP). Double-stranded oligodeoxyribonucleotides identical to the TATA-containing promoter regions of the GRIN1, ASCL3, and NOS1 genes (reference and minor alleles) and recombinant human TBP were employed to study in vitro (by an electrophoretic mobility shift assay) kinetic characteristics of the formation of TBP–TATA complexes and their affinity. It was found, for example, that allele A of rs1402667001 in the GRIN1 promoter increases TBP–TATA affinity 1.4-fold, whereas allele C in the TATA box of the ASCL3 promoter decreases the affinity 1.4-fold. The lifetime of the complexes in both cases decreased by ~20 % due to changes in the rates of association and dissociation of the complexes (ka and kd, respectively). Our experimental results are consistent with the literature showing GRIN1 underexpression in schizophrenic disorders as well as an increased risk of cervical, bladder, and kidney cancers and lymphoma during ASCL3 underexpression. The effect of allele A of the –27G>A SNP (rs1195040887) in the NOS1 promoter is suggestive of an increased risk of ischemic damage to the brain in carriers. A comparison of experimental TBP–TATA affinity values (KD) of wild-type and minor alleles with predicted ones showed that the data correlate well (linear correlation coefficient r = 0.94, p < 0.01).

Group Theory of Syntactical Freedom in DNA Transcription and Genome Decoding

Transcription factors (TFs) are proteins that recognize specific DNA fragments in order to decode the genome and ensure its optimal functioning. TFs work at the local and global scales by specifying cell type, cell growth and death, cell migration, organization and timely tasks. We investigate the structure of DNA-binding motifs with the theory of finitely generated groups. The DNA ‘word’ in the binding domain—the motif—may be seen as the generator of a finitely generated group Fdna on four letters, the bases A, T, G and C. It is shown that, most of the time, the DNA-binding motifs have subgroup structures close to free groups of rank three or less, a property that we call ‘syntactical freedom’. Such a property is associated with the aperiodicity of the motif when it is seen as a substitution sequence. Examples are provided for the major families of TFs, such as leucine zipper factors, zinc finger factors, homeo-domain factors, etc. We also discuss the exceptions to the existence of such DNA syntactical rules and their functional roles. This includes the TATA box in the promoter region of some genes, the single-nucleotide markers (SNP) and the motifs of some genes of ubiquitous roles in transcription and regulation.

A bioinformatic search for correspondence between differentially expressed genes of domestic versus wild animals and orthologous human genes altering reproductive potential

One of the greatest achievements of genetics in the 20th century is D.K. Belyaev’s discovery of destabilizing selection during the domestication of animals and that this selection affects only gene expression regulation (not gene structure) and inf luences systems of neuroendocrine control of ontogenesis in a stressful environment. Among the experimental data generalized by Belyaev’s discovery, there are also f indings about accelerated extinc tion of testes’ hormonal function and disrupted seasonality of reproduction of domesticated foxes in comparison with their wild congeners. To date, Belyaev’s discovery has already been repeatedly conf irmed, for example, by independent observations during deer domestication, during the use of rats as laboratory animals, after the reintroduction of endangered species such as Przewalski’s horse, and during the creation of a Siberian reserve population of the Siberian grouse when it had reached an endangered status in natural habitats. A genome-wide comparison among humans, several domestic animals, and some of their wild congeners has given rise to the concept of self-domestication syndrome, which includes autism spectrum disorders. In our previous study, we created a bioinformatic model of human self-domestication syndrome using differentially expressed genes (DEGs; of domestic animals versus their wild congeners) orthologous to the human genes (mainly, nervous-system genes) whose changes in expression affect reproductive potential, i. e., growth of the number of humans in the absence of restrictions caused by limiting factors. Here, we applied this model to 68 human genes whose changes in expression alter the reproductive health of women and men and to 3080 DEGs of domestic versus wild animals. As a result, in domestic animals, we identif ied 16 and 4 DEGs, the expression changes of which are codirected with changes in the expression of the human orthologous genes decreasing and increasing human reproductive potential, respectively. The wild animals had 9 and 11 such DEGs, respectively. This difference between domestic and wild animals was signif icant according to Pearson’s χ2 test (p < 0.05) and Fisher’s exact test (p < 0.05). We discuss the results from the standpoint of restoration of endangered animal species whose natural habitats are subject to an anthropogenic impact.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

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 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).

[Candidate SNP-markers altering TBP binding affinity for promoters of the Y-linked genes CDY2A, SHOX, and ZFY are lowering many indexes of reproductive potential in men]

Репродуктивный потенциал – уровень физического и психического состояния организма, позволяю- щий при достижении социальной зрелости воспроизводить здоровое потомство. В узком биомедицинском смыс- ле определение включает комплекс функциональных показателей репродуктивной системы, но в более широком смысле его можно рассматривать как совокупность физиологических, поведенческих, адаптивных, ментальных, ан- тропометрических и генетических характеристик особи, способствующих размножению. Целью настоящей работы было расширить область применимости созданного ранее Web-сервиса SNP_TATA_Z-tester для поиска кандидатных маркеров однонуклеотидного полиморфизма (SNP) на хромосоме Y человека, связанных с мужским репродуктив- ным потенциалом (МРП). В поиске кандидатных SNP-маркеров для МРП мы сосредоточились на генах хромосомы Y человека. Изучены 35 SNP в промоторах генов CDY2A, SHOX и ZFY, представляющих все три типа генов хромосомы Y человека: уникальный, псевдоаутосомный и паралог гена хромосомы X человека соответственно. Предсказаны 11 кандидатных SNP-маркеров ослабления МРП из-за изменения сродства TATA-связывающего белка (TBP) к этим промоторам. Выборочно верифицированы in vitro величины сродства «TBP-промотор», предсказанные в этой рабо- те. Установлена достоверная корреляция (r = 0.94, p < 0.005) между ними и результатами измерения in vitro сродства ТВР человека к олигонуклеотидам, идентичным сайтам ТВР-связывания исследуемых промоторов. Проведя поиск в базе данных PubMed по ключевым словам, мы нашли клиническое описание патологических состояний человека, соответствующих изменению экспрессии генов, несущих предсказанные нами кандидатные SNP-маркеры. Среди них оказались такие патологии, как нарушение сперматогенеза (ZFY: rs1388535808 и rs996955491), задержка поло- вого созревания (CDY2A: rs200670724), нарушения эмбриогенеза (SHOX: rs28378830) и непропорционально низкий рост с деформациями Маделунга (SHOX: rs1452787381). Они свидетельствуют, что в случае SNP-промоторов генов хромосомы Y человека следует ожидать изменений широкого круга показателей МРП, выходящих далеко за рамки генетического контроля собственно мужской репродуктивной функции.

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.

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.

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.

Сandidate SNP-markers of rheumatoid arthritis that can significantly alter the affinity of the TATA-binding protein for human gene promoters

Rheumatoid polyarthritis (RA) is an autoimmune disease with autoantibodies, including antibodies to citrullant antigens and proinflammatory cytokines, such as TNF-α and IL-6, which are involved in the induction of chronic synovitis, bone erosion, followed by deformity. Immunopathogenesis is based on the mechanisms of the breakdown of immune tolerance to its own antigens, which is characterized by an increase in the activity of T-effector cells, causing RA symptomatology. At the same time, against the background of such increased activity of effector lymphocytes, a decrease in the activity of a number of regulatory cells, including regulatory T-cells (Treg) and myeloid suppressor cells, is recorded. There is reason to say that it is the change in the activity of suppressor cells that is the leading element in RA pathogenesis. That is why only periods of weakening (remission) of RA are spoken of. According to the more powerful female immune system compared to the male one, the risk of developing RA in women is thrice as high, this risk decreases during breastfeeding and grows during pregnancy as well as after menopause in proportion to the level of sex hormones. It is believed that 50 % of the risk of developing RA depends on the conditions and lifestyle, while the remaining 50 % is dependent on genetic predisposition. That is why, RA fits the main idea of postgenomic predictive-preventive personalized medicine that is to give a chance to those who would like to reduce his/her risk of diseases by bringing his/her conditions and lifestyle in line with the data on his/her genome sequenced. This is very important, since doctors consider RA as one of the most frequent causes of disability. Using the Web service SNP_TATA_Z-tester (http://beehive.bionet.nsc.ru/cgi-bin/mgs/tatascan_fox/start.pl), 227 variants of single nucleotide polymorphism (SNP) of the human gene promoters were studied. As a result, 43 candidate SNP markers for RA that can alter the affinity of the TATA-binding protein (TBP) for the promoters of these genes were predicted.

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).