BMC Genomics最新文献

The recent emergence of Elizabethkingia anophelis as a human pathogen is a major concern for global public health. This organism has the potential to cause severe infections and has inherent antimicrobial resistance. The potential for widespread outbreaks and rapid global spread highlights the critical importance of understanding the biology and transmission dynamics of this infectious agent. We performed a large-scale analysis of available 540 E. anophelis, including one novel strain isolated from raw milk and sequenced in this study. Pan-genome analysis revealed an open and diverse pan-genome in this species, characterized by the presence of many accessory genes. This suggests that the species has a high level of adaptability and can thrive in a variety of environments. Phylogenetic analysis has also revealed a complex population structure, with limited source-lineage correlation. We identified diverse antimicrobial resistance factors, including core-genome and accessory ones often associated with mobile genetic elements within specific lineages. Mobilome analysis revealed a dynamic landscape primarily composed of genetic islands, integrative and conjugative elements, prophage elements, and small portion of plasmids emphasizing a complex mechanism of horizontal gene transfer. Our study underscores the adaptability of E. anophelis, characterized by a diverse range of antimicrobial resistance genes, putative virulence factors, and genes enhancing fitness. This adaptability is also supported by the organism's ability to acquire genetic material through horizontal gene transfer, primarily facilitated by mobile genetic elements such as integrative and conjugative elements (ICEs). The potential for rapid evolution of this emerging pathogen poses a significant challenge to public health efforts.

The quality of Recombination signal sequences (RSSs), location, and genetics of mammalian V, D, and J genes synergistically affect the recombination frequency of genes; however, the specific regulatory mechanism and efficiency have not been elucidated. By taking advantage of single-cell RNA-sequencing (scRNA-seq) and high-throughput sequencing (HTS) to investigate V(D)J rearrangement characteristics in the CDR3 repertoire, we found that the distal and proximal V genes (or J genes) "to D" gene were involved in rearrangement significantly more frequently than the middle V genes (or J genes) in the TRB locus among various species, including Primates (human and rhesus monkey), Rodentia (BALB/c, C57BL/6, and Kunming mice), Artiodactyla (buffalo), and Chiroptera (Rhinolophus affinis). The RSS quality of the V and J genes affected their frequency in rearrangement to varying degrees, especially when the V-RSSs with recombination signal information content (RIC) score < -45 significantly reduced the recombination frequency of the V gene. The V and J genes that were "away from D" had the dual advantages of recombinant structural accessibility and relatively high-quality RSSs, which promoted their preferential utilization in rearrangement. The quality of J-RSSs formed during mammalian evolution was apparently greater than that of V-RSSs, and the D-J distance was obviously shorter than that of V-D, which may be one of the reasons for guaranteeing that the "D-to-J preceding V-to-DJ rule" occurred when rearranged. This study provides a novel perspective on the mechanism and efficiency of V-D-J rearrangement in the mammalian TRB locus, as well as the biased utilization characteristics and application of V and J genes in the initial CDR3 repertoire.

Background: Ophiocordyceps sinensis (O. sinensis) is the dominant bacterium in the asexual stage of Chinese cordyceps, and its growth usually suffers from water stress. Thus, simulating its ecological growth environment is crucial for artificial cultivation. This study aimed to reveal the mechanism underlying the water stress tolerance of Ophiocordyceps sinensis (O. sinensis) by combining metabolomic and transcriptome analyses to identify crucial pathways related to differentially expressed genes (DEGs) and metabolites (DEMs) involved in the response to water stress.

Results: Gene coexpression analysis revealed that many genes related to 'betalain biosynthesis', 'tyrosine metabolism', 'linoleic acid metabolism', 'fructose and mannose metabolism', and 'starch and sucrose metabolism' were highly upregulated after 20d-water stress. Metabolomic analysis revealed that many metabolites regulated by these genes in these metabolic pathways were markedly decreased. On the one hand, we surmised that carbohydrate metabolism and the β-oxidation pathway worked cooperatively to generate enough acyl-CoA and then entered the TCA cycle to provide energy when exposed to water stress. On the other hand, the betalain biosynthesis and tyrosine metabolism pathway might play crucial roles in response to water stress in O. sinensis by enhancing cell osmotic potential and producing osmoregulatory substances (betaine) and antioxidant pigments (eumelanin).

Conclusions: Overall, our findings provide important information for further exploration of the mechanism underlying the water stress tolerance of O. sinensis for the industrialization of artificial cultivation of Chinese cordyceps.

Background: The deep-sea cold seep zone is characterized by high pressure, low temperature, darkness, and oligotrophy. Vesicomyidae clams are the dominant species within this environment, often forming symbiotic relationships with chemosynthetic microbes. Understanding the mechanisms by which Vesicomyidae clams adapt to the cold seep environment is significant. Acetylation modification of lysine is known to play a crucial role in various metabolic processes. Consequently, investigating the role of lysine acetylation in the adaptation of Vesicomyidae clams to deep-sea environments is worthwhile. So, a comparative study of lysine acetylation in cold seep clam Archivesica marissinica and shallow water shellfish Ruditapes philippinarum was conducted.

Results: A total of 539 acetylated proteins were identified with 1634 acetylation sites. Conservative motif enrichment analysis revealed that the motifs -KacR-, -KacT-, and -KacF- were the most conserved. Subsequent gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were conducted on significantly differentially expressed acetylated proteins. The GO enrichment analysis indicated that acetylated proteins are crucial in various biological processes, including cellular response to stimulation, and other cellular processes ( p < 0.05 and false discovery rate (FDR) < 0.25). The results of KEGG enrichment analysis indicated that acetylated proteins are involved in various cellular processes, including tight junction, motor proteins, gap junction, phagosome, cGMP-PKG signaling pathways, endocytosis, glycolysis/gluconeogenesis, among others (p < 0.05 and FDR < 0.25). Notably, a high abundance of lysine acetylation was observed in the glycolysis/glycogenesis pathways, and the acetylation of glyceraldehyde 3-phosphate dehydrogenase might facilitate ATP production. Subsequent investigation into acetylation modifications associated with deep-sea adaptation revealed the specific identification of key acetylated proteins. Among these, the adaptation of cold seep clam hemoglobin and heat shock protein to high hydrostatic pressure and low temperature might involve an increase in acetylation levels. Acetylation of arginine kinase might be related to ATP production and interaction with symbiotic bacteria. Myosin heavy chain (Ama01085) has the most acetylation sites and might improve the actomyosin system stability through acetylation. Further validation is required for the acetylation modification from Vesicomyidae clams.

Conclusion: A novel comparative analysis was undertaken to investigate the acetylation of lysine in Vesicomyidae clams, yielding novel insights into the regulatory role of lysine acetylation in deep-sea organisms. The findings present many potential proteins for further exploration of acetylation functions in cold seep clams and other deep-sea mollusks.

Background: Gram-negative bacteria are the main bacterial pathogens infecting Chinese giant salamanders (Andrias davidianus; CGS) in captivity and the wild, causing substantial economic losses in the CGS industry. However, the molecular mechanisms underlying pathogenesis following infection remain unclear.

Results: Spleen-derived macrophages from healthy CGS were isolated, cultured, and identified using density gradient centrifugation and immunofluorescence. A macrophage transcriptome database was established 0, 6, and 12 h post lipopolysaccharide stimulation using RNA-sequencing. In the final database 76,743 unigenes and 4,698 differentially expressed genes (DEGs) were functionally annotated. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment results showed that DEGs were concentrated in toll-like receptor-nuclear factor kappa B-related immune pathways. Ten DEGs were validated 12 h after lipopolysaccharide (LPS) stimulation. Although the common LPS recognition receptor toll-like receptor 4 was not activated and the key adaptor protein MyD88 showed no significant response, we observed significant up-regulation of the following adaptors: toll/interleukin-1 receptor domain-containing adaptor inducing interferon-β, tumour necrosis factor receptor-associated factor 6, and transforming growth factor-β activated kinase 1, which are located downstream of the non-classical MyD88 pathway.

Conclusions: In contrast to that in other species, macrophage activation in CGS could depend on the non-classical MyD88 pathway in response to bacterial infection. Our study provides insights into the molecular mechanisms regulating CGS antibacterial responses, with implications for disease prevention and understanding immune evolution in amphibians.

Background: Hypersaline habitats, as extreme environments, are a great source of well-adapted organisms with unique properties as they have evolved various strategies to cope with these extreme conditions. Bioinformatics and genomic mining may shed light on evolutionary relationships among them. Therefore, the aim of this study was to assess the biodiversity and especially the strategies evolved within the Idiomarina genus, with the primary focus on the taxonomy and genomic adaptations of two novel strains affiliated with Idiomarina genus isolated from unique environment - brines of two Early Miocene salt deposits.

Results: Both analyzed species belonged to the Idiomarina loihiensis cluster with similarity levels of 16S rRNA gene sequences as high as 99.5% and showed a significant genome size reduction, known characteristic of Idiomarina genomes, though within the genome of Sol25 strain the lowest extent of the carbohydrate utilization genes reduction was observed t among the Idiomarina species. Moreover, the comparative genome analyses indicated that despite both strains being isolated from geographically and geologically similar environments (brines from at least 12 Ma), the species showed higher relatedness to other Idiomarina species than to each other.

Conclusion: The present findings highlighted the importance of genomic data in resolving taxonomic uncertainties and understanding of adaptation strategies of extremophiles. Geographic isolation likely contributed to population divergence of the Idiomarina genus, and the recent study offered insights into biogeographic patterns and allopatric speciation of this bacterial group.

Methylobacterium sp. XJLW converts formaldehyde into methanol and formic acid via a Cannizzaro reaction in response to environmental formaldehyde stress. Methanol is further assimilated without formaldehyde or formic acid formation, whereas formic acid accumulates without undergoing further metabolism. Synthetic biology-based biotransformation of methanol to generate additional products can potentially achieve carbon neutrality. However, practical applications are hampered by limitations such as formaldehyde tolerance. In this study, we aimed to explore the specific mechanism of strain XJLW in response to formaldehyde stress. Thus, a transcriptomic analysis of XJLW under formaldehyde treatment was performed, revealing changes in the expression of specific genes related to one-carbon metabolism. Central metabolic genes were downregulated, whereas metabolic bypass genes were upregulated to maintain methanol assimilation in XJLW's response to formaldehyde treatment. In total, 100 genes potentially related to methyl transfer were identified. The function of only one gene, RS27765, was similar to that of glyA, which encodes a methyltransferase involved in one-carbon metabolism. The double-mutant strain, lacking RS27765 and glyA, lost its ability to grow in methanol, whereas the single-mutant strain, lacking only one of these genes, still grew in methanol. Co-expression of RS27765 and RS31205 (YscQ/HrcQ type III secretion apparatus protein) enabled Escherichia coli BL21 (DE3) to effectively degrade methanol. Using protein sequence analysis and molecular docking, we proposed a model wherein RS27765 is necessary for cell growth by using methanol generated via formaldehyde cannizzaro reaction. This process enables direct assimilation of methanol without producing formaldehyde and formic acid as intermediate metabolites. The RS27765 gene cluster, in conjunction with metabolic bypass genes, constitutes a novel auxiliary pathway facilitating formaldehyde stress tolerance in the strain.

Understanding the syntenic relationships among genomes is crucial to elucidate the genomic mechanisms that drive the evolution of species. The nematode Caenorhabditis is a good model for studying genomic evolution due to the well-established biology of Caenorhabditis elegans and the availability of > 50 genomes in the genus. However, effective alignment of more than ten species in Caenorhabditis has not been conducted before, and there is currently no tool to visualize the synteny of more than two species. In this study, we used Progressive Cactus, a recently developed multigenome aligner, to align the genomes of eleven Caenorhabditis species. Through the progressive alignment, we reconstructed nine ancestral genomes, analyzed the mutational types that cause genomic rearrangement during speciation, and found that insertion and duplication are the major driving forces for genome expansion. Dioecious species appear to expand their genomes more than androdioecious species. We then built an online interactive app called WormSynteny to visualize the syntenic relationship among the eleven species. Users can search the alignment dataset using C. elegans query sequences, construct synteny plots at different genomic scales, and use a set of options to control alignment output and plot presentation. We showcased the use of WormSynteny to visualize the syntenic conservation of one-to-one orthologues among species, tandem and dispersed gene duplication in C. elegans, and the evolution of exon and intron structures. Importantly, the integration of orthogroup information with synteny linkage in WormSynteny allows the easy visualization of conserved genomic blocks and disruptive rearrangement. In conclusion, WormSynteny provides immediate access to the syntenic relationships among the most widely used Caenorhabditis species and can facilitate numerous comparative genomics studies. This pilot study with eleven species also serves as a proof-of-concept to a more comprehensive larger-scale analysis using hundreds of nematode genomes, which is expected to reveal mechanisms that drive genomic evolution in the Nematoda phylum. Finally, the WormSynteny software provides a generalizable solution for visualizing the output of Progressive Cactus with interactive graphics, which would be useful for a broad community of genome researchers.

Background: The ubiquity of sex across eukaryotes, given its high costs, strongly suggests it is evolutionarily advantageous. Asexual lineages can avoid, for example, the risks and energetic costs of recombination, but suffer short-term reductions in adaptive potential and long-term damage to genome integrity. Despite these costs, lichenized fungi have frequently evolved asexual reproduction, likely because it allows the retention of symbiotic algae across generations. The lichenized fungal genus Lepraria is thought to be exclusively asexual, while its sister genus Stereocaulon completes a sexual reproductive cycle. A comparison of sister sexual and asexual clades should shed light on the evolution of asexuality in lichens in general, as well as the apparent long-term maintenance of asexuality in Lepraria, specifically.

Results: In this study, we assembled and annotated representative long-read genomes from the putatively asexual Lepraria genus and its sexual sister genus Stereocaulon, and added short-read assemblies from an additional 22 individuals across both genera. Comparative genomic analyses revealed that both genera were heterothallic, with intact mating-type loci of both idiomorphs present across each genus. Additionally, we identified and assessed 29 genes involved in meiosis and mitosis and 45 genes that contribute to formation of fungal sexual reproductive structures (ascomata). All genes were present and appeared functional in nearly all Lepraria, and we failed to identify a general pattern of relaxation of selection on these genes across the Lepraria lineage. Together, these results suggest that Lepraria may be capable of sexual reproduction, including mate recognition, meiosis, and production of ascomata.

Conclusions: Despite apparent maintenance of machinery essential for fungal sex, over 200 years of careful observations by lichenologists have produced no evidence of canonical sexual reproduction in Lepraria. We suggest that Lepraria may have instead evolved a form of parasexual reproduction, perhaps by repurposing MAT and meiosis-specific genes. This may, in turn, allow these lichenized fungi to avoid long-term consequences of asexuality, while maintaining the benefit of an unbroken bond with their algal symbionts.

Background: Buffalo is a globally important livestock species, but its reproductive performance is relatively low than cattles. At present, dominant follicle development specific process and mechanistic role of follicular growth related genes in water buffaloes are not well understood. Therefore, we comprehensively performed transcriptomics of granulosa cells and oocytes from different-sized follicles in water buffalo to identify key candidate genes that influence follicle development and diameter, and further explored the potential regulatory mechanisms of granulosa cells and oocytes in the process of water buffalo follicle development.

Results: In this study, we found918 granulosa cell transcripts and 1401 oocyte transcripts were correlated in follicles of different diameters, and the expression differences were significant. Subsequent enrichment analysis of the co-expressed differentially expressed transcripts identified several genes targeted by long non-coding RNAs (lncRNAs) and associated with follicular development. Notably, the upregulation of BUB1 regulated by MSTRG.41325.4 and interactive action of SMAD2 and SMAD7 might have key regulatory role in follicular development. Additionally, we also detected key differentially expressed genes that potentially influence follicular hormone metabolism and growth, like ID2, CHRD, TGIF2 and MAD2L1, and constructed an interaction network between lncRNA transcripts and mRNAs.

Conclusions: In summary, this study preliminarily revealed the differences in gene expression patterns among buffalo follicles of different sizes and their potential molecular regulatory mechanisms. It provides a new perspective for exploring the mechanism of buffalo follicular dominance and improving buffalo reproductive performance.