Next Article in Journal
Pro-Inflammatory and Neurotrophic Factor Responses of Cells Derived from Degenerative Human Intervertebral Discs to the Opportunistic Pathogen Cutibacterium acnes
Previous Article in Journal
Correction: Yan, L., et al. Combined Taurine, Epigallocatechin Gallate and Genistein Therapy Reduces HSC-T6 Cell Proliferation and Modulates the Expression of Fibrogenic Factors. Int. J. Mol. Sci. 2013, 14, 20543–20554
Previous Article in Special Issue
Transcriptomic Changes Related to Cellular Processes with Particular Emphasis on Cell Activation in Lysosomal Storage Diseases from the Group of Mucopolysaccharidoses
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 5
10.3390/ijms22052346
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Gennady Vasiliev
1,2,
Irina Chadaeva
2,
Dmitry Rasskazov
2,
Petr Ponomarenko
2,
Ekaterina Sharypova
2,
Irina Drachkova
2,
Anton Bogomolov
2,
Ludmila Savinkova
2,
Mikhail Ponomarenko
2,*,
Nikolay Kolchanov
2,
Alexander Osadchuk
2,
Dmitry Oshchepkov
2 and
Ludmila Osadchuk
2
1
Novosibirsk State University, 630090 Novosibirsk, Russia
2
Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(5), 2346; https://doi.org/10.3390/ijms22052346
Submission received: 21 November 2020 / Revised: 13 February 2021 / Accepted: 23 February 2021 / Published: 26 February 2021
(This article belongs to the Special Issue Utilizing RNA-Seq and Genome Sequencing to Uncover Complexities of Genetic Disorders)
Browse Figure
Versions Notes

Abstract

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

1. Introduction

As Theodore Dobzhansky noted, “man is genetically specialized to be unspecialized” [1]. People have settled down, indeed, in all landscapes and all climates of all continents [2,3]. To this end, humans did not adapt anthropometrically to the natural environment like animals do but rather adapted it (as an artificial anthropogenic environment) to their own lives, including the domestication of other organisms to satisfy their own needs and even whims [4,5]. Of note, many behavioral and anatomical differences between domestic and wild animals can correspond to human disease symptoms. For example, the diminished weakened heart of domestic ducks versus the large powerful heart of wild ducks reflects settledness versus spring-and-autumn migrations [6], and this relation may correspond to myocardial infarction as the most common cause of death in humans, as rated by the World Health Organization (WHO) [7]. Currently, translational medicine researchers [8] already use domestic animals, including transgenic ones [9], in preclinical trials of drugs intended for human treatments. Thus, bioinformatics models of human diseases on the basis of differentially expressed genes (DEGs) of domestic versus wild animals (from animal genome-wide transcriptomes determined in vivo) are relevant and useful in the post-genomic era of life sciences.
According to Ernst Mayer’s concept [10], the cornerstone for the existence of a biological species is a set of populations related to one another due to reproduction, as follows: each previous generation (parents) must leave behind the next generation (offspring). In 1987, the WHO classified human reproduction biomedical risk factors (in decreasing order): female infertility (38%), seeming infertility in biologically parental couples (27%), male infertility (20%), and everything else (15%) [11]; this ranking has remained true to date [12]. That is why the Centers for Disease Control and Prevention (USA) are uniformly systematizing and coordinating world-wide studies on female infertility, whereas male infertility studies are still scattered across diverse local databases that collect anecdotal information only on those individual patients who are anxious about their own reproductive disorders [12].
At the same time, phylogenetic reconstruction of the coevolution of chromosomes X and Y of a wide range of species (opossums, cattle, rats, mice, marmosets, rhesus monkeys, chimpanzees, and humans) on the basis of data from whole-genome high-throughput sequencing [13] indicates that in addition to reproductive function, these chromosomes determine gender dimorphism in both predisposition to diseases and in organism viability as a whole. With this in mind, at either a populational or individual scale, reproductive potential is the most common conventional indicator of the probability to survive, produce offspring, and bring it to a reproductive state under the best conditions [14,15]. As for human biomedicine, in the narrowest sense, reproductive potential means only a set of infertility risk factors (e.g., [16]), whereas its broadest interpretation includes a wide variety of physiological, mental, behavioral, social, anthropometric, and genetic indicators contributing to humankind’s reproduction as a whole [17].
According to Konrad Lorenz’s ideas [18], if there is a need to protect one’s own life, offspring, home, territory, or access to limited resources, then reproductive behavior can also be accompanied by aggressiveness. In social animals (including humans), the aggressiveness of an individual determines his/her socio-hierarchical rank, which affects the quality of his/her life and his/her lifespan [19]. Moreover, in many human diseases, aggressiveness can be a risk factor (e.g., [20]), a symptom (e.g., [21]), or a complication of treatment (e.g., [22]). On a whole-genome scale, many animals are already subjected to aggressiveness studies, including the trout [23], chickens [24], foxes [25,26], dogs [27,28], pigs [29,30], rats [29,31], mice [32], rabbits, and guinea pigs [29]; however, this research still does not take into account the reference human genome [33], human variome [34], or the classification of diseases and health problems ICD-11 [35].
Nevertheless, now among the most challenging problems is how to restore many critically endangered animal species within the unclear natural limitations of population restoration by captive breeding [36]. One possible step towards solving this urgent problem could be 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, in order to determine how the anthropogenic environment, compared to natural wildlife, can alter the animal gene expression on the genome-wide scale and thus affect reproductive potential.
In our previous works, we performed genome-wide analyses of single-nucleotide polymorphisms (SNPs) of TATA-binding protein (TBP)-binding sites in the promoters of genes associated with human reproductive health [37,38,39,40,41]. As a result of these bioinformatic estimates, we proposed that human reproductive potential diminishes during self-domestication [42]. This is why, here, we, for the first time, proposed to utilize domestic versus wild animal transcriptomes determined in vivo as an animal bioinformatics model of human diseases, and, then, within this model framework, verified our above-mentioned hypothesis in silico on experimental datasets of the DEGs of domestic versus wild animals. Finally, we discussed the results obtained with respect to how it could be possible to domesticate new economically valuable animals without reducing their reproductive potential as well as why the return from anthropogenic habitat conditions to natural ones (in the wild) is actually contributing to the restoration of endangered animal species.

2. Results and Discussion

In this work, we verified our whole-genome sequence-based a priori hypothesis (derived in silico [41]) that human reproductive potential can diminish during self-domestication [42], as depicted in Figure 1 and described in the Section 3.

2.1. The Bioinformatics Model Developed

As the bioinformatics model development, first of all, we compiled the set of all the 275 human genes whose effects on human reproductive potential have been a priori estimated by means of SNPs within their 70 bp proximal promoters in our previous articles on this matter [37,38,39,40,41] and, next, updated literary sources in line with the current state of PubMed [43], as shown in Figure 1 (Step-1) and presented in Table S1 (hereinafter, see Supplementary Materials). Then, by means of PubMed [43], we compiled independent publicly available experimental in vivo RNA-Seq data sets on domestic versus wild animals [26,28,29], as depicted in Figure 1 (Step-2) and presented in Table 1. For the RNA-Seq data analysis, here, we, for the first time, resumed the use of one of the eldest [44] commonly accepted (e.g., see [45]) and widely used (e.g., see [46]) concepts of “divergence from the nearest common ancestor (NCA)” within the framework of the microevolution theory in relation to domesticated and wild animals, whose differentially expressed genes (DEGs) were compared with orthologous human genes under this study, as depicted in Figure 1. After that, we found 28 pairs of orthologous genes in humans and animals among all genes analyzed in this work, as depicted by a Venn diagram in Figure 1 (Step-3). The largest part (13) of the 28 DEGs (46%) characterizes the guinea pig (Cavia porcellus) in comparison with the cavy (C. aperea) [29]. Additionally, the remaining 15 DEGs (54%) characterize domestic rabbits (Oryctolagus cuniculus) [29], dogs (Canis familiaris) [28,29], and foxes (Vulpes vulpes) [26] versus their wild congeners (including the wolf, C. lupus). Next, we compared the effects of co-directional changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their NCA, as shown in Figure 1 (Step-4). On this basis, we formatted Fisher’s binary correlation tables 2 × 2, depicted in Figure 1 (Step-5). Finally, we used the standard software STATISTICA (StatsoftTM, Tulsa, OK, USA), the pathway “Statistics” → “Nonparametric” → “table 2 × 2” of which was provided by the Fisher’s exact test, Pearson χ2 test and binomial distribution test to verify our findings (Figure 1: Step-6).
Each of these DEGs is discussed in the following.

2.2. DEGs of the Guinea Pig Versus Cavy and How Their Human Orthologous Genes Change Reproductive Potential

Human gene CETP encodes cholesteryl ester transfer protein, whose deficit lowers the risk of myocardial infarction [47] and improves human reproductive potential consistently with cavy gene Cetp downregulation during this wild animal’s microevolution from the above NCA [29]. Conversely, CETP overexpression is a risk factor of hypercholesterolemia during late pregnancy [48] and reduces human reproductive potential; the direction of expression change is consistent with the guinea pig gene Cetp overexpression in microevolution [29], as shown in Table 2.
Human genes CHRNA3 and CHRNA6 correspond to cholinergic receptor nicotinic subunits α3 and α6, deficiencies of which elevate human reproductive potential owing to improvements in finding opposite sex congeners [49] and in maternal behavior [51], respectively, in agreement with expression changes of cavy orthologous genes during microevolution [29]. Conversely, CHRNA3 overexpression and CHRNA6 overexpression worsen human reproductive potential through worse effects of nicotine compounds on primordial oocytes [50] and via a higher risk of social defeats [52], in agreement with the direction of expression change of the guinea pig orthologous genes in microevolution [29] (Table 2).
Human gene DRD5 codes for dopamine receptor D5, whose deficit reduces conditioned fear [53] and improves reproductive potential (Table 2) co-directedly with guinea pig Drd5 underexpression during microevolution [29]. As an alternative, DRD5 upregulation raises risks of mental disorders [54] and worsens human reproductive potential, co-directedly with the cavy gene Cyp17a1 decrease in expression during microevolution [29].
Human gene FLT4 produces feline McDonough sarcoma-like tyrosine kinase 4, a deficit of which suppresses melanoma metastasis [55] and enhances human reproductive potential consistently with the cavy gene Flt4 expression change as compared with the NCA of the guinea pig and cavy [29]. On the other hand, FLT4 excess exacerbates post-traumatic inflammatory neovascularization in humans [56] and worsens human reproductive potential co-directedly with the guinea pig orthologous gene during microevolution [29] (Table 2).
Human genes GFRA3 and GFRA4 correspond to GDNF family receptors α3 and α4, low expression of which worsens reproductive potential via accelerated neurodegeneration [57] and premature adolescent bone formation [59], respectively, in agreement with the direction of expression change of the corresponding guinea pig and cavy orthologous genes during their microevolution (Table 2). Conversely, both GFRA3 and GFRA4 excesses can increase human reproductive potential because of improvements in neural regeneration [58] and neuronal survival [60], respectively.
Human gene HTR3A encodes serotonin receptor 3A, whose deficit elevates the risk of sudden cardiac death in pregnant women [61] consistently with the guinea pig orthologous gene in microevolution [29] (Table 2), whereas HTR3A overexpression improves mood and behavior [62].
Human gene IL1B encodes interleukin 1β, the decreased expression of which prevents bone deformations in infections [63] and improves reproductive potential in agreement with cavy Il1b downregulation in microevolution [29]. Contrariwise, IL1B overexpression enhances circadian pain hypersensitivity [64], worsening reproductive potential in line with Il1b excess in guinea pig microevolution [29], as shown in Table 2.
Human gene NR5A1 encodes steroidogenic factor 1, whose underexpression and overexpression correspond to impaired [65] and improved [66] gonadal development, respectively, as well as to expression changes of the guinea pig and cavy orthologous genes during their divergence from their NCA [29] (Table 2).
Human gene PDGFRL (platelet-derived growth factor receptor–like protein) underexpression reduces tumor mutation burden [67], thus raising reproductive potential consistently with the cavy ortholog expression alteration during microevolution [29], whereas a PDGFRL excess leads to hypertensive behavior and myocardial hypertrophy [68] (Table 2).
Human gene PDYN encodes prodynorphin, underexpression of which reduces reproductive potential through obesity-related subfertility [69] in agreement with the direction of expression change of the cavy ortholog downregulation in the course of microevolution [29], but a PDYN excess prevents conditioned fear behavior [70] (Table 2).
Human gene SLC6A4 produces serotonin transporter 1, whose deficiency improves small-intestine function [71], thus elevating reproductive potential consistently with the cavy orthologous gene expression change during microevolution [29], whereas SLC6A4 overexpression worsens depression, anxiety, and inertia [72] (Table 2).
Table 3 summarizes the results of the comparison of the above orthologous genes of humans and guinea pigs, namely, 11 and two of these guinea pig DEGs were found to correspond to human gene-markers of worsened and improved reproductive potential, and the same is true for three and ten DEGs in the cavy. This means that the DEGs in guinea pigs significantly correspond to the human orthologous genes reducing human reproductive potential according to three independent tests, namely, Pearson’s χ2 test (p < 0.01), Fisher’s exact test (p < 0.05), and binomial distribution (p < 0.05), in contrast to cavy DEGs, which almost equally fit human orthologous genes worsening and improving this human trait (p > 0.05, binomial distribution), thus enforcing the wild-type norm.
To verify that this phenomenon is not species-specific among animals, we considered the DEGs of domestic rabbits, dogs, and foxes as compared with their wild congeners in relation to human gene-markers of reproductive potential changes, as revealed in this work (Table 4).

2.3. DEGs of Domestic Versus Wild Animals and How the Human Orthologous Genes Alter Reproductive Potential

Human gene F7 encodes coagulation factor VII, and its underexpression correlates with life-threatening bleeding [73], thus reducing reproductive potential in agreement with the direction of expression change of domestic rabbit F7 relative to the NCA of domestic and wild rabbits [29], whereas recombinant F7 is a drug saving life and fertility during intractable obstetric bleeding in women [74] (Table 4).
Human gene PDGFRA encodes platelet-derived growth factor receptor α, both a deficit and excess of which worsen reproductive potential through skeletal defects in newborns [75] and predisposition to infertility after infections [76] (Table 4).
Human gene GABARAPL2 produces GABA type A receptor-associated protein-like 2, whose downregulation retards wound healing [77] and reduces reproductive potential consistently with the dog orthologous gene’s expression change during both dog and wolf divergence from their NCA [28]. Conversely, a GABARAPL2 excess improves tooth injury healing [78], as shown in Table 4.
Human gene GH1 encodes growth hormone 1, whose underexpression raises the risks of morbidity and mortality from cardiovascular diseases [79] and reduces reproductive potential, in agreement with the wolf Gh1 underexpression in microevolution [28] (Table 4). On the contrary, recombinant GH1 is used as a drug to prolong the reproductive age in women [80].
Human gene HBB codes for hemoglobin subunit β, whose deficit (thalassemia) worsens women’s reproductive health [81], in agreement with the direction of expression change of the dog orthologous gene—deficit during microevolution [28] (Table 4). Conversely, in traditional Chinese medicine, the Jian-Pi-Yi-Sheng decoction (JPYS) is employed to raise the HBB level for treating anemia in chronic kidney diseases [82]
Human gene NRP2 produces neuropilin 2, the downregulation of which improves survival after radiochemotherapy [83], thereby raising reproductive potential consistently with wolf orthologous gene downregulation in microevolution [28] (Table 4). Conversely, an NRP2 excess causes post-traumatic vascular neointimal hyperplasia [84] (Table 4).
Human gene TAC3 codes for tachykinin precursor 3, whose deficit and excess correspond to high and low risks of socially induced subfertility [85] as well as to expression changes of the Tac3 gene during the divergence of dogs and wolfs from their NCA [28] (Table 4).
Human gene TGFB3 encodes transforming growth factor β3, whose deficit and excess worsen reproductive potential in men [86] and women [87], respectively (Table 4).
Human gene ESR2 corresponds to estrogen receptor 2; both deficiency and overabundance in adolescents worsen spermatogenesis in adult males [88] (Table 4).
Human gene GRIN3A codes for glutamate ionotropic receptor NMDA type subunit 3A whose underexpression prevents cocaine addiction [89], thus improving reproductive potential consistently with the wild fox Grin3a deficiency in microevolution [26] (Table 4). Oppositely, a GRIN3A excess increases the risk of inattentive behavior [90] and, therefore, reduces reproductive potential that fits a Grin3a excess during the tame fox microevolution [26], as shown in Table 4.
Human gene HTR3B encodes serotonin receptor 3B, the downregulation of which reduces anger-resolutive behavior [91], thereby reducing reproductive potential in agreement with the expression alteration of the tame fox orthologous gene during microevolution [26] (Table 4), whereas HTR3B overexpression decreases the risk of pulmonary embolism [92].
Human gene IL6ST codes for interleukin 6 signal transducer; both deficit and excess worsen reproductive potential through increased risk of mortality during sepsis [93] and sensitivity to fatigue [94], respectively, as shown in Table 4.
Human gene IL9R encrypts interleukin 9 receptor; both downregulation and overabundance impair reproductive potential via impaired trophoblast implantation in preeclampsia [95] and increased risk of life-threatening anaphylaxis [96].
Human gene NPY codes for neuropeptide Y, whose deficit [97] and excess [98] cause subfertility, as presented in Table 4.
Human gene TGFB2 produces transforming growth factor β2; its downregulation and upregulation reduce reproductive potential through increased risk of perinatal mortality [99] and impaired wound healing [100], respectively (Table 4).
Table 5 sums up the findings of the comparative analysis of the above orthologous genes from humans, rabbits, dogs, wolfs, and foxes. For example, 14 and one of these domestic animal DEGs were found to correspond to human gene-markers of reduced and elevated reproductive potential, and the same is true for eight and seven DEGs in the wild animals. Accordingly, once again, we observed that the DEGs in domestic animals reliably correspond to their human orthologous genes impairing human reproductive potential, according to Pearson’s χ2 test (p < 0.05), Fisher’s exact test (p < 0.05), and binomial distribution (p < 0.001). On the other hand, DEGs of wild animals correlate equally to human orthologous genes, which weaken and enhance human reproductive potential (p > 0.5, binomial distribution), which corresponds to the wild-type norm, as indicated in Table 5.

2.4. DEGs in Domestic Animals Reliably Correspond to Their Human Orthologs Reducing Reproductive Potential

Generalizing Table 3 and Table 5, we found that 25 and three DEGs in domestic animals as compared with 11 and 17 DEGs in wild animals correspond to the human orthologous gene-markers of reduced and increased reproductive potential (Table 6). Therefore, DEGs in domestic animals reliably correspond to their human orthologous genes that diminish reproductive potential, judging by Pearson’s χ2 test (p < 0.001), Fisher’s exact test (p < 0.05), and binomial distribution (p < 0.0001). By contrast, DEGs of wild animals are equally fitting human orthologous genes decreasing and increasing reproductive potential (in terms of binomial distribution, p > 0.1), thereby possibly enforcing the wild-type norm.
This is why it is important to underscore that our bioinformatic hypothesis examined here was confirmed in vivo in animal models of human diseases by means of independent experimental RNA-Seq data from domestic and wild animals and was formulated in accordance with the concept of human reproductive potential as the most comprehensive indicator of chances for successful survival, production of offspring, and bringing these offspring to reproductive age [10,14,15]. As for the actual realization of these chances in practice, we found literary evidence both in favor and against one-to-one correspondence between domestication and reduced reproductive potential in animals. First of all, the observed correspondence between the dog-versus-wolf DEGs and clinically proven markers of reduced human reproductive potential (Table 4) is consistent with the correspondence between some deviations characteristic of autism spectrum disorders in humans and both physiological and behavioral differences of dogs from wolves [42]. Moreover, tame foxes (Table 4) as compared to wild ones have a worse female endocrine system [101], reduced sexual activity in first-year males [102,103], accelerated extinction of testicle hormonal function [104], impaired reproductive seasonality [105], and reduced gonad mass together with heterochrony of pituitary–spermatic complex development in the embryo [106]. Additionally, compared to aggressive male rats, tame ones show delayed puberty [107]. Additionally, the results of this work are supported by the finding that a return of Saiga antelope (Saiga tatarica) [108] and Amur tigers (Panthera tigris altaica) [109] from anthropogenic habitat conditions to natural ones (in the wild) has already successfully helped to restore these endangered species.
Curiously, as a counterargument to these five observations, there is a good example of domestic pigs, which outperform wild boars on both sperm quality and spermatogenesis [110]. This phenomenon may be due to artificial selection aimed at improving their fertility for greater meat production. This means that if during the domestication of new economically valuable animal species (e.g., the musk deer Moschus berezovskii [111]), artificial selection for improved target traits is supplemented with selection for higher fertility, then this approach may compensate the risks of the domestication process, e.g., reduced reproductive potential, as reported by many authors elsewhere [42,101,102,103,104,105,106,107].

3. Materials and Methods

3.1. Human Genes under Study

Here, we studied 275 human genes, which are described in Table S1 (see Supplementary Materials) according to the results of our in silico analysis of the effect of the SNPs (located in proximal promoters) on human reproductive health [37,38,39,40,41]; the literature supporting these data was updated according to the current state of the PubMed database [43], as depicted in Figure 1 (Step-1).

3.2. DEGs of Domestic Animals Compared to Their Wild Congeners

In this work, we used publicly available independent experimental RNA-Seq datasets on transcriptomes of domestic versus wild animals [26,28,29]. Although, here, we compare the DEG of each animal with one human orthologous gene, we nevertheless limited our analysis to only those DEGs that were statistically significant according to Fisher’s Z-test, with those corrections on multiple comparisons (PADJ < 0.05), which are publicly available within the RNA-Seq data studied here, as published by their authors [26,28,29]. As a result, a total of 1740 DEGs were analyzed, namely, within the frontal cortex, there were 883 DEGs of guinea pigs versus cavies, 30 DEGs of pigs versus boars, 20 DEGs of tame versus aggressive rats, 17 DEGs of domesticated versus wild rabbits, and 13 DEGs of dogs versus wolves [29], as well as 450 DEGs in the blood of dogs versus wolves [28] and 327 DEGs in the pituitary gland of tame versus aggressive foxes [26] (Table 1 and Figure 1: Step-2).
The only novelty of this work is that for the RNA-Seq data analysis, we, for the first time, resumed the use of one of the most ancient [44] conventional (e.g., see [45]) and actual (e.g., see [46]) microevolutional concepts of “divergence from the nearest common ancestor (NCA)” in the case of domestic and wild animals, whose differentially expressed genes (DEGs) were compared with orthologous human genes under this study (Figure 1: Step-3). This allows us to compare the effects of co-directional changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their NCA, as shown in Table 2 and Table 4, as well as in Figure 1 (Step-4).

3.3. Statistical Analysis

Using Table 2 and Table 4, we made the standard statistical tables 2 × 2 (Figure 1: Step-5), which are the input data for the standard software STATISTICA (StatsoftTM, Tulsa, OK, USA), where its pathway “Statistics” → “Nonparametric” → “table 2 × 2” led to the Fisher’s exact test, Pearson χ2 test, and binomial distribution test to verify the significance of our results (Figure 1: Step-6).

3.4. The Knowledge Base on Domestic Animals’ DEGs Whose Human Orthologous Genes Can Change Reproductive Potential

We Excel-compatibly formatted all the identified associations between DEGs of domestic versus wild animals and the effects of the human orthologous genes on reproductive potential as a textual flat file and finally converted it into the PetDEGsDB knowledge base format in Web environment MariaDB 10.2.12 (MariaDB Corp AB, Espoo, Finland). The PetDEGsDB knowledge base made by this work is publically available on the web page www.sysbio.ru/domestic-wild (accessed on 26 February 2021).

4. Conclusions

In this work, we, for the first time, proposed to utilize differences between domestic and wild animal transcriptomes as a bioinformatics model of human diseases. Within this model framework, we successfully confirmed (Table 2 and Table 4) our a priori hypothesis (derived in silico) that human reproductive potential can diminish during self-domestication [42]. This hypothesis was formulated on the basis of our genome-wide sequence-based analysis of SNPs within 70 bp proximal promoters of human Y-linked genes [41]. As we presented in Table 2, Table 3, Table 4, Table 5 and Table 6, DEGs of domestic versus wild animals of the same species indicate genetic differences among all intergroup differences that are statistically significant in comparison with the variability within the species in question. (Some variability is necessary for the existence of a species.) This property makes these DEGs a promising tool for microevolution research (Table 2 and Table 4).
Finally, in this work, we fully analyzed the genes of human neurotransmitter and neurotrophinergic systems as well as genes on the human Y chromosome. Therefore, in the future, an expansion of the list of such genes to, for example, the human endocrine and immune systems may increase the completeness of the reproductive potential analysis. With this in mind, in the future, a similar extended bioinformatic analysis of RNA-Seq data on the multifactorial human diseases will be interesting, because it will estimate the effects of predisposition and resistance to such diseases on human reproductive potential.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/22/5/2346/s1.

Author Contributions

Conceptualization and supervision, N.K.; methodology, A.O., L.O. and A.B.; software, A.B.; validation, I.C., E.S., I.D. and L.S.; investigation, D.R. and P.P.; resources, D.O.; data curation, G.V.; writing—original draft preparation, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

Manuscript writing was supported by Russian government project #0259-2021-0009 (for M.P.). The data curation was supported by project #28.12487.2018/12.1 from the Ministry of Science and Higher Education of the Russian Federation (for G.V.). The conceptualization and supervision of the investigation were supported by the Russian Federal Science & Technology Program for the Development of Genetic Technologies (for N.K., D.R., P.P., I.C., E.S., I.D., L.S. and D.O.). The methodology and software development were supported by project #19-15-00075 from the Russian Science Foundation (for A.O., L.O. and A.B.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We are grateful to Nikolai Shevchuk (Shevchuk Editing Co., Brooklyn, NY, USA) for useful discussions of the work and assistance in adapting it to the requirements of the English language for scientific publications.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DEGdifferentially expressed gene
NCAnearest common ancestor
SNPsingle-nucleotide polymorphism
WHOWorld Health Organization

References

  1. Dobzhansky, T. Genetic entities in hominid evolution. In Classification and Human Evolution; Washburn, S.L., Ed.; Aldine Publishing Co.: Chicago, IL, USA, 1963; pp. 347–362. [Google Scholar]
  2. Darwin, C.R. A Naturalist’s Voyage Round the World: The Voyage of the Beagle; John Murray Publishing: London, UK, 1913. [Google Scholar]
  3. Henry, J.P. Genetics and origin of Homo sapiens. Med. Sci. 2019, 35, 39–45. [Google Scholar]
  4. Darwin, C.R. On the Origin of Species by Means of Natural Selection, or the Preservation of Favoured Races in the Struggle for Life; John Murray Publishing: London, UK, 1859. [Google Scholar]
  5. Hammond, T.T.; Ortiz-Jimenez, C.A.; Smith, J.E. Anthropogenic change alters ecological relationships via interactive changes in stress physiology and behavior within and among organisms. Integr. Comp. Biol. 2020, 60, 57–69. [Google Scholar] [CrossRef] [PubMed]
  6. Charuta, A.; Mankowska-Pliszka, H.; Bartyzel, B.J.; Wysocki, J. Size of heart of the domestic Pekin duck (Anas platyrhynchos f. domestica) and wild duck (Anas platyrhynchos, Linnaeus, 1758). Acta Sci. Pol. Med. Vet. 2005, 4, 11–19. [Google Scholar]
  7. Barquera, S.; Pedroza-Tobias, A.; Medina, C.; Hernandez-Barrera, L.; Bibbins-Domingo, K.; Lozano, R.; Moran, A.E. Global overview of the epidemiology of atherosclerotic cardiovascular disease. Arch. Med. Res. 2015, 46, 328–338. [Google Scholar] [CrossRef]
  8. McGonigle, P.; Ruggeri, B. Animal models of human disease: Challenges in enabling translation. Biochem. Pharmacol. 2014, 87, 1621–1671. [Google Scholar] [CrossRef]
  9. Lutz, C.; Maher, L.; Lee, C.; Kang, W. COVID-19 preclinical models: Human angiotensin-converting enzyme 2 transgenic mice. Hum. Genom. 2020, 14, 20. [Google Scholar] [CrossRef] [PubMed]
  10. Mayr, E. Systematics and the Origin of Species, from the Viewpoint of a Zoologist; Harvard University Press: Cambridge, MA, USA, 1942. [Google Scholar]
  11. World Health Organization. Towards more objectivity in diagnosis and management of male infertility: Results of a World Health Organization multicentre study. Int. J. Androl. 1987, 7, 1–53. [Google Scholar]
  12. Hanson, B.M.; Hotaling, J.M. Limitations and opportunities in male fertility databases. Transl. Urol. 2018, 7, S292–S302. [Google Scholar] [CrossRef] [PubMed]
  13. Bellott, D.W.; Hughes, J.F.; Skaletsky, H.; Brown, L.G.; Pyntikova, T.; Cho, T.-J.; Koutseva, N.; Zaghlul, S.; Graves, T.; Rock, S.; et al. Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators. Nature 2014, 508, 494–499. [Google Scholar] [CrossRef]
  14. Chapman, R.N. Animal Ecology with Special Reference to Insects; McGraw-Hill Book Co.: New York, NY, USA, 1931. [Google Scholar]
  15. Pianka, E.R. Natural selection of optimal reproductive tactics. Am. Zool. 1976, 16, 775–784. [Google Scholar] [CrossRef]
  16. du Fosse, N.A.; van der Hoorn, M.P.; van Lith, J.M.M.; le Cessie, S.; Lashley, E.E.L.O. Advanced paternal age is associated with an increased risk of spontaneous miscarriage: A systematic review and meta-analysis. Hum. Reprod. Update 2020, 26, 650–669. [Google Scholar] [CrossRef] [PubMed]
  17. Lam, D. How the world survived the population bomb: Lessons from 50 years of extraordinary demographic history. Demography 2011, 48, 1231–1262. [Google Scholar] [CrossRef] [Green Version]
  18. Lorenz, K. On Aggression; Psychology Press: Hove, UK, 2002. [Google Scholar]
  19. Michopoulos, V.; Higgins, M.; Toufexis, D.; Wilson, M.E. Social subordination produces distinct stress-related phenotypes in female rhesus monkeys. Psychoneuroendocrinology 2012, 37, 1071–1085. [Google Scholar] [CrossRef] [Green Version]
  20. Louise, S.; Warrington, N.M.; McCaskie, P.A.; Oddy, W.H.; Zubrick, S.R.; Hands, B.; Mori, T.A.; Briollais, L.; Silburn, S.; Palmer, L.J.; et al. Associations between aggressive behaviour scores and cardiovascular risk factors in childhood. Pediatr. Obes. 2012, 7, 319–328. [Google Scholar] [CrossRef] [PubMed]
  21. Johnson, S.R.; Winkleby, M.A.; Boyce, W.T.; McLaughlin, R.; Broadwin, R.; Goldman, L. The association between hemoglobin and behavior problems in a sample of low-income Hispanic preschool children. J. Dev. Behav. Pediatr. 1992, 13, 209–214. [Google Scholar] [CrossRef]
  22. Ji, N.Y.; Findling, R.L. Pharmacotherapy for mental health problems in people with intellectual disability. Curr. Opin. Psychiatry 2016, 29, 103–125. [Google Scholar] [CrossRef] [Green Version]
  23. Christie, M.R.; Marine, M.L.; Fox, S.E.; French, R.A.; Blouin, M.S. A single generation of domestication heritably alters the expression of hundreds of genes. Nat. Commun. 2016, 7, 10676. [Google Scholar] [CrossRef] [PubMed]
  24. Natt, D.; Rubin, C.J.; Wright, D.; Johnsson, M.; Belteky, J.; Andersson, L.; Jensen, P. Heritable genome-wide variation of gene expression and promoter methylation between wild and domesticated chickens. BMC Genom. 2012, 13, 59. [Google Scholar] [CrossRef] [Green Version]
  25. Kukekova, A.V.; Johnson, J.L.; Xiang, X.; Feng, S.; Liu, S.; Rando, H.M.; Kharlamova, A.V.; Herbeck, Y.; Serdyukova, N.A.; Xiong, Z.; et al. Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviors. Nat. Ecol. Evol. 2018, 2, 1479–1491. [Google Scholar] [CrossRef]
  26. Hekman, J.P.; Johnson, J.L.; Edwards, W.; Vladimirova, A.V.; Gulevich, R.G.; Ford, A.L.; Kharlamova, A.V.; Herbeck, Y.; Acland, G.M.; Raetzman, L.T.; et al. Anterior pituitary transcriptome suggests differences in ACTH release in tame and aggressive foxes. G3 2018, 8, 859–873. [Google Scholar] [CrossRef] [Green Version]
  27. Zapata, I.; Serpell, J.A.; Alvarez, C.E. Genetic mapping of canine fear and aggression. BMC Genom. 2016, 17, 572. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, X.; Zhang, H.; Shang, J.; Liu, G.; Xia, T.; Zhao, C.; Sun, G.; Dou, H. Comparative analysis of the blood transcriptomes between wolves and dogs. Anim. Genet. 2018, 49, 291–302. [Google Scholar] [CrossRef]
  29. Albert, F.W.; Somel, M.; Carneiro, M.; Aximu-Petri, A.; Halbwax, M.; Thalmann, O.; Blanco-Aguiar, J.A.; Plyusnina, I.Z.; Trut, L.; Villafuerte, R.; et al. A comparison of brain gene expression levels in domesticated and wild animals. PLoS Genet. 2012, 8, e1002962. [Google Scholar] [CrossRef] [Green Version]
  30. Zhou, Z.Y.; Li, A.M.; Adeola, A.C.; Liu, Y.H.; Irwin, D.M.; Xie, H.B.; Zhang, Y.P. Genome-wide identification of long intergenic noncoding RNA genes and their potential association with domestication in pigs. Genome Biol. Evol. 2014, 6, 138713–138792. [Google Scholar] [CrossRef] [Green Version]
  31. Alttoa, A.; Koiv, K.; Hinsley, T.A.; Brass, A.; Harro, J. Differential gene expression in a rat model of depression based on persistent differences in exploratory activity. Eur. Neuropsychopharmacol. 2010, 20, 288–300. [Google Scholar] [CrossRef] [PubMed]
  32. Kudryavtseva, N.N.; Smagin, D.A.; Kovalenko, I.L.; Galyamina, A.G.; Vishnivetskaya, G.B.; Babenko, V.N.; Orlov, Y.L. Serotonergic genes in the development of anxiety/depression-like state and pathology of aggressive behavior in male mice: RNA-seq data. Mol. Biol. 2017, 51, 288–300. [Google Scholar] [CrossRef]
  33. Telenti, A.; Pierce, L.; Biggs, W.; di Iulio, J.; Wong, E.; Fabani, M.; Kirkness, E.; Moustafa, A.; Shah, N.; Xie, C.; et al. Deep sequencing of 10,000 human genomes. Proc. Natl. Acad. Sci. USA 2016, 113, 11901–11906. [Google Scholar] [CrossRef] [Green Version]
  34. Sherry, S.; Ward, M.; Kholodov, M.; Baker, J.; Phan, L.; Smigielski, E.; Sirotkin, K. dbSNP: The NCBI database of genetic variation. Nucleic Acids Res. 2001, 29, 308–311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Pocai, B. The ICD-11 has been adopted by the World Health Assembly. World Psychiatry 2019, 18, 371–372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Esmaeili, S.; Hemami, M.R.; Goheen, J.R. Human dimensions of wildlife conservation in Iran: Assessment of human-wildlife conflict in restoring a wide-ranging endangered species. PLoS ONE 2019, 14, e0220702. [Google Scholar] [CrossRef] [Green Version]
  37. Chadaeva, I.; Ponomarenko, M.; Rasskazov, D.; Sharypova, E.; Kashina, E.; Matveeva, M.; Arshinova, T.; Ponomarenko, P.; Arkova, O.; Bondar, N.; et al. 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. BMC Genom. 2016, 17, 995. [Google Scholar] [CrossRef] [Green Version]
  38. Chadaeva, I.; Rasskazov, D.; Sharypova, E.; Savinkova, L.; Ponomarenko, P.; Ponomarenko, M. Candidate SNP markers of social dominance, which may affect the affinity of the TATA-binding protein for human gene promoters. Russ. J. Genet. Appl. Res. 2017, 7, 523–537. [Google Scholar] [CrossRef]
  39. Chadaeva, I.; Ponomarenko, P.; Rasskazov, D.; Sharypova, E.; Kashina, E.; Zhechev, D.; Drachkova, I.; Arkova, O.; Savinkova, L.; Ponomarenko, M.; et al. Candidate SNP markers of reproductive potential are predicted by a significant change in the affinity of TATA-binding protein for human gene promoters. BMC Genom. 2018, 19, 15–37. [Google Scholar] [CrossRef] [Green Version]
  40. Chadaeva, I.; Ponomarenko, P.; Rasskazov, D.; Sharypova, E.; Kashina, E.; Kleshchev, M.; Ponomarenko, M.; Naumenko, V.; Savinkova, L.; Kolchanov, N.; et al. 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. Front. Genet. 2019, 10, 73. [Google Scholar] [CrossRef] [Green Version]
  41. Ponomarenko, M.; Kleshchev, M.; Ponomarenko, P.; Chadaeva, I.; Sharypova, E.; Rasskazov, D.; Kolmykov, S.; Drachkova, I.; Vasiliev, G.; Gutorova, N.; et al. Disruptive natural selection by male reproductive potential prevents underexpression of protein-coding genes on the human Y chromosome as a self-domestication syndrome. BMC Genet. 2020, 21, 89. [Google Scholar] [CrossRef]
  42. Theofanopoulou, C.; Gastaldon, S.; O’Rourke, T.; Samuels, B.D.; Martins, P.T.; Delogu, F.; Alamri, S.; Boeckx, C. Self-domestication in Homo sapiens: Insights from comparative genomics. PLoS ONE 2017, 12, e0185306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Lu, Z. PubMed and beyond: A survey of web tools for searching biomedical literature. Database 2011, 2011, baq036. [Google Scholar] [CrossRef] [Green Version]
  44. Samet, H. A top-down quadtree traversal algorithm. IEEE Trans. Pattern Anal. Mach. Intell. 1985, 7, 94–98. [Google Scholar] [CrossRef] [PubMed]
  45. Sun, G.L.; Shen, W.; Wen, J.F. Triosephosphate isomerase genes in two trophic modes of euglenoids (euglenophyceae) and their phylogenetic analysis. J. Eukaryot. Microbiol. 2008, 55, 170–177. [Google Scholar] [CrossRef] [PubMed]
  46. Morozova, O.V.; Alekseeva, A.E.; Sashina, T.A.; Brusnigina, N.F.; Epifanova, N.V.; Kashnikov, A.U.; Zverev, V.V.; Novikova, N.A. Phylodynamics of G4P[8] and G2P[4] strains of rotavirus A isolated in Russia in 2017 based on full-genome analyses. Virus Genes 2020, 56, 537–545. [Google Scholar] [CrossRef] [PubMed]
  47. Plengpanich, W.; Le Goff, W.; Poolsuk, S.; Julia, Z.; Guerin, M.; Khovidhunkit, W. CETP deficiency due to a novel mutation in the CETP gene promoter and its effect on cholesterol efflux and selective uptake into hepatocytes. Atherosclerosis 2011, 216, 370–373. [Google Scholar] [CrossRef]
  48. Silliman, K.; Tall, A.R.; Kretchmer, N.; Forte, T.M. Unusual high-density lipoprotein subclass distribution during late pregnancy. Metabolism 1993, 42, 1592–1599. [Google Scholar] [CrossRef]
  49. Louder, M.I.M.; Hauber, M.E.; Balakrishnan, C.N. Early social experience alters transcriptomic responses to species-specific song stimuli in female songbirds. Behav. Brain Res. 2018, 347, 69–76. [Google Scholar] [CrossRef]
  50. Grondahl, M.L.; Borup, R.; Vikesa, J.; Ernst, E.; Andersen, C.Y.; Lykke-Hartmann, K. The dormant and the fully competent oocyte: Comparing the transcriptome of human oocytes from primordial follicles and in metaphase II. Mol. Hum. Reprod. 2013, 19, 600–617. [Google Scholar] [CrossRef] [Green Version]
  51. Mann, P.E. Gene expression profiling during pregnancy in rat brain tissue. Brain Sci. 2014, 4, 125–135. [Google Scholar] [CrossRef] [PubMed]
  52. Tomas-Roig, J.; Piscitelli, F.; Gil, V.; Quintana, E.; Ramio-Torrenta, L.L.; Del Rio, J.A.; Moore, T.P.; Agbemenyah, H.; Salinas, G.; Pommerenke, C.; et al. Effects of repeated long-term psychosocial stress and acute cannabinoid exposure on mouse corticostriatal circuitries: Implications for neuropsychiatric disorders. CNS Neurosci. Ther. 2018, 24, 528–538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Inoue, T.; Izumi, T.; Maki, Y.; Muraki, I.; Koyama, T. Effect of the dopamine D(1/5) antagonist SCH 23390 on the acquisition of conditioned fear. Pharmacol. Biochem. Behav. 2000, 66, 573–578. [Google Scholar] [CrossRef]
  54. Schneider, T.; Ilott, N.; Brolese, G.; Bizarro, L.; Asherson, P.J.; Stolerman, I.P. Prenatal exposure to nicotine impairs performance of the 5-choice serial reaction time task in adult rats. Neuropsychopharmacology 2011, 36, 1114–1125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lee, J.Y.; Hong, S.H.; Shin, M.; Heo, H.R.; Jang, I.H. Blockade of FLT4 suppresses metastasis of melanoma cells by impaired lymphatic vessels. Biochem. Biophys. Res. Commun. 2016, 478, 733–738. [Google Scholar] [CrossRef]
  56. Leedom, A.J.; Sullivan, A.B.; Dong, B.; Lau, D.; Gronert, K. Endogenous LXA4 circuits are determinants of pathological angiogenesis in response to chronic injury. Am. J. Pathol. 2010, 176, 74–84. [Google Scholar] [CrossRef]
  57. Wang, S.; Albers, K.M. Behavioral and cellular level changes in the aging somatosensory system. Ann. N. Y. Acad. Sci. 2009, 1170, 745–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Jankowski, M.P.; Rau, K.K.; Soneji, D.J.; Anderson, C.E.; Koerber, H.R. Enhanced artemin/GFRα3 levels regulate mechanically insensitive, heat-sensitive C-fiber recruitment after axotomy and regeneration. J. Neurosci. 2010, 30, 16272–16283. [Google Scholar] [CrossRef] [Green Version]
  59. Lindfors, P.H.; Lindahl, M.; Rossi, J.; Saarma, M.; Airaksinen, M.S. Ablation of persephin receptor glial cell line-derived neurotrophic factor family receptor alpha4 impairs thyroid calcitonin production in young mice. Endocrinology 2006, 147, 2237–2244. [Google Scholar] [CrossRef]
  60. Yang, J.; Lindahl, M.; Lindholm, P.; Virtanen, H.; Coffey, E.; Runeberg-Roos, P.; Saarma, M. PSPN/GFRalpha4 has a significantly weaker capacity than GDNF/GFRalpha1 to recruit RET to rafts, but promotes neuronal survival and neurite outgrowth. FEBS Lett. 2004, 569, 267–271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Park, H.; Oh, C.M.; Park, J.; Park, H.; Cui, S.; Kim, H.S.; Namkung, J.; Park, S.K.; Pak, H.N.; Lee, M.H.; et al. Deletion of the serotonin receptor type 3a in mice leads to sudden cardiac death during pregnancy. Circ. J. 2015, 79, 1807–1815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Kondo, M.; Nakamura, Y.; Ishida, Y.; Shimada, S. The 5-HT3 receptor is essential for exercise-induced hippocampal neurogenesis and antidepressant effects. Mol. Psychiatry 2015, 20, 1428–1437. [Google Scholar] [CrossRef] [PubMed]
  63. Sasaki, Y.; Otsuka, K.; Arimochi, H.; Tsukumo, S.I.; Yasutomo, K. Distinct roles of IL-1β and IL-18 in NLRC4-induced autoinflammation. Front. Immunol. 2020, 11, 591713. [Google Scholar] [CrossRef]
  64. Olkkonen, J.; Kouri, V.P.; Hynninen, J.; Konttinen, Y.T.; Mandelin, J. Differentially expressed in chondrocytes 2 (DEC2) increases the expression of IL-1β and is abundantly present in synovial membrane in rheumatoid arthritis. PLoS ONE 2015, 10, e0145279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Nagaraja, M.R.; Gubbala, S.P.; Delphine Silvia, C.R.W.; Amanchy, R. Molecular diagnostics of disorders of sexual development: An Indian survey and systems biology perspective. Syst. Biol. Reprod. Med. 2019, 65, 105–120. [Google Scholar] [CrossRef]
  66. Jameson, J.L. Of mice and men: The tale of steroidogenic factor-1. J. Clin. Endocrinol. Metab. 2004, 89, 5927–5929. [Google Scholar] [CrossRef] [Green Version]
  67. Li, L.; Chen, X.; Hao, L.; Chen, Q.; Liu, H.; Zhou, Q. Exploration of immune-related genes in high and low tumor mutation burden groups of chromophobe renal cell carcinoma. Biosci. Rep. 2020, 40, bsr20201491. [Google Scholar] [CrossRef] [PubMed]
  68. Coan, P.M.; Hummel, O.; Garcia Diaz, A.; Barrier, M.; Alfazema, N.; Norsworthy, P.J.; Pravenec, M.; Petretto, E.; Hubner, N.; Aitman, T.J. Genetic, physiological and comparative genomic studies of hypertension and insulin resistance in the spontaneously hypertensive rat. Dis. Model Mech. 2017, 10, 297–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Nakao, K.; Iwata, K.; Takeshita, T.; Ozawa, H. Expression of hypothalamic kisspeptin, neurokinin B, and dynorphin A neurons attenuates in female Zucker fatty rats. Neurosci. Lett. 2018, 665, 135–139. [Google Scholar] [CrossRef]
  70. Szklarczyk, K.; Korostynski, M.; Golda, S.; Solecki, W.; Przewlocki, R. Genotype-dependent consequences of traumatic stress in four inbred mouse strains. Genes Brain Behav. 2012, 11, 977–985. [Google Scholar] [CrossRef] [PubMed]
  71. Greig, C.J.; Zhang, L.; Cowles, R.A. Potentiated serotonin signaling in serotonin re-uptake transporter knockout mice increases enterocyte mass and small intestinal absorptive function. Physiol. Rep. 2019, 7, e14278. [Google Scholar] [CrossRef] [Green Version]
  72. Morgan, J.A.; Singhal, G.; Corrigan, F.; Jaehne, E.J.; Jawahar, M.C.; Breen, J.; Pederson, S.; Baune, B.T. Ceasing exercise induces depression-like, anxiety-like, and impaired cognitive-like behaviours and altered hippocampal gene expression. Brain Res. Bull. 2019, 148, 118–130. [Google Scholar] [CrossRef]
  73. Senol, B.K.; Zulfikar, B. Clinical problems and surgical interventions in inherited factor VII deficiency. Turk Pediatri Ars. 2020, 55, 184–190. [Google Scholar]
  74. Burad, J.; Bhakta, P.; Sharma, J. Timely ‘off-label’ use of recombinant activated factor VII (NovoSeven(®)) can help in avoiding hysterectomy in intractable obstetric bleeding complicated with disseminated intravascular coagulation: A case report and review of the literature. Indian J. Anaesth. 2012, 56, 69–71. [Google Scholar]
  75. Qian, C.; Wong, C.W.Y.; Wu, Z.; He, Q.; Xia, H.; Tam, P.K.H.; Wong, K.K.Y.; Lui, V.C.H. Stage specific requirement of platelet-derived growth factor receptor-α in embryonic development. PLoS ONE 2017, 12, e0184473. [Google Scholar] [CrossRef] [Green Version]
  76. Foster, D.S.; Marshall, C.D.; Gulati, G.S.; Chinta, M.S.; Nguyen, A.; Salhotra, A.; Jones, R.E.; Burcham, A.; Lerbs, T.; Cui, L.; et al. Elucidating the fundamental fibrotic processes driving abdominal adhesion formation. Nat. Commun. 2020, 11, 4061. [Google Scholar] [CrossRef]
  77. Brigger, D.; Torbett, B.E.; Chen, J.; Fey, M.F.; Tschan, M.P. Inhibition of GATE-16 attenuates ATRA-induced neutrophil differentiation of APL cells and interferes with autophagosome formation. Biochem. Biophys. Res. Commun. 2013, 438, 283–288. [Google Scholar] [CrossRef] [Green Version]
  78. Cho, H.S.; Park, S.Y.; Kim, S.M.; Kim, W.J.; Jung, J.Y. Autophagy-related protein MAP1LC3C plays a crucial role in odontogenic differentiation of human dental pulp cells. Tissue Eng. Regen. Med. 2020. [Google Scholar] [CrossRef] [PubMed]
  79. Jorgensen, J.O.L.; Juul, A. Therapy of endocrine disease: Growth hormone replacement therapy in adults: 30 years of personal clinical experience. Eur. J. Endocrinol. 2018, 179, R47–R56. [Google Scholar] [CrossRef] [PubMed]
  80. Regan, S.L.P.; Knight, P.G.; Yovich, J.L.; Arfuso, F.; Dharmarajan, A. Growth hormone during in vitro fertilization in older women modulates the density of receptors in granulosa cells, with improved pregnancy outcomes. Fertil. Steril. 2018, 110, 1298–1310. [Google Scholar] [CrossRef]
  81. Takhviji, V.; Zibara, K.; Azarkeivan, A.; Mehrvar, N.; Mehrvar, N.; Mezginejad, F.; Khosravi, A. Fertility and pregnancy in Iranian thalassemia patients: An update on transfusion complications. Transfus. Med. 2020, 30, 352–360. [Google Scholar] [CrossRef] [PubMed]
  82. Wang, F.; Yu, H.; Huang, S.; Zheng, L.; Zheng, P.; Zhang, S.; Li, S.; Chen, J. Jian-Pi-Yi-Shen regulates EPO and iron recycling protein expressions in anemic rats with chronic kidney disease: Accumulation of hypoxia inducible factor-2α via ERK signaling. Evid. Based Complement. Alternat. Med. 2020, 2020, 8894257. [Google Scholar] [CrossRef] [PubMed]
  83. Schulz, A.; Gorodetska, I.; Behrendt, R.; Fuessel, S.; Erdmann, K.; Foerster, S.; Datta, K.; Mayr, T.; Dubrovska, A.; Muders, M.H. Linking NRP2 with EMT and chemoradioresistance in bladder cancer. Front. Oncol. 2020, 9, 1461. [Google Scholar] [CrossRef]
  84. Pellet-Many, C.; Mehta, V.; Fields, L.; Mahmoud, M.; Lowe, V.; Evans, I.; Ruivo, J.; Zachary, I. Neuropilins 1 and 2 mediate neointimal hyperplasia and re-endothelialization following arterial injury. Cardiovasc. Res. 2015, 108, 288–298. [Google Scholar] [CrossRef] [Green Version]
  85. Yang, J.J.; Caligioni, C.S.; Chan, Y.M.; Seminara, S.B. Uncovering novel reproductive defects in neurokinin B receptor null mice: Closing the gap between mice and men. Endocrinology 2012, 153, 1498–1508. [Google Scholar] [CrossRef]
  86. Nilsson, L.L.; Hornstrup, M.B.; Perin, T.L.; Lindhard, A.; Funck, T.; Bjerrum, P.J.; Mule, H.T.; Scheike, T.; Nielsen, H.S.; Hviid, T.V.F. Soluble HLA-G and TGF-β in couples attending assisted reproduction—A possible role of TGF-β isoforms in semen? J. Reprod. Immunol. 2020, 137, 102857. [Google Scholar] [CrossRef]
  87. Ciebiera, M.; Wlodarczyk, M.; Wrzosek, M.; Męczekalski, B.; Nowicka, G.; Lukaszuk, K.; Ciebiera, M.; Słabuszewska-Jozwiak, A.; Jakiel, G. Role of transforming growth factor β in uterine fibroid biology. Int. J. Mol. Sci. 2017, 18, 2435. [Google Scholar] [CrossRef] [Green Version]
  88. Ivanski, F.; de Oliveira, V.M.; de Oliveira, I.M.; de Araujo Ramos, A.T.; de Oliveira Tonete, S.T.; de Oliveira Hykavei, G.; Bargi-Souza, P.; Schiessel, D.L.; Martino-Andrade, A.J.; Romano, M.A.; et al. Prepubertal acrylamide exposure causes dose-response decreases in spermatic production and functionality with modulation of genes involved in the spermatogenesis in rats. Toxicology 2020, 436, 152428. [Google Scholar] [CrossRef]
  89. Creed, M.; Kaufling, J.; Fois, G.R.; Jalabert, M.; Yuan, T.; Luscher, C.; Georges, F.; Bellone, C. Cocaine exposure enhances the activity of ventral tegmental area dopamine neurons via calcium-impermeable NMDARs. J. Neurosci. 2016, 36, 10759–10768. [Google Scholar] [CrossRef] [Green Version]
  90. Dela Pena, I.; Bang, M.; Lee, J.; de la Pena, J.B.; Kim, B.N.; Han, D.H.; Noh, M.; Shin, C.Y.; Cheong, J.H. Common prefrontal cortical gene expression profiles between adolescent SHR/NCrl and WKY/NCrl rats which showed inattention behavior. Behav. Brain Res. 2015, 291, 268–276. [Google Scholar] [CrossRef]
  91. Guo, Y.; Zhang, H.; Gao, J.; Wei, S.; Song, C.; Sun, P.; Qiao, M. Study of genes associated with the ‘anger-in’ and ‘anger-out’ emotions of humans using a rat model. Exp. Ther. Med. 2015, 9, 1448–1454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Jin, Y.; Wang, L.; Duan, Q.; Gong, Z.; Yang, F.; Song, Y. Differential expression of 5-HT-related genes in symptomatic pulmonary embolism patients. Int. J. Clin. Exp. Med. 2015, 8, 512–518. [Google Scholar]
  93. Sackett, S.D.; Otto, T.; Mohs, A.; Sander, L.E.; Strauch, S.; Streetz, K.L.; Kroy, D.C.; Trautwein, C. Myeloid cells require gp130 signaling for protective anti-inflammatory functions during sepsis. FASEB J. 2019, 33, 6035–6044. [Google Scholar] [CrossRef] [PubMed]
  94. Gray, S.R.; Robinson, M.; Nimmo, M.A. Response of plasma IL-6 and its soluble receptors during submaximal exercise to fatigue in sedentary middle-aged men. Cell Stress Chaperones 2008, 13, 247–251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Sun, Y.; Liu, S.; Hu, R.; Zhou, Q.; Li, X. Decreased placental IL9 and IL9R in preeclampsia impair trophoblast cell proliferation, invasion, and angiogenesis. Hypertens. Pregnancy 2020, 39, 228–235. [Google Scholar] [CrossRef] [PubMed]
  96. Osterfeld, H.; Ahrens, R.; Strait, R.; Finkelman, F.D.; Renauld, J.C.; Hogan, S.P. Differential roles for the IL-9/IL-9 receptor alpha-chain pathway in systemic and oral antigen-induced anaphylaxis. J. Allergy Clin. Immunol. 2010, 125, 469–476.e2. [Google Scholar] [CrossRef] [Green Version]
  97. Hahm, S.; Mizuno, T.M.; Wu, T.J.; Wisor, J.P.; Priest, C.A.; Kozak, C.A.; Boozer, C.N.; Peng, B.; McEvoy, R.C.; Good, P.; et al. Targeted deletion of the Vgf gene indicates that the encoded secretory peptide precursor plays a novel role in the regulation of energy balance. Neuron 1999, 23, 537–548. [Google Scholar] [CrossRef] [Green Version]
  98. Sainsbury, A.; Schwarzer, C.; Couzens, M.; Herzog, H. Y2 receptor deletion attenuates the type 2 diabetic syndrome of ob/ob mice. Diabetes 2002, 51, 3420–3427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Sanford, L.P.; Ormsby, I.; Gittenberger-de Groot, A.C.; Sariola, H.; Friedman, R.; Boivin, G.P.; Cardell, E.L.; Doetschman, T. TGFbeta2 knockout mice have multiple developmental defects that are non-overlapping with other TGFbeta knockout phenotypes. Development 1997, 124, 2659–2670. [Google Scholar] [PubMed]
  100. Kim, J.H.; Lee, Y.W.; Park, Y.M.; Park, K.A.; Park, S.H.; Lee, W.T.; Lee, J.E. Agmatine-reduced collagen scar area accompanied with surface righting reflex recovery after complete transection spinal cord injury. Spine 2011, 36, 2130–2138. [Google Scholar] [CrossRef]
  101. Osadchuk, L.V. Endocrine gonadal function in silver fox under domestication. Scientifur 1992, 16, 116–121. [Google Scholar]
  102. Osadchuk, L.V. Some peculiarities in reproduction in silver fox males under domestication. Scientifur 1992, 16, 285–288. [Google Scholar]
  103. Osadchuk, L.V. Black-silver fox Vulpes vulpes male reproductive potential after longitudinal selection on domestic behavior. Zh. Evol. Biokhim. Fiziol. 2006, 42, 146–152. [Google Scholar]
  104. Osadchuk, L.V.; Krass, P.M.; Trut, L.N.; Belyaev, D.K. Effects of selection for behavior on the endocrine function of the gonads in male silver-black foxes. Dokl. Akad. Nauk USSR 1978, 240, 1255–1258. [Google Scholar]
  105. Osadchuk, L.V.; Krass, P.M.; Trut, L.N.; Ivanova, L.N. Gonadal endocrine function in male silver foxes with different hereditary determined forms of defensive behavior. Proc. Sib. Branch Acad. Nauk USSR 1978, 79–86. [Google Scholar]
  106. Osadchuk, L.V. Biosynthesis of testosterone in the gonads in silver fox embryos after long-term selection for domesticated behavior. Genetika 1998, 34, 941–946. [Google Scholar] [PubMed]
  107. Prasolova, L.A.; Gerbek, Y.E.; Gulevich, R.G.; Shikhevich, S.G.; Konoshenko, M.Y.; Kozhemyakina, R.V.; Oskina, I.N.; Plyusnina, I.Z. The effects of prolonged selection for behavior on the stress response and activity of the reproductive system of male grey mice (Rattus norvegicus). Russ. J. Genet. 2014, 50, 846–852. [Google Scholar] [CrossRef]
  108. Liu, X.; Mawolo, J.B.; Du, X.; Zhou, Y.; Wang, H.; Liu, F.; He, Z.; Marela, H.A. Investigation of biochemical and physiological parameters of the newborn Saiga antelope (Saiga tatarica) in Gansu Province, China. PLoS ONE 2019, 14, e0224822. [Google Scholar] [CrossRef] [PubMed]
  109. Xiao, W.; Feng, L.; Mou, P.; Miquelle, D.G.; Hebblewhite, M.; Goldberg, J.F.; Robinson, H.S.; Zhao, X.; Zhou, B.; Wang, T.; et al. Estimating abundance and density of Amur tigers along the Sino-Russian border. Integr. Zool. 2016, 11, 322–332. [Google Scholar] [CrossRef]
  110. Almeida, F.F.; Leal, M.C.; Franca, L.R. Testis morphometry, duration of spermatogenesis, and spermatogenic efficiency in the wild boar (Sus scrofa scrofa). Biol. Reprod. 2006, 75, 792–799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  111. Cai, Y.; Yang, J.; Wang, J.; Yang, Y.; Fu, W.; Zheng, C.; Cheng, J.; Zeng, Y.; Zhang, Y.; Xu, L.; et al. Changes in the population genetic structure of captive forest musk deer (Moschus berezovskii) with the increasing number of generation under closed breeding conditions. Animals 2020, 10, 255. [Google Scholar] [CrossRef] [Green Version]
Ijms 22 02346 g001 550
Figure 1. An algorithmic flow chart of the 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. Legend: Δ, gene expression change: deficit (↓), excess (↑); ♂♀, reproductive potential change: worse (←), improved (→); NCA, nearest common ancestor; PADJ, significance (Fisher’s Z-test with those corrections on multiple comparisons, which are available within the RNA-Seq data, as published by their authors [26,28,29]).
Figure 1. An algorithmic flow chart of the 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. Legend: Δ, gene expression change: deficit (↓), excess (↑); ♂♀, reproductive potential change: worse (←), improved (→); NCA, nearest common ancestor; PADJ, significance (Fisher’s Z-test with those corrections on multiple comparisons, which are available within the RNA-Seq data, as published by their authors [26,28,29]).
Ijms 22 02346 g001
Table
Table 1. The investigated genome-wide RNA-Seq transcriptomes of domestic animals with their wild congeners, publicly available in database PubMed.
Table 1. The investigated genome-wide RNA-Seq transcriptomes of domestic animals with their wild congeners, publicly available in database PubMed.
#Domestic AnimalsWild AnimalsTissueNumber of DEGs[Ref]
1tame foxes (Vulpes vulpes)
6 males:		aggressive foxes (V. vulpes)
6 males:		pituitary327[26]
2dogs (Canis familiaris):
1 female and 1 male		wolves (C. lupus):
2 females and 1 male		blood450[28]
3dogs (C. familiaris):
2 females and 3 males		wolves (C. lupus):
2 females and 1 male		frontal cortex13[29]
4pigs (Sus scrofa):
5 females		boars (S. scrofa):
5 females		frontal cortex30[29]
5guinea pigs (Cavia porcellus):
3 females and 3 males		cavy (C. aperea):
3 females and 3 males		frontal cortex883[29]
6domesticated rabbits
(Oryctolagus cuniculus domesticus):
3 females and 3 males		wild rabbits
(Oryctolagus cuniculus):
3 females and 3 males		frontal cortex17[29]
7tame rats (Rattus norvegicus):
3 females and 3 males		aggressive rats (R. norvegicus):
3 females and 3 males		frontal cortex20[29]
Σ6 domesticated animal species:
17 females and 19 males		6 wild animal species:
18 females and 17 males		3 tissues17403 Refs
Note: DEG, differentially expressed gene; Ref, reference.
Table
Table 2. Comparing the effects of changes in the expression of orthologous genes on human reproductive potential and on the divergence of the guinea pig and cavy from their nearest common ancestor (NCA) [29].
Table 2. Comparing the effects of changes in the expression of orthologous genes on human reproductive potential and on the divergence of the guinea pig and cavy from their nearest common ancestor (NCA) [29].
HumansAnimals
GeneEffect of Gene Expression Change (Δ) on Reproductive Potential (♂♀)DEGRNA-SeqΔ during Diver-Gence from NCA
Deficit (↓)♂♀Excess (↑)♂♀logPADJJDeficit (↓)Excess (↑)
CETPlesser risk of myo-cardial infarction [47]→hypercholesterolemia in pregnancy [48]←Cetp2.110−3wilddomestic
CHRNA3improved finding oppo-site sex congeners [49]→greater nicotine effects on oocytes [50]←Chrna30.90.05wilddomestic
CHRNA6improved maternal behavior [51]→higher risk of social defeats [52]←Chrna60.90.05wilddomestic
DRD5reduces conditioned fear [53]→higher risk of mental disorders [54]←Drd5–1.910−2domesticwild
FLT4suppressed melanoma metastasis [55]→worse post-injury neo- vascularization [56]←Flt40.80.05wilddomestic
GFRA3accelerated neuro-degeneration [57]←improved neural regeneration [58]→Gfra3–1.00.05domesticwild
GFRA4premature adolescent bone formation [59]←improved neuronal survival [60]→Gfra41.50.05wilddomestic
HTR3Ahigher risk of death during pregnancy [61]←improved mood and behavior [62]→Htr3a–2.910−14domesticwild
IL1Bless bone deformation in infections [63]→circadian pain hypersensitivity [64]←Il1b2.310−2wilddomestic
NR5A1higher risk of gonadal dysgenesis in men [65]←improved gonadal development [66]→Nr5a1–2.20.05domesticwild
PDGFRLreduced tumor mutation burden [67]→myocardial hypertrophy [68]←Pdgfrl1.310−8wilddomestic
PDYNobesity-related subfertility [69]←prevented conditioned fear behavior [70]→Pdyn0.910−2wilddomestic
SLC6A4improved small-intestinal function [71]→worse depression, anxiety, inertia [72]←Slc6a42.910−2wilddomestic
Note: see the caption of Figure 1; log, the ratio of a domestic-animal gene expression level to that in wild animals (log2 units); PADJ, significance (Fisher’s Z-test with those corrections on multiple comparisons, which are available within the RNA-Seq data, as published by their authors [26,28,29]).
Table
Table 3. Correlations between the effects of the co-directed changes in the expression of orthologous genes on human reproductive potential and on the divergence of the guinea pig and cavy from their NCA.
Table 3. Correlations between the effects of the co-directed changes in the expression of orthologous genes on human reproductive potential and on the divergence of the guinea pig and cavy from their NCA.
HumansChange in Reproductive Potential (♂♀)Binomial Distributionχ2-TestFisher’s Test
Animals Worse (←)Improved (→)χ2pValuep
change during divergence from NCAdomestic1120.059.9110−20.0030.05
wild3100.06
Note: see Figure 1 caption.
Table
Table 4. Comparing the effects of changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their nearest common ancestor (NCA).
Table 4. Comparing the effects of changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their nearest common ancestor (NCA).
HumansAnimals
GeneEffect of Gene Expression Change (Δ) on Reproductive Potential (♂♀)DEGRNA-SeqΔ during Diver-Gence from NCA
Deficit (↓)♂♀Excess (↑)♂♀logPADJJDeficit (↓)Excess (↑)
domesticated versus wild rabbits [29]
F7higher risk of bleeding [73]←recombinant F7 treats obstetric bleeding [74]→F7–2.70.05domesticwild
dog versus wolf (frontal cortex) [29]
PDGFRAskeletal defects in newborns [75]←higher risk of infertility [76]←Pdgfra1.510−3wilddomestic
dog versus wolf (blood) [28]
GABARAPL2impaired wound healing [77]←improved tooth injury healing [78]→Gabarapl2–2.010−2domesticwild
GH1higher risk of mortality [79]←prolonged reproductive age in women [80]→Gh12.810−2wilddomestic
HBBworse reproductive health in women [81]←relieved anemia in kidney diseases [82]→Hbbl–5.910−8domesticwild
NRP2better survival after radiochemotherapy [83]→vascular neointimal hyperplasia [84]←Nrp21.80.05wilddomestic
TAC3higher risk of subfertility [85]←lesser socially induced subfertility [85]→Tac3–4.710−8domesticwild
TGFB3lowers semen quality and infertility [86]←higher risk of female infertility [87]←Tgfb33.30.05wilddomestic
tame versus aggressive foxes [26]
ESR2impaired spermatogenesis [88]←impaired spermatogenesis [88]←Esr2–0.30.05domesticwild
GRIN3Aprevents cocaine addiction [89]→higher risk of inatten-tive behavior [90]←Grin3a0.510−2wilddomestic
HTR3Breduced anger-reso-lutive behavior [91]←lesser risk of pulmo-nary embolism [92]→Htr3b–0.50.05domesticwild
IL6SThigher risk of mortality during sepsis [93]←increased sensitivity to fatigue [94]←Il6st0.30.05wilddomestic
IL9Rimpaired trophoblast implantation [95]←higher risk of anaphylaxis [96]←Il9r0.40.05wilddomestic
NPYhigher risk of infertility [97]←higher risk of obesity and subfertility [98]←Npy0.410−2wilddomestic
TGFB2higher risk of perinatal mortality [99]←impaired wound healing [100]←Tgfb20.510−2wilddomestic
Note: see the caption of Figure 1 and the footnote of Table 2.
Table
Table 5. Correlations between the effects of the co-directed changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their NCA (without guinea pig and cavy, the correlations for which are outlined in Table 2 and Table 3).
Table 5. Correlations between the effects of the co-directed changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their NCA (without guinea pig and cavy, the correlations for which are outlined in Table 2 and Table 3).
HumansChange in Reproductive Potential (♂♀)Binomial Distributionχ2-TestFisher’s Test
Animals Worse (←)Improved (→)χ2pValuep
change during divergence from NCAdomestic14110−36.140.050.040.05
wild870.5
Note: see Figure 1 caption.
Table
Table 6. Correlations between the effects of the co-directed changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their NCA.
Table 6. Correlations between the effects of the co-directed changes in the expression of orthologous genes on human reproductive potential and on the divergence of domestic and wild animals from their NCA.
HumansChange in Reproductive Potential (♂♀)Binomial Distributionχ2-TestFisher’s Test
Animals Worse (←)Improved (→)χ2pValuep
change during divergence from NCAdomestic25310−415.210−310−40.05
wild11170.1
Note: see Figure 1 caption.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vasiliev, G.; Chadaeva, I.; Rasskazov, D.; Ponomarenko, P.; Sharypova, E.; Drachkova, I.; Bogomolov, A.; Savinkova, L.; Ponomarenko, M.; Kolchanov, N.; et al. 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. Int. J. Mol. Sci. 2021, 22, 2346. https://doi.org/10.3390/ijms22052346

AMA Style

Vasiliev G, Chadaeva I, Rasskazov D, Ponomarenko P, Sharypova E, Drachkova I, Bogomolov A, Savinkova L, Ponomarenko M, Kolchanov N, et al. 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. International Journal of Molecular Sciences. 2021; 22(5):2346. https://doi.org/10.3390/ijms22052346

Chicago/Turabian Style

Vasiliev, Gennady, Irina Chadaeva, Dmitry Rasskazov, Petr Ponomarenko, Ekaterina Sharypova, Irina Drachkova, Anton Bogomolov, Ludmila Savinkova, Mikhail Ponomarenko, Nikolay Kolchanov, and et al. 2021. "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" International Journal of Molecular Sciences 22, no. 5: 2346. https://doi.org/10.3390/ijms22052346

APA Style

Vasiliev, G., Chadaeva, I., Rasskazov, D., Ponomarenko, P., Sharypova, E., Drachkova, I., Bogomolov, A., Savinkova, L., Ponomarenko, M., Kolchanov, N., Osadchuk, A., Oshchepkov, D., & Osadchuk, L. (2021). 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. International Journal of Molecular Sciences, 22(5), 2346. https://doi.org/10.3390/ijms22052346

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop