Transcriptome comparison in the pituitary-adrenal axis between Beagle and Chinese Field dogs after chronic stress exposure

Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent

Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1-13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.

Genes Participating in the Ensheathment of Neurons Are Affected by Postnatal Stress and Maternal Immune Activation in the Pituitary Gland

Immune challenges during gestation are associated with neurodevelopmental disorders and can interact with stress later in life. The pituitary gland participates in endocrine- and immune-related processes that influence development, growth, and reproduction and can modulate physiological and behavioral responses to challenges. The objective of this study was to investigate the effect of stressors at different time points on the molecular mechanisms of the pituitary gland and detect sex differences. RNA sequencing was used to profile the pituitary glands of female and male pigs exposed to weaning stress and virally induced maternal immune activation (MIA), relative to unchallenged groups. Significant effects (FDR-adjusted p-value < 0.05) of MIA and weaning stress were detected in 1829 and 1014 genes, respectively. Of these, 1090 genes presented significant interactions between stressors and sex. The gene ontology biological process of the ensheathment of neurons (GO:0007272), substance abuse, and immuno-related pathways, including the measles disease (ssc05162), encompasses many genes with profiles impacted by MIA and weaning stress. A gene network analysis highlighted the under-expression of myelin protein zero (Mpz) and inhibitors of DNA binding 4 (Id4) among the non-stressed males exposed to MIA, relative to the control and non-MIA males exposed to weaning stress, relative to non-stressed pigs. The detection of changes in the molecular mechanisms of the pituitary gland could advance our understanding of disruptions in the formation of the myelin sheath and the transmission of neuron-to-neuron signals in behavioral disorders associated with maternal immune activation and stress.

An age stratified analysis of the biomarkers in patients with colorectal cancer

Colorectal cancer (CRC), a common malignant tumor of the digestive tract, has a high incidence and mortality rate. Several recent studies have found that aging is associated with the increasing risk of cancer. Nevertheless, the expression status and function of age-related genes in CRC is still not well understood. In the study, we comprehensively analyzed the gene expression data of CRC patients from The Cancer Genome Atlas (TCGA) database. Age-related differential expression genes (age-related DEGs) in tumor tissues compared with normal tissues of CRC were further identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of age-related DEGs were performed by clusterProfiler of R. Afterwards, we used the STRING database to map the protein-protein interaction network of DEGs. We constructed prognostic model through univariate and multivariate COX regression analyses, and further evaluated their predictive power. The prognostic gene signature-related functional pathways were explored by gene set enrichment analysis (GSEA). The weighted gene co-expression network analysis (WGCNA) was used to identify key module associated with two prognostic gene signatures. Finally, we used the Metascape to perform functional enrichment analysis of genes in the key module. A total of 279 age-related DEGs were identified from the TCGA database. GO and KEGG enrichment analysis showed that the age-related DEGs were enriched in the Modulation of chemical synaptic transmission and Neuroactive ligand-receptor interaction. Moreover, we established a novel age-related gene signature (DLX2 and PCOLCE2) for overall survival in CRC, which was further predicted in both the training and validation sets. The results of GSEA demonstrated that numerous disease-related pathways were enriched in the high-risk group. We identified 43 genes related to the DLX2 and PCOLCE2 by the WGCNA co-expression network. We also found that these 43 genes were enriched in the cancer-related pathways. To sum up, the study identified an age-related gene signature for predicting the prognosis of CRC patients, which is conducive to the identification of novel prognostic molecular markers.

A Genome-Wide Detection of Copy Number Variations Using SNP Genotyping Arrays in Braque Français Type Pyrénées Dogs

Copy number variants (CNVs) are an important source of genetic variation complementary to single nucleotide polymorphisms (SNPs). Only few studies have been conducted in dogs on CNVs derived from high-density SNP array data, and many canine breeds still remain uncharacterized, e.g., the Braque Français, type Pyrénées breed (BRA). Therefore, in an effort to more comprehensively investigate the canine genome for CNVs, we used a high-density SNP array (170 K) to discover CNVs in BRA. The CNV regions (CNVRs) were identified through the merging of two different CNVRs datasets, obtained separately from SNP data using the PennCNV and SVS software. A total of 45 stringent CNVRs, ranging from 3.5 kb to 458,716 kb in length were detected in 26 dog samples. Results overlapped moderately in comparison with previous studies on CNVs in dogs, leading to the identification of 16 novel CNVRs. A total of 159 genes were annotated in the CNVRs detected with stringent quality criteria in particular high classification stringency and false discovery rate correction. The gene ontology enrichment analysis provided information on biological processes and cellular components related to muscle structure development and muscle cell differentiation. Considering that BRA is a breed used for speed in hunting and retrieval, for the ability to find feathered game, and for pointing, we can hypothesize that selection for such hunting behavior could have driven, at least in part, the presence of these genes into the CNVRs.

Sex-biased transcriptomic response of the reproductive axis to stress

Stress is a well-known cause of reproductive dysfunction in many species, including birds, rodents, and humans, though males and females may respond differently. A powerful way to investigate how stress affects reproduction is by examining its effects on a biological system essential for regulating reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. Often this is done by observing how a stressor affects the amount of glucocorticoids, such as cortisol or corticosterone, circulating in the blood and their relationship with a handful of known HPG-producing reproductive hormones, like testosterone and estradiol. Until now, we have lacked a full understanding of how stress affects all genomic activity of the HPG axis and how this might differ between the sexes. We leveraged a highly replicated and sex-balanced experimental approach to test how male and female rock doves (Columba livia) respond to restraint stress at the level of their transcriptome. Females exhibit increased genomic responsiveness to stress at all levels of their HPG axis as compared to males, and these responsive genes are mostly unique to females. Reasons for this may be due to fluctuations in the female endocrine environment over the reproductive cycle and/or their evolutionary history, including parental investment and the potential for maternal effects. Direct links between genome to phenome cause and effect cannot be ascertained at this stage; however, the data we report provide a vital genomic foundation on which sex-specific reproductive dysfunction and adaptation in the face of stress can be further experimentally studied, as well as novel gene targets for genetic intervention and therapy investigations.

Comparison of miRNA expression profiles in pituitary-adrenal axis between Beagle and Chinese Field dogs after chronic stress exposure

MicoRNAs (miRNAs), usually as gene regulators, participate in various biological processes, including stress responses. The hypothalamus–pituitary–adrenal axis (HPA axis) is an important pathway in regulating stress response. Although the mechanism that HPA axis regulates stress response has been basically revealed, the knowledge that miRNAs regulate stress response within HPA axis, still remains poor. The object of this study was to investigate the miRNAs in the pituitary and adrenal cortex that regulate chronic stress response with high-throughput sequencing. The pituitary and adrenal cortex of beagles and Chinese Field dogs (CFD) from a stress exposure group (including beagle pituitary 1 (BP1), CFD pituitary 1 (CFDP1), beagle adrenal cortex 1 (BAC1), CFD adrenal cortex 1 (CFDAC1)) and a control group (including beagle pituitary 2 (BP2), CFD pituitary 2 (CFDP2), beagle adrenal cortex 2 (BAC2), CFD adrenal cortex 2 (CFDAC2)), were selected for miRNA-seq comparisons. Comparisons, that were made in pituitary (including BP1 vs. BP2, CFDP1 vs. CFDP2, BP1 vs. CFDP1 and BP2 vs. CFDP2) and adrenal cortex (including BAC1 vs. BAC2, CFDAC1 vs. CFDAC2, BAC1 vs. CFDAC1 and BAC2 vs. CFDAC2), showed that a total of 39 and 18 common differentially expressed miRNAs (DE-miRNAs) (Total read counts > 1,000, Fold change > 2 & p-value < 0.001), that shared in at least two pituitary comparisons and at least two adrenal cortex comparisons, were detected separately. These identified DE-miRNAs were predicted for target genes, thus resulting in 3,959 and 4,010 target genes in pituitary and adrenal cortex, respectively. Further, 105 and 10 differentially expressed genes (DEGs) (Fold change > 2 & p-value < 0.05) from those target genes in pituitary and adrenal cortex were obtained separately, in combination with our previous corresponding transcriptome study. Meanwhile, in line with that miRNAs usually negatively regulated their target genes and the dual luciferase reporter assay, we finally identified cfa-miR-205 might play an important role by upregulating MMD in pituitary and hippocampus, thus enhancing the immune response, under chronic stress exposure. Our results shed light on the miRNA expression profiles in the pituitary and adrenal cortex with and without chronic stress exposure, and provide a new insight into miR-205 with its feasible role in regulating chronic stress in the pituitary and hippocampus through targeting MMD.