Search

ABSTRACTAlthough vaginitis is closely related to vaginal microecology in females, the precise composition and functional potential of different types of vaginitis remain unclear. Here, metagenomic sequencing was applied to analyze the vaginal flora in patients with various forms of vaginitis, including cases with a clue cell proportion ranging from 1% to 20% (Clue1_20), bacterial vaginitis (BV), vulvovaginal candidiasis (VVC), and BV combined with VVC (VVC_BV). Our results identified Prevotella as an important biomarker between BV and Clue1_20. Moreover, a gradual decrease was observed in the relative abundance of shikimic acid metabolism associated with bacteria producing indole as well as a decline in the abundance of Gardnerella vaginalis in patients with BV, Clue1_20, and healthy women. Interestingly, the vaginal flora of patients in the VVC_BV group exhibited structural similarities to that of the VVC group, and its potentially functional characteristics resembled those of the BV and VVC groups. Finally, Lactobacillus crispatus was found in high abundance in healthy samples, greatly contributing to the stability of the vaginal environment. For the further study of L. crispatus, we isolated five strains of L. crispatus from healthy samples and evaluated their capacity to inhibit G. vaginalis biofilms and produce lactic acid in vitro to select the potential probiotic candidate for improving vaginitis in future clinical studies. Overall, we successfully identified bacterial biomarkers of different vaginitis and characterized the dynamic shifts in vaginal flora between patients with BV and healthy females. This research advances our understanding and holds great promise in enhancing clinical approaches for the treatment of vaginitis.IMPORTANCEVaginitis is one of the most common gynecological diseases, mostly caused by infections of pathogens such as Candida albicans and Gardnerella vaginalis. In recent years, it has been found that the stability of the vaginal flora plays an important role in vaginitis. Furthermore, the abundant Lactobacillus-producing rich lactic acid in the vagina provides a healthy acidic environment such as Lactobacillus crispatus. The metabolites of Lactobacillus can inhibit the colonization of pathogens. Here, we collected the vaginal samples of patients with bacterial vaginitis (BV), vulvovaginal candidiasis (VVC), and BV combined with VVC to discover the differences and relationships among the different kinds of vaginitis by metagenomic sequencing. Furthermore, because of the importance of L. crispatus in promoting vaginal health, we isolated multiple strains from vaginal samples of healthy females and chose the most promising strain with potential probiotic benefits to provide clinical implications for treatment strategies.

Red pandas evolved from carnivores to herbivores and are unique within Carnivora. Red pandas and carnivorous mammals consume milk during the suckling period, while they consume bamboo and meat during the adult period, respectively. Red pandas and carnivorous mammal ferrets have a close phylogenetic relationship. To further investigate the molecular mechanisms of dietary changes and nutrient utilization in red pandas from suckling to adult, comparative analysis of the whole transcriptome was performed on stomach tissues from red pandas and ferrets during the suckling and adult periods. The main results are as follows: (1) we identified ncRNAs for the first time in stomach tissues of both species, and found significant expression changes of 109 lncRNAs and 106 miRNAs in red pandas and 756 lncRNAs and 109 miRNAs in ferrets between the two periods; (2) up-regulated genes related to amino acid transport regulated by lncRNA-miRNA-mRNA networks may efficiently utilize limited bamboo amino acids in adult red pandas, while up-regulated genes related to amino acid degradation regulated by lncRNAs may maintain the balance of amino acid metabolism due to larger daily intakes in adult ferrets; and (3) some up-regulated genes related to lipid digestion may contribute to the utilization of rich nutrients in milk for the rapid growth and development of suckling red pandas, while up-regulated genes associated with linoleic acid metabolism regulated by lncRNA-miRNA-mRNA networks may promote cholesterol decomposition to reduce health risks for carnivorous adult ferrets. Collectively, our study offers evidence of gene expression adaptation and ncRNA regulation in response to specific dietary changes and nutrient utilization in red pandas during suckling and adult periods.

Abstract Rhesus macaques (Macaca mulatta, RMs) are widely used in sexual maturation studies due to their high genetic and physiological similarity to humans. However, judging sexual maturity in captive RMs based on blood physiological indicators, female menstruation, and male ejaculation behavior can be inaccurate. Here, we explored changes in RMs before and after sexual maturation based on multi-omics analysis and identified markers for determining sexual maturity. We found that differentially expressed microbiota, metabolites, and genes before and after sexual maturation showed many potential correlations. Specifically, genes involved in spermatogenesis (TSSK2, HSP90AA1, SOX5, SPAG16, and SPATC1) were up-regulated in male macaques, and significant changes in gene (CD36), metabolites (cholesterol, 7-ketolithocholic acid, and 12-ketolithocholic acid), and microbiota (Lactobacillus) related to cholesterol metabolism were also found, suggesting the sexually mature males have stronger sperm fertility and cholesterol metabolism compared to sexually immature males. In female macaques, most differences before and after sexual maturity were related to tryptophan metabolism, including changes in IDO1, IDO2, IFNGR2, IL1Β, IL10, L-tryptophan, kynurenic acid (KA), indole-3-acetic acid (IAA), indoleacetaldehyde, and Bifidobacteria, indicating that sexually mature females exhibit stronger neuromodulation and intestinal immunity than sexually immature females. Cholesterol metabolism-related changes (CD36, 7-ketolithocholic acid, 12-ketolithocholic acid) were also observed in female and male macaques. Exploring differences before and after sexual maturation through multi-omics, we identified potential biomarkers of sexual maturity in RMs, including Lactobacillus (for males) and Bifidobacterium (for females) valuable for RM breeding and sexual maturation research.

Abstract Macaques (genus Macaca) are the most widely distributed non‐human primates, and their evolutionary history, gene expression profiles, and genetic differences have been extensively studied. However, the DNA methylomes of macaque species are not available in public databases, which hampers understanding of epigenetic differences among macaque species. Epigenetic modifications can potentially affect development, physiology, behavior, and evolution. Here, we investigated the methylation patterns of the Tibetan macaque (M. thibetana; TM), Chinese rhesus macaque (M. mulatta lasiota; CR), and crab‐eating macaque (M. fascicularis; CE) through whole‐genome bisulfite sequencing from peripheral blood. We compared genome‐wide methylation site information for the three species. We identified 12,128 (CR vs. CE), 59,165 (CR vs. TM), and 39,751 (CE vs. TM) differentially methylated regions (DMRs) in the three macaques. Furthermore, we obtained the differentially expressed genes (DEGs) among the three macaque species. The differences between CR and CE were smaller at both the methylome and transcriptome levels than compared with TM (CR vs. TM and CE vs. TM). We also found a change in the density of single nucleotide mutations in DMRs relative to their flanking regions, indicating a potential mechanism through which genomic alterations may modulate methylation landscapes, thereby influencing the transcriptome. Functional enrichment analyses showed the DMR‐related genes were enriched in developmental processes and neurological functions, such as the growth hormone‐related pathway, insulin secretion pathway, thyroid hormone synthesis pathway, morphine addiction, and GABAergic synapses. These differences may be associated with variations in physiology and habitat among the macaques. Our study provides one of the first genome‐wide comparisons of genetic, gene expression, and epigenetic variations across different macaques. Our results should facilitate further research on comparative genomic and genetic differences in macaque species.

IntroductionIn an era increasingly defined by the challenge of antibiotic resistance, this study offers groundbreaking insights into the antibacterial properties of two distinct Lactiplantibacillus plantarum strains, TE0907 and TE1809, hailing from the unique ecosystem of Bufo gargarizans. It uniquely focuses on elucidating the intricate components and mechanisms that empower these strains with their notable antibacterial capabilities.MethodsThe research employs a multi-omics approach, including agar diffusion tests to assess antibacterial efficacy and adhesion assays with HT-29 cells to understand the preliminary mechanisms. Additionally, gas chromatography-mass spectrometry (GC-MS) is employed to analyze the production of organic acids, notably acetic acid, and whole-genome sequencing is utilized to identify genes linked to the biosynthesis of antibiotics and bacteriocin-coding domains.ResultsThe comparative analysis highlighted the exceptional antibacterial efficacy of strains TE0907 and TE1809, with mean inhibitory zones measured at 14.97 and 15.98 mm, respectively. A pivotal discovery was the significant synthesis of acetic acid in both strains, demonstrated by a robust correlation coefficient (cor ≥ 0.943), linking its abundance to their antimicrobial efficiency. Genomic exploration uncovered a diverse range of elements involved in the biosynthesis of antibiotics similar to tetracycline and vancomycin and potential regions encoding bacteriocins, including Enterolysin and Plantaricin.ConclusionThis research illuminates the remarkable antibacterial efficacy and mechanisms intrinsic to L. plantarum strains TE0907 and TE1809, sourced from B. gargarizans. The findings underscore the strains' extensive biochemical and enzymatic armamentarium, offering valuable insights into their role in antagonizing enteric pathogens. These results lay down a comprehensive analytical foundation for the potential clinical deployment of these strains in safeguarding animal gut health, thereby enriching our understanding of the role of probiotic bacteria in the realm of antimicrobial interventions.

Background Up the reproductive tract, there are large differences in the composition of vaginal microbes. Throughout the menstrual cycle, the structure of the vaginal microbiome shifts. Few studies have examined both in combination. Our study was designed to explore trends in the microbiome of different parts of the vagina in healthy women over the menstrual cycle. Methods We performed metagenomic sequencing to characterize the microbiome differences between the cervical orifice and mid-vagina throughout the menstrual cycle. Results Our results showed the vaginal microbiome of healthy women in the cervical orifice and the mid-vagina was similar during the periovulatory and luteal phases, with Lactobacillus being the dominant bacteria. In the follicular phase, Acinetobacter was detected in the cervical orifice. From the follicular phase to the luteal phase, the community state types (all five community status types were defined as CSTs) in samples No. 10 and No. 11 changed from CST III to CST I. In addition, the composition of the vaginal microbiome in healthy women from different regions of China was significantly different. We also detected viruses including Human alphaherpesvirus 1 (HSV-1) during periovulatory phase. Conclusion This study is valuable for understanding whether the microbial composition of the vagina is consistent in different parts of the menstrual cycle.

ABSTRACTRhesus macaque (RM, Macaca mulatta), as an important model animal, commonly suffers from chronic diarrheal disease, challenging the breeding of RMs. Gut microbiomes play key roles in maintaining intestinal health of RMs. However, it is still unclear about more features of gut microbiome as responsible for intestinal health of RMs. In this study, we performed de novo assembly of metagenome-assembled genomes (MAGs) based on fecal metagenomes from chronic diarrheal RMs and asymptomatic individuals. In total of 731 non-redundant MAGs with at least 80% completeness were reconstructed in this study. More than 97% MAGs were novel genomes compared with more than 250,000 reference genomes. MAGs of Campylobacter and Helicobacteraceae from RM guts mainly carried flagella-associated virulence genes and chemotaxis-associated virulence genes, which might mediate motility and adhesion of bacteria. Comparing to MAGs of Campylobacter from humans, distributions and functions of these MAGs of Campylobacter from RMs exhibited significant differences. Most members of Bacteroidota, Spirochaetota, Helicobacteraceae, Lactobacillaceae and Anaerovibrio significantly decreased in guts of chronic diarrhea RMs. More than 92% MAGs in this study were not contained in 2,985 MAGs previously reported from other 22 non-human primates (NHPs), expanding the microbial diversity in guts of NHPs. The distributions and functions of gut microbiome were prominently influenced by host phylogeny of NHPs. Our results could help to more clearly understand about the diversity and function of RMs gut microbiome.

Ticks and tick-borne diseases have negative impacts on the health of wild animals including endangered and vulnerable species. The giant panda (Ailuropoda melanoleuca), a vulnerable and iconic flagship species, is threatened by tick infestation as well. Not only can ticks cause anemia and immunosuppression in the giant panda, but also bacterial and viral diseases. However, previous studies regarding tick infestation on giant pandas were limited in scope as case reports from sick or dead animals. In this study, an investigation focusing on the tick infestation of a reintroduced giant panda at the Daxiangling Reintroduction Base in Sichuan, China was conducted. Ticks were routinely collected and identified from the ears of the giant panda from March to September in 2021. A linear model was used to test the correlation between tick abundance and climate factors. All ticks were identified as Ixodes ovatus. Tick abundance was significantly different among months. Results from the linear model showed temperature positively correlated to tick abundance, while air pressure had a negative correlation with tick abundance. To the best of our knowledge, this study is the first reported investigation of tick species and abundance on a healthy giant panda living in the natural environment, and provides important information for the conservation of giant pandas and other species sharing the same habitat.

Abstract There are only six isolated living giant panda populations, and a comprehensive understanding of their genetic health status is crucial for the conservation of this vulnerable species. Liangshan Mountains is one of the main distribution areas of living giant pandas and is outside the newly established Giant panda national park. In this study, 971 giant panda fecal samples were collected in the heartland of Liangshan Mountains (Mabian Dafengding Nature Reserve: MB; Meigu Dafengding Nature Reserve: MG; and Heizhugou Nature Reserve: HZG). Microsatellite markers and mitochondrial D‐loop sequences were used to estimate population size and genetic diversity. We identified 92 individuals (MB: 27, MG: 22, HZG: 43) from the three reserves. Our results showed that: (1) genetic diversity of three giant panda populations was moderate; (2) several loci deviated significantly from the Hardy–Weinberg equilibrium and almost all these deviated loci showed significant heterozygote deficiencies and inbreeding; (3) three giant panda populations have substantial genetic differentiation with the most differentiation between MB and the two other populations; and (4) a large amount of giant panda feces outside the three reserves were found, implying the existence of protection gap. These results indicated that under stochastic events, the giant panda populations in Liangshan Mountains are at risk of genetic decline or extinction and urgent need of human management. This study revealed that high attention should be paid to the protection of these giant panda populations outside the Giant panda national park, to ensure their survival in their distribution areas.

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

Abstract The Heilongjiang brown frog (Rana amurensis) is widely used in traditional Chinese medicine. In particular, the oviduct and skin have been developed into various health products. However, limited numbers of complete genomes of amphibian species have been reported, excluding the Heilongjiang brown frog. Here, the transcriptomes of 45 samples from the liver, spleen, heart, ovaries, thigh muscles, skin, oviduct, stomach and intestine of five Heilongjiang brown frog were reassembled and analyzed. A total of 1,085,532 unigenes with an average length of 676.6 bp and N50 of 722 bp were obtained. Comparative transcriptomics of different tissues detected tissue-specific expression. There were 3248 differentially expressed genes (DEGs) in the ovary, and the number of unique DEGs between the ovary and spleen was the largest. The results of DEGs enrichment showed there were many pathways and items related to protein synthesis and metabolism in the oviduct. The DEGs of the skin were enriched with many bacterial defense items, indicating that there were a large number of antimicrobial peptides in the skin. Thus, these were suitable as biological sources for the development and extraction of antimicrobial peptides. Through the assembly of transcriptome sequencing data and functional annotation of the Heilongjiang brown frog genome, this study provides reference materials for further exploring and utilizing functional gene resources of frogs and lays a foundation for medical research and the development of new products.

Abstract Background Spodoptera litura is an important polyphagous pest that causes significant damage to the agricultural sector. We performed RNA-seq of 15 S. litura individuals from larval (fifth and sixth instar larvae), chrysalis, and adult developmental stages. We also compared the S. litura transcriptome data with Spodoptera frugiperda across the same developmental stages, which was sequenced in our previous study. Results A total of 101,885 differentially expressed transcripts (DETs) were identified in S. litura. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that S. litura may undergo active xenobiotic and detoxifying metabolism during its larval and adult stages, which may explain difficulties with current population control measures. We also found that DETs of single-copy orthologous genes between S. litura and S. frugiperda were involved in basic metabolism and development. However, energy and metabolic processes genes had a higher expression in S. litura, whereas nervous and olfactory function genes had a higher expression in S. frugiperda. Metagenomics analysis in larval S. litura and S. frugiperda revealed that microbiota participate in the detoxification and metabolism processes, but the relative abundance of detoxification-related microbiota was more abundant in S. frugiperda. Transcriptome results also confirmed the detoxification-related pathway of S. frugiperda was more abundant than in S. litura. Conclusions Significant changes at transcriptional level were identified during the different development stages of S. litura. Importantly, we also identified detoxification associated genes and gut microbiota between S. litura and S. frugiperda at different developmental stages, which will be valuable in revealing possible mechanisms of detoxification and development in these two lepidopterans.

Abstract Background The giant panda (Ailuropoda melanoleuca) is a threatened species endemic to China. Alopecia, characterized by thinning and broken hair, mostly occurs in breeding males. Alopecia significantly affects the health and public image of the giant panda and the cause of alopecia is unclear. Results Here, we researched gene expression profiles of four alopecia giant pandas and seven healthy giant pandas. All pandas were approximately ten years old and their blood samples collected during the breeding season. A total of 458 up-regulated DEGs and 211 down-regulated DEGs were identified. KEGG pathway enrichment identified that upregulated genes were enriched in the Notch signaling pathway and downregulated genes were enriched in ribosome, oxidative phosphorylation, and thermogenesis pathways. We obtained 28 hair growth-related DEGs, and identified three hub genes NOTCH1, SMAD3, and TGFB1 in PPI analysis. Five hair growth-related signaling pathways were identified with abnormal expression, these were Notch, Wnt, TGF-β, Mapk, and PI3K-Akt. The overexpression of NOTCH1 delays inner root sheath differentiation and results in hair shaft abnormalities. The delayed hair regression was associated with a significant decrease in the expression levels of TGFB1. Conclusions Our data confirmed the abnormal expression of several hair-related genes and pathways and identified alopecia candidate genes in the giant panda. Results of this study provide theoretical basis for the establishment of prevention and treatment strategies for giant pandas with alopecia.

Abstract Background Trueperella pyogenes and Pseudomonas aeruginosa are two important bacterial pathogens closely relating to the occurrence and development of forest musk deer respiratory purulent disease. Although T. pyogenes is the causative agent of the disease, the subsequently invaded P. aeruginosa will predominate the infection by producing a substantial amount of quorum-sensing (QS)-controlled virulence factors, and co-infection of them usually creates serious difficulties for veterinary treatment. In order to find a potential compound that targets both T. pyogenes and P. aeruginosa, the antibacterial and anti-virulence capacities of 55 compounds, which have similar core structure to the signal molecules of P. aeruginosa QS system, were tested in this study by performing a series of in vitro screening experiments. Results We identified that furazolidone could significantly reduce the cell densities of T. pyogenes in mono-culture or in the co-culture with P. aeruginosa. Although the growth of P. aeruginosa could also be moderately inhibited by furazolidone, the results of phenotypic identification and transcriptomic analysis further revealed that sub-inhibitory furazolidone had remarkable inhibitory effect on the biofilm production, motility, and QS system of P. aeruginosa. Moreover, furazolidone could efficiently protect Caenorhabditis elegans models from P. aeruginosa infection under both fast-killing and slow-killing conditions. Conclusions This study reports the antibacterial and anti-virulence abilities of furazolidone on T. pyogenes and P. aeruginosa, and provides a promising strategy and molecular basis for the development of novel anti-infectious drugs to dealing with forest musk deer purulent disease, or other diseases caused by T. pyogenes and P. aeruginosa co-infection.

Abstract Background Chronic diarrhea is a common disease causing morbidity and mortality of captive rhesus macaques (RMs, Macaca mulatta). Chronic diarrhea in RMs is typically characterized by long-term diarrhea and a weak response to antibiotic treatment. Diarrhea is also a common disease in humans and can cause death. However, the etiology of about half of diarrheal cases of humans is still unclear. Therefore, we performed shotgun metagenomic sequencing to characterize the differences in the gut microbiome and resistome of chronic diarrhea RMs and asymptomatic individuals. Results Our results showed Lactobacillus spp. (mainly L. johnsonii, L. reuteri and L. amylovorus) were significantly depleted in chronic diarrhea RM guts compared to asymptomatic individuals (5.2 vs 42.4%). Functional annotation of genes suggested these Lactobacillus spp. carried genes involved in the adhesion of intestinal epithelial cells and production of bacteriocin. Chronic diarrhea RM guts also had a significantly greater abundance of many other gut bacteria, including mucin-degrading bacteria and opportunistic pathogens. The metabolic pathways of chronic diarrhea RM gut microbiome were enriched in aerobactin biosynthesis, while the metabolic pathways of asymptomatic RM gut microbiome were enriched in the production of short-chain fatty acids (SCFAs). Chronic diarrhea RM guts had a significantly greater abundance of antibiotic resistance genes (ARGs), such as ermF, aph(3’)-IIIa, ermB, and floR. The strains isolated from feces and tissue fluid of chronic diarrhea RMs had higher resistance rates to the majority of tested antibiotics, but not cephamycin and carbapenem antibiotics. Gut microbial composition comparisons showed that several captive nonhuman primate (NHP) guts were more similar to the guts of humans with a non-westernized diet than humans with a westernized diet. Chronic diarrhea RM gut microbiome was strikingly similar to rural-living humans with diarrhea and humans with a non-westernized diet than asymptomatic RMs. Conclusions Our results suggested chronic diarrhea significantly altered the composition and metabolic pathways of the RM gut microbiome. The frequent use of antibiotics caused antibiotic resistance in chronic diarrhea RM gut microbiome with serious consequences for individual treatment and survival. The findings of this study will help us to improve the effective prevention and treatment of diarrhea in RMs. Video Abstract

DNA methylation modification can regulate gene expression without changing the genome sequence, which helps organisms to rapidly adapt to new environments. However, few studies have been reported in non-model mammals. Giant panda (Ailuropoda melanoleuca) is a flagship species for global biodiversity conservation. Wildness and reintroduction of giant pandas are the important content of giant pandas’ protection. However, it is unclear how wildness training affects the epigenetics of giant pandas, and we lack the means to assess the adaptive capacity of wildness training giant pandas. We comparatively analyzed genome-level methylation differences in captive giant pandas with and without wildness training to determine whether methylation modification played a role in the adaptive response of wildness training pandas. The whole genome DNA methylation sequencing results showed that genomic cytosine methylation ratio of all samples was 5.35%–5.49%, and the methylation ratio of the CpG site was the highest. Differential methylation analysis identified 544 differentially methylated genes (DMGs). The results of KEGG pathway enrichment of DMGs showed that VAV3, PLCG2, TEC and PTPRC participated in multiple immune-related pathways, and may participate in the immune response of wildness training giant pandas by regulating adaptive immune cells. A large number of DMGs enriched in GO terms may also be related to the regulation of immune activation during wildness training of giant pandas. Promoter differentially methylation analysis identified 1,199 genes with differential methylation at promoter regions. Genes with low methylation level at promoter regions and high expression such as, CCL5, P2Y13, GZMA, ANP32A, VWF, MYOZ1, NME7, MRPS31 and TPM1 were important in environmental adaptation for wildness training giant pandas. The methylation and expression patterns of these genes indicated that wildness training giant pandas have strong immunity, blood coagulation, athletic abilities and disease resistance. The adaptive response of giant pandas undergoing wildness training may be regulated by their negatively related promoter methylation. We are the first to describe the DNA methylation profile of giant panda blood tissue and our results indicated methylation modification is involved in the adaptation of captive giant pandas when undergoing wildness training. Our study also provided potential monitoring indicators for the successful reintroduction of valuable and threatened animals to the wild.

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.

The red panda is an endangered forest species distributed on the edge of the Qinghai Tibet Plateau. The species has been conserved in ex-situ in many countries and its survival is threatened by many diseases. Its immune system is vulnerable to age-associated alterations, which accumulate and result in a progressive deterioration that leads to an increased incidence of diseases. We identified 2,219 differentially expressed genes (DEGs) between geriatric (11–16 years) and adult individuals (4–8 years), and 1690 DEGs between adults and juveniles (1 year). The gene expression and functional annotation results showed that the innate immunity of red pandas increases significantly in geriatric individuals, whereas its change remains unclear when comparing adults and juveniles. We found that the adaptive immunity of red pandas first increased and then decreased with age. We identified CXCR3, BLNK, and CCR4 as the hub genes in the age-related protein–protein interaction network, which showed their central role in age-related immune changes. Many DNA repair genes were down-regulated in geriatric red pandas, suggesting that the DNA repair ability of the blood tissue in geriatric red pandas is significantly reduced. The significantly up-regulated TLR5 in geriatric individuals also suggests the possibility of enhancing the vaccination immune response by incorporating flagellin, which could be used to address decreased vaccine responses caused by age-related declines in immune system function. This work provides an insight into gene expression changes associated with aging and paves the way for effective disease prevention and treatment strategies for red pandas in the future.

Abstract Size changes in brain and brain regions along altitudinal gradients provide insight into the trade‐off between energetic expenditure and cognitive capacity. We investigated the brain size variations of the Asiatic Toad (Bufo gargarizans) across altitudes from 700 m to 3,200 m. A total of 325 individuals from 11 sites and two transects were sampled. To reduce confounding factors, all sampling sites within each transect were within a maximum distance of 85 km and an altitudinal difference close to 2,000 m. Brains were dissected, and five regions were both measured directly and with 3D CT scan. There is a significant negative correlation between the relative whole‐brain volume (to snout‐vent length) and altitude. Furthermore, the relative volumes (to whole‐brain volume) of optic tectum and cerebellum also decrease along the altitudinal gradients, while the telencephalon increases its relative volume along the gradients. Therefore, our results are mostly consistent with the expensive brain hypothesis and the functional constraint hypothesis. We suggest that most current hypotheses are not mutually exclusive and data supporting one hypothesis are often partially consistent with others. More studies on mechanisms are needed to explain the brain size evolution in natural populations.

The giant panda (Ailuropoda melanoleuca) is the epitome of a flagship species for wildlife conservation and also an ideal model of adaptive evolution. As an obligate bamboo feeder, the giant panda relies on the olfaction for food recognition. The number of olfactory receptor (OR) genes and the rate of pseudogenes are the main factors affecting the olfactory ability of animals. In this study, we used the chromosome-level genome of the giant panda to identify OR genes and compared the genome sequences of OR genes with five other Ursidae species (spectacled bear (Tremarctos ornatus), American black bear (Ursus americanus), brown bear (Ursus arctos), polar bear (Ursus maritimus) and Asian black bear (Ursus thibetanus)). The giant panda had 639 OR genes, including 408 functional genes, 94 partial OR genes and 137 pseudogenes. Among them, 222 OR genes were detected and distributed on 18 chromosomes, and chromosome 8 had the most OR genes. A total of 448, 617, 582, 521 and 792 OR genes were identified in the spectacled bear, American black bear, brown bear, polar bear and Asian black bear, respectively. Clustering analysis based on the OR protein sequences of the six species showed that the OR genes distributed in 69 families and 438 subfamilies based on sequence similarity, and the six mammals shared 72 OR gene subfamilies, while the giant panda had 31 unique OR gene subfamilies (containing 35 genes). Among the 35 genes, there are 10 genes clustered into 8 clusters with 10 known human OR genes (OR8J3, OR51I1, OR10AC1, OR1S2, OR1S1, OR51S1, OR4M1, OR4M2, OR51T1 and OR5W2). However, the kind of odor molecules can be recognized by the 10 known human OR genes separately, which needs further research. The phylogenetic tree showed that 345 (about 84.56%) functional OR genes were clustered as Class-II, while only 63 (about 15.44%) functional OR genes were clustered as Class-I, which required further and more in-depth research. The potential odor specificity of some giant panda OR genes was identified through the similarity to human protein sequences. Sequences similar to OR2B1, OR10G3, OR11H6 and OR11H7P were giant panda-specific lacking, which may be related to the transformation and specialization from carnivore to herbivore of the giant panda. Since our reference to flavoring agents comes from human research, the possible flavoring agents from giant panda-specific OR genes need further investigation. Moreover, the conserved motifs of OR genes were highly conserved in Ursidae species. This systematic study of OR genes in the giant panda will provide a solid foundation for further research on the olfactory function and variation of the giant panda.

Abstract Background The giant panda (Ailuropoda melanoleuca) is a threatened endemic Chinese species and a flagship species of national and global conservation concern. Life history theory proposes that reproduction and immunity can be mutually constraining and interrelated. Knowledge of immunity changes of male giant pandas during the breeding season is limited. Results Here, we researched peripheral blood gene expression profiles associated with immunity. Thirteen captive giant pandas, ranging from 9 to 11 years old, were divided into two groups based on their reproductive status. We identified 318 up-regulated DEGs and 43 down-regulated DEGs, which were enriched in 87 GO terms and 6 KEGG pathways. Additionally, we obtained 45 immune-related genes with altered expression, mostly up-regulated, and identified four hub genes HSPA4, SUGT1, SOD1, and IL1B in PPI analysis. These 45 genes were related to pattern recognition receptors, autophagy, peroxisome, proteasome, natural killer cell, antigen processing and presentation. SUGT1 and IL1B were related to pattern recognition receptors. HSP90AA1 was the most up-regulated gene and is a member of heat shock protein 90 family. HSP90 contributes to the translocation of extracellular antigen. KLRD1 encodes CD94, whose complex is an inhibitor of the cytotoxic activity of NK cells, was down-regulated. IGIP, which has the capability of inducing IgA production by B cells, was down-regulated, suggesting low concentration of IgA in male giant pandas. Our results suggest that most immune-related genes were up-regulated and more related to innate immune than adaptive immune. Conclusions Our results indicated that breeding male giant pandas presented an immunoenhancement in innate immunity, enhanced antigen presentation and processing in cellular immunity compared to non-breeding males. The humoral immunity of male giant pandas may show a tendency to decrease during the breeding season. This study will provide a foundation for further studies of immunity and reproduction in male giant pandas.

Abstract The blood pheasant (Ithaginis cruentus), the only species in the genus Ithaginis, lives in an extremely inhospitable high‐altitude environment, coping with hypoxia and ultraviolet (UV) radiation. To further investigate the phylogeny of Phasianidae species based on complete genomes and understand the molecular genetic mechanisms of the high‐altitude adaptation of the blood pheasant, we de novo assembled and annotated the complete genome of the blood pheasant. The blood pheasant genome size is 1.04 Gb with scaffold N50 of 10.88 Mb. We identified 109.92 Mb (10.62%) repetitive elements, 279,037 perfect microsatellites, and 17,209 protein‐coding genes. The phylogenetic tree of Phasianidae based on whole genomes revealed three highly supported major clades with the blood pheasant included in the “erectile clade.” Comparative genomics analysis showed that many genes were positively selected in the blood pheasant, which was associated with response to hypoxia and/or UV radiation. More importantly, among these positively selected genes (PSGs) which were related to high‐altitude adaptation, sixteen PSGs had blood pheasant‐specific missense mutations. Our data and analysis lay solid foundation to the study of Phasianidae phylogeny and provided new insights into the potential adaptation mechanisms to the high altitude employed by the blood pheasant.

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.

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

Abstract The gut microbiomes of the host are large and complex communities, which helps to maintain homeostasis, improves digestive efficiency, and promotes the development of the immune system. The small mammals distributed in Sichuan Province are the most popular species for biodiversity research in Southwest China. However, the effects of different diets on the structure and function of the gut microbial community of these small mammals are poorly understood. In this study, whole‐metagenome shotgun sequencing has been used to analyze the composition and functional structures of the gut microbiota of seven small mammals in Laojunshan National Nature Reserve, Sichuan Province, China. Taxonomic classification revealed that the most abundant phyla in the gut of seven small mammals were Bacteroides, Proteobacteria, and Firmicutes. Moreover, Hafnia, Lactobacillus, and Yersinia were the most abundant genus in the gut microbiomes of these seven species. At the functional level, we annotated a series of KEGG functional pathways, six Cazy categories, and 46,163 AROs in the gut microbiomes of the seven species. Comparative analysis found that the difference in the gut microbiomes between the Soricidea and Muridae concentrated on the increase in the F/B (Firmicutes/Bacteroides) ratio in the Soricidea group, probably driven by the high‐fat and ‐calorie digestive requirements due to their insectivorous diet. The comparative functional profiling revealed that functions related to metabolism and carbohydrates were significantly more abundant in Muridae group, which may be attributed to their high carbohydrate digestion requirements caused by their herbivorous diet. These data suggested that different diets in the host may play an important role in shaping the gut microbiota, and lay the foundation for teasing apart the influences of heritable and environmental factors on the evolution of gut microbial communities.

Microsatellites (SSRs) are widely distributed in the genomes of organisms and are an important genetic basis for genome evolution and phenotypic adaptation. Although the distribution patterns of microsatellites have been investigated in many phylogenetic lineages, they remain unclear within the morphologically and physiologically diverse avian clades. Here, based on high-quality chromosome-level genomes, we examined the microsatellite distribution patterns for 53 birds from 16 orders. The results demonstrated that each type of SSR had the same ratio between taxa. For example, the frequency of imperfect SSRs (I-SSRs) was 69.90–84.61%, while perfect SSRs (P-SSRs) were 14.86–28.13% and compound SSRs (C-SSRs) were 0.39–2.24%. Mononucleotide SSRs were dominant for perfect SSRs (32.66–76.48%) in most bird species (98.11%), and A(n) was the most abundant repeat motifs of P-SSRs in all birds (5.42–68.22%). Our study further confirmed that the abundance and diversity of microsatellites were less effected by evolutionary history but its length. The number of P-SSRs decreased with increasing repeat times, and longer P-SSRs motifs had a higher variability coefficient of the repeat copy number and lower diversity, indicating that longer motifs tended to have more stable preferences in avian genomes. We also found that P-SSRs were mainly distributed at the gene ends, and the functional annotation for these genes demonstrated that they were related to signal transduction and cellular process. In conclusion, our research provided avian SSR distribution patterns, which will help to explore the genetic basis for phenotypic diversity in birds.

Sperm whales are a marine mammal famous for the aromatic substance, the ambergris, produced from its colon. Little is known about the biological processes of ambergris production, and this study aims to investigate the genetic mechanism of ambergris production in the sperm whale based on its chromosome-level genome. Comparative genomics analyses found 1207 expanded gene families and 321 positive selected genes (PSGs) in the sperm whale, and functional enrichment analyses suggested revelatory pathways and terms related to the metabolism of steroids, terpenoids, and aldosterone, as well as microbiota interaction and immune network in the intestine. Furthermore, two sperm-whale-specific missense mutations (Tyr393His and Leu567Val) were detected in the PSG LIPE, which has been reported to play vital roles in lipid and cholesterol metabolism. In total, 46 CYP genes and 22 HSD genes were annotated, and then mapped to sperm whale chromosomes. Furthermore, phylogenetic analysis of CYP genes in six mammals found that CYP2E1, CYP51A and CYP8 subfamilies exhibited relative expansion in the sperm whale. Our results could help understand the genetic mechanism of ambergris production, and further reveal the convergent evolution pattern among animals that produce similar odorants.

Lepidoptera (moths and butterflies) are widely distributed in the world, but high-level phylogeny in Lepidoptera remains uncertain. More mitochondrial genome (mitogenome) data can help to conduct comprehensive analysis and construct a robust phylogenetic tree. Here, we sequenced and annotated 17 complete moth mitogenomes and made comparative analysis with other moths. The gene order of trnM-trnI-trnQ in 17 moths was different from trnI-trnQ-trnM of ancestral insects. The number, type, and order of genes were consistent with reported moths. The length of newly sequenced complete mitogenomes ranged from 14,231 bp of Rhagastis albomarginatus to 15,756 bp of Numenes albofascia. These moth mitogenomes were typically with high A+T contents varied from 76.0% to 81.7% and exhibited negative GC skews. Among 13 protein coding genes (PCGs), some unusual initiations and terminations were found in part of newly sequenced moth mitogenomes. Three conserved gene-overlapping regions and one conserved intergenic region were detected among 17 mitogenomes. The phylogenetic relationship of major superfamilies in Macroheterocera was as follows: (Bombycoidea + Lasiocampoidea) + ((Drepanoidea + Geometroidea) + Noctuoidea)), which was different from previous studies. Moreover, the topology of Noctuoidea as (Notodontidae + (Erebidae + Noctuidae)) was supported by high Bayesian posterior probabilities (BPP = 1.0) and bootstrapping values (BSV = 100). This study greatly enriched the mitogenome database of moth and strengthened the high-level phylogenetic relationships of Lepidoptera.

Colour plays a key role in animal social communication including as an indicator of individual quality. Using spectrophotometry, we examined colour variation in the throat and venter of the crocodile lizard (Shinisaurus crocodilurus), an endangered species native to southern China and northern Vietnam. We detected two broad colour variants, individuals with and without red, for each body region and each sex. A cluster analysis of spectral colour measurements (hue, chroma, luminance) revealed discrete throat and ventral morphs when measured in a single snapshot in time. However, photographic evidence revealed that the amount of red relative to body size increased as they got older. Individuals with red were equally likely to be male or female and throat colour was unrelated to ventral colour. Therefore, it is premature to claim that crocodile lizards have discrete colour morphs. We used visual modelling to show that the throat and venter were easily discriminable to a lizard visual system, suggesting they function in social communication. We also asked whether colour variation signalled individual quality. Females with red throats had greater bite force while males with red throats were older. In addition, females with red venters had larger heads. We also detected differences in morphology linked to colour. Females with red throats had slender bodies and longer tails, while individuals lacking red on their throats were stouter and had shorter tails. Finally, throat and ventral colour were unrelated to reproductive output (litter size and mass) in females. Males with greater ventral luminance contrast sired offspring from litters with greater litter mass (including stillborns), while males with greater ventral chromatic contrast sired offspring whose collective live mass (excluding stillborns) was greater. Males with greater luminance contrast also sired more live offspring (excluding stillborns). Collectively, these results suggest that male ventral colour signals individual quality in males. Conservation initiatives should take colour variation into account when planning future captive breeding and release programs for this endangered species.

Moniezia expansa (M. expansa) parasitizes the small intestine of sheep and causes inhibited growth and development or even death. Being globally distributed, it causes considerable economic losses to the animal husbandry industry. Here, using Illumina, PacBio and BioNano techniques, we obtain a high-quality genome assembly of M. expansa, which has a total length of 142 Mb, a scaffold N50 length of 7.27 Mb and 8,104 coding genes. M. expansa has a very high body fat content and a specific type of fatty acid metabolism. It cannot synthesize any lipids due to the loss of some key genes involved in fatty acid synthesis, and it may can metabolize most lipids via the relatively complete fatty acid β-oxidation pathway. The M. expansa genome encodes multiple lipid transporters and lipid binding proteins that enable the utilization of lipids in the host intestinal fluid. Although many of its systems are degraded (with the loss of homeobox genes), its reproductive system is well developed. PL10, AGO, Nanos and Pumilio compose a reproductive stem cell regulatory network. The results suggest that the high body lipid content of M. expansa provides an energy source supporting the high fecundity of this parasite. Our study provides insight into host interaction, adaptation, nutrient acquisition, strobilization, and reproduction in this parasite and this is also the first genome published in Anoplocephalidae.

Abstract The American cockroach (Periplaneta americana) is a globally invasive pest that can cause significant economic loss and threaten human health. Although it is abundant and lives in close proximity to humans, few studies have investigated the genetic diversity of P. americana. Our study analyzed 1,053 P. americana and other Periplaneta species' samples from different locations in China and the United States. A traditional tree‐based method using 17 unique mitochondrial COI haplotypes of P. americana and 20 haplotypes of the other Periplaneta species accurately identified P. americana with a barcoding threshold of 5.1%. To identify the population genetic structure of P. americana, we investigated wingless gene and pooled them with obtained mtDNA data for a combined analysis. Although the genetic diversity of the USA group was relatively higher than the China group, the number of haplotypes and alleles of both groups was small. The analysis of molecular variance (AMOVA), intraspecific phylogeny, and haplotype networks indicated that P. americana had very little global genetic differentiation. The weak geographic genetic structure might reflect the human‐mediated dispersal of P. americana. Despite no apparent phylogeographic assignment of mtDNA and nuclear lineages was observed in both BI trees, the integrated COI sequence data identified four distinct P. americana haplotype groups, showing four ancient maternal lineages of P. americana in China and the United States.

Abstract Background Microbial communities are susceptible to the public goods dilemma, whereby individuals can gain an advantage within a group by utilizing, but not sharing the cost of producing, public goods. In bacteria, the development of quorum sensing (QS) can establish a cooperation system in a population by coordinating the production of costly and sharable extracellular products (public goods). Cooperators with intact QS system and robust ability in producing public goods are vulnerable to being undermined by QS-deficient defectors that escape from QS but benefit from the cooperation of others. Although microorganisms have evolved several mechanisms to resist cheating invasion in the public goods game, it is not clear why cooperators frequently coexist with defectors and how they form a relatively stable equilibrium during evolution. Results We show that in Pseudomonas aeruginosa, QS-directed social cooperation can select a conditional defection strategy prior to the emergence of QS-mutant defectors, depending on resource availability. Conditional defectors represent a QS-inactive state of wild type (cooperator) individual and can invade QS-activated cooperators by adopting a cheating strategy, and then revert to cooperating when there are abundant nutrient supplies irrespective of the exploitation of QS-mutant defector. Our mathematical modeling further demonstrates that the incorporation of conditional defection strategy into the framework of iterated public goods game with sound punishment mechanism can lead to the coexistence of cooperator, conditional defector, and defector in a rock-paper-scissors dynamics. Conclusions These findings highlight the importance of behavioral heterogeneity in stabilizing the population structure and provide a potential reasonable explanation for the maintenance and evolution of cooperation in microbial communities.

Abstract Background The Hainan Partridge (Arborophila ardens, Phasianidae, Galliformes) is an endemic species of Hainan Island, China, and it is classified as globally vulnerable species. There are at least 16 species in genus Arborophila and no genome sequence is available. Methods The whole genome of Hainan Partridge was de novo sequenced (with shotgun approach on the Illumina 2000 platform) and assembled. Results The genome size of Arborophila ardens is about 1.05 Gb with a high N50 scaffold length of 8.28 Mb and it is the first high quality genome announced in Arborophila genus. About 9.19% of the genome was identified as repeat sequences and about 5.88 million heterozygous SNPs were detected. A total of 17,376 protein-coding genes were predicted and their functions were annotated. The genome comparison between Hainan Partridge and Red Junglefowl (Gallus gallus) demonstrated a conserved genome structure. The phylogenetic analysis indicated that the Hainan Partridge possessed a basal phylogenetic position in Phasianidae and it was most likely derived from a common ancestor approximately 36.8 million years ago (Mya). We found that the Hainan Partridge population had experienced bottleneck and its effective population decreased from about 1,040,000 individuals 1.5 Mya to about 200,000 individuals 0.2 Mya, and then recovered to about 460,000 individuals. The number of 1:1 orthologous genes that were predicted to have undergone positive selection in the Hainan Partridge was 504 and some environmental adaptation related categories, such as response to ultraviolet radiation were represented in GO distribution analysis. Conclusions We announced the first high quality genome in Arborophila genus and it will be a valuable genomic resource for the further studies such as evolution, adaption, conservation, not only on Hainan Partridge but also on Arborophila or Phasianidae species.

Trueperella pyogenes (T. pyogenes) is a crucial opportunistic pathogen normally causing mastitis, abscesses and pneumonia in economically important ruminants. Although only one commercial vaccine of T. pyogenes is currently obtainable, its immunoprotective effect is limited. Pyolysin (PLO) is the most predominant virulence factor highly expressed in T. pyogenes and is an excellent target for the development of novel vaccines against T. pyogenes. In this study, we designed a heterologous prime-boost vaccination scheme combining a DNA vaccine pVAX1-PLO and a subunit vaccine His-PLO to maximize host responses in mice. Humoral and cellular immune responses and protective effects were evaluated in mice to compare the immunogenicity induced by different immunization schemes. Compared to the PBS-control group, in vivo immunization results showed that better immune responses of mice immunized with the pVAX1-PLO plasmids and His-PLO proteins were induced. The residual bacterial burdens from the liver and peritoneal fluid were remarkably decreased in the immunized mice compared with the PBS group. Notably, the heterologous prime-boost vaccination groups significantly enhanced host humoral and cellular immune responses and protected mice from different virulent T. pyogenes strains infection. Conclusively, this study provides a favorable strategy for the further development of next-generation vaccines against T. pyogenes infections.

Mitophagy has recently been implicated in bacterial infection but the underlying mechanism remains largely unknown. Here, we uncover a role of microRNA-302/367 cluster in regulating mitophagy and its associated host response against bacterial infection. We demonstrate that miR-302/367 cluster expression was significantly increased after Pseudomonas aeruginosa infection. Enhanced expression of miR-302/367 cluster accelerated the mitophagic response in macrophages, thus increasing clearance of invading P. aeruginosa by regulating production of reactive oxygen species (ROS), while application of miR-302/367 cluster inhibitors decreased bacterial clearance. Blocking mitophagy with siRNA against mitophagy receptor prohibitin 2 (PHB2) reduced the effect of miR-302/367 cluster on induction of mitophagy and its-associated P. aeruginosa elimination. In addition, we found that miR-302/367 cluster also increased bacterial clearance in mouse model. Mechanistically, we illustrate that miR-302/367 cluster binds to the 3′-untranslated region of nuclear factor kappa B (NF-κB), a negative regulator of mitophagy, accelerated the process of mitophagy by inhibiting NF-κB. Furthermore, inhibition of NF-κB in macrophages attenuated the ROS or cytokines production and may reduce cell injury by P. aeruginosa infection to maintain cellular homeostasis. Collectively, our findings elucidate that miR-302/367 cluster functions as potent regulators in mitophagy-mediated P. aeruginosa elimination and pinpoint an unexpected functional link between miRNAs and mitophagy.

Abstract Background Trueperella pyogenes is a commensal and opportunistic pathogen that normally causes mastitis, liver abscesses and pneumonia of economically important livestock. To develop efficacious and potent vaccine against T. pyogenes, chimeric gene DNA vaccines were constructed and encapsulated in chitosan nanoparticles (pPCFN-CpG-CS-NPs). Results The pPCFN-CpG-CS-NPs consists of the plo, cbpA, fimA, and nanH gene of T. pyogenes and CpG ODN1826. It was produced with good morphology, high stability, a mean diameter of 93.58 nm, and a zeta potential of + 5.27 mV. Additionally, chitosan encapsulation was confirmed to protect the DNA plasmid from DNase I digestion. The immunofluorescence assay indicated that the four-chimeric gene could synchronously express in HEK293T cells and maintain good bioactivity. Compared to the mice immunized with the control plasmid, in vivo immunization showed that mice immunized with the pPCFN-CpG-CS-NPs had better immune responses, and release of the plasmid DNA was prolonged. Importantly, immunization with pPCFN-CpG-CS-NPs could significantly protect mice from highly virulent T. pyogenes TP7 infection. Conclusions This study indicates that chitosan-DNA nanoparticles are potent immunization candidates against T. pyogenes infection and provides strategies for the further development of novel vaccines encapsulated in chitosan nanoparticles.

Abstract Background Insect mitochondrial genomes (mitogenomes) exhibit high diversity in some lineages. The gene rearrangement and large intergenic spacer (IGS) have been reported in several Coleopteran species, although very little is known about mitogenomes of Meloidae. Results We determined complete or nearly complete mitogenomes of seven meloid species. The circular genomes encode 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and two ribosomal RNAs (rRNAs), and contain a control region, with gene arrangement identical to the ancestral type for insects. The evolutionary rates of all PCGs indicate that their evolution is based on purifying selection. The comparison of tRNA secondary structures indicates diverse substitution patterns in Meloidae. Remarkably, all mitogenomes of the three studied Hycleus species contain two large intergenic spacers (IGSs). IGS1 is located between trnW and trnC, including a 9 bp consensus motif. IGS2 is located between trnS2 (UCN) and nad1, containing discontinuous repeats of a pentanucleotide motif and two 18-bp repeat units in both ends. To date, IGS2 is found only in genera Hycleus across all published Coleopteran mitogenomes. The duplication/random loss model and slipped-strand mispairing are proposed as evolutionary mechanisms for the two IGSs (IGS1, IGS2). The phylogenetic analyses using MrBayes, RAxML, and PhyloBayes methods based on nucleotide and amino acid datasets of 13 PCGs from all published mitogenomes of Tenebrionoids, consistently recover the monophylies of Meloidae and Tenebrionidae. Within Meloidae, the genus Lytta clusters with Epicauta rather than with Mylabris. Although data collected thus far could not resolve the phylogenetic relationships within Meloidae, this study will assist in future mapping of the Meloidae phylogeny. Conclusions This study presents mitogenomes of seven meloid beetles. New mitogenomes retain the genomic architecture of the Coleopteran ancestor, but contain two IGSs in the three studied Hycleus species. Comparative analyses of two IGSs suggest that their evolutionary mechanisms are duplication/random loss model and slipped-strand mispairing.

In this paper, we determined and described the complete mitochondrial genome of Robin Accentor (Prunella rubeculoides), the first complete mitogenome. The complete genome of P. rubeculoides was 16,796 bp in length and contained 13 protein-coding genes, 22 transfer RNA genes, two ribosome RNA genes, and one non-coding control region. The overall base composition of the mitochondrial DNA was 29.7% for A, 23.7% for T, 15.6% for G, and 31.0% for C, with a GC content of 46.6%. This information of P. rubeculoides mitogenome is significance for phylogenetic studies of the family Prunellidea.

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals, and its conservation has significant ecosystem and cultural service value. Cytosine DNA methylation (5mC) is a stable epigenetic modification to the genome and has multiple functions such as gene regulation. However, DNA methylome of giant panda and its function have not been reported as of yet. Bisulfite sequencing was performed on a 4-day-old male giant panda’s brain, liver and pancreatic tissues. We found that the whole genome methylation level was about 0.05% based on reads normalization and mitochondrial DNA was not methylated. Three tissues showed similar methylation tendency in the protein-coding genes of their genomes, but the brain genome had a higher count of methylated genes. We obtained 467 and 1,013 different methylation regions (DMR) genes in brain vs. pancreas and liver, while only 260 DMR genes were obtained in liver vs pancreas. Some lncRNA were also DMR genes, indicating that methylation may affect biological processes by regulating other epigenetic factors. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis indicated that low methylated promoter, high methylated promoter and DMR genes were enriched at some important and tissue-specific items and pathways, like neurogenesis, metabolism and immunity. DNA methylation may drive or maintain tissue specificity and organic functions and it could be a crucial regulating factor for the development of newborn cubs. Our study offers the first insight into giant panda’s DNA methylome, laying a foundation for further exploration of the giant panda’s epigenetics.

Microorganisms evolve kinds of elaborate interaction models that can form relatively stable communities in a wide range of ecosystems. It is recognized that the spatial genetic structure of microbes in surface-attached environments lays a good foundation for the persistence of polymicrobial communities in adverse conditions. However, the interacting dynamics of microbes in facilitating the formation and stabilization of community structure still remains elusive. In this study, we identify a hawk-dove game like interspecific relationship between the two Gram-negative opportunistic pathogens Pseudomonas aeruginosa and Klebsiella pneumoniae, which naturally coexist in insect gut and can cocolonize human tissues. Specifically, although P. aeruginosa had significant competitive advantage over cocultured K. pneumoniae on solid medium with rich nutrient factors, K. pneumoniae could resist the suppression of P. aeruginosa by enhancing the expression of membrane transporters induced by the extracellular metabolites of P. aeruginosa. By contrast, under the condition that K. pneumoniae had a growth advantage but P. aeruginosa met a metabolic burden in producing quorum-sensing-controlled extracellular products, the frequency of K. pneumoniae would be slightly higher than P. aeruginosa during the coexistence because K. pneumoniae was also capable of exploiting the extracellular metabolite from P. aeruginosa. In addition, P. aeruginosa quorum-sensing variant could reap benefits from K. pneumoniae in turn and reach a relatively stable two species equilibrium. These findings provide an explanation for the formation and maintenance of polymicrobial communities in different spatially structured environments, and thus may contribute to understanding the complex interspecific interactions of microbes in local communities and shed new light on the development of social microbiology.

It is reported that a wide range of bacterial infections are polymicrobial, and the members in a local microcommunity can influence the growth of neighbors through physical and chemical interactions. Pseudomonas aeruginosa is an important opportunistic pathogen that normally causes a variety of acute and chronic infections, and clinical evidences suggest that P. aeruginosa can be frequently coisolated with other pathogens from the patients with chronic infections. However, the interspecific interaction and the coexisting mechanism of P. aeruginosa with coinfecting bacterial species during evolution still remain largely unclear. In this study, the relationships of P. aeruginosa with other Gram-positive (Staphylococcus aureus) and Gram-negative (Klebsiella pneumoniae) are investigated by using a series of on-plate proximity assay, in vitro coevolution assay, and RNA-sequencing. We find that although the development of a quorum-sensing system contributes P. aeruginosa a significant growth advantage to compete with S. aureus and K. pneumoniae, the quorum-sensing regulation of P. aeruginosa will be decreased during evolution and thus provides a basis for the formation of interspecific coexistence. The results of comparative transcriptomic analyses suggest that the persistent survival of S. aureus in the microcommunity has no significant effect on the intracellular transcriptional pattern of P. aeruginosa, while a more detailed competition happens between P. aeruginosa and K. pneumoniae. Specifically, the population of P. aeruginosa with decreased quorum-sensing regulation can still restrict the proportion increase of K. pneumoniae by enhancing the type VI secretion system-elicited cell aggressivity during further coevolution. These findings provide a general explanation for the formation of a dynamic stable microcommunity consisting of more than two bacterial species, and may contribute to the development of population biology and clinical therapy.

The blue-winged minla (Minla cyanouroptera) belongs to Timaliidae, which is a near threatened species in the IUCN Red List. In this study, the complete mitochondrial genome of Minla cyanouroptera was sequenced and characterized. The genome is 17,862 bases in length (GenBank accession no. MK940810). The overall base composition of M. cyanouroptera mitogenome is 30.9% for A, 31.5% for C, 13.8% for G, and 23.8% for T. The nucleotide sequence data of 13 protein-coding genes of 12 Passeriformes species were used for phylogenetic analyses. Trees constructed by using Bayesian and maximum-likelihood methods demonstrated that Minla ignotincta was closest to M. cyanouroptera.

Mitochondrial genome sequences are valuable resources for systematics and conservation biology studies. In this paper, we present the complete mitogenome of Aquila nipalensis which was 18,450 bp in length. The gene content and arrangement were typical for avian mtDNA. The overall A + T content of was 54.1%, and the AT skew was calculated as 0.12 for the complete mitogenome of A. nipalensis. The maximum-likelihood (ML) tree based on the concatenated 12 protein-coding genes (PCGs) revealed the basal phylogenetic position of A. nipalensis in Aquila.

In our research, the complete mitochondrial genome of Culicicapa ceylonensise was determined and described, which was the first complete mitogenome reported in the Culicicapa. The complete mitochondrial genome sequence of C. ceylonensise was 16,851 bp in length and contains 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and one control region (D-loop). The overall base composition of the mitochondrial DNA was 29.8% for A, 23.5% for T, 31.8% for C, and 14.9% for G, and the percentage of GC content was 46.7%. The complete mitochondrial genome information of C. ceylonensis will contribute to the phylogenetic studies of the Passeriformes.

Horornis fortipes has been listed on the IUCN red list of endangered species - Least Concern (LC). In this study, we sequenced the complete mitochondrial genome of H. fortipes by the polymerase chain reaction (PCR) method. The complete mitochondrial genome of H. fortipes is a circular molecule of 16,908 bp in total length, and it contains 13 typical vertebrate protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes (12srRNA, 16srRNA), and 2 control regions. This is the fist complete mitochondrial genome reported in Cettiidae. Phylogenetic analysis showed that Cettiidae diverged earlier than other family except Pycnonotidae, Acrocephalidae, and Pnoepygidae among 15 families within Passeriformes.

Historically, people made a mistake on species identification on musk deer from Shanxi population. In order to further confirm the species of this population, we sequenced its mitochondrial genome and conducted the species identification. Moreover, the genetic relationships among the sister species of musk deer were also investigated. The results showed that the species of musk deer from the Shanxi population was Moschus berezovskii, and M. anhuiensis was an independent species that had a very close relationship with M. berezovskii. Additionally, M. anhuiensis and M. berezovskii might have experienced a recent history of reproductive isolation.

The complete mitochondrial genome of Otus sunia was 17,835 bp in length and its overall base composition was estimated to be 31.0% for A, 23.3% for T, 32.4% for C, and 13.3% for G. It consisted of the typical structure of 13 PCGs, 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and 1 control region. The phylogenetic analysis demonstrated that O. sunia was more closely related to O. scops than O. bakkamoena.

All the genes of Paradoxornis gularis and Niltava davidi complete mitochondrial genomes encoded on the H-strand, apart from one PCG (ND6) and eight tRNAs (tRNAGln, tRNAAla, tRNAAsn, tRNACys, tRNATyr, tRNASer(UCN), tRNAPro, and tRNAGlu). Strikingly large tandem repeats were found in the control region (CR) of P. gularis mitogenome. The phylogenetic analysis demonstrated that P. gularis possessed a basal phylogenetic position within the genus Paradoxornis, and N. davidi was more closely related to Cyanoptila cyanomelana than other Muscicapidae species.

Complete mitochondrial genomes (mitogenomes) of two cockroach species, Periplaneta australasiae and Neostylopyga rhombifolia, 15,605 bp and 15,711 bp in length, respectively, were determined. As reported for other cockroach mitogenomes, the two mitogenomes possessed typical ancestral insect mitogenome gene composition and arrangement. Only several small intergenic spacers were found: one, which was common in all sequenced cockroach mitogenomes except for the genus Cryptocercus, was between tRNA-Ser (UCN) and ND1 and contained a 7bp highly conserved motif (WACTTAA). Three different types of short tandem repeats in the N. rhombifolia control region (CR) were observed. The homologous alignments of these tandem repeats with other six cockroach mitogenome CRs revealed a low similarity. Three conserved sequence blocks (CSB) were detected in both cockroach mitochondrial CRs. CSB1 was specific for blattinine mitogenomes and was highly conserved with 95% similarity, speculating that this block was a possible molecular synapomorphy for this subfamily. CSB3 located nearby downstream of CSB1 and has more variations within blattinine mitogenomes compared with CSB1. The CSB3 was capable of forming stable stem-loop structure with a small T-stretch in the loop portion. We assessed the influence of four datasets and two inference methods on topology within Orthopteroidea. All genes excluding the third codon positions of PCGs could generate more stable topology, and higher posterior probabilities than bootstrap values were presented at some branch nodes. The phylogenetic analysis with different datasets and analytical methods supported the monophyly of Dictyoptera and supported strongly the proposal that Isoptera should be classified as a family (Termitidae) of the Blattaria. Specifically, Shelfordella lateralis was inserted in the clade Periplaneta. Considering the K2P genetic distance, morphological characters, and the phylogenetic trees, we suggested that S. lateralis should be placed in the genus Periplaneta.

The lineage of the Bar-backed Partridge (Arborophila brunneopectus) was investigated to determine the phylogenetic relationships within Arborophila as the species is centrally distributed within an area covered by the distributions of 22 South-east Asian hill partridge species. The complete mitochondrial genome (mitogenome) of A. brunneopectus was determined and compared with four other hill partridge species mitogenomes. NADH subunit genes are radical in hill partridge mitogenomes and contain the most potential positive selective sites around where variable sites are abundant. Together with 44 other mitogenomes of closely related species, we reconstructed highly resolved phylogenetic trees using maximum likelihood (ML) and Bayesian inference (BI) analyses and calculated the divergence and dispersal history of Arborophila using combined datasets composed of their 13-protein coding sequences. Arborophila is reportedly be the oldest group in Phasianidae whose ancestors probably originated in Asia. A. rufipectus shares a closer relationship with A. ardens and A. brunneopectus compared to A. gingica and A. rufogularis, and such relationships were supported and profiled by NADH dehydrogenase subunit 5 (ND5). The intragenus divergence of all five Arborophila species occurred in the Miocene (16.84~5.69 Mya) when there were periods of climate cooling. We propose that these cooling events in the Miocene forced hill partridges from higher to lower altitudes, which led to geographic isolation and speciation. We demonstrated that the apparently deleterious +1 frameshift mutation in NADH dehydrogenase subunit 3 (ND3) found in all Arborophila is an ancient trait that has been eliminated in some younger lineages, such as Passeriformes. It is unclear of the biological advantages of this elimination for the relevant taxa and this requires further investigation.