Search

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.

Aging is a gradual process of natural change that occurs after reaching sexual maturity. It is also a known risk factor for many chronic diseases. Recent research has shown that senolytics can extend the lifespans and health spans of model organisms, and they have also been demonstrated effective in treating age-related diseases. In this study, we conducted a high-throughput screening of 156 drugs that targeted the PI3K/AKT/mTOR pathway to identify potential senolytic medications. Among these drugs, PF-04691502 was selected for further investigation to understand its molecular mechanism of action. Our findings indicate that PF-04691502, a dual inhibitor of PI3K/AKT and mTOR, specifically eliminates senescent cells. It reduces the expression levels of key markers of cellular senescence, such as SA-β-Gal, senescence-associated secretory phenotypes (SASPs) and p16INK4a. Additionally, PF-04691502 inhibits the phosphorylation of S6K and AKT, leading to the apoptosis of senescent cells. These results suggest that PF-04691502 holds promise as a new senolytic drug. This paper provides important insights into the potential application of PF-04691502 in the study of cell senescence.

The geochemical characteristics of a 2.1 m BBW25 core, collected from the Beibu Gulf, have been investigated in terms of the major and trace elements, organic matter, and CaCO3 and AMS 14C dating by XRF, ICP-OES, ICP-MS, and more. We have found through previous research that there are issues with unclear delineation of sedimentary evolution environments and inexact responses between chemical weathering intensity and major paleoclimate events in the Beibu Gulf. The AMS 14C dating results indicate that the sedimentary age at the bottom was 19.24 ky b.p. CaCO3, δ13C, C/N, and Sr/Ba indexes show a sedimentary environment change from terrestrial to marine environments and a “jump” of ~4000 years in continent–ocean changes. The evolution of the sedimentary environment of Beibu Gulf was divided into three environments and five sub-environments. The changes in chemical weathering intensity indicators recorded by the CIX and the Fe/Al ratio respond well to the East Asian monsoon cycle, the meltwater events, and the alternation of cold and warm events. This study explains the chemical weathering intensity and sedimentary environment in the BBW25 core by geochemical characteristics and further reveals the paleoenvironmental characteristics and possible driving mechanisms over the past ~20,000 years.

Abstract There have been many reports on the genetic mechanism in rhesus macaques (RMs) for environmental adaptation to high altitudes, but the synergistic involvement of gut microbiota in this adaptation remains unclear. Here we performed fecal metagenomic and metabolomic studies on samples from high- and low-altitude populations to assess the synergistic role of gut microbiota in the adaptation of RMs to high-altitude environments. Microbiota taxonomic annotation yielded 7471 microbiota species. There were 37 bacterial species whose abundance was significantly enriched in the high-altitude populations, 16 of which were previously reported to be related to the host’s dietary digestion and energy metabolism. Further functional gene enrichment found a stronger potential for gut microbiota to synthesize energy substrate acetyl-CoA using CO2 and energy substrate pyruvate using oxaloacetate, as well as a stronger potential to transform acetyl-CoA to energy substrate acetate in high-altitude populations. Interestingly, there were no apparent differences between low-altitude and high-altitude populations in terms of genes enriched in the main pathways by which the microbiota consumed the three energy substrates, and none of the three energy substrates were detected in the fecal metabolites. These results strongly suggest that gut microbiota plays an important energy compensatory role that helps RMs to adapt to high-altitude environments. Further functional enrichment after metabolite source analysis indicated the abundance of metabolites related to the degradation of exogenous toxins was also significantly higher in high-altitude populations, which suggested a contributory role of gut microbiota to the degradation of exogenous toxins in wild RMs adapted to high-altitude environments.

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.

Enterocytozoon bieneusi is a zoonotic pathogen with a wide range of animal host. There are only few reports of E. bieneusi infection in wild Chinese rhesus macaques. Here, we determined the prevalence of E. bieneusi in nine wild rhesus macaque populations and assessed their zoonotic potential by performed genotype of ITS gene. A total of 324 fecal samples of rhesus macaque were collected in nine geographical populations from five Chinese provinces (Sichuan, Chongqing, Qinghai, Tibet and Hainan). 38 samples (11.7%) were found to be infected with E. bieneusi, and 11 genotypes were identified including three known genotypes (D, EbpC and SCC-2) and eight novel genotypes named Mul6∼13. Genotype D (63.2%) was the most prevalent, being observed in seven populations except XZ-2 and QH, and other genotypes were identified only in a single area. According to the phylogenetic analysis, Mul6∼9, Mul11∼13 and zoonotic genotype D were clustered into Group 1, indicating that these genotypes may be potentially zoonotic. Among nine populations, population SC-3 had the highest infection rate (26.3%), and the lowest was the wild QH population without infection, but the difference of infection rate among the nine populations is not significant. It is concluded that, rhesus macaque populations are generally infected E. bieneusi in many areas of China, and there may be a risk of cross infection with E. bieneusi in some areas found having zoonotic genotypes, and these areas should be paid more attention to prevent.

ABSTRACT Trichosporon asahii is an emerging opportunistic pathogen that causes potentially fatal disseminated trichosporonosis. The global prevalence of coronavirus disease 2019 (COVID-19) poses an increasing fungal infection burden caused by T. asahii. Allicin is the main biologically active component with broad-spectrum antimicrobial activity in garlic. In this study, we performed an in-depth analysis of the antifungal characteristics of allicin against T. asahii based on physiological, cytological, and transcriptomic assessments. In vitro, allicin inhibited the growth of T. asahii planktonic cells and biofilm cells significantly. In vivo, allicin improved the mean survival time of mice with systemic trichosporonosis and reduced tissue fungal burden. Electron microscopy observations clearly demonstrated damage to T. asahii cell morphology and ultrastructure caused by allicin. Furthermore, allicin increased intracellular reactive oxygen species (ROS) accumulation, leading to oxidative stress damage in T. asahii cells. Transcriptome analysis showed that allicin treatment disturbed the biosynthesis of cell membrane and cell wall, glucose catabolism, and oxidative stress. The overexpression of multiple antioxidant enzymes and transporters may also place an additional burden on cells, causing them to collapse. Our findings shed new light on the potential of allicin as an alternative treatment strategy for trichosporonosis. IMPORTANCE Systemic infection caused by T. asahii has recently been recognized as an important cause of mortality in hospitalized COVID-19 patients. Invasive trichosporonosis remains a significant challenge for clinicians, due to the limited therapeutic options. The present work suggests that allicin holds great potential as a therapeutic candidate for T. asahii infection. Allicin demonstrated potent in vitro antifungal activity and potential in vivo protective effects. In addition, transcriptome sequencing provided valuable insights into the antifungal effects of allicin.

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.

IntroductionFull-thickness skin wound healing remains a serious undertaking for patients. While stem cell-derived exosomes have been proposed as a potential therapeutic approach, the underlying mechanism of action has yet to be fully elucidated. The current study aimed to investigate the impact of exosomes derived from human umbilical cord mesenchymal stem cells (hucMSC-Exosomes) on the single-cell transcriptome of neutrophils and macrophages in the context of wound healing.MethodsUtilizing single-cell RNA sequencing, the transcriptomic diversity of neutrophils and macrophages was analyzed in order to predict the cellular fate of these immune cells under the influence of hucMSC-Exosomes and to identify alterations of ligand-receptor interactions that may influence the wound microenvironment. The validity of the findings obtained from this analysis was subsequently corroborated by immunofluorescence, ELISA, and qRT-PCR. Neutrophil origins were characterized based on RNA velocity profiles.ResultsThe expression of RETNLG and SLC2A3 was associated with migrating neutrophils, while BCL2A1B was linked to proliferating neutrophils. The hucMSC-Exosomes group exhibited significantly higher levels of M1 macrophages (215 vs 76, p < 0.00001), M2 macrophages (1231 vs 670, p < 0.00001), and neutrophils (930 vs 157, p < 0.00001) when compared to control group. Additionally, it was observed that hucMSC-Exosomes elicit alterations in the differentiation trajectories of macrophages towards more anti-inflammatory phenotypes, concomitant with changes in ligand-receptor interactions, thereby facilitating healing.DiscussionThis study has revealed the transcriptomic heterogeneity of neutrophils and macrophages in the context of skin wound repair following hucMSC-Exosomes interventions, providing a deeper understanding of cellular responses to hucMSC-Exosomes, a rising target of wound healing intervention.

As Drosophila is an extensively used genetic model system, understanding of its regulatory networks has great significance in revealing the genetic mechanisms of ageing and human diseases. Competing endogenous RNA (ceRNA)-mediated regulation is an important mechanism by which circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) regulate ageing and age-related diseases. However, extensive analyses of the multiomics (circRNA/miRNA/mRNA and lncRNA/miRNA/mRNA) characteristics of adult Drosophila during ageing have not been reported. Here, differentially expressed circRNAs and microRNAs (miRNAs) between 7 and 42-day-old flies were screened and identified. Then, the differentially expressed mRNAs, circRNAs, miRNAs, and lncRNAs between the 7- and 42-day old flies were analysed to identify age-related circRNA/miRNA/mRNA and lncRNA/miRNA/mRNA networks in ageing Drosophila. Several key ceRNA networks were identified, such as the dme_circ_0009500/dme_miR-289-5p/CG31064, dme_circ_0009500/dme_miR-289-5p/frizzled, dme_circ_0009500/dme_miR-985-3p/Abl, and XLOC_027736/dme_miR-985-3p/Abl XLOC_189909/dme_miR-985-3p/Abl networks. Furthermore, real-time quantitative PCR (qPCR) was used to verify the expression level of those genes. Those results suggest that the discovery of these ceRNA networks in ageing adult Drosophila provide new information for research on human ageing and age-related diseases.

Animal gut microbiota plays an indispensable role in host adaptation to different altitude environments. At present, little is known about the mechanism of animal gut microbiota in host adaptation to high altitude environments. Here, we selected wild macaques, humans, and dogs with different levels of kinship and intimate relationships in high altitude and low altitude environments, and analyzed the response of their gut microbiota to the host diet and altitude environments. Alpha diversity analysis found that at high altitude, the gut microbiota diversity of wild macaques with more complex diet in the wild environments is much higher than that of humans and dogs with simpler diet (p

In-depth analysis of the characteristics of water and sedimentary changes at the confluence of the upper reaches of Weigan River and Heizi River is the prerequisite and basis for studying the future water and sedimentary situation in the Kizil Reservoir. This article analyzes the trend, mutation point and periodicity of the water and sedimentary sequence with multi-methods. The double cumulative curve method and the attenuation analysis method were used to analyze the impact of precipitation and human activities on runoff and sedimentary transport, and other influencing factors were qualitatively analyzed also. The results showed that: (1) In the past 60 years, Kizil Reservoir's runoff and sedimentary transport have shown an upward trend; the mutation point of the reservoirs' runoff was in 1994; the annual runoff had a periodicity of 7, 13, and 28 years; the annual sedimentary transport had a periodicity of 4, 12, and 26 years. (2) The impact of precipitation and human activities on runoff were 485.37% and −385.37%, respectively; the impact of precipitation on sedimentary transport was first high and then low, while the impact of human activities on sedimentary transport was first low and then high. Solar activities, ENSO, underlying surface changes and glacier snow melt also contributes to water-sediment changes to some extent. HIGHLIGHTS Various analytical methods were used to analyze the changing trends and mutation points, and comprehensively determined the results.; The ensemble empirical mode decomposition (EEMD) method was used to accurately analyze the periodicity of runoff and sedimentary transport.; The quantitative and qualitative systematic analysis of the factors which affect water and sedimentary changes were carried out.;

Chrysolophus amherstiae and Chrysolophus pictus are two related species of Chrysolophus. Understanding the differences in their mitochondrial genome structure is of great significance for studying their phylogenetic relationship. In this study, the full mitochondrial genome of C. amherstiae was sequenced and annotated and analyzed with C. pictus to reveal their structural differences. At the same time, the phylogenetic trees were constructed based on the BI and ML methods with other species of the same family to analyze the phylogenetic relationship of Phasianidae. The results showed that the basic characteristics of mitochondrial genomes of C. amherstiae and C. pictus were similar. By comparing the secondary structures of rRNA and the relative synonymous codon usage of protein-coding genes, the subtle differences of mitochondrial structures between the two were further demonstrated.

Forkhead box H1 (FoxH1) is a sexually dimorphic gene in Oreochromis niloticus, Oplegnathus fasciatus, and Acanthopagrus latus, indicating that it is essential for gonadal development. In the present study, the molecular characteristics and potential function of FoxH1 and the activation of the cyp19a1a promoter in vitro were evaluated in Monopterus albus. The levels of foxh1 in the ovaries were three times higher than those in the testes and were regulated by gonadotropins (Follicle-Stimulating Hormone and Human Chorionic Gonadotropin). FoxH1 colocalized with Cyp19a1a in the oocytes and granulosa cells of middle and late vitellogenic follicles. In addition, three FoxH1 binding sites were identified in the proximal promoter of cyp19a1a, namely, FH1 (−871/−860), FH2 (−535/−524), and FH3 (−218/−207). FoxH1 overexpression significantly attenuated the activity of the cyp19a1a promoter in CHO cells, and FH1/2 mutation increased promoter activity. Taken together, these results suggest that FoxH1 may act as an important regulator in the ovarian development of M. albus by repressing cyp19a1a promoter activity, which provides a foundation for the study of FoxH1 function in bony fish reproductive processes.

To explore the differences in the growth characteristics and population dynamics of Schizothorax wangchiachii populations in the Jinsha River (JSR) and the Yalong River (YLR), samples were collected in the upper reaches of the JSR (n = 230) from 2019 to 2020 and the middle reaches of the YLR (n = 187) from 2017 to 2018. In the JSR and YLR populations, the age range was 11 and 12 years old, respectively, and the best growth equation was the Von Bertalanffy equation. The comparative analysis of the two populations showed that the growth coefficient, initial sexual maturity age and age at first capture of the YLR population were greater than those of the JSR population. Comparing the mortality rates of the two groups, we found that the YLR population had the higher female mortality rate (0.658 years−1) and the lower male mortality rate (0.453 years−1). Our assessment of the three natural mortality rates showed that the Fcur of both male and female populations was greater than F25%, indicating that both populations were in an overexploited state. Therefore, we suggest considering the two groups as separate protection units and implementing management measures such as ecological regulation, restoration of tributary habitat and strengthening of fishing ban monitoring to protect their resources.

The extent of reservoir sedimentation is an important index related to the functional operation of reservoirs. Therefore, it is vital to accurately conduct sedimentary assessment. In this paper, the analytic hierarchy process was used to determine subjective weights, gray correlation analysis and entropy weight method were used to determine objective weights. The combination weights obtained using optimized combination weighting method based on genetic algorithm were more suitable for the comprehensive analysis of the impact of reservoir sedimentation. This was then used to constructed a multi-level fuzzy comprehensive evaluation model based on improved cloud model. A reservoir was selected as the study object, and its sedimentation impact level was evaluated: the numerical characteristics of the stratus cloud of the comment on the impact of the reservoir sedimentation were (0.6372, 0.0664, 0.0795). The results showed that the reservoir sedimentation is considered as severe influence, and the sedimentation has become a major problem that needed to be solved urgently. The results of this paper could provide insight for reservoir research domestically and abroad. Furthermore, it could also enable managers to more accurately grasp the severity of reservoir sedimentation. HIGHLIGHTS The combination of the cloud model and genetic algorithm in computing was used to study and evaluate the influence of reservoir sedimentation. The calculation result was proved to be much accurate than that of the traditional method, and is more in line with the actual purpose of the project.; The multi-level fuzzy comprehensive evaluation model was constructed based on the improved cloud model for the sole purpose of analyzing the influence of sedimentation on reservoir function.; The improved combination weighting method based on genetic algorithm was used to solve the combination weight issues.;

Slow lorises are small arboreal and nocturnal primates. Due to the illegal trade, a large number of slow lorises were confiscated into wildlife sanctuaries or rescue centers. The re-release has been considered a preferable approach for alleviating the captive pressure, but inappropriate and long-term confinement make it difficult to achieve this goal. In this study, we investigated and compared the fecal and oral microbiome of Bengal slow lorises (Nycticebus bengalensis) under long-term captivity (LC) and short-term captivity (SC) groups based on 16s rRNA high-throughput gene sequencing. The oral microbiome displayed higher Chao1 richness but lower Shannon and Simpson indices than the fecal microbiome. The Bengal slow lorises under long-term captivity had abundant pathogenic genera in both gut and oral microbiomes, such as Desulfovibrio, Actinomyces, Capnocytophaga, Neisseria, and Fusobacterium, while some specific bacterial taxa associated with intestinal balance were more enriched in the SC group. Due to the plant gum scarcity in the diet, both groups had a low abundance of Bifidobacterium. Function profile prediction indicated that the LC group was enriched with genetic information processing and metabolism pathways due to the stable food intake. The increased membrane transport and xenobiotic metabolism and degradation functions in the SC group could be explained by the function of the host microbiome in facilitating adaptation to changing environments and diets. The results demonstrated that the oral microbiome had the potential to be used as a regular surveillance tool. Also, current captive management should be improved to ensure reintroduction success.

Abstracts: Sex plasticity is easily influenced by external conditions in teleosts. In the present study, the morphogenesis of gonad development was described, and the sex ratio and gonadal development were also analysed after induction by exogenous E2. The primordial gonad was first visible at 9 dph, and PGCs entered mitosis at 90 dph. The time of transition from single oogonia to clusters of oogonia was from 95 dph until 126 dph, and the formation of chromatin nucleolar oocytes appeared at 212 dph. Compared with the ovary, the putative testis had fewer germ cells and a smaller gonad, which started to differentiate at 186 dph and formed spermatocytes at 237 dph. Exogenous E2 treatment promoted the expression of Cyp19a1a in differentiating gonads, and the female rate of all E2-treated groups was significantly higher than control groups (P

tRNA-derived small RNAs (tsRNAs) are derived from tRNA and include tRNA halves (tiRNAs) and tRNA fragments (tRFs). tsRNAs have been implicated in a variety of important biological functions, such as cell growth, transcriptional regulation, and apoptosis. Emerging evidence has shown that Ago1-guided and Ago2-guided tsRNAs are expressed at 3 and 30 days in Drosophila and that tRF biogenesis in fruit flies affects tRNA processing and tRNA methylation. However, a wide analysis of tsRNA patterns in different ages of Drosophila have not been reported via the small RNA sequencing method. In the present study, tsRNAs of young (7 days) and old (42 days) Drosophila were sequenced and their expression characteristics were analysed. Then, a specific tRF (named tRF-Trp-CCA-014) was determined and was found to be conserved in fruit flies, mice, and humans. The expression patterns of tRF-Trp-CCA-014 in different tissues and stages of fruit flies and mice, and mouse NIH/3T3 cells were detected. Furthermore, mouse embryonic fibroblast NIH/3T3 cells were used as a model to analyse the function and targets of tRF-Trp-CCA-014. The RNA-seq data of six groups (Mimics, Mimic NC, Inhibitors, Inhibitor NC, Aging (adriamycin), and Control (Normal)) in mouse NIH3T3 cells were analysed. The results showed that the number of tsRNAs at 42 days (417) was more than at 7 days (288); thus, it was enriched with age. tRFs-1 were the most enriched, followed by 5′-tRFs and 3′-tRFs. Twenty-one differentially expressed tsRNAs were identified between 7 days and 42 days. Then, the conserved tRF tRF-Trp-CCA-014 was identified and found to accumulate in aged fruit flies and aged mouse NIH3T3 cells. RNA-seq data showed that most differentially expressed genes were involved in the immune system, cancer: overview, and signal translation. Furthermore, tRF-Trp-CCA-014 was found to bind to the 3′UTR of H3C4 in a dual-luciferase reporter gene assay. tRF-Trp-CCA-014 and H3C4 were detected in the cytoplasm of aged NIH3T3 cells by RNA in situ hybridization. These results suggest that the H3C4 gene is the target of tRF-Trp-CCA-014. This study will advance the current understanding of tRF roles and their implication in Drosophila and mouse studies.

Enterocytozoon bieneusi is an obligate intracellular parasite and the most common pathogen of microsporidiosis in humans and animals. In this study, a total of 198 fecal samples were collected from 11 captive populations of Rhesus macaque in Chinese zoos, to investigate the prevalence and analyze the zoonotic potential of E. bieneusi by genotype of the Internal Transcribed Spacer (ITS) gene on the rRNA of E. bieneusi via nested PCR. Results showed that the average infection rate of E. bieneusi in the 11 populations was 13.6%, and the highest infection rate was 56.5% in the population of Xinjiang Tianshan Zoo. Seven genotypes were identified including 2 known genotypes (D and CM1) and 5 novel genotypes (Mul1, Mul2, Mul3, Mul4 and Mul5). Phylogenetic analysis revealed that the novel genotypes Mul2, Mul3, Mul4 and Mul5 belonged to Group 1 showed the zoonotic potential．These findings extend the distribution of E. bieneusi genotypes and provide baseline data for controlling E. bieneusi infection.

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.

Abstract Reintroduction is a key approach in the conservation of endangered species. In recent decades, many reintroduction projects have been conducted for conservation purposes, but the rate of success has been low. Given the important role of gut microbiota in health and diseases, we questioned whether gut microbiota would play a crucial role in giant panda's wild‐training process. The wild procedure is when captive‐born babies live with their mothers in a wilderness enclosure and learn wilderness survival skills from their mothers. During the wild‐training process, the baby pandas undergo wilderness survival tests and regular physical examinations. Based on their performance through these tests, the top subjects (age 2–3 years old) are released into the wild while the others are translocated to captivity. After release, we tracked one released panda (Zhangxiang) and collected its fecal samples for 5 months (January 16, 2013 to March 29 2014). Here, we analyzed the Illumina HiSeq sequencing data (V4 region of 16S rRNA gene) from captive pandas (n = 24), wild‐training baby pandas (n = 8) of which 6 were released and 2 were unreleased, wild‐training mother pandas (n = 8), one released panda (Zhangxiang), and wild giant pandas (n = 18). Our results showed that the gut microbiota of wild‐training pandas is significantly different from that of wild pandas but similar to that of captive ones. The gut microbiota of the released panda Zhangxiang gradually changed to become similar to those of wild pandas after release. In addition, we identified several bacteria that were enriched in the released baby pandas before release, compared with the unreleased baby pandas. These bacteria include several known gut‐health related beneficial taxa such as Roseburia, Coprococcus, Sutterella, Dorea, and Ruminococcus. Therefore, our results suggest that certain members of the gut microbiota may be important in panda reintroduction.

A balanced, diverse gut microbiota is vital for animal health. The microbial population is shaped by multiple factors including genetic background and environment, but other determinants remain controversial. Numerous studies suggest that the dominant factor is genetic background while others emphasize the environmental factors. Here, we bred asexual hybridization queens (AHQs) of honeybees through nutritional crossbreeding (laid in Apis mellifera colony but bred in Apis cerana colony), sequenced their gut microbiome, and compared it with normally bred sister queens to determine the primary factor shaping the gut microbiota. Our results showed that the dominant genera in the gut microbiota of AHQs were Brevundimonas, Bombella, and Lactobacillus, and its microbial community was more related to A. mellifera queens. The AHQs had a moderate number of different bacterial species and diversity, but total bacterial numbers were low. There were more significant taxa identified in the comparison between AHQ and A. cerana queen according to LEfSe analysis results. The only genetic-specific taxon we figured out was Brevundimonas. The growth of core bacterial abundance showed different characteristics among different queen groups in the first week after emerging. Collectively, this study suggested that the genetic background played a more dominant role than environmental factors in shaping the gut microbiota of honeybee queen and the microbiota of midgut was more sensitive than that of rectum to this impact.

Bmpr2 plays a central role in the regulation of reproductive development in mammals, but its role during ovarian development in fish is still unclear. To ascertain the function of bmpr2 in ovarian development in the ricefield eel, we isolated and characterized the bmpr2 cDNA sequence; the localization of Bmpr2 protein was determined by immunohistochemical staining; and the expression patterns of bmpr2 in ovarian tissue incubated with FSH and hCG in vitro were analyzed. The full-length bmpr2 cDNA was 3311 bp, with 1061 amino acids encoded. Compared to other tissues, bmpr2 was abundantly expressed in the ovary and highly expressed in the early yolk accumulation (EV) stages of the ovary. In addition, a positive signal for Bmpr2 was detected in the cytoplasm of oocytes in primary growth (PG) and EV stages. In vitro, the expression level of gdf9, the ligand of bmpr2, in the 10 ng/mL FSH treatment group was significantly higher after incubation for 4 h than after incubation for different durations. However, bmpr2 expression in the 10 ng/mL FSH treatment group at 2 h, 4 h and 10 h was significantly lower. Importantly, the expression level of bmpr2 and gdf9 in the 100 IU/mL hCG group had similar changes that were significantly decreased at 4 h and 10 h. In summary, Bmpr2 might play a pivotal role in ovarian growth in the ricefield eel, and these results provide a better understanding of the function of bmpr2 in ovarian development and the basic data for further exploration of the regulatory mechanism of gdf9 in oocyte development.

To understand the characteristics of the fish community structure and biomass particle-size spectrum in the upper reaches of the Jinsha River, fish and environmental surveys were conducted in 21 segments of the upper reaches of the Jinsha River in September 2019 and June 2020. A total of 4062 fish belonging to 2 orders, 5 families, 18 genera, and 28 species were collected. Among them, Cyprinidae fish were the most abundant (14 species), accounting for 50.00%. The Shannon index and Pielou evenness index values varied from 0.402–1.770 and 0.254–0.680, respectively. The dominant species of fish were Triplophysa stenura, Schizothorax wangchiachii, and Schizopygopsis malacanthus. Redundancy analysis (RDA) was used to analyse the relationship between the fish community and environmental factors. Velocity, altitude, and dissolved oxygen were the main influencing factors of fish community structure differences in the upper reaches of the Jinsha River. The abundance/biomass curves showed that the fish communities in the upper reaches of the Jinsha River were moderately or severely disturbed. The standardized biomass particle-size spectrum of fish showed that the degree of disturbance of fish in tributaries was much lower than that in the main stream. Compared with the historical data, the fish community structure in the Jinsha River has changed significantly, with the number of exotic species increasing, and the individual fish showing miniaturization and younger ages. It is suggested that habitat conservation strategies be adopted in the upper tributaries of the Jinsha River to provide a reference for the restoration of fishery resources and the conservation of fish diversity in the Yangtze River.

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

Blastocystis sp. is a common intestinal anaerobic parasite infecting non-human primates and many other animals. This taxon threatens the health of NHPs due to its high genetic diversity, impeding efforts to improve confined management and subsequent conservation practices. This study collected 100 and 154 fecal samples from captive macaques, gibbons, and slow lorises in the summer and winter, respectively. The Blastocystis infection, its gene subtypes, and its zoonotic potential based on small subunit ribosomal RNA (SSU rRNA) were analyzed. The prevalence of Blastocystis in the three primate genera was 57.79% (89/154) in the summer (2021) and 29.00% (29/100) in the winter (2020). Four zoonotic subtypes—ST1, ST2, ST3, and ST4—were identified. ST2 was the most prevalent subtype, suggesting that these animals may serve as reservoirs for pathogens of human Blastocystis infections. The macaques showed a more significant variation in Blastocystis infection between seasons than gibbons and slow lorises. The slow lorises in small cages and enclosure areas were potentially more infected by Blastocystis in the summer, indicating that inappropriate captive management may have detrimental effects on their health.

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy “LTPF” tRNA gene cluster compared to the novel “TLPF” order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

Apoptosis plays a key role in the effective removal of excessive and defective germ cells, which is essential for sequential hermaphroditism and sex change in vertebrates. The ricefield eel, Monopterus albus is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. Previous studies have demonstrated that apoptosis is involved in sex change in M. albus. However, the apoptotic signaling pathway is unclear. In the current study, we explored the underlying mechanism of apoptosis during gonadal development and focused on the role of the mitochondrial apoptosis signaling pathway in sex change in M. albus. Flow cytometry was performed to detect apoptosis in gonads at five sexual stages and ovary tissues exposed to hydrogen peroxide (H2O2) in vitro. Then the expression patterns of key genes and proteins in the mitochondrial pathway, death receptor pathway and endoplasmic reticulum (ER) pathway were examined. The results showed that the apoptosis rate was significantly increased in the early intersexual stage and then decreased with the natural sex change from female to male. Quantitative real-time PCR revealed that bax, tnfr1, and calpain were mainly expressed in the five stages. ELISA demonstrated that the relative content of cytochrome-c (cyt-c) in the mitochondrial pathway was significantly higher than that of caspase8 and caspase12, with a peak in the early intersexual stage, while the levels of caspase8 and caspase12 peaked in the late intersexual stage. Interestingly, the Pearson’s coefficient between cyt-c and the apoptosis rate was 0.705, which suggests that these factors are closely related during the gonadal development of M. albus. Furthermore, the cyt-c signal was found to be increased in the intersexual stage by immunohistochemistry. After incubation with H2O2, the mRNA expression of mitochondrial pathway molecules such as bax, apaf-1, and caspase3 increased in ovary tissues. In conclusion, the present results suggest that the mitochondrial apoptotic pathway may play a more important role than the other apoptotic pathways in sex change in M. albus.

Schizothorax kozlovi Nikolsky is a unique cold−water fish in the upper reaches of the Yangtze River in China and has high economic value. In our study, genetic diversity and population structure analyses were performed on seven wild populations (originating from the Jinsha River, Yalong River, and Wujiang River) in the upper reaches of the Yangtze River by genotyping by sequencing (GBS). The results indicated that a total of 303,970 single−nucleotide polymorphisms (SNPs) were identified from the seven wild populations. Lower genetic diversity was exhibited among the intrapopulations of the three tributaries, and the Wujiang River population had significant genetic differentiation when compared to the Jinsha River and Yalong River populations. Furthermore, the selected SNPs were enriched in cellular processes, environmental adaptation, signal transduction, and related metabolic processes between the Wujiang population and the other two populations. The above results indicate that the populations of S. kozlovi have different degrees of tolerance and selection pressure in response to temperature and altitude. The Wujiang intrapopulation has greater genetic diversity and differentiation than the Jinsha River and Yalong River populations, which demonstrates that the Jinsha and Yalong populations require more attention and resources for their protection. The results of this study will increase our understanding of the diversity of S. kozlovi in the upper reaches of the Yangtze River and provide a basis for the conservation and utilization of wild resources.

The Asiatic softshell turtle, also known as the black-rayed softshell turtle (Amyda cartilaginea; Accession no: MT039230), is found in northeastern India (Mizoram), Brunei Darussalam, Indonesia, Malaysia, Singapore, Myanmar, Laos, Vietnam, Cambodia, and Thailand. This turtle is thought to have been introduced into the Sunda Islands, Sulawesi, and Yunnan, China, through the Malay Peninsula to Sumatra, Java, and Borneo. Herein, we determined the complete mitochondrial genome of A. cartilaginea for the first time using next-generation sequencing (NGS). The assembled mitogenome was 16,763 bp in length and encoded 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes (12S rRNA and 16S rRNA), and one control region (CR). The PCGs based maximum-likelihood phylogeny discriminated A. cartilaginea from other Testudines and clusters within family Trionychidae with the sister taxa of Nilssonia nigricans.

Abstract In this paper, a variant of Mehrotra-type predictor–corrector algorithm is proposed for P∗(κ) $P_{*}(\kappa )$ linear complementarity problems. In this algorithm, a safeguard step is used to avoid small step sizes and a new corrector direction is adopted. The algorithm has polynomial iteration complexity and the iteration bound is O((14κ+11)(1+4κ)(1+2κ)nlog(x0)Ts0ε) $O ((14\kappa +11)\sqrt{(1+4 \kappa )(1+2\kappa )}~n\log \frac{(x^{0})^{T}s^{0}}{\varepsilon } )$. Some numerical results are reported as well.

Bengal slow lorises (Nycticebus bengalensis) are threatened by illegal trade. Subsequently, numerous wild-born individuals are rescued and transferred to rescue centers. Metabonomic analysis of intestinal microbiomes has increasingly played a vital role in evaluating the effects of dietary alteration on the captive status of endangered non-human primates. A synthetic analysis was done to test the differences in gut microbes and fecal metabolites between two dietary groups of Bengal slow lorises across 8 weeks. Dietary interventions led to intra-group convergence and inter-group variation in the composition of intestinal flora, metabolites, and short-chain fatty acids (SCFAs). The control diet, consisting of gums and honey, significantly increased the abundance of some potential probiotics, such as Bifidobacterium and Roseburia, and the concentration of some anti-disease related metabolites. The decrease in some amino acid metabolites in the original group fed without gums was attributed to poor body condition. Some distinct SCFAs found in the control group indicated the dietary alteration herein was fat-restricted but fiber deficient. Cognizant of this, plant exudates and fiber-enriched food supplies should be considered an optimal approach for dietary improvement of the confiscated and captive Bengal slow lorises.

Previous researches have highlighted the finding that the availability of food resources plays an important role in reshaping the gut microbiota of wild animals as they to adapt to the surrounding habitat. The gut microbial characteristics were examined in Tibetan macaques (Macaca thibetana) from eastern China (Huangshan, Anhui), but little is known about other geographical populations. Herein, we collected 40 fecal samples of Tibetan macaques from the Labahe Natural Reserve in western China over four seasons and analyzed the seasonal variation of gut microbiota using 16S rRNA gene V3–V4 sequencing. The analysis showed that gut microbial diversity and composition varied by season. Autumn exhibited the lowest alpha diversity and remarkably separated from the other three seasons for beta diversity. There were statistically significant differences in the composition of gut microbiota during different seasons. Autumn presented the lowest Firmicutes/Bacteroidetes ratio while winter had the highest. The family Ruminococcaceae and Lactobacillaceae within Firmicutes accounted for a high proportion in each season, but the abundance of Lactobacillaceae was significantly lower in winter. During autumn, Preotella, which can utilize carbohydrates increased significantly. Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) prediction revealed that metabolism of cofactors of vitamin pathways and glycan biosynthesis pathways were highly enriched in autumn. Our results revealed that the gut microbiota potentially play an important role in the adaptation of wild Tibetan macaques to the seasonal fluctuations of food availability, which may contribute to the conservation of this endangered species.

MicroRNA (miRNA) is known to be an important regulator of muscle growth and development. The regulation of microRNA on the skeletal muscle phenotype of animals is mainly achieved by regulating the proliferation and differentiation of myoblasts. In this study, we sequenced a total of 60 samples from 15 developing stages of the pectoral muscle and five other tissues at 300 days of Tibetan chicken. We characterized the expression patterns of miRNAs across muscle developmental stages, and found that the chicken growth and development stage was divided into early-embryonic and late-embryonic as well as postnatal stages. We identified 81 and 21 DE-miRNAs by comparing the miRNA profiles of pectoral muscle of three broad periods and different tissues, respectively; and 271 miRNAs showed time-course patterns. Their potential targets were predicted and used for functional enrichment to understand their regulatory functions. Significantly, GgmiRNA-454 is a time-dependent and tissue-differential expression miRNA. In order to elucidate the role of gga-miRNA-454 in the differentiation of myoblasts, we cultured chicken myoblasts in vitro. The results show that although gga-miRNA-454-3p initiates increase and thereafter decrease during the chicken myoblasts differentiation, it had no effect on primary myoblasts proliferation. Furthermore, we confirm that gga-miRNA-454 inhibits myoblast differentiation by targeting the myotube-associated protein SBF2.

The liver is the major organ of lipid biosynthesis in the chicken. In laying hens, the liver synthesizes most of the yolk precursors and transports them to developing follicles to produce eggs. However, a systematic investigation of the long non-coding RNA (lncRNA) and mRNA transcriptome in liver across developmental stages is needed. Here, we constructed 12 RNA libraries from liver tissue during four developmental stages: juvenile (day 60), sexual maturity (day 133), peak laying (day 220), and broodiness (day 400). A total of 16,930 putative lncRNAs and 18,260 mRNAs were identified. More than half (53.70%) of the lncRNAs were intergenic lncRNAs. The temporal expression pattern showed that lncRNAs were more restricted than mRNAs. We identified numerous differentially expressed lncRNAs and mRNAs by pairwise comparison between the four developmental stages and found that VTG2, RBP, and a novel protein-coding gene were differentially expressed in all stages. Time-series analysis showed that the modules with upregulated genes were involved in lipid metabolism processes. Co-expression networks suggested functional relatedness between mRNAs and lncRNAs; the DE-lncRNAs were mainly involved in lipid biosynthesis and metabolism processes. We showed that the liver transcriptome varies across different developmental stages. Our results improve our understanding of the molecular mechanisms underlying liver development in chickens.

We sequenced and annotated the complete mitogenome sequence of Oreolalax omeimontis (17,675 bp long). The mitogenome encoded 13 protein-coding genes (PCG), 2 ribosomal RNA (rRNA) genes, 23 transfer RNA (tRNA) genes, and a control region (GenBank accession number MN803321). The phylogenetic tree conforms the close relationship of O. omeimontis and O. multipunctatus, and a monophyletic clade of genus Oreolalax.

We newly obtained the near complete mitochondrial genome of Oreolalax schmidti (18,481bp) by using polymerase chain reaction (PCR) in this study. It includes 13 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 23 transfer RNA (tRNA) genes. The phylogenetic tree indicates that the Oreolalax schmidti is closely related to the O. lichuanensis. This report will help the further studies of Oreolalax species classification and source protection.

A complete mitogenome sequence of the blue-crested lizard (Calotes mystaceus) was determined in this study. The 16,506 bp genome consists 13 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes, and a control region. The phylogenetic tree reveals that the Calotes mystaceus is closely related to the C. versicolor. This report provides the basic data for further studies of Calotes species classification and phylogeny.

Red-and-white giant flying squirrel (Petaurista alborufus) is a widely distributed but rarely studied species. Here, we sequenced the complete mitochondrial genome about it, which is a circular genome of 16,511 bp and have a typical structure with mammals. The base composition shows higher composition of A\T. There are 61 variations and eight base deletions between two sequences when compared with the reported sequence (JQ743657). In this study, we obtain a new mitochondrial genome sequence of red-and-white giant flying squirrel, which provide novel molecular information that could prove fundamental to our understanding of diversity research and phylogenetic structure.

Cynops cyanurus chuxiongensis is distributed in Yunnan plateau. In this paper, we sequenced and determined the complete mitochondrial genome of C. c. chuxiongensis. The assembled mitogenome is 16,465 bp in length and encoding 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes (12S rRNA and 16S rRNA), and one control region (D-loop). The phylogenetic trees indicated C. c. chuxiongensis (KY418068) has the closest relationship with C. orientalis + C. orphicus, and clustered within the clade of genus Cynops.

Due to the multiple causes, the population of roe deer has declined significantly. In this study, we analyzed the complete mitogenome of Capreolus pygargus bedfordi, whose genome was 16,357 bp long. There were 13 protein-coding genes (PCG), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA genes, and one control region. Nine PCGs started with ATG, while NAD2, NAD3, and NAD5 genes commenced with ATA, and ND4L began with GTG. ND6 and eight tRNA genes were encoded on the L-strand. These results provide newer molecular information, which contribute to its molecular and phylogenetic studies, and genetic diversity conservation.

Ferret-badger (Melogale moschata) is a widely distributed and rarely studied species. Therefore, we sequenced its complete mitochondrial genome, which is circular in shape with 16,511 bp length and have a structure typical to mammals. There is 98.67% similarity when compared with the reported sequence (HM106328). It contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control regions (D-loop). This study can provide novel molecular information and genetic diversity of ferret-badger.

The complete mitochondrial of Cynops cyanurus cyanurus was determined by using Sanger sequencing and the mitogenome size of C. c. cyanurus was 16,465 bp. It contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs, 1 control region (D-loop) and a non-coding region (NC). In addition, the phylogenetic tree shows that C. c. cyanurus has the closest relationship with C. c. chuxiongensis. Hoping this study can improve our understanding of the species evolutionary relationship of Cynops.

The complete mitochondrial genome of the Oreolalax major (17,786 bp long) was obtained in this study. It includes 13 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 23 transfer RNA (tRNA) genes (GenBank accession number MN803320). The phylogenetic tree indicates that the O. major is closely related to the O. xiangchengensis.