Age-related gene expression and DNA methylation changes in rhesus macaque

Age-associated cortical similarity networks correlate with cell type-specific transcriptional signatures

Human brain structure shows heterogeneous patterns of change across adults aging and is associated with cognition. However, the relationship between cortical structural changes during aging and gene transcription signatures remains unclear. Here, using structural magnetic resonance imaging data of two separate cohorts of healthy participants from the Cambridge Centre for Aging and Neuroscience (n = 454, 18-87 years) and Dallas Lifespan Brain Study (n = 304, 20-89 years) and a transcriptome dataset, we investigated the link between cortical morphometric similarity network and brain-wide gene transcription. In two cohorts, we found reproducible morphometric similarity network change patterns of decreased morphological similarity with age in cognitive related areas (mainly located in superior frontal and temporal cortices), and increased morphological similarity in sensorimotor related areas (postcentral and lateral occipital cortices). Changes in morphometric similarity network showed significant spatial correlation with the expression of age-related genes that enriched to synaptic-related biological processes, synaptic abnormalities likely accounting for cognitive decline. Transcription changes in astrocytes, microglia, and neuronal cells interpreted most of the age-related morphometric similarity network changes, which suggest potential intervention and therapeutic targets for cognitive decline. Taken together, by linking gene transcription signatures to cortical morphometric similarity network, our findings might provide molecular and cellular substrates for cortical structural changes related to cognitive decline across adults aging.

Gene expressions between obligate bamboo-eating pandas and non-herbivorous mammals reveal converged specialized bamboo diet adaptation

BACKGROUND

It is inevitable to change the function or expression of genes during the environmental adaption of species. Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to Carnivora and have developed similar adaptations to the same dietary switch to bamboos at the morphological and genomic levels. However, the genetic adaptation at the gene expression level is unclear. Therefore, we aimed to examine the gene expression patterns of giant and red panda convergent specialized bamboo-diets. We examined differences in liver and pancreas transcriptomes between the two panda species and other non-herbivorous species.

RESULTS

The clustering and PCA plots suggested that the specialized bamboo diet may drive similar expression shifts in these two species of pandas. Therefore, we focused on shared liver and pancreas DEGs (differentially expressed genes) in the giant and red panda relative to other non-herbivorous species. Genetic convergence occurred at multiple levels spanning carbohydrate metabolism, lipid metabolism, and lysine degradation. The shared adaptive convergence DEGs in both organs probably be an evolutionary response to the high carbohydrate, low lipid and lysine bamboo diet. Convergent expression of those nutrient metabolism-related genes in both pandas was an intricate process and subjected to multi-level regulation, including DNA methylation and transcription factor. A large number of lysine degradation and lipid metabolism related genes were hypermethylated in promoter regions in the red panda. Most genes related to carbohydrate metabolism had reduced DNA methylation with increased mRNA expression in giant pandas. Unlike the red panda, the core gene of the lysine degradation pathway (AASS) doesn't exhibit hypermethylation modification in the giant panda, and dual-luciferase reporter assay showed that transcription factor, NR3C1, functions as a transcriptional activator in AASS transcription through the binding to AASS promoter region.

CONCLUSIONS

Our results revealed the adaptive expressions and regulations of the metabolism-related genes responding to the unique nutrients in bamboo food and provided data accumulation and research hints for the future revelation of complex mechanism of two pandas underlying convergent adaptation to a specialized bamboo diet.

Skeletal muscle transcriptome is affected by age in severely burned mice

Severe burn results in muscle wasting affecting quality of life in both children and adults. Biologic metabolic profiles are noticeably distinctive in childhood. We posit that muscle gene expression profiles are differentially regulated in response to severe burns in young animals. Twelve C57BL6 male mice, including young (5 weeks-old) and adults (11 weeks-old), received either scald burn, or sham procedure. Mouse muscle tissue was harvested 24 h later for Next Generation Sequence analysis. Our results showed 662 downregulated and 450 upregulated genes in gastrocnemius of young mice compared to adults without injury. After injury, we found 74/75 downregulated genes and 107/128 upregulated genes in both burned groups compared to respective uninjured age groups. VEGFA-VEGFR2, focal adhesion, and nuclear receptor meta-pathways were the top 3 gene pathways undergoing a differential change in response to age. Of note, the proteasome degradation pathway showed the most similar changes in both adult and young burned animals. This study demonstrates the characteristic profile of gene expression in skeletal muscle in young and adult burned mice. Prominent age effects were revealed in transcriptional levels with increased alterations of genes, miRNAs, pathways, and interactions.

Genome-Wide DNA Methylation Profile Indicates Potential Epigenetic Regulation of Aging in the Rhesus Macaque Thymus

We analyzed whole-genome bisulfite sequencing (WGBS) and RNA sequencing data of two young (1 year old) and two adult (9 years old) rhesus macaques (Macaca mulatta) to characterize the genomic DNA methylation profile of the thymus and explore the molecular mechanism of age-related changes in the thymus. Combining the two-omics data, we identified correlations between DNA methylation and gene expression and found that DNA methylation played an essential role in the functional changes of the aging thymus, especially in immunity and coagulation. The hypomethylation levels of C3 and C5AR2 and the hypermethylation level of C7 may lead to the high expressions of these genes in adult rhesus macaque thymuses, thus activating the classical complement pathway and the alternative pathway and enhancing their innate immune function. Adult thymuses had an enhanced coagulation pathway, which may have resulted from the hypomethylation and upregulated expressions of seven coagulation-promoting factor genes (F13A1, CLEC4D, CLEC4E, FCN3, PDGFRA, FGF2 and FGF7) and the hypomethylation and low expression of CPB2 to inhibit the degradation of blood clots. Furthermore, the functional decline in differentiation, activation and maturation of T cells in adult thymuses was also closely related to the changes in methylation levels and gene expression levels of T cell development genes (CD3G, GAD2, ADAMDEC1 and LCK) and the thymogenic hormone gene TMPO. A comparison of the age-related methylated genes among four mammal species revealed that most of the epigenetic clocks were species-specific. Furthermore, based on the genomic landscape of allele-specific DNA methylation, we identified several age-related clustered sequence-dependent allele-specific DNA methylated (cS-ASM) genes. Overall, these DNA methylation patterns may also help to assist with understanding the mechanisms of the aging thymus with the epigenome.

Multi-omics analysis reveals the host-microbe interactions in aged rhesus macaques

Aging is a complex multifactorial process that greatly affects animal health. Multi-omics analysis is widely applied in evolutionary biology and biomedical research. However, whether multi-omics can provide sufficient information to reveal comprehensive changes in aged non-human primates remains unclear. Here, we explored changes in host-microbe interactions with aging in Chinese rhesus macaques (Macaca mulatta lasiota, CRs) using multi-omics analysis. Results showed marked changes in the oral and gut microbiomes between young and aged CRs, including significantly reduced probiotic abundance and increased pathogenic bacterial abundance in aged CRs. Notably, the abundance of Lactobacillus, which can metabolize tryptophan to produce aryl hydrocarbon receptor (AhR) ligands, was decreased in aged CRs. Consistently, metabolomics detected a decrease in the plasma levels of AhR ligands. In addition, free fatty acid, acyl carnitine, heparin, 2-(4-hydroxyphenyl) propionic acid, and docosahexaenoic acid ethyl ester levels were increased in aged CRs, which may contribute to abnormal fatty acid metabolism and cardiovascular disease. Transcriptome analysis identified changes in the expression of genes associated with tryptophan metabolism and inflammation. In conclusion, many potential links among different omics were found, suggesting that aged CRs face multiple metabolic problems, immunological disorders, and oral and gut diseases. We determined that tryptophan metabolism is critical for the physiological health of aged CRs. Our findings demonstrate the value of multi-omics analyses in revealing host-microbe interactions in non-human primates and suggest that similar approaches could be applied in evolutionary and ecological research of other species.

Comparative Transcriptomics and Methylomics Reveal Adaptive Responses of Digestive and Metabolic Genes to Dietary Shift in Giant and Red Pandas

Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to the order Carnivora, but have changed their dietary habits to eating bamboo exclusively. The convergent evolution characteristics of their morphology, genome and gut flora have been found in the two pandas. However, the research on the convergent adaptation of their digestion and metabolism to the bamboo diet, mediated by the dietary shift of the two pandas at the gene-expression and epigenetic regulation levels, is still lacking. We therefore used RNA sequencing among five species (two pandas and three non-herbivore mammals) and bisulfite sequencing among three species (two pandas and a carnivore ferret) to sequence key digestion and metabolism tissues (stomach and small intestine). Our results provide evidence that the convergent differentially expressed genes (related to carbohydrate utilization, bile secretion, Lys and Arg metabolism, vitamin B12 utilization and cyanide detoxification) of the two pandas are adaptive responses to the bamboo diet containing low lipids, low Lys and Arg, low vitamin B12 and high cyanide. We also profiled the genome-wide methylome maps of giant panda, red panda and ferret, and the results indicated that the promoter methylation of the two pandas may regulate digestive and metabolic genes to adapt to sudden environmental changes, and then, transmit genetic information to future generations to evolve into bamboo eaters. Taken together, our study provides new insights into the molecular mechanisms of the dietary shift and the adaptation to a strict bamboo diet in both pandas using comparative transcriptomics and methylomics.

Epigenomic profiling indicates a role for DNA methylation in the postnatal liver and pancreas development of giant pandas

Giant pandas are unique within Carnivora with a strict bamboo diet. Here, the epigenomic profiles of giant panda liver and pancreas tissues collected from three important feeding stages were investigated using BS-seq. Few differences in DNA methylation profiles were exhibited between no feeding and suckling groups in both tissues. However, we observed a tendency toward a global loss of DNA methylation in the gene-body and promoter region of metabolism-related genes from newborn to adult. Correlation analysis revealed a significant negative correlation between the changes in methylation levels within gene promoters and gene expression. The majority of genes related to nutrition metabolism had lost DNA methylation with increased mRNA expression in adult giant pandas. The few galactose metabolism and unsaturated fatty acid metabolism related genes that were hypomethylated and highly-expressed at early stages of giant panda development may meet the nutritional requirement of this species' highly altricial neonates.

Obesity-Associated Differentially Methylated Regions in Colon Cancer

Obesity with adiposity is a common disorder in modern days, influenced by environmental factors such as eating and lifestyle habits and affecting the epigenetics of adipose-based gene regulations and metabolic pathways in colorectal cancer (CRC). We compared epigenetic changes of differentially methylated regions (DMR) of genes in colon tissues of 225 colon cancer cases (154 non-obese and 71 obese) and 15 healthy non-obese controls by accessing The Cancer Genome Atlas (TCGA) data. We applied machine-learning-based analytics including generalized regression (GR) as a confirmatory validation model to identify the factors that could contribute to DMRs impacting colon cancer to enhance prediction accuracy. We found that age was a significant predictor in obese cancer patients, both alone (p = 0.003) and interacting with hypomethylated DMRs of ZBTB46, a tumor suppressor gene (p = 0.008). DMRs of three additional genes: HIST1H3I (p = 0.001), an oncogene with a hypomethylated DMR in the promoter region; SRGAP2C (p = 0.006), a tumor suppressor gene with a hypermethylated DMR in the promoter region; and NFATC4 (p = 0.006), an adipocyte differentiating oncogene with a hypermethylated DMR in an intron region, are also significant predictors of cancer in obese patients, independent of age. The genes affected by these DMR could be potential novel biomarkers of colon cancer in obese patients for cancer prevention and progression.

Comparative transcriptomes of three different skin sites for the Asiatic toad (Bufo gargarizans)

Toads release toxic dry secretions from glands in their skin. Toxin possesses a wide range of biological effects, but little is known about its specific gene expression pattern and regulatory mechanisms. The Asiatic toad (Bufo gargarizans) is widely used to produce toxin. Here, we explored the gene expression of 30 tissue samples from three different skin sites (parotoid gland, dorsal skin, and abdomen skin) of B. gargarizans. After de novo assembly, 783,130 unigenes with an average length of 489 bp (N50 = 556 bp) were obtained. A total of 9,248 significant differentially expressed genes (DEGs) were detected. There were 8,819 DEGs between the parotoid gland and abdomen skin and 1,299 DEGs between the dorsal skin and abdomen skin, while only 1,283 DEGs were obtained between the parotoid gland and dorsal skin. Through enrichment analysis, it was found that the detected differential gene expressions corresponded to the different functions of different skin sites. Our key findings were the genetic expression of toxin secretion, the protection function of skin, and the related genes such as HSD3B, Cyp2c, and CAT, LGALS9. In conclusion, we provide useful transcript resources to study the gene expression and gene function of B. gargarizans and other amphibians. The detected DEGs between different sites of the skin provided better insights into the genetic mechanisms of toxin secretion and the protection function of skin for amphibians.

Transcriptomic landscape of persistent diarrhoea in rhesus macaques and comparison with humans and mouse models with inflammatory bowel disease

Diarrhoea is a widespread disease in captive rhesus macaques (Macaca mulatta) and a small proportion of individuals may experience persistent diarrhoea. Persistent diarrhoea can lead to a compromised immune system, intestinal inflammation and malnutrition. We analyzed the blood transcriptomes of 10 persistent diarrhoeal and 12 healthy rhesus macaques to investigate the gene expression differences between the two groups. We identified 330 DEGs between persistent diarrhoeal and healthy rhesus macaques. The 211 up-regulated DEGs in the diarrhoeal group were mainly enriched in immune-related and interleukin-related categories. Among them, three interleukin (IL) 18 related DEGs (IL18, IL18R1, and IL18BP) played important roles in actively regulating pro-inflammatory responses. Interestingly, the up- and down-regulated DEGs were both enriched in the same immune-related categories. Thus, we applied a new method to examine the distribution of DEGs in all child categories. We found that interleukin and T cell related categories were mainly occupied by up-regulated DEGs, while immunoglobulin production and B cell related categories were enriched by down-regulated DEGs. We also compared rhesus macaque DEGs with the DEGs of inflammatory bowel disease (IBD) humans and IBD mouse models and found that 30-40% of macaque DEGs were shared with IBD humans and mouse models. In conclusion, our results showed that there were significant immune differences between persistent diarrhoeal rhesus macaques and healthy macaques, which was similar to the expression differences in IBD patients and mouse models.

Sex-specific gene expression in the blood of four primates

Blood is an important non-reproductive tissue, but little is known about the sex-specific gene expressions in the blood. Therefore, we investigated sex-specific gene expression differences in the blood tissues of four primates, rhesus macaques (Macaca mulatta), Tibetan macaques (M. thibetana), yellow baboons (Papio cynocephalus), and humans. We identified seven sex-specific differentially expressed genes (SDEGs) in each non-human primate and 31 SDEGs in humans. The four primates had only one common SDEG, MAP7D2. In humans, immune-related SDEGs were identified as up-regulated, but also down-regulated in females. We also found that most of the X-Y gene pairs had similar expression levels between species, except pair EIF1AY/EIF1AX. The expression level of X-Y gene pairs of rhesus and Tibetan macaques showed no significant differential expression levels, while humans had six significant XY-biased and three XX-biased X-Y gene pairs. Our observed sex differences in blood should increase understanding of sex differences in primate blood tissue.

Mining of RNA Methylation-Related Genes and Elucidation of Their Molecular Biology in Gallbladder Carcinoma

Gallbladder carcinoma (GBC), which has high invasion and metastasis risks, remains the most common biliary tract malignancy. Surgical resection for GBC is the only effective treatment, but most patients miss the opportunity for curative surgery because of a lack of timely diagnosis. The aim of this study was to identify and verify early candidate diagnostic and prognostic RNA methylation related genes for GBC via integrated transcriptome bioinformatics analysis. Lists of GBC-related genes and methylation-related genes were collected from public databases to screen differentially expressed genes (DEGs) by using the limma package and the RobustRankAggreg (RRA) package. The core genes were collected with batch effects corrected by the RRA algorithm through protein interaction network analysis, signaling pathway enrichment analysis and gene ranking. Four modules obtained from four public microarray datasets were found to be related to GBC, and FGA, F2, HAO1, CFH, PIPOX, ITIH4, GNMT, MAT1A, MTHFD1, HPX, CTH, EPHX2, HSD17B6, AKR1C4, CFHR3, ENNP1, and NAT2 were revealed to be potential hub genes involved in methylation-related pathways and bile metabolism-related pathways. Among these, FGA, CFH, F2, HPX, and PIPOX were predicted to be methylated genes in GBC, but POPIX had no modification sites for RNA methylation. Furthermore, survival analysis of TCGA (the Cancer Genome Atlas) database showed that six genes among the hub genes, FGA, CFH, ENPP1, CFHR3, ITIH4, and NAT2, were highly expressed and significantly correlated with worse prognosis. Gene correlation analysis revealed that the FGA was positively correlated with the ENPP1, NAT2, and CFHR3, while CFH was positively correlated with the NAT2, CFHR3, and FGA. In addition, the results of immunohistochemistry (IHC) showed that the expressions of FGA, F2, CFH, PIPOX, ITIH4, GNMT, MAT1A, MTHFD1, HPX, CFHR3, NAT2, and ENPP1 were higher in GBC tissues than that in control tissues. In conclusion, two genes, FGA and CFH, were identified as RNA methylation-related genes also involved in bile metabolism in GBC, which may be novel biomarkers to early diagnose and evaluate prognosis for GBC.