Gut microbiota is involved in male reproductive function: a review.

Globally, ~8%-12% of couples confront infertility issues, male-related issues being accountable for 50%. This review focuses on the influence of gut microbiota and their metabolites on the male reproductive system from five perspectives: sperm quality, testicular structure, sex hormones, sexual behavior, and probiotic supplementation. To improve sperm quality, gut microbiota can secrete metabolites by themselves or regulate host metabolites. Endotoxemia is a key factor in testicular structure damage that causes orchitis and disrupts the blood-testis barrier (BTB). In addition, the gut microbiota can regulate sex hormone levels by participating in the synthesis of sex hormone-related enzymes directly and participating in the enterohepatic circulation of sex hormones, and affect the hypothalamic-pituitary-testis (HPT) axis. They can also activate areas of the brain that control sexual arousal and behavior through metabolites. Probiotic supplementation can improve male reproductive function. Therefore, the gut microbiota may affect male reproductive function and behavior; however, further research is needed to better understand the mechanisms underlying microbiota-mediated male infertility.

Evaluation of immunoregulation and immunoprotective efficacy of Glaesserella parasuis histidine kinase QseC.

QseC is a membrane sensor kinase that enables bacteria to perceive autoinducers -3, adrenaline, and norepinephrine to initiate downstream gene transcription. In this study, we found that the QseC protein of Glaesserella parasuis can serve as an effective antigen to activate the host's immune response. Therefore, we investigated the immunogenicity and host protective effect of this protein. ELISA and indirect immunofluorescence results showed that QseC protein can induce high titer levels of humoral immunity in mice and regularly generate specific serum antibodies. We used MTS reagents to detect lymphocyte proliferation levels and found that QseC protein can cause splenic lymphocyte proliferation with memory and specificity. Further immunological analysis of the spleen cell supernatant revealed significant upregulation of levels of IL-1ß, IL-4 and IFN-γ in the QseC + adjuvant group. In the mouse challenge experiment, it was found that QseC + adjuvant can provide effective protection. The results of this study demonstrate that QseC protein provides effective protection in a mouse model and has the potential to serve as a candidate antigen for a novel subunit vaccine for further research.

Gut microbiota affects obesity susceptibility in mice through gut metabolites.

Introduction: It is well-known that different populations and animals, even experimental animals with the same rearing conditions, differ in their susceptibility to obesity. The disparity in gut microbiota could potentially account for the variation in susceptibility to obesity. However, the precise impact of gut microbiota on gut metabolites and its subsequent influence on susceptibility to obesity remains uncertain. Methods: In this study, we established obesity-prone (OP) and obesity-resistant (OR) mouse models by High Fat Diet (HFD). Fecal contents of cecum were examined using 16S rDNA sequencing and untargeted metabolomics. Correlation analysis and MIMOSA2 analysis were used to explore the association between gut microbiota and intestinal metabolites. Results: After a HFD, gut microbiota and gut metabolic profiles were significantly different between OP and OR mice. Gut microbiota after a HFD may lead to changes in eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA), a variety of branched fatty acid esters of hydroxy fatty acids (FAHFAs) and a variety of phospholipids to promote obesity. The bacteria g_Akkermansia (Greengene ID: 175696) may contribute to the difference in obesity susceptibility through the synthesis of glycerophosphoryl diester phosphodiesterase (glpQ) to promote choline production and the synthesis of valyl-tRNA synthetase (VARS) which promotes L-Valine degradation. In addition, gut microbiota may affect obesity and obesity susceptibility through histidine metabolism, linoleic acid metabolism and protein digestion and absorption pathways.

Genome analysis of the mpox (formerly monkeypox) virus and characterization of core/variable regions.

Since smallpox was eradicated in 1980, the monkeypox virus (MPXV) has emerged as the most threatening orthopoxvirus in the world. In this study, we conducted a comprehensive analysis of the currently published complete genome sequences of the monkeypox virus. The core/variable regions were identified through core-pan analysis of MPXV. Besides single-nucleotide polymorphisms, our study also revealed that specific genes, multi-copy genes, repeat sequences, and recombination fragments are primarily distributed in the variable region. This result suggests that variable regions are not only more susceptible to single-base mutations, but also to events such as gene loss or gain, as well as recombination. Taken together, our results demonstrate the genomic characteristics of the core/variable regions of MPXV, and contribute to our understanding of the evolution of MPXV.

Comparative genomic analysis of alloherpesviruses: Exploring an available genus/species demarcation proposal and method.

The family Alloherpesviridae contains herpesviruses of fish and amphibians. Due to the significant economic losses to aquaculture that herpesviruses can cause, the primary areas of research interest are concerning their pathogenesis and prevention. Despite alloherpesvirus genomic sequences becoming more widely accessible, methods regarding their genus/species classification are still relatively unexplored. In the present study, the phylogenetic relationships between 40 completely sequenced alloherpesviruses were illustrated by the viral proteomic tree (ViPTree), which was divided into three monophyletic groups, namely Cyprinivirus, Ictalurivirus and Batrachovirus. Additionally, average nucleotide identity (ANI) and average amino acid identity (AAI) analyses were performed across all available sequences and clearly displayed species boundaries with the threshold value of ANI/AAI set at 90%. Subsequently, core-pan analysis uncovered 809 orthogroups and 11 core genes shared by all 40 alloherpesvirus genome sequences. For the former, a 15 percent identity depicts a clear genus boundary; for the latter, 8 of them may be qualified for phylogenetic analysis based on amino acid or nucleic acid sequences after being verified using maximum likelihood (ML) or neighbor-joining (NJ) phylogenetic trees. Finally, although the dot plot analysis was valid for the members within Ictalurivirus, it was unsuccessful for Cyprinivirus and Batrachovirus. Taken together, the comparison of individual methodologies provides a wide range of alternatives for alloherpesviruses classification under various circumstances.

Comprehensive analysis of transcriptome and metabolome analysis reveal new targets of Glaesserella parasuis glucose-specific enzyme IIBC (PtsG).

The ptsG (hpIIBCGlc) gene, belonging to the glucose-specific phosphotransferase system, encodes the bacterial glucose-specific enzyme IIBC. In this study, the effects of a deletion of the ptsG gene were investigated by metabolome and transcriptome analyses. At the transcriptional level, we identified 970 differentially expressed genes between ΔptsG and sc1401 (Padj<0.05) and 2072 co-expressed genes. Among these genes, those involved in methane metabolism, amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, pyruvate metabolism, phosphotransferase system (PTS), biotin metabolism, Two-component system and Terpenoid backbone biosynthesis showed significant changes in the ΔptsG mutant strain. Metabolome analysis revealed that a total of 310 metabolites were identified, including 20 different metabolites (p < 0.05). Among them, 15 metabolites were upregulated and 5 were downregulated in ΔptsG mutant strain. Statistical analysis revealed there were 115 individual metabolites having correlation, of which 89 were positive and 26 negative. These metabolites include amino acids, phosphates, amines, esters, nucleotides, benzoic acid and adenosine, among which amino acids and phosphate metabolites dominate. However, not all of these changes were attributable to changes in mRNA levels and must also be caused by post-transcriptional regulatory processes. The knowledge gained from this lays the foundation for further study on the role of ptsG in the pathogenic process of Glaesserella parasuis (G.parasuis).

Corrigendum: The Microbiota and It's Correlation With Metabolites in the Gut of Mice With Nonalcoholic Fatty Liver Disease.

[This corrects the article DOI: 10.3389/fcimb.2022.870785.].

Comparative genomic analysis reveals new evidence of genus boundary for family Iridoviridae and explores qualified hallmark genes.

Members of the family Iridoviridae (iridovirids) are globally distributed and trigger adverse economic and ecological impacts on aquaculture and wildlife. Iridovirids taxonomy has previously been studied based on a limited number of genomes, but this is not suitable for the current and future virological studies as more iridovirids are emerging. In our study, 57 representative iridovirids genomes were selected from a total of 179 whole genomes available on NCBI. Then 18 core genes were screened out for members of the family Iridoviridae. Average amino acid sequence identity (AAI) analysis indicated that a cut-off value of 70% is more suitable for the current iridovirids genome database than ICTV-defined 50% threshold to better clarify viral genus boundaries. In addition, more subgroups were divided at genus level with the AAI threshold of 70%. This observation was further confirmed by genomic synteny analysis, codon usage preference analysis, genome GC content and length analysis, and phylogenic analysis. According to the pairwise comparison analysis of core genes, 9 hallmark genes were screened out to conduct preliminary identification and investigation at the genus level of iridovirids in a more convenient and economical manner.

ANI analysis of poxvirus genomes reveals its potential application to viral species rank demarcation.

Average nucleotide identity (ANI) is a prominent approach for rapidly classifying archaea and bacteria by recruiting both whole genomic sequences and draft assemblies. To evaluate the feasibility of ANI in virus taxon demarcation, 685 poxviruses were assessed. Prior to the analysis, the fragment length and threshold of the ANI value were optimized as 200 bp and 98 per cent, respectively. After ANI analysis and network visualization, the resulting sixty-one species (ANI species rank) were clustered and largely consistent with the groupings found in National Center for Biotechnology Information Virus [within the International Committee on Taxonomy of Viruses (ICTV) Master Species List]. The species identities of thirty-four other poxviruses (excluded by the ICTV Master Species List) were also identified. Subsequent phylogenetic analysis and Guanine-Cytosine (GC) content comparison done were found to support the ANI analysis. Finally, the BLAST identity of concatenated sequences from previously identified core genes showed 91.8 per cent congruence with ANI analysis at the species rank, thus showing potential as a marker gene for poxviruses classification. Collectively, our results reveal that the ANI analysis may serve as a novel and efficient method for poxviruses demarcation.

The Microbiota and It's Correlation With Metabolites in the Gut of Mice With Nonalcoholic Fatty Liver Disease.

In recent years, nonalcoholic fatty liver disease (NAFLD) has become the most common liver disease in the world. As an important model animal, the characteristics of gut microbiota alteration in mice with NAFLD have been studied but the changes in metabolite abundance in NAFLD mice and how the gut microbiota affects these intestinal metabolites remain unclear. In this experiment, a mouse model for NAFLD was established by a high-fat diet. The use of 16S rDNA technology showed that while there were no significant changes in the alpha diversity in the cecum of NAFLD mice, the beta diversity changed significantly. The abundance of Blautia, Unidentified-Lachnospiraceae, Romboutsia, Faecalibaculum, and Ileibacterium increased significantly in NAFLD mice, while Allobaculum and Enterorhabdus decreased significantly. Amino acids, lipids, bile acids and nucleotide metabolites were among the 167 significantly different metabolites selected. The metabolic pathways of amino acids, SFAs, and bile acids were significantly enhanced, while the metabolic pathways of PUFAs, vitamins, and nucleotides were significantly inhibited. Through correlation and MIMOSA2 analysis, it is suggested that gut microbiota does not affect the changes of lipids and bile acids but can reduce thiamine, pyridoxine, and promote L-phenylalanine and tyramine production. The findings of this study will help us to better understand the relationship between gut microbiota and metabolites in NAFLD.

Impact of intracellular response regulator QseB in quorum sensing regulatory network in a clinical isolate SC1401 of Glaesserella parasuis.

Two component systems (TCS) mediate specific responses to different conditions and/or pressures. In the quorum sensing Glaesserella parasuis (QSE) BC TCS, qseB, as a response regulator, is closely related to the transcriptional regulation of multiple downstream genes. In this study, the effects of qseB gene deletion, which encodes the response regulator of population density sensing in G. parasuis, were studied through biological characteristics and metabolomic analysis. Based on previous research, we further explored the virulence of ΔqseB mutant strains through cell morphology, adhesion and invasion. The ΔqseB mutant and parent strains were sequenced by metabolome and combined with the previous transcriptome sequencing results for joint analysis. This study aims to clarify the regulatory effect of QseB on the virulence of G. parasuis and lay the foundation for revealing the pathogenic mechanism of G. parasuis. We detected 476 different metabolites, of which 30 metabolites (6.3%) had a significant difference in abundance between SC1401 and ΔqseB (p < 0.05). We conducted a comparative analysis of pathway enrichment on the transcriptome and metabolome, and found that the two omics participate in seven metabolic pathways together. The top 10 KEGG pathways with the largest number of genes and metabolites identified in this experiment are ABC transporters, Biosynthesis of secondary metabolites, Cysteine and methionine metabolism, Purine metabolism, Pyrimidine metabolism, Metabolic pathways, and Nicotinate and nicotinamide metabolism. Analysis of metabolome sequencing results showed that differential metabolites were also enriched in metabolic pathways, such as Purine metabolism, cGMP-PKG signaling pathway and cAMP signaling pathway, which were not found in transcriptome sequencing data. The internal coloration of the mutant strain ΔqseB was uneven, and the adhesion and invasion ability of PAM cell lines were significantly reduced. We speculate that QseB may affect the adhesion and invasion ability of Glaesserella parasuis by influencing substance transport and signal transduction.

Deletion of two-component system QseBC weakened virulence of Glaesserella parasuis in a murine acute infection model and adhesion to host cells.

The widespread two-component system (TCS), QseBC, involves vital virulence regulators in Enterobacteriaceae and Pasteurellaceae. Here we studied the function of QseBC in Glaesserella parasuis. A ΔqseBC mutant was constructed using a Glaesserella parasuis serovar 11 clinical strain SC1401 by natural transformation. Immunofluorescence was used to evaluate cellular adhesion, the levels of inflammation and apoptosis. The ability of ΔqseBC and ΔqseC mutant strains to adhere to PAM and MLE-12 cells was significantly reduced. Additionally, by focusing on the clinical signs, H&E, and IFA for inflammation and apoptosis, we found that the ΔqseBC mutant weakened virulence in the murine models. Together, these findings suggest that QseBC plays an important role in the virulence of Glaesserella parasuis.

Landscape of keratinocytes transcriptome alterations in response to Trichophyton mentagrophytes infection.

Dermatophytosis is an intractable superficial fungal infection of keratinized structures, with approximately 20% incidence in humans. Alterations of keratinocytes in the pathogenesis of dermatophytosis at the transcriptome level remain unclear. To understand and characterize such responses, keratinocytes were infected with Trichophyton mentagrophytes. After infection with 1 × 105 conidia/mL T. mentagrophytes for 24 h, the adherence of fungal hyphae to keratinocytes and the damage caused to cell morphology and structure were observed by light microscopy and transmission electron microscopy, respectively. Levels of pro-inflammatory cytokines IL-1α, IL-1ß, TNFα, and IL-8 significantly increased after infection. RNA-seq and bioinformatic analyses revealed that 766 genes were significantly whereas 2207 genes were repressed in the T. mentagrophyte-infected cells. Some of the differentially expressed genes (DEGs) were related to inflammation, immune responses, wound healing, metabolism, and oxidative stress. GO and KEGG pathway enrichment analyses revealed that DEGs and pathways involved in inflammatory response, immune response, and pathogen-induced dysfunction were significantly enriched in the infected cells. Furthermore, gene set enrichment analysis revealed that higher expression gene sets were mainly involved in immune responses, whereas lower expression gene sets were related to cell component organization or biogenesis and transporter activity. Furthermore, protein-protein interaction network and function analyses revealed that JUN, TP53, FOS, MYC, and HSP90AA1 play a key role in immune responses. Overall, our study systematically uncovered the transcriptome-level response of keratinocytes to T. mentagrophyte and provided insights into dermatophytosis treatment.

Genome-wide identification of RING finger genes in flax (Linum usitatissimum) and analyses of their evolution.

BACKGROUND: Flax (Linum usitatissimum) is an important crop for its seed oil and stem fiber. Really Interesting New Gene (RING) finger genes play essential roles in growth, development, and biotic and abiotic stress responses in plants. However, little is known about these genes in flax. METHODS: Here, we performed a systematic genome-wide analysis to identify RING finger genes in flax. RESULTS: We identified 587 RING domains in 574 proteins and classified them into RING-H2 (292), RING-HCa (181), RING-HCb (23), RING-v (53), RING-C2 (31), RING-D (2), RING-S/T (3), and RING-G (2). These proteins were further divided into 45 groups according to domain organization. These genes were located in 15 chromosomes and clustered into three clades according to their phylogenetic relationships. A total of 312 segmental duplicated gene pairs were inferred from 411 RING finger genes, indicating a major contribution of segmental duplications to the RING finger gene family expansion. The non-synonymous/synonymous substitution ratio of the segmentally duplicated gene pairs was less than 1, suggesting that the gene family was under negative selection since duplication. Further, most RING genes in flax were differentially expressed during seed development or in the shoot apex. This study provides useful information for further functional analysis of RING finger genes in flax and to develop gene-derived molecular markers in flax breeding.

Genomic analysis of Poxviridae and exploring qualified gene sequences for phylogenetics.

The members of the Poxviridae family are globally distributed all over the world and can cause infectious diseases. Although genome sequences are publicly available for representative isolates of all genera, studies on the criteria for genome-based classification within the Poxviridae family have rarely been reported. In our study, 60 Poxviridae genomes were re-annotated using Prokka. By using BLAST filtration and MCScanX, synteny and similarity of whole genomic amino acid sequences were visualized. According to the analysis pattern, the Chordopoxvirinae and Entomopoxvirinae subfamilies can be subdivided into five and two categories respectively, which is consistent with the phylogenetic tree constructed based on whole genomic amino acid sequences and Poxvirus core genes. Finally, four genes (Early transcription factor, DNA-directed RNA polymerase, RNA polymerase-associated transcription-specificity factor and DNA-dependent RNA polymerase) were selected from Poxvirus core genes by substitution saturation analysis and phylogenetic tree verification. Phylogenetic trees constructed based on single gene and concatenated sequences of the four selected genes showed that the classification of subgroups was consistent with the phylogenetic trees based on genome. Conclusion: a new method based on the similarity of whole genomic amino acid sequences was proposed for Poxviridae taxon demarcation, and the use of the four selected qualified genes will help make phylogenic identification of newly discovered Poxviridae isolates more convenient and accurate.

The Insights of Genomic Synteny and Codon Usage Preference on Genera Demarcation of Iridoviridae Family.

The members of the family Iridoviridae are large, double-stranded DNA viruses that infect various hosts, including both vertebrates and invertebrates. Although great progress has been made in genomic and phylogenetic analyses, the adequacy of the existing criteria for classification within the Iridoviridae family remains unknown. In this study, we redetermined 23 Iridoviridae core genes by re-annotation, core-pan analysis and local BLASTN search. The phylogenetic tree based on the 23 re-annotated core genes (Maximum Likelihood, ML-Tree) and amino acid sequences (composition vector, CV-Tree) were found to be consistent with previous reports. Furthermore, the information provided by synteny analysis and codon usage preference (relative synonymous codon usage, correspondence analysis, ENC-plot and Neutrality plot) also supports the phylogenetic relationship. Collectively, our results will be conducive to understanding the genera demarcation within the Iridoviridae family based on genomic synteny and component (codon usage preference) and contribute to the existing taxonomy methods for the Iridoviridae family.

Genomic analysis of Ranavirus and exploring alternative genes for phylogenetics.

Ranaviruses can infect both captive and wild cold-blooded vertebrates, leading to significant economic and environmental losses. With the cases of ranavirus infection increasing, many ranavirus genomic sequences were published, but little is known about ranavirus taxonomy on a whole-genome level. In this study, 44 ranaviruses core genes were identified in 32 ranaviruses genome sequences by using PanX. The neighbour-joining phylogenetic trees (NJ-tree) based on 44 ranaviruses core genes and 24 iridoviridae core genes and composition vector phylogenetic tree (CV-Tree) based on whole genome were constructed. The three of phylogenetic trees showed that 32 ranavirus isolates can be divided into 4 different subgroups including SGIV-like, EHNV-like, FV3-like and CMTV-like, and subgroups taxonomic position of three phylogenetic trees were consistent. However, the phylogenetic position of ToRV could not be determined if it belongs to FV3-like or CMTV-like group. Subsequently, we carried out dot plot analysis and confirmed that ToRV should belong to CMTV-like group. Based on dot plot analysis and phylogenetic trees, the taxonomic classification of ranaviruses was confirmed. Finally, four genes which are suitable for the construction of phylogenetic tree were selected from ranavirus core genes by recombination analysis, substitution saturation analysis and single-gene phylogenetic analysis. Phylogenetic tree based on concatenated sequences of the four selected genes showed that the classification of subgroups was identical with three of the phylogenetic trees. Conclusion: Our results confirmed taxonomic identification of ranaviruses; the four selected genes used in phylogenetic analysis will make taxonomic identification more convenient and accurate.

The complete chloroplast genome of Syringa oblata (Oleaceae).

Syringa oblata Lindl. is a popular ornamental shrub with aroma compounds. Here, we sequenced and assembled the complete chloroplast genome of S. oblata. The complete chloroplast genome of S. oblata is 155,648 bp in length, containing a pair of inverted repeated (IRa and IRb) region of 25,732 bp that are separated by a large single copy (LSC) region of 86,247 bp, and a small single copy (SSC) region of 17,937 bp. A total of 132 functional genes were annotated, including 88 protein-coding genes, 36 tRNA genes, and eight rRNA genes. The Neighbour-joining phylogenetic tree based on complete chloroplast genomes suggested that S. oblata is most closely related to S. vulgaris.

Genome-wide identification and evolution of HECT genes in wheat.

BACKGROUND: As an important class of E3 ubiquitin ligases in the ubiquitin proteasome pathway, proteins containing homologous E6-AP carboxyl terminus (HECT) domains are crucial for growth, development, metabolism, and abiotic and biotic stress responses in plants. However, little is known about HECT genes in wheat (Triticum aestivum L.), one of the most important global crops. METHODS: Using a genome-wide analysis of high-quality wheat genome sequences, we identified 25 HECT genes classified into six groups based on the phylogenetic relationship among wheat, rice, and Arabidopsis thaliana. RESULTS: The predicted HECT genes were distributed evenly in 17 of 21 chromosomes of the three wheat subgenomes. Twenty-one of these genes were hypothesized to be segmental duplication genes, indicating that segmental duplication was significantly associated with the expansion of the wheat HECT gene family. The Ka/Ks ratios of the segmental duplication of these genes were less than 1, suggesting purifying selection within the gene family. The expression profile analysis revealed that the 25 wheat HECT genes were differentially expressed in 15 tissues, and genes in Group II, IV, and VI (UPL8, UPL6, UPL3) were highly expressed in roots, stems, and spikes. This study contributes to further the functional analysis of the HECT gene family in wheat.