Orthologous Groups Research Articles

Whole-genome analysis showed the promotion of genetic diversity and coevolution in Staphylococcus aureus lytic bacteriophages and their hosts mediated by prophages via worldwide recombination events.

Prophages as a part of Staphylococcus aureus genome contribute to the genetic diversity as well as survival strategies of their host. Some S. aureus prophages also have an imminent risk of host cell lysis and become a lytic phage. Nonetheless, interactions among S. aureus prophages, lytic phages, and their hosts, as well as the genetic diversity of S. aureus prophages, remain unclear. We identified 579 intact and 1,389 incomplete prophages in the genomes of 493 S. aureus isolates obtained from the NCBI database. The structural diversity and gene content of intact and incomplete prophages were investigated and compared with 188 lytic phages. Mosaic structure comparison, ortholog group clustering, phylogenetic analysis, and recombination network analysis were performed to estimate genetic relatedness among S. aureus intact prophages, incomplete prophages, and lytic phages. The intact and incomplete prophages harbored 148 and 522 distinct mosaic structures, respectively. The major difference between lytic phages and prophages was the lack of functional modules and genes. Compared to the lytic phages, both the S. aureus intact and incomplete prophages harbored multiple antimicrobial resistance (AMR) and virulence factor (VF) genes. Several functional modules of lytic phages 3_AJ_2017 and 23MRA shared more than 99% nucleotide sequence identity with S. aureus intact (ST20130943_p1 and UTSW_ MRSA_55_ip3) and incomplete prophages (SA3_LAU_ip3 and MRSA_FKTN_ip4); other modules showed little nucleotide sequence similarity. Ortholog and phylogenetic analyses revealed a common gene pool shared between the prophages and lytic Siphoviridae phages. Moreover, most shared sequences existed within intact (43428/137294, 31.6%) and incomplete prophages (41248/137294, 30.0%). Therefore, the maintenance or loss of functional modules in intact and incomplete prophages is key to balance the costs and benefits of large prophages harboring various AMR and VF genes in the bacterial host. The shared identical functional modules between S. aureus lytic phages and prophages are likely to result in the exchange, acquisition, and loss of functional modules, and therefore contribute to their genetic diversity. Moreover, constant recombination events within prophages globally were responsible for the coevolution of lytic phages and their bacterial hosts.

Genome-wide analysis of bromodomain gene family in Arabidopsis and rice.

The bromodomain-containing proteins (BRD-proteins) belongs to family of 'epigenetic mark readers', integral to epigenetic regulation. The BRD-members contain a conserved 'bromodomain' (BRD/BRD-fold: interacts with acetylated-lysine in histones), and several additional domains, making them structurally/functionally diverse. Like animals, plants also contain multiple Brd-homologs, however the extent of their diversity and impact of molecular events (genomic duplications, alternative splicing, AS) therein, is relatively less explored. The present genome-wide analysis of Brd-gene families of Arabidopsis thaliana and Oryza sativa showed extensive diversity in structure of genes/proteins, regulatory elements, expression pattern, domains/motifs, and the bromodomain (w.r.t. length, sequence, location) among the Brd-members. Orthology analysis identified thirteen ortholog groups (OGs), three paralog groups (PGs) and four singleton members (STs). While more than 40% Brd-genes were affected by genomic duplication events in both plants, AS-events affected 60% A. thaliana and 41% O. sativa genes. These molecular events affected various regions (promoters, untranslated regions, exons) of different Brd-members with potential impact on expression and/or structure-function characteristics. RNA-Seq data analysis indicated differences in tissue-specificity and stress response of Brd-members. Analysis by RT-qPCR revealed differential abundance and salt stress response of duplicate A. thaliana and O. sativa Brd-genes. Further analysis of AtBrd gene, AtBrdPG1b showed salinity-induced modulation of splicing pattern. Bromodomain (BRD)-region based phylogenetic analysis placed the A. thaliana and O. sativa homologs into clusters/sub-clusters, mostly consistent with ortholog/paralog groups. The bromodomain-region displayed several conserved signatures in key BRD-fold elements (α-helices, loops), along with variations (1-20 sites) and indels among the BRD-duplicates. Homology modeling and superposition identified structural variations in BRD-folds of divergent and duplicate BRD-members, which might affect their interaction with the chromatin histones, and associated functions. The study also showed contribution of various duplication events in Brd-gene family expansion among diverse plants, including several monocot and dicot plant species.

Pan-Chromosome and Comparative Analysis of Agrobacterium fabrum Reveal Important Traits Concerning the Genetic Diversity, Evolutionary Dynamics, and Niche Adaptation of the Species.

Agrobacterium fabrum has been critical for the development of plant genetic engineering and agricultural biotechnology due to its ability to transform eukaryotic cells. However, the gene composition, evolutionary dynamics, and niche adaptation of this species is still unknown. Therefore, we established a comparative genomic analysis based on a pan-chromosome data set to evaluate the genetic diversity of A. fabrum. Here, 25 A. fabrum genomes were selected for analysis by core genome phylogeny combined with the average nucleotide identity (ANI), amino acid identity (AAI), and in silico DNA-DNA hybridization (DDH) values. An open pan-genome of A. fabrum exhibits genetic diversity with variable accessorial genes as evidenced by a consensus pan-genome of 12 representative genomes. The genomic plasticity of A. fabrum is apparent in its putative sequences for mobile genetic elements (MGEs), limited horizontal gene transfer barriers, and potentially horizontally transferred genes. The evolutionary constraints and functional enrichment in the pan-chromosome were measured by the Clusters of Orthologous Groups (COG) categories using eggNOG-mapper software, and the nonsynonymous/synonymous rate ratio (dN/dS) was determined using HYPHY software. Comparative analysis revealed significant differences in the functional enrichment and the degree of purifying selection between the core genome and non-core genome. We demonstrate that the core gene families undergo stronger purifying selection but have a significant bias to contain one or more positively selected sites. Furthermore, although they shared similar genetic diversity, we observed significant differences between chromosome 1 (Chr I) and the chromid in their functional features and evolutionary constraints. We demonstrate that putative genetic elements responsible for plant infection, ecological adaptation, and speciation represent the core genome, highlighting their importance in the adaptation of A. fabrum to plant-related niches. Our pan-chromosome analysis of A. fabrum provides comprehensive insights into the genetic properties, evolutionary patterns, and niche adaptation of the species. IMPORTANCE Agrobacterium spp. live in diverse plant-associated niches such as soil, the rhizosphere, and vegetation, which are challenged by multiple stressors such as diverse energy sources, plant defenses, and microbial competition. They have evolved the ability to utilize diverse resources, escape plant defenses, and defeat competitors. However, the underlying genetic diversity and evolutionary dynamics of Agrobacterium spp. remain unexplored. We examined the phylogeny and pan-genome of A. fabrum to define intraspecies evolutionary relationships. Our results indicate an open pan-genome and numerous MGEs and horizontally transferred genes among A. fabrum genomes, reflecting the flexibility of the chromosomes and the potential for genetic exchange. Furthermore, we observed significant differences in the functional features and evolutionary constraints between the core and accessory genomes and between Chr I and the chromid, respectively.

Phylogenetic placement of whole genome duplications in yeasts through quantitative analysis of hierarchical orthologous groups

Phylogenetic placement of whole genome duplications in yeasts through quantitative analysis of hierarchical orthologous groups

Intestinal Microbiota of Anser fabalis Wintering in Two Lakes in the Middle and Lower Yangtze River Floodplain.

The intestinal microbiota of migratory birds participate in the life activities of the host and are affected by external environmental factors. The difference in habitat environment provides diversity in external environmental selection pressure for the same overwintering waterfowl, which may be reflected in their intestinal microbiota. Caizi lake and Shengjin Lake in the Middle and Lower Yangtze River Floodplain are the main habitats for migratory waterfowl in winter, especially the Anser fabalis (A. fabalis). It is important to explore the changes in intestinal microbiota composition and function of A. fabalis in the early overwintering period to clarify the effect of habitat size and protection status on intestinal microbiota. In this study, the composition and structural characteristics of the intestinal microbiota of A. fabalis in Shengjin Lake (SL) and Caizi Lake (CL) were preliminarily explored in order to obtain data for the migratory birds. In both SL and CL groups, 16S rRNA amplicon sequencing analysis showed that Firmicutes was the dominant bacterial phylum, but the relative abundance showed significant differences. Lactobacillus was the most abundant genus in both SL and CL groups. At the species level, the abundance of L. aviaries was the highest, with a relative abundance in both SL and CL groups of more than 34%. When comparing the average relative abundance of the 15 most abundant genera, it was found that Subdoligranulum, Exiguobacterium, and Terrisporobacter had higher abundances in the intestinal microbiota of CL A. fabalis, while Streptococcus and Rothia had higher abundances in the intestinal microbiota of SL A. fabalis. There was only a positive correlation between Bacteroidota and Proteobacteria in the intestinal microbiota flora of SL A. fabalis, and the species were closely related. At the same time, there were positive and negative correlations between Firmicutes and Actinomycetes. However, CL is mainly associated with a positive correlation between Firmicutes and Actinomycetes, and there are also a small number of connections between Firmicutes. PICRUSt1 prediction analysis revealed that the Clusters of Orthologous Groups (COG) functions of SL and CL involve energy production and transformation, amino acid transport and metabolism, carbohydrate transport and metabolism, and transcription. Understanding the changes in intestinal microbiota in Aves during the overwintering period is of great importance to explore the adaptation mechanism of migratory Aves to the overwintering environment. This work provides basic data for an A. fabalis intestinal microbiota study.

OrthoPhy: A Program to Construct Ortholog Data Sets Using Taxonomic Information.

Species phylogenetic trees represent the evolutionary processes of organisms, and they are fundamental in evolutionary research. Therefore, new methods have been developed to obtain more reliable species phylogenetic trees. A highly reliable method is the construction of an ortholog data set based on sequence information of genes, which is then used to infer the species phylogenetic tree. However, although methods for constructing an ortholog data set for species phylogenetic analysis have been developed, they cannot remove some paralogs, which is necessary for reliable species phylogenetic inference. To address the limitations of current methods, we developed OrthoPhy, a program that excludes paralogs and constructs highly accurate ortholog data sets using taxonomic information dividing analyzed species into monophyletic groups. OrthoPhy can remove paralogs, detecting inconsistencies between taxonomic information and phylogenetic trees of candidate ortholog groups clustered by sequence similarity. Performance tests using evolutionary simulated sequences and real sequences of 40 bacteria revealed that the precision of ortholog inference by OrthoPhy is higher than that of existing programs. Additionally, the phylogenetic analysis of species was more accurate when performed using ortholog data sets constructed by OrthoPhy than that performed using data sets constructed by existing programs. Furthermore, we performed a benchmark test of the Quest for Orthologs using real sequence data and found that the concordance rate between the phylogenetic trees of orthologs inferred by OrthoPhy and those of species was higher than the rates obtained by other ortholog inference programs. Therefore, ortholog data sets constructed using OrthoPhy enabled a more accurate phylogenetic analysis of species than those constructed using the existing programs, and OrthoPhy can be used for the phylogenetic analysis of species even for distantly related species that have experienced many evolutionary events.

Genome-Wide Identification of G Protein-Coupled Receptors in Ciliated Eukaryotes

G protein-coupled receptors (GPCRs) are the largest family of transmembrane receptors and play important roles in many physiological processes. As a representative group of protozoa, ciliates represent the highest stage of eukaryotic cell differentiation and evolution in terms of their reproductive mode, two-state karyotype, and extremely diverse cytogenesis patterns. GPCRs have been poorly reported in ciliates. In this study, we identified 492 GPCRs in 24 ciliates. Using the existing classification system for animals, GPCRs in ciliates can be assigned to four families, including families A, B, E, and F. Most (377 members) belong to family A. The number of GPCRs is extremely different in different ciliates; the Heterotrichea ciliates usually have more GPCRs than other ciliates. Parasitic or symbiotic ciliates usually have only a few GPCRs. Gene/genome duplication events seem to play important roles in the expansion of the GPCR superfamily in ciliates. GPCRs in ciliates displayed seven typical domain organizations. GPCRs in an ortholog group are common and conserved in all ciliates. The gene expression analysis of the members in this conserved ortholog group in the model ciliate, Tetrahymena thermophila, suggested that these GPCRs play important roles in the life cycle of ciliates. In summary, this study provides the first comprehensive genome-wide identification of GPCRs in ciliates, improving our understanding of the evolution and function of GPCR in ciliates.

Sphingobium lignivorans sp. nov., isolated from river sediment downstream of a paper mill.

A bacterial isolate, B1D3AT, was isolated from river sediment collected from the Hiwassee River near Calhoun, TN, by enrichment culturing with a model 5-5' lignin dimer, dehydrodivanillate, as its sole carbon source. B1D3AT was also shown to utilize several model lignin-derived monomers and dimers as sole carbon sources in a variety of minimal media. Cells were Gram-stain-negative, aerobic, motile, rod-shaped and formed yellow/cream-coloured colonies on rich agar. Optimal growth occurred at 30 °C, pH 7-8, and in the absence of NaCl. The major fatty acids of B1D3AT were C18 : 1 ω7c and C17 : 1 ω6c. The predominant hydroxy fatty acids were C14 : 0 2-OH and C15 : 0 2-OH. The polar lipid profile consisted of a mixture of phosphatidylethanolamine, phosphatidylglycerol, diphosphatidylglycerol, phosphatidyldimethylethanolamine and sphingoglycolipid. B1D3AT contained spermidine as the only major polyamine. The major isoprenoid quinone was Q-10 with minor amounts of Q-9 and Q-11. The genomic DNA G+C content of B1D3AT was 65.6 mol%. Phylogenetic analyses based on 16S rRNA gene sequences and coding sequences of 49 core, universal genes defined by Clusters of Orthologous Groups gene families indicated that B1D3AT was a member of the genus Sphingobium. B1D3AT was most closely related to Sphingobium sp. SYK-6, with a 100 % 16S rRNA gene sequence similarity. B1D3AT showed 78.1-89.9 % average nucleotide identity and 19.5-22.2% digital DNA-DNA hybridization identity with other type strains from the genus Sphingobium. On the basis of phenotypic and genotypic properties and phylogenetic inference, strain B1D3AT should be classified as representing a novel species of the genus Sphingobium, for which the name Sphingobium lignivorans sp. nov. is proposed. The type strain is strain B1D3AT (ATCC TSD-279T=DSM 111877T).

InParanoiDB 9: Ortholog Groups for Protein Domains and Full-Length Proteins

Prediction of orthologs is an important bioinformatics pursuit that is frequently used for inferring protein function and evolutionary analyses. The InParanoid database is a well known resource of ortholog predictions between a wide variety of organisms. Although orthologs have historically been inferred at the level of full-length protein sequences, many proteins consist of several independent protein domains that may be orthologous to domains in other proteins in a way that differs from the full-length protein case. To be able to capture all types of orthologous relations, conventional full-length protein orthologs can be complemented with orthologs inferred at the domain level. We here present InParanoiDB 9, covering 640 species and providing orthologs for both protein domains and full-length proteins. InParanoiDB 9 was built using the faster InParanoid-DIAMOND algorithm for orthology analysis, as well as Domainoid and Pfam to infer orthologous domains. InParanoiDB 9 is based on proteomes from 447 eukaryotes, 158 bacteria and 35 archaea, and includes over one billion predicted ortholog groups. A new website has been built for the database, providing multiple search options as well as visualization of groups of orthologs and orthologous domains. This release constitutes a major upgrade of the InParanoid database in terms of the number of species as well as the new capability to operate on the domain level. InParanoiDB 9 is available at https://inparanoidb.sbc.su.se/.

Genome-Wide Identification and Evolutionary Analyses of SrfA Operon Genes in Bacillus.

A variety of secondary metabolites contributing to plant growth are synthesized by bacterial nonribosomal peptide synthases (NRPSs). Among them, the NRPS biosynthesis of surfactin is regulated by the SrfA operon. To explore the molecular mechanism for the diversity of surfactins produced by bacteria within the genus Bacillus, we performed a genome-wide identification study focused on three critical genes of the SrfA operon-SrfAA, SrfAB and SrfAC-from 999 Bacillus genomes (belonging to 47 species). Gene family clustering indicated the three genes can be divided into 66 orthologous groups (gene families), of which a majority comprised members of multiple genes (e.g., OG0000009 had members of all three SrfAA, SrfAB and SrfAC genes), indicating high sequence similarity among the three genes. Phylogenetic analyses also found that none of the three genes formed monophyletic groups, but were usually arranged in a mixed manner, suggesting the close evolutionary relationship among the three genes. Considering the module structure of the three genes, we propose that self-duplication, especially tandem duplications, might have contributed to the initial establishment of the entire SrfA operon, and further gene fusion and recombination as well as accumulated mutations might have continuously shaped the different functional roles of SrfAA, SrfAB and SrfAC. Overall, this study provides novel insight into metabolic gene clusters and operon evolution in bacteria.

Assessment of Genomic and Metabolic Characteristics of Cholesterol-Reducing and GABA Producer Limosilactobacillus fermentum AGA52 Isolated from Lactic Acid Fermented Shalgam Based on "In Silico" and "In Vitro" Approaches.

This study aimed to characterize the genomic and metabolic properties of a novel Lb. fermentum strain AGA52 which was isolated from a lactic acid fermented beverage called "shalgam." The genome size of AGA52 was 2,001,184bp, which is predicted to carry 2024 genes, including 50 tRNAs, 3 rRNAs, 3 ncRNAs, 15 CRISPR repeats, 14 CRISPR spacers, and 1 CRISPR array. The genome has a GC content of 51.82% including 95 predicted pseudogenes, 56 complete or partial transposases, and 2 intact prophages. The similarity of the clusters of orthologous groups (COG) was analyzed by comparison with the other Lb. fermentum strains. The detected resistome on the genome of AGA52 was found to be intrinsic originated. Besides, it has been determined that AGA52 has an obligate heterofermentative carbohydrate metabolism due to the absence of the 1-phosphofructokinase (pfK) enzyme. Furthermore, the strain is found to have a better antioxidant capacity and to be tolerant to gastrointestinal simulated conditions. It was also observed that the AGA52 has antimicrobial activity against Yersinia enterocolitica ATCC9610, Bacillus cereus ATCC33019, Salmonella enterica sv. Typhimurium, Escherichia coli O157:h7 ATCC43897, Listeria monocytogenes ATCC7644, Klebsiella pneumoniae ATCC13883, and Proteus vulgaris ATCC8427. Additionally, AGA52 exhibited 42.74 ± 4.82% adherence to HT29 cells. Cholesterol assimilation (33.9 ± 0.005%) and GABA production capacities were also confirmed by "in silico" and "in vitro." Overall, the investigation of genomic and metabolic features of the AGA52 revealed that is a potential psychobiotic and probiotic dietary supplement candidate and can bring functional benefits to the host.

Genome Characterization and Development of Real-Time PCR Assays for Ditylenchus Dipsaci and D. Weischeri.

The stem and bulb nematode Ditylenchus dipsaci is a destructive nematode pest on many crops and is internationally quarantined in many countries, whereas Ditylenchus weischeri, only known to infect a weed plant (Cirsium arvense), is an unregulated nematode species with no known economic importance. In this study, we used comparative genomics to identify multiple gene regions and developed novel real-time PCR assays for the detection of D. dipsaci and D. weischeri. We sequenced the genomes of two mixed-stage nematode populations of D. dipsaci and two mixed-stage nematode populations of D. weischeri. The assembled genomes of D. dipsaci were 228.2 Mb and 239.5 Mb, and the genomes of D. weischeri were 177.0 Mb and 196.3 Mb. Depending on the species, 21,403-27,365 gene models were predicted. Using orthologous group analysis, single-copy and species-specific genes were identified. Primers and probes were designed targeting two species-specific genes in each species. The assays detected as low as 12 pg of DNA from the target species, or as few as five nematodes, with a Cq of 31 cycles or less. Our study provides genome data for two additional D. dipsaci isolates and two D. weischeri isolates, and four new and validated molecular assays to be used for rapid detection and identification of the two species.

Antifatigue and microbiome reshaping effects of yak bone collagen peptides on Balb/c mice

As the main by-product of yak (Bos grunniens) processing, yak bone was discarded or simply processed for low-value production. The development of yak bone collagen peptides (YBP) has shown great potential. In the study, the 1 g/kg· BW peptides were supplied (the group of YBP) or not supplie the control group) in the Balb/c mice diet for 18 weeks. Peptide supplementation prolonged the swimming time of mice by 25% and 37% for twice tests, respectively, accompanied by the increased content of glucose (Glu), lactate dehydrogenase (LDH), superoxide dismutase (SOD), glutathione (GSH), and decreased content of lactic acid (LA), blood ammonia (BA) in blood. YBP intervention promoted the relative abundance of Bacteroides, Blautia, Muribaculum, etc. The proportion of amino acid and carbohydrate transport and metabolism increased among all the predicted functions from the orthologous groups of proteins (COG) database. Besides that, after YBP intervention, many up-regulated proteins in the liver were involved in the protein and carbohydrate metabolism-related pathways, and the proteins in lipid metabolism-related pathways were down-regulated. This study firstly clarified the antifatigue and gut microbiome reshaping effects of YBP on Balb/c mice and provided the theoretical basis for future development.

Comparative Transcriptomics Reveal Metabolic Rather than Genetic Control of Divergent Antioxidant Metabolism in the Primary Root Elongation Zone of Water-Stressed Cotton and Maize.

Under water stress, the primary root elongation zones of cotton and maize exhibit both conserved and divergent metabolic responses, including variations in sulfur and antioxidant metabolism. To explore the relative importance of metabolic and genetic controls of these responses for each species, and the extent to which responses are mediated by similar gene expression networks within the framework of ortholog groups, comparative transcriptomics analyses were conducted under conditions of equivalent tissue water stress. Ortholog analysis revealed that 86% of the transcriptome response to water stress was phylogenetically unrelated between cotton and maize. Elevated transcript abundances for genes involved in abscisic acid (ABA) biosynthesis and signaling, as well as key enzymes that enable osmotic adjustment, were conserved between the species. In contrast, antioxidant responses, at least with regard to glutathione metabolism and anti-oxidative enzymes, did not exhibit such a transcript abundance adaptive signature. In particular, previously characterized differential responses of the glutathione and sulfur metabolic pathways between cotton and maize were not evident in the transcriptomic responses. The findings indicate that the antioxidant response in both species results from a metabolic acclimation to water stress, and thus represents an example of water stress-related metabolic plasticity.

Label-free LC-MS/MS proteomics analyses reveal CLIC1 as a predictive biomarker for bladder cancer staging and prognosis.

Bladder cancer (BC) is a significant carcinoma of the urinary system that has a high incidence of morbidity and death owing to the challenges in accurately identifying people with early-stage BC and the lack of effective treatment options for those with advanced BC. Thus, there is a need to define new markers of prognosis and prediction. In this study, we have performed a comprehensive proteomics experiment by label-free quantitative proteomics to compare the proteome changes in the serum of normal people and bladder cancer patients-the successful quantification of 2064 Quantifiable proteins in total. A quantitative analysis was conducted to determine the extent of changes in protein species' relative intensity and reproducibility. There were 43 upregulated proteins and 36 downregulated proteins discovered in non-muscle invasive bladder cancer and normal individuals. Sixty-four of these proteins were elevated, and 51 were downregulated in muscle-invasive and non-muscle-invasive bladder cancer, respectively. Functional roles of differentially expressed proteins were annotated using Gene Ontology (GO) and Clusters of Orthologous Groups of Proteins (COG). To analyze the functions and pathways enriched by differentially expressed proteins, GO enrichment analysis, protein domain analysis, and KEGG pathway analysis were performed. The proteome differences were examined and visualized using radar plots, heat maps, bubble plots, and Venn diagrams. As a result of combining the Venn diagram with protein-protein interactions (PPIs), Chloride intracellular channel 1 (CLIC1) was identified as the primary protein. Using the Gene Set Cancer Analysis (GSCA) website, the influence of CLIC1 on immune infiltration was analyzed. A negative correlation between CD8 naive and CLIC1 levels was found. For validation, immunohistochemical (IHC), qPCR, and western blotting (WB) were performed.Further, we found that CLIC1 was associated with a poor prognosis of bladder cancer in survival analysis. Our research screened CLIC1 as a tumor-promoting protein in bladder cancer for the first time using serum mass spectrometry. And CLIC1 associated with tumor stage, and immune infiltrate. The prognostic biomarker and therapeutic target CLIC1 may be new for bladder cancer patients.

Screening and isolation of cold-adapted cellulose degrading bacterium: A candidate for straw degradation and De novo genome sequencing analysis.

Degradation of crop straw in natural environment has been a bottleneck. There has been a recent increase in the exploration of cold-adapted microorganisms as they can solve the problem of corn straw degradation under low temperatures and offer new alternatives for the sustainable development of agriculture. The study was conducted in low-temperature (10°C) and high-efficiency cellulose-degrading bacteria were screened using carboxymethyl cellulose (CMC) selection medium and subjected to genome sequencing by the third-generation Pacbio Sequl and the second-generation Illumina Novaseq platform, and their cellulase activity was detected by 3,5-dinitrosalicylic acid (DNS) method. The results showed that the low-temperature (10°C) and high-efficiency cellulose-degrading bacterium Bacillus subtilis K1 was 4,060,823 bp in genome size, containing 4,213 genes, with 3,665, 3,656, 2,755, 3,240, 1,261, 3,336 and 4,003 genes annotated in the non-redundant protein sequence database (NR), Pfam, clusters of orthologous groups of proteins (COGs), Genome Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Annotation databases, respectively. In addition, a large number of lignocellulose degradation-related genes were annotated in the genome. The cellulose activity of B. subtilis K1 was higher, exhibiting the highest activity of endo-β-glucanase (24.69 U/ml), exo-β-glucanase (1.72 U/ml) and β-glucosaccharase (1.14 U/ml). It was found that through adding cold-adapted cellulose-degrading bacteriaK1 in the corn straw composting under 6°C (ambient temperature), the average temperature of straw composting was 58.7°C, and higher 86.7% as compared to control. The HA/FA was higher 94.02% than the control and the lignocellulose degradation rate was lower 18.01-41.39% than the control. The results provide a theoretical basis for clarifying the degradation potential of cold-adapted cellulose-degrading bacteria and improving the cellulose degradation efficiency.

Exploring Evolutionary Relationships within Neodermata Using Putative Orthologous Groups of Proteins, with Emphasis on Peptidases

The phylogenetic relationships within Neodermata were examined based on putative orthologous groups of proteins (OGPs) from 11 species of Monogenea, Trematoda, and Cestoda. The dataset included OGPs from BUSCO and OMA. Additionally, peptidases were identified and evaluated as phylogenetic markers. Phylogenies were inferred using the maximum likelihood method. A network analysis and a hierarchical grouping analysis of the principal components (HCPC) of orthologous groups of peptidases were performed. The phylogenetic analyses showed the monopisthocotylean monogeneans as the sister-group of cestodes, and the polyopisthocotylean monogeneans as the sister-group of trematodes. However, the sister-group relationship between Monopisthocotylea and Cestoda was not statistically well supported. The network analysis and HCPC also showed a cluster formed by polyopisthocotyleans and trematodes. The present study supports the non-monophyly of Monogenea. An analysis of mutation rates indicated that secreted peptidases and inhibitors, and those with multiple copies, are under positive selection pressure, which could explain the expansion of some families such as C01, C19, I02, and S01. Whilst not definitive, our study presents another point of view in the discussion of the evolution of Neodermata, and we hope that our data drive further discussion and debate on this intriguing topic.

Altered gut microbiota composition in children and their caregivers infected with the SARS-CoV-2 Omicron variant.

Gut microbiota alterations have been implicated in the pathogenesis of coronavirus disease 2019 (COVID-19). This study aimed to explore gut microbiota changes in a prospective cohort of COVID-19 children and their asymptomatic caregivers infected with the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) Omicron variant. A total of 186 participants, including 59 COVID-19 children, 50 asymptomatic adult caregivers, 52 healthy children (HC), and 25 healthy adults (HA), were recruited between 15 April and 31 May 2022. The gut microbiota composition was determined by 16S rRNA gene sequencing in fecal samples collected from the participants. Gut microbiota functional profiling was performed by using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) software. The gut microbiota analysis of beta diversity revealed that the fecal microbial community of COVID-19 children remained far distantly related to HC. The relative abundances of the phyla Actinobacteria and Firmicutes were decreased, whereas Bacteroidetes, Proteobacteria, and Verrucomicrobiota were increased in COVID-19 children. Feces from COVID-19 children exhibited notably lower abundances of the genera Blautia, Bifidobacterium, Fusicatenibacter, Streptococcus, and Romboutsia and higher abundances of the genera Prevotella, Lachnoclostridium, Escherichia-Shigella, and Bacteroides than those from HC. The enterotype distributions of COVID-19 children were characterized by a high prevalence of enterotype Bacteroides. Similar changes in gut microbiota compositions were observed in asymptomatic caregivers. Furthermore, the microbial metabolic activities of KEGG (Kyoto Encyclopedia of Genes and Genomes) and COG (cluster of orthologous groups of proteins) pathways were perturbed in feces from subjects infected with the SARS-CoV-2 Omicron variant. Our data reveal altered gut microbiota compositions in both COVID-19 children and their asymptomatic caregivers infected with the SARS-CoV-2 Omicron variant, which further implicates the critical role of gut microbiota in COVID-19 pathogenesis.

Transcriptome Analysis of Rhododendron liliiflorum H. Lév. Flower Colour Differences

Rhododendron liliiflorum H. Lév., with white outer edges and yellow inner edges of petals, is an ornamental flower that originated in China. In this study, we analysed the white (W) and yellow (Y) parts of R. liliiflorum flowers by RNA sequencing. Then, unigene assembly, unigene annotation, and classification of Eukaryotic Orthologous Groups (KOGs) were performed. Gene ontology (GO) classification and pathway enrichment analysis for unigenes were also conducted. A total of 219,221 transcripts and 180,677 unigenes of R. liliiflorum were obtained from 48.52 Gb of clean reads. Differentially expressed gene (DEG) analysis indicated that 2310 unigenes were upregulated and 3062 were downregulated in W vs. Y. Thirty-six of these DEGs were involved in the flavonoid biosynthesis pathway. Pathway enrichment analysis showed that DEGs were significantly enriched in phenylpropanoid, flavonoid, and isoflavone biosynthesis. The expression of dihydroflavonol-4-reductase (DFR) and chalcone synthase (CHS) may affect differences in R. liliiflorum flower colour. The findings on flavonoid biosynthesis and other related genes in this study will provide guidance for exploring the mechanism of flower colour formation in Rhododendron.

Comparative genomics reveals cellobiose hydrolysis mechanism of Ruminiclostridium thermocellum M3, a cellulosic saccharification bacterium.

The cellulosome of Ruminiclostridium thermocellum was one of the most efficient cellulase systems in nature. However, the product of cellulose degradation by R. thermocellum is cellobiose, which leads to the feedback inhibition of cellulosome, and it limits the R. thermocellum application in the field of cellulosic biomass consolidated bioprocessing (CBP) industry. In a previous study, R. thermocellum M3, which can hydrolyze cellulosic feedstocks into monosaccharides, was isolated from horse manure. In this study, the complete genome of R. thermocellum M3 was sequenced and assembled. The genome of R. thermocellum M3 was compared with the other R. thermocellum to reveal the mechanism of cellulosic saccharification by R. thermocellum M3. In addition, we predicted the key genes for the elimination of feedback inhibition of cellobiose in R. thermocellum. The results indicated that the whole genome sequence of R. thermocellum M3 consisted of 3.6 Mb of chromosomes with a 38.9% of GC%. To be specific, eight gene islands and 271 carbohydrate-active enzyme-encoded proteins were detected. Moreover, the results of gene function annotation showed that 2,071, 2,120, and 1,246 genes were annotated into the Clusters of Orthologous Groups (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively, and most of the genes were involved in carbohydrate metabolism and enzymatic catalysis. Different from other R. thermocellum, strain M3 has three proteins related to β-glucosidase, and the cellobiose hydrolysis was enhanced by the synergy of gene BglA and BglX. Meanwhile, the GH42 family, CBM36 family, and AA8 family might participate in cellobiose degradation.