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.

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.

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.

CisCross: A gene list enrichment analysis to predict upstream regulators in Arabidopsis thaliana.

Having DNA-binding profiles for a sufficient number of genome-encoded transcription factors (TFs) opens up the perspectives for systematic evaluation of the upstream regulators for the gene lists. Plant Cistrome database, a large collection of TF binding profiles detected using the DAP-seq method, made it possible for Arabidopsis. Here we re-processed raw DAP-seq data with MACS2, the most popular peak caller that leads among other ones according to quality metrics. In the benchmarking study, we confirmed that the improved collection of TF binding profiles supported a more precise gene list enrichment procedure, and resulted in a more relevant ranking of potential upstream regulators. Moreover, we consistently recovered the TF binding profiles that were missing in the previous collection of DAP-seq peak sets. We developed the CisCross web service (https://plamorph.sysbio.ru/ciscross/) that gives more flexibility in the analysis of potential upstream TF regulators for Arabidopsis thaliana genes.

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.

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 Catalog of Human Genes Associated With Pathozoospermia and Functional Characteristics of These Genes.

Genetic causes of the global decline in male fertility are among the hot spots of scientific research in reproductive genetics. The most common way to evaluate male fertility in clinical trials is to determine semen quality. Lower semen quality is very often accompanied by subfertility or infertility, occurs in many diseases and can be caused by many factors, including genetic ones. The following forms of lowered semen quality (pathozoospermia) are known: azoospermia, oligozoospermia, asthenozoospermia, teratozoospermia, and some combined forms. To systematize information about the genetic basis of impaired spermatogenesis, we created a catalog of human genes associated with lowered semen quality (HGAPat) and analyzed their functional characteristics. The catalog comprises data on 126 human genes. Each entry of the catalog describes an association between an allelic variant of the gene and a particular form of lowered semen quality, extracted from the experimental study. Most genes included into the catalog are located on autosomes and are associated with such pathologies as non-obstructive azoospermia, oligozoospermia or asthenozoospermia. Slightly less than half of the included genes (43%) are expressed in the testes in a tissue-specific manner. Functional annotation of genes from the catalog showed that spermatogenic failure can be associated with mutations in genes that control biological processes essential for spermiogenesis (regulating DNA metabolism, cell division, formation of cellular structures, which provide cell movement) as well as with mutations in genes that control cellular responses to unfavorable conditions (stress factors, including oxidative stress and exposure to toxins).

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

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.

Origin and Evolution of the Neo-Sex Chromosomes in Pamphagidae Grasshoppers through Chromosome Fusion and Following Heteromorphization.

In most phylogenetic lineages, the evolution of sex chromosomes is accompanied by their heteromorphization and degradation of one of them. The neo-sex chromosomes are useful model for studying early stages of these processes. Recently two lineages of the neo-sex chromosomes on different stages of heteromorphization was discovered in Pamphagidae family. The neo-sex chromosome heteromorphization was analyzed by generation of DNA probes derived from the neo-Xs and neo-Ys followed with chromosome painting in nineteen species of Pamphagidae family. The homologous regions of the neo-sex chromosomes were determined in closely related species with the painting procedure and image analysis with application of the Visualization of the Specific Signal in Silico software package. Results of these analyses and distribution of C-positive regions in the neo-sex chromosomes revealed details of the heteromorphization of the neo-sex chromosomes in species from both phylogenetic lineages of Pamphagidae grasshoppers. The hypothetical mechanism of the neo-Y degradation was suggested. It includes expansion of different repeats from the proximal neo-Y chromosome region by inversions, spreading them towards distal region. Amplification of these repeats leads to formation of C-positive regions and elimination of the C-negative regions located between them.

Spatial organization of fibroblast and spermatocyte nuclei with different B-chromosome content in Korean field mouse, Apodemus peninsulae (Rodentia, Muridae).

Korean field mouse (Apodemus peninsulae) shows a wide variation in the number of B chromosomes composed of constitutive heterochromatin. For this reason, it provides a good model to study the influence of the number of centromeres and amount of heterochromatin on spatial organization of interphase nuclei. We analyzed the three-dimensional organization of fibroblast and spermatocyte nuclei of the field mice carrying a different number of B chromosomes using laser scanning microscopy and 3D fluorescence in situ hybridization. We detected a co-localization of the B chromosomes with constitutive heterochromatin of the chromosomes of the basic set. We showed a non-random distribution of B chromosomes in the spermatocyte nuclei. Unpaired B chromosomes showed a tendency to occur in the compartment formed by the unpaired part of the XY bivalent.