Comparative analysis of the blood transcriptomes between wolves and dogs

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 (log₂ 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 log₂ 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.

Transcriptome analysis provides new insights into cold adaptation of corsac fox ( Vulpes Corsac )

Vulpesare widely distributed throughout the world and have undergone drastic physiological and phenotypic changes in response to their environment. However, little is known about the underlying genetic causes of these traits, especially Vulpes corsac. In this study, RNA‐Seq was used to obtain a comprehensive dataset for multiple pooled tissues of corsac fox, and selection analysis of orthologous genes was performed to identify the genes that may be influenced by the low‐temperature environment. More than 6.32 Gb clean reads were obtained and assembled into a total of 173,353 unigenes with an average length of 557 bp for corsac fox. Selective pressure analysis showed that 16 positively selected genes (PSGs) were identified in corsac fox, red fox, and arctic fox. Enrichment analysis of PSGs showed that the LRP11 gene was enriched in several pathways related to the low‐temperature response and might play a key role in response to environmental stimuli of foxes. In addition, several positively selected genes were related to DNA damage repair (ELP2 and CHAF1A), innate immunity (ARRDC4 and S100A12), and the respiratory chain (NDUFA5), and these positively selected genes might play a role in adaptation to harsh wild fox environments. The results of common orthologous gene analysis showed that gene flow or convergent evolution might be an important factor in promoting regional differentiation of foxes. Our study provides a valuable transcriptomic resource for the evolutionary history of the corsac fox and the adaptations to the extreme environments.

Poly(A) RNA sequencing reveals age-related differences in the prefrontal cortex of dogs

Dogs may possess a unique translational potential to investigate neural aging and dementia because they are prone to age-related cognitive decline, including an Alzheimer’s disease–like pathological condition. Yet very little is known about the molecular mechanisms underlying canine cognitive decline. The goal of the current study was to explore the transcriptomic differences between young and old dogs’ frontal cortex, which is a brain region often affected by various forms of age-related dementia in humans. RNA isolates from the frontal cortical brain area of 13 pet dogs, which represented 7 different breeds and crossbreds, were analyzed. The dogs were euthanized for medical reasons, and their bodies had been donated by their owners for scientific purposes. The poly(A) tail RNA subfraction of the total transcriptome was targeted in the sequencing analysis. Cluster analyses, differential gene expression analyses, and gene ontology analyses were carried out to assess which genes and genetic regulatory mechanisms were mostly affected by aging. Age was the most prominent factor in the clustering of the animals, indicating the presence of distinct gene expression patterns related to aging in a genetically variable population. A total of 3436 genes were found to be differentially expressed between the age groups, many of which were linked to neural function, immune system, and protein synthesis. These findings are in accordance with previous human brain aging RNA sequencing studies. Some genes were found to behave more similarly to humans than to rodents, further supporting the applicability of dogs in translational aging research.

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

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

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

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

Identification of key molecular markers of acute coronary syndrome using peripheral blood transcriptome sequencing analysis and mRNA-lncRNA co-expression network construction

Acute coronary syndrome (ACS) is a term used to describe major cardiovascular diseases, and treatment of in-stent restenosis in patients with ACS remains a major clinical challenge. Further investigation into molecular markers of ACS may aid early diagnosis, and the treatment of ACS and post-treatment recurrence. In the present study, total RNA was extracted from the peripheral blood samples of 3 patients with ACS, 3 patients with percutaneous coronary intervention (PCI)_non-restenosis, 3 patients with PCI_restenosis and 3 healthy controls. Subsequently, RNA library construction and high-throughput sequencing were performed. DESeq2 package in R was used to screen genes that were differentially expressed between the different samples. Moreover, the intersection of the differentially expressed mRNAs (DEmRNAs) and differentially expressed long noncoding RNAs (DElncRNAs) obtained. GeneCodis4.0 was used to perform function enrichment for DEmRNAs, and lncRNA-mRNA co-expression network was constructed. The GSE60993 dataset was utilized for diagnostic analysis, and the aforementioned investigations were verified using in vitro studies. Results of the present study revealed a large number of DEmRNAs and DElncRNAs in the different groups. We selected genes in the top 10 of differential expression and also involved in the co-expression of lncRNA-mRNA for diagnostic analysis in the GSE60993 dataset. The area under curve (AUC) of PDZK1IP1 (0.747), PROK2 (0.769) and LAMP3 (0.725) were all >0.7. These results indicated that the identified mRNAs and lncRNAs may act as potential clinical biomarkers, and more specifically, PDZK1IP1, PROK2 and LAMP3 may act as potential biomarkers for the diagnosis of ACS.

Comparative transcriptome provides insights into the selection adaptation between wild and farmed foxes

The silver fox and blue fox are economically important fur species and were domesticated by humans from their wild counterparts, the arctic fox and red fox, respectively. Farmed foxes show obvious differences from their wild counterparts, including differences in physiology, body size, energy metabolism, and immunity. However, the molecular mechanisms underlying these differences are presently unclear. In this study, we built transcriptome libraries from multiple pooled tissues for each species of farmed fox, used RNA-seq to obtain a comprehensive dataset, and performed selection analysis and sequence-level analyses of orthologous genes to identify the genes that may be influenced by human domestication. More than 153.3, 248.0, 81.6, and 65.8 million clean reads were obtained and assembled into a total of 118,577, 401,520, 79,900, and 186,988 unigenes with an average length range from 521 to 667 bp for AF, BF, RF, and SF, respectively. Selective pressure analysis showed that 11 and 14 positively selected genes were identified, respectively, in the two groups (AF vs. BF and RF vs. SF). Several of these genes were associated with natural immunity (CFI and LRRFIP1), protein synthesis (GOLGA4, CEP19 and SLC35A2), and DNA damage repair (MDC1). Further functional enrichment analyses demonstrated that two positively selected genes (ACO1 and ACAD10) were involved in metabolic process (GO:0008152, p-value = .032), representing a significant enrichment. Sequence analysis of 117 orthologous genes shared by the two groups showed that the LEMD2, RRBP1, and IGBP1 genes might be affected by artificial selection in farmed foxes, with mutation sites located within sequences that are otherwise highly conserved across most mammals. Our results provide a valuable transcriptomic resource for future genetic studies and improvement in the assisted breeding of foxes and other farmed animals.

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

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

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

Gene expression profile of peripheral blood mononuclear cells in mild to moderate obesity in dogs

Background

Molecular mechanisms and early diagnosis on the development of mild to moderate of canine obesity are not understood although recent dog obesity is a widespread problem. To understand the differences between normal weight and mild to moderate obesity, the purpose of this study is to investigate the gene expression profiles of peripheral blood mononuclear cells (PBMC) in dogs.

Methods

This study comprised a sample of 12 privately-owned Miniature Dachshund, which were divided into two groups (obese and control) based on body condition scores (BCS). Serum biochemical parameters and PBMC gene expression profiles were compared between groups.

Results

A statistically significant between group differences was recorded for body weight (BW), BCS, serum Insulin and triglyceride (TG) levels (p < 0.05). RNA-seq revealed the upregulated 154 genes and the downregulated 198 genes in obese dogs at more than 3.5-fold change compared with control animals. Hemoglobin subunits alpha- and beta-like were detected in the downregulated genes. RT-PCR analysis showed downregulation of FOLH1, ALAS2 and LOC100855540 genes, and upregulation of BCL2L15 gene, suggesting that the metabolic difference between normal and mild to moderate obesity was involved in the hemoglobin metabolism.

Conclusions

This study revealed significant differences in the gene expression of BCL2L15, FOLH1, ALAS2, and hemoglobin subunits such as LOC100855540 between normal weight and mild to moderate obese dogs, which indicate that these genes may prevent the obesity in dogs and be potentially useful for diagnosis of mild to moderate obesity.

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

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions

The underlying neurological events accompanying dog domestication remain elusive. To reconstruct the domestication process in an experimental setting, silver foxes (Vulpes vulpes) have been deliberately bred for tame vs aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. The hypothalamus is an essential part of the hypothalamic-pituitary-adrenal axis and regulates the fight-or-flight response, and thus, we hypothesized that selective breeding for tameness/aggressiveness has shaped the hypothalamic transcriptomic profile. RNA-seq analysis identified 70 differentially expressed genes (DEGs). Seven of these genes, DKKL1, FBLN7, NPL, PRIMPOL, PTGRN, SHCBP1L and SKIV2L, showed the same direction expression differences in the hypothalamus, basal forebrain and prefrontal cortex. The genes differentially expressed across the three tissues are involved in cell division, differentiation, adhesion and carbohydrate processing, suggesting an association of these processes with selective breeding. Additionally, 159 transcripts from the hypothalamus demonstrated differences in the abundance of alternative spliced forms between the tame and aggressive foxes. Weighted gene coexpression network analyses also suggested that gene modules in hypothalamus were significantly associated with tame vs aggressive behavior. Pathways associated with these modules include signal transduction, interleukin signaling, cytokine-cytokine receptor interaction and peptide ligand-binding receptors (eg, G-protein coupled receptor [GPCR] ligand binding). Current studies show the selection for tameness vs aggressiveness in foxes is associated with unique hypothalamic gene profiles partly shared with other brain regions and highlight DEGs involved in biological processes such as development, differentiation and immunological responses. The role of these processes in fox and dog domestication remains to be determined.

Comparative Analysis of the Liver and Spleen Transcriptomes between Holstein and Yunnan Humped Cattle

Simple Summary

Cattle are important agricultural animals that provide essential sources of meat, milk, fertilizer for crops, clothing, and animal traction for human use, and the demand for these products has increased in recent years. There are existing differences in disease resistance between different cattle breeds. However, the genetic basis underlying disease resistance differences is poorly understood and requires further investigation. In this study, many immune- and disease-relevant genes and pathways were identified between Holstein and Yunnan humped cattle using RNA-sequencing. The novel findings regarding the genetic basis underlying disease resistance differences between zebu cattle and taurine cattle will provide a scientific basis and key technical support for disease-resistant breeding of domestic cattle, and thus have important social and economic significance.

Abstract

Previous studies have shown that Yunnan humped cattle have higher disease resistance than pure taurine cattle, such as Holsteins. However, there exists limited information about the molecular genetic basis underlying disease resistance differences between them. The objective of this study was to compare differentially expressed genes (DEGs) in the liver and spleen tissues of Holstein and Yunnan humped cattle through comparative transcriptome analysis, using RNA-sequencing. In total, 1564 (647 up- and 917 down-regulated genes) and 1530 (716 up- and 814 down-regulated genes) DEGs were obtained in the liver and spleen tissues of Holstein and Yunnan humped cattle comparison groups, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the DEGs were mainly associated with the RIG-I signaling pathway, immune responses, major histocompatibility complex (MHC) class I protein complex and complement activation, human T-cell lymphotropic virus type-I (HTLV-I) infection. Some genes related to immune function, such as C1QB, CD55, MASP2, C4BPA, MAVS, NOD2, and CD46, were up-regulated in Yunnan humped cattle, while C2, SERPING1, SERPINE1, TIRAP, TLR2, and TLR6 were down-regulated. The expression levels of 11 selected DEGs, analyzed by quantitative reverse-transcription polymerase chain reaction (RT-qPCR), were consistent with the deep sequencing results by RNA-sequencing. Our results will provide a scientific basis and key technical support for disease-resistant breeding of domestic cattle.