Search

Abstract Background N6-Methyladenosine (m6A) methylation, a common form of RNA modification, play an important role in the pathogenesis of various diseases and in the ontogeny of organisms. Nevertheless, the precise function of m6A methylation in photoaging remains unknown. Objectives This study aims to investigate the biological role and underlying mechanism of m6A methylation in photoaging. Methods m6A dot blot, Real-time quantitative PCR (RT-qPCR), western blot and immunohistochemical (IHC) assays were employed to detect the m6A level and specific m6A methylase in ultraviolet ray (UVR)-induced photoaging tissue. The profile of m6A-tagged mRNA was identified by methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA-seq analysis. Finally, we investigated the regulatory mechanism of KIAA1429 by MeRIP-qPCR, RNA knockdown and immunofluorescence assay. Results m6A levels were increased in photoaging and were closely associated with the upregulation of KIAA1429 expression. 1331 differentially m6A methylated genes were identified in the UVR group compared with the control group, of which 1192 (90%) were hypermethylated. Gene ontology analysis showed that genes with m6A hypermethylation and mRNA downregulation were mainly involved in extracellular matrix metabolism and collagen metabolism-related processes. Furthermore, KIAA1429 knockdown abolished the downregulation of TGF-bRII and upregulation of MMP1 in UVR-irradiated human dermal fibroblasts (HDFs). Mechanically, we identified MFAP4 as a target of KIAA1429-mediated m6A modification and KIAA1429 might suppress collagen synthesis through an m6A-MFAP4-mediated process. Conclusions The increased expression of KIAA1429 hinders collagen synthesis during UVR-induced photoaging, suggesting that KIAA1429 represents a potential candidate for targeted therapy to mitigate UVR-driven photoaging.

Abstract Background Cyp19a1a is a key enzyme in the pathway that converts androgens into estrogen and is regulated by TGF-β signaling. Smad4 and FoxH1 are downstream effectors of TGF-β signaling and may play important roles in ovarian development in M. albus. Methods We investigated the expression pattern of the Smad4 and FoxH1 using qRT‒PCR and immunofluorescence, then tested the changes of smad4 and foxh1 by qRT‒PCR after ovary incubation with FSH in vitro, and analysed the regulation of cyp19a1a transcription by Smad4 and FoxH1 by dual-luciferase reporter assays. Results We found that Smad4 encoded a putative protein of 449 amino acids and harbored the three conserved domains typical of this protein family. Smad4 and foxh1 exhibited similar expression patterns during ovarian development and after FSH incubation, with Pearson’s coefficients of 0.873 and 0.63–0.81, respectively. Furthermore, Smad4, FoxH1 and Cyp19a1a colocalized in the granulosa cells and theca cells of ovaries during the mid-to-late vitellogenic stage. Smad4 repressed cyp19a1a activity via SBE1 (− 1372/−1364) and SBE2 (− 415/−407) in the cyp19a1a promoter, whereas mutating SBE1 or SBE2 restored cyp19a1a promoter activity. Co-overexpression of Smad4 and FoxH1 significantly reduced cyp19a1a promoter activity. Conclusions This study provides new insights into the potential functions of transcription factors Smad4 and FoxH1 in ovarian development and the transcriptional regulation mechanism of cyp19a1a in M. albus, which will reveal Smad4/FoxH1-mediated TGF-β signaling in reproduction and the regulation of the cyp19a1a. Plain English summary Aromatase, encoded by cyp19a1a, is involved in ovarian development and plays an important role in the quality of eggs, as well the sex ratio, of the teleost fish, M. albus. The research on the transcriptional regulation of cyp19a1a has contributed to the understanding of its role in ovarian development. In previous study, it was shown that FoxH1 inhibits cyp19a1a transcription. In the present study, Smad4 was confirmed as a cyp19a1a transcriptional repressor and Smad4 may also coordinate with FoxH1 to repress cyp19a1a transcription. At present, we provide a new perspective for the transcriptional regulation of cyp19a1a by transcription factors Smad4 and FoxH1 in teleost fish ovary. In the future, the regulatory networks of Smad4 and FoxH1 will be further studied and the gene editing technology will be applied to screen specific regulatory factors of cyp191a1a gene, so as to alter the female cycle and modulate the sex ratio of the eggs production.

Abstract Background The ricefield eel Monopterus albus undergoes a natural sex change from female to male during its life cycle, and previous studies have shown the potential mechanisms of this transition at the transcriptional and protein levels. However, the changes in protein levels have not been fully explored, especially in the intersexual stage. Results In the present study, the protein expression patterns in the gonadal tissues from five different periods, the ovary (OV), early intersexual stage gonad (IE), middle intersexual stage gonad (IM), late intersexual stage gonad (IL), and testis (TE), were determined by untargeted proteomics sequencing. A total of 5125 proteins and 394 differentially expressed proteins (DEPs) were detected in the gonadal tissues. Of the 394 DEPs, there were 136 between the OV and IE groups, 20 between the IM and IE groups, 179 between the IL and IM groups, and 59 between the TE and IL groups. Three candidate proteins, insulin-like growth factor 2 mRNA-binding protein 3 isoform X1 (Igf2bp3), triosephosphate isomerase (Tpi), and Cu-Zn superoxide dismutase isoform X1 [(Cu-Zn) Sod1], were validated by western blotting to verify the reliability of the data. Furthermore, metal metabolite-related proteins were enriched in the IL vs. IM groups and TE vs. IL groups, which had close relationships with sex change, including Cu2+-, Ca2+-, Zn2+- and Fe2+/Fe3+-related proteins. Analysis of the combined transcriptome data revealed consistent protein/mRNA expression trends for two metal metabolite-related proteins/genes [LOC109953912 and calcium Binding Protein 39 Like (cab39l)]. Notably, we detected significantly higher levels of Cu2+ during the sex change process, suggesting that Cu2+ is a male-related metal metabolite that may have an important function in male reproductive development. Conclusions In summary, we analyzed the protein profiles of ricefield eel gonadal tissues in five sexual stages (OV, IE, IM, IL, and TE) and verified the plausibility of the data. After preforming the functional enrichment of metal metabolite-related DEPs, we detected the contents of the metal metabolites Zn2+, Cu2+, Ca2+, and Fe2+/Fe3+ at these five stages and screened for (Cu-Zn) Sod1 and Mmp-9 as possible key proteins in the sex reversal process.

Abstract Background Schizothorax o’connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o’connori using the whole-genome resequencing method. Results The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o’connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism. Conclusion Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o’connori. There was no obvious genetic differentiation at the genome level between the S. o’connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o’connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Abstract Background Blastocystis is a common intestinal protist with a wide range of hosts. Thus far, 38 subtypes have been identified. In recent years, wild animals have been confronted with habitat fragmentation as well as an increasing risk of zoonotic disease transmission due to human disturbance. Only limited data are available on Blastocystis infection and subtype distribution in wild rhesus macaques in China. The aim of the present study was to investigate the prevalence and genetic diversity of Blastocystis in nine wild rhesus macaque populations in China. Methods A total of 276 faecal samples were collected from five high-altitude populations (high-altitude group [HAG]; 2800–4100 m a.s.l.) and four low-altitude populations (low-altitude group [LAG]; 5–1000 m a.s.l) of rhesus macaques. PCR-based analysis, using a new primer pair for the amplification of a 1690-bp sequence of the small subunit ribosomal RNA (SSU rRNA) gene, was used for prevalence and genetic diversity analysis. Result Analysis of faecal samples revealed that Blastocystis infection was common in rhesus macaques, with an infection positivity rate of 80.1% (n = 221/276 samples). There was no significant difference (P = 0.121) in positivity rate between the LAG (84.3%) and HAG (76.8%). Overall, 33 haplotypes were obtained and classified into four subtypes (STs), of which three were potentially zoonotic subtypes (ST1, 29.7%; ST2, 16.7%; ST3, 31.9%) and one that was first identified in this study and named ST39 (12.0%). The STs were distributed differently among the rhesus macaque populations, except for ST3, which was found in all populations. Phylogenetic analyses revealed two major divergent clades of ST3 for the HAG and LAG. Genetic diversity analysis showed a high genetic diversity of ST3 (haplotype diversity: 0.846; nucleotide diversity: 0.014) in the rhesus macaques, but a high genetic differentiation (FST > 0.25) and a low gene flow (Nm = 0.09) between the HAG and LAG. Conclusion Our study, which is the first investigation on Blastocystis infection in multiple wild rhesus macaque populations in China, indicates a potential risk of zoonotic transmission of Blastocystis in the study areas. Blastocystis ST3 showed high genetic diversity in wild rhesus macaques and significant genetic differentiation between the HAG and LAG. Our results provide fundamental information on the genetic diversity and prevalence of Blastocystis in wild rhesus macaque populations. Graphical Abstract

Abstract Background The expression and biological functions of circular RNAs (circRNAs) in reproductive organs have been extensively reported. However, it is still unclear whether circRNAs are involved in sex change. To this end, RNA sequencing (RNA-seq) was performed in gonads at 5 sexual stages (ovary, early intersexual stage gonad, middle intersexual stage gonad, late intersexual stage gonad, and testis) of ricefield eel, and the expression profiles and potential functions of circRNAs were studied. Results Seven hundred twenty-one circRNAs were identified, and the expression levels of 10 circRNAs were verified by quantitative real-time PCR (qRT–PCR) and found to be in accordance with the RNA-seq data, suggesting that the RNA-seq data were reliable. Then, the sequence length, category, sequence composition and the relationship between the parent genes of the circRNAs were explored. A total of 147 circRNAs were differentially expressed in the sex change process, and GO and KEGG analyses revealed that some differentially expressed (such as novel_circ_0000659, novel_circ_0004005 and novel_circ_0005865) circRNAs were closely involved in sex change. Furthermore, expression pattern analysis demonstrated that both circSnd1 and foxl2 were downregulated in the process of sex change, which was contrary to mal-miR-135b. Finally, dual-luciferase reporter assay and RNA immunoprecipitation showed that circSnd1 and foxl2 can combine with mal-miR-135b and mal-miR-135c. These data revealed that circSnd1 regulates foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. Conclusion Our results are the first to demonstrate that circRNAs have potential effects on sex change in ricefield eel; and circSnd1 could regulate foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. These data will be useful for enhancing our understanding of sequential hermaphroditism and sex change in ricefield eel or other teleosts.

Abstract Background An increasing number of long noncoding RNAs (lncRNAs) have been found to play important roles in sex differentiation and gonad development by regulating gene expression at the epigenetic, transcriptional and posttranscriptional levels. The ricefield eel, Monopterus albus, is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. However, the roles of lncRNA in the sex change is unclear. Results Herein, we performed RNA sequencing to analyse lncRNA expression patterns in five different stages of M. albus development to investigate the roles of lncRNAs in the sex change process. A total of 12,746 lncRNAs (1503 known lncRNAs and 11,243 new lncRNAs) and 2901 differentially expressed lncRNAs (DE-lncRNAs) were identified in the gonads. The target genes of the DE-lncRNAs included foxo1, foxm1, smad3, foxr1, camk4, ar and tgfb3, which were mainly enriched in signalling pathways related to gonadal development, such as the insulin signalling pathway, MAPK signalling pathway, and calcium signalling pathway. We selected 5 highly expressed DE-lncRNAs (LOC109952131, LOC109953466, LOC109954337, LOC109954360 and LOC109958454) for full length amplification and expression pattern verification. They were all expressed at higher levels in ovaries and intersex gonads than in testes, and exhibited specific time-dependent expression in ovarian tissue incubated with follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG). The results of quantitative real-time PCR (qRT-PCR) analysis and a dual-luciferase assay showed that znf207, as the gene targeted by LOC109958454, was expressed in multiple tissues and gonadal developmental stages of M. albus, and its expression was also inhibited by the hormones FSH and hCG. Conclusions These results provide new insights into the role of lncRNAs in gonad development, especially regarding natural sex changes in fish, which will be useful for enhancing our understanding of sequential hermaphroditism and sex changes in the ricefield eel (M. albus) and other teleosts.

Abstract Background The mammal intestinal microbiota is involved in various physiological processes and plays a key role in host environment adaption. However, for non-human primates (NHPs), little is known about their gut microbial community in high-altitude environments and even less about their adaption to such habitats. We characterised the gut microbial community of rhesus macaques from multiple high-altitude environments and compared it to those of low-altitude populations. Results We collected faecal samples of rhesus macaques from four high-altitude populations (above 3000 m) and three low-altitude populations (below 500 m). By calculating the alpha diversity index, we found that high-altitude populations exhibited a higher diversity. Statistical analysis of beta diversity indicated significant differences between high- and low-altitude populations. Significant differences were also detected at the phylum and family levels. At the phylum level, the high-altitude gut microbial community was dominated by Firmicutes (63.42%), while at low altitudes, it was dominated by Bacteroidetes (47.4%). At the family level, the high-altitude population was dominated by Ruminococcaceae (36.2%), while the low-altitude one was dominated by Prevotellaceae (39.6%). Some families, such as Christensenellaceae and Rikenellaceae, were consistently higher abundant in all high-altitude populations. We analysed the overlap of operational taxonomic units (OTUs) in high-altitude populations and determined their core OTUs (shared by all four high-altitude populations). However, when compared with the low-altitude core OTUs, only 65% were shared, suggesting a divergence in core OTUs. Function prediction indicated a significant difference in gene copy number of 35 level-2 pathways between high- and low-altitude populations; 29 of them were higher in high altitudes, especially in membrane transport and carbohydrate metabolism. Conclusions The gut microbial community of high-altitude rhesus macaques was significantly distinct from that of low-altitude populations in terms of diversity, composition and function. High-altitude populations were dominated by Firmicutes and Ruminococcace, while in low-altitude populations, Bacteroidetes and Prevotellaceae were dominant. The difference in gut microbiota between these two populations may be caused by differences in host diet, environmental temperature and oxygen pressure. These differentiated gut microbial microorganisms may play a critical role in the adaptive evolution of rhesus macaques to high-altitude environments.