Search

Abstract Background Hydrogen peroxide (H2O2), a novel water treatment agent, can be used for disinfection, water quality adjustment, and disease prevention, while excessive H2O2 can injure farm animals, even leading to death. Hydrogen peroxide is a recommended disinfectant and bactericide for treating gill diseases and vibriosis in the greenfin horse-faced filefish Thamnaconus septentrionalis. However, its cumulative effect, toxic molecular mechanism and relevant signal transduction/metabolic networks in marine fishes are largely unknown. Results We employed a multi-omics approach to investigate the detrimental effects of 50 mg/L H2O2 exposure (2 h/d) on filefish for 2 d, 4 d, and 6 d. Transcriptome sequencing showed that differentially expressed genes (DEGs) were mainly classified into functions such as signal transduction, nervous system, liver and bile acid metabolism, energy metabolism, cell adhesion and communication, inflammation and immune response. Metabolomic analysis found that the significantly changed metabolites (SCMs) were involved in phenylalanine metabolism, inflammatory mediator regulation, linoleic acid metabolism, and necroptosis. The main SCMs were cholic acid, carnitine C12:1, dimethylmalonic acid, glutamic acid, L-lactic acid, shikimic acid, 2-methylsuccinic acid, and others. Moreover, H2O2-induced oxidative stress also disturbs the balance of the gut microbiota, altering the microbial composition and affecting digestive processes. Conclusions Integrated multiomics analysis revealed that H2O2-induced detrimental impacts include mucosal damage, inflammatory and immune responses, altered energy metabolism, and gut microbiota disorders. These findings offer novel insights into the harmful effects and signal transduction/metabolic pathways triggered by H2O2 exposure in marine fishes.

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

Abstract Background Two-line hybrid wheat technology system is one way to harness wheat heterosis both domestically and internationally. Seed vigor is a crucial parameter for assessing seed quality, as enhanced seed vigor can lead to yield increments of over 20% to a certain extent. MicroRNAs (miRNAs) were known to participate in the development and vigor of seed in plants, but its impact on seed vigor in two-line hybrid wheat remains poorly elucidated. Results The hybrid (BS1453/11GF5135) wheat exhibited superiority in seed vigor and anti-aging capacity, compared to its male parent (11GF5135, MP) and female parent (BS1453, FP). We identified four miRNAs associated with seed vigor, all of which are novel miRNAs. The majority of targets of miRNAs were related to ubiquitin ligases, kinases, sucrose synthases and hydrolases, involving in starch and sucrose metabolism, hydrolysis, catalysis, plant hormone signal transduction, and other pathways, which played crucial roles in seed development. Additionally, we also found miR531 was differentially expressed in both male parent and hybrid, and its target gene was a component of the E1 subunit of α-ketoate dehydrogenase complex, which interacted with dihydrolipoamide acetyltransferase (E2) and dihydrolipoyl dehydrogenase (E3). Finally, We established a presumptive interaction model to speculate the relationship of miR531 and seed vigor. Conclusions This study analyzed the seed vigor of two-line hybrid wheat, and screened seed vigor-related miRNAs. Meanwhile speculated the genetic relationship of hybrid and parents, in terms of miRNAs. Consequently, the present study provides new insights into the miRNA-mediated gene and protein interaction network that regulates seed vigor. These findings hold significance for enhancing the yield and quality of two-line hybrid wheat, facilitating its future applications.

Abstract Background Garuga Roxb. is a genus endemic to southwest China and other tropical regions in Southeast Asia facing risk of extinction due to the loss of tropical forests and changes in land use. Conducting a genome survey of G. forrestii contribute to a deeper understanding and conservation of the genus. Results This study utilized genome survey of G. forrestii generated approximately 54.56 GB of sequence data, with approximately 112 × coverage. K-mer analysis indicated a genome size of approximately 0.48 GB, smaller than 0.52GB estimated by flow cytometry. The heterozygosity is of about 0.54%, and a repeat rate of around 51.54%. All the shotgun data were assembled into 339,729 scaffolds, with an N50 of 17,344 bp. The average content of guanine and cytosine was approximately 35.16%. A total of 330,999 SSRs were detected, with mononucleotide repeats being the most abundant at 70.16%, followed by dinucleotide repeats at 20.40%. We conducted a preliminary ploidy assessment using Smudgeplot and observed a clear bimodal distribution in G. forrestii at 1/2 relative coverage depth and total coverage depth (2n), suggesting a potential diploid genome structure. A pseudo chromosome of G. forrestii and a gemone of Boswellia sacra were used as reference genome to perform a primer population resequencing analysis within three Garuga species. Principal component analysis (PCA) indicated three distinct groups, but genome wide phylogenetics represented conflicting both between the dataset of different reference genomes and between maternal and nuclear genome. Conclusion In summary, the genome of G. forrestii is small, and the phylogenetic relationships within the Garuga genus are complex. The genetic data presented in this study holds significant value for comprehensive whole-genome analyses, the evaluation of population genetic diversity, investigations into adaptive evolution, the advancement of artificial breeding efforts, and the support of species conservation and restoration initiatives. Ultimately, this research contributes to reinforcing the conservation and management of natural ecosystems, promoting biodiversity conservation, and advancing sustainable development.

Abstract Background The 4 serotypes of dengue virus (DENV1-4) can each cause potentially deadly dengue disease, and are spreading globally from tropical and subtropical areas to more temperate ones. Nepal provides a microcosm of this global phenomenon, having met each of these grim benchmarks. To better understand DENV transmission dynamics and spread into new areas, we chose to study dengue in Nepal and, in so doing, to build the onsite infrastructure needed to manage future, larger studies. Methods and results During the 2022 dengue season, we enrolled 384 patients presenting at a hospital in Kathmandu with dengue-like symptoms; 79% of the study participants had active or recent DENV infection (NS1 antigen and IgM). To identify circulating serotypes, we screened serum from 50 of the NS1+ participants by RT-PCR and identified DENV1, 2, and 3 – with DENV1 and 3 codominant. We also performed whole-genome sequencing of DENV, for the first time in Nepal, using our new on-site capacity. Sequencing analysis demonstrated the DENV1 and 3 genomes clustered with sequences reported from India in 2019, and the DENV2 genome clustered with a sequence reported from China in 2018. Conclusion These findings highlight DENV’s geographic expansion from neighboring countries, identify China and India as the likely origin of the 2022 DENV cases in Nepal, and demonstrate the feasibility of building onsite capacity for more rapid genomic surveillance of circulating DENV. These ongoing efforts promise to protect populations in Nepal and beyond by informing the development and deployment of DENV drugs and vaccines in real time.

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

Abstract Background Grapes are highly valued for their nutritional and economic benefits, and have been widely studied for their biological attributes such as fruit development, quality formation, and stress resistance. One significant threat to grape quality is gray mold, caused by Botrytis cinerea, which can infect the flowers, fruits, leaves, and stems. The quantitative real-time PCR (qRT-PCR), known for its high sensitivity and quantitative accuracy, is an essential tool for analyzing gene expression related to the pathogenesis of gray mold, thereby providing deeper insights into the disease. Result In this study, we aim to identify stable internal reference genes crucial for accurate gene expression analysis via qRT-PCR. Utilizing transcriptome data from grapes under various disease stresses, we identified twelve candidate reference genes with consistently high expression levels. The stability of these genes was assessed through delta-CT, geNorm, NormFinder, BestKeeper, and RefFinder analyses after establishing the cycling thresholds (Ct) in different grape varieties treated with Botrytis cinerea. Conclusions Our findings reveal that VIT-17s0000g02750 and VIT-06s0004g04280 exhibit stable expression and are suitable as new reference genes. This foundational work supports further research into the molecular mechanisms of grape biological processes.

Abstract Background Multiheme cytochromes c (MHC) provide prokaryotes with a broad metabolic versatility that contributes to their role in the biogeochemical cycling of the elements and in energy production in bioelectrochemical systems. However, MHC have only been isolated and studied in detail from a limited number of species. Among these, Desulfuromonadia spp. are particularly MHC-rich. To obtain a broad view of the diversity of MHC, we employed bioinformatic tools to study the cytochromome encoded in the genomes of the Desulfuromonadia class. Results We found that the distribution of the MHC families follows a different pattern between the two orders of the Desulfuromonadia class and that there is great diversity in the number of heme-binding motifs in MHC. However, the vast majority of MHC have up to 12 heme-binding motifs. MHC predicted to be extracellular are the least conserved and show high diversity, whereas inner membrane MHC are well conserved and show lower diversity. Although the most prevalent MHC have homologues already characterized, nearly half of the MHC families in the Desulforomonadia class have no known characterized homologues. AlphaFold2 was employed to predict their 3D structures. This provides an atlas of novel MHC, including examples with high beta-sheet content and nanowire MHC with unprecedented high numbers of putative heme cofactors per polypeptide. Conclusions This work illuminates for the first time the universe of experimentally uncharacterized cytochromes that are likely to contribute to the metabolic versatility and to the fitness of Desulfuromonadia in diverse environmental conditions and to drive biotechnological applications of these organisms.

Abstract Background The tea plant Camellia sinensis (L.) O. Kuntze is a perennial crop, invaded by diversity of insect pest species, and pink tea mite is one of the most devastating pests for sustainable tea production. However, molecular mechanism of defense responses against pink tea mites in tea is still unknown. In this study, metabolomics and transcriptome profiles of susceptible and resistant tea varieties were compared before and after pink tea mite infestation. Results Metabolomics analysis revealed that abundance levels of polyphenol-related compounds changed significantly before and after infestation. At the transcript level, nearly 8 GB of clean reads were obtained from each sequenced library, and a comparison of infested plants of resistant and susceptible tea varieties revealed 9402 genes with significant differential expression. An array of genes enriched in plant pathogen interaction and biosynthetic pathways of phenylpropanoids showed significant differential regulation in response to pink tea mite invasion. In particular, the functional network linkage of disease resistant proteins, phenylalanine ammonia lyase, flavanone -3-hydroxylase, hydroxycinnamoyl-CoA shikimate transferase, brassinosteroid-6-oxidase 1, and gibberellin 2 beta-dioxygenase induced dynamic defense signals to suppress prolonged pink tea mite attacks. Further integrated analyses identified a complex network of transcripts and metabolites interlinked with precursors of various flavonoids that are likely modulate resistance against to pink tea mite. Conclusions Our results characterized the profiles of insect induced metabolic and transcript reprogramming and identified a defense regulatory network that can potentially be used to fend off pink tea mites damage.

Abstract Background With the improvement of living standards, the quality of chicken has become a significant concern. Chinese Dagu Chicken (dual-purpose type) and Arbor Acres plus broiler (AA+ broiler) (meat-type) were selected as the research subjects in this study, the meat quality of the breast and leg muscles were measured. However, the molecular mechanism(s) underlying regulation of muscle development are not yet fully elucidated. Therefore, finding molecular markers or major genes that regulate muscle quality has become a crucial breakthrough in chicken breeding. Unraveling the molecular mechanism behind meat traits in chicken and other domestic fowl is facilitated by identifying the key genes associated with these developmental events. Here, a comparative transcriptomic analysis of chicken meat was conducted on breast muscles (BM) and leg muscles (LM) in AA+ broilers (AA) and Dagu chicken (DG) to explore the differences in their meat traits employing RNA-seq. Results Twelve cDNA libraries of BM and LM from AA and DG were constructed from four experimental groups, yielding 14,464 genes. Among them, Dagu chicken breast muscles (DGB) vs AA+ broilers breast muscles (AAB) showed 415 upregulated genes and 449 downregulated genes, Dagu chicken leg muscles (DGL) vs AA+ broilers leg muscles (AAL) exhibited 237 upregulated genes and 278 downregulated genes, DGL vs DGB demonstrated 391 upregulated genes and 594 downregulated genes, and AAL vs AAB displayed 122 upregulated genes and 154 downregulated genes. 13 genes, including nine upregulated genes (COX5A, COX7C, NDUFV1, UQCRFS1, UQCR11, BRT-1, FGF14, TMOD1, MYOZ2) and four downregulated genes (MYBPC3, MYO7B, MTMR7, and TNNC1), were found to be associated with the oxidative phosphorylation signaling pathway. Further analysis revealed that the differentially expressed genes (DEGs) from muscle were enriched in various pathways, such as metabolic pathways, oxidative phosphorylation, carbon metabolism, glycolysis, extracellular matrix-receptor interaction, biosynthesis of amino acids, focal adhesion, vascular smooth muscle contraction, and cardiac muscle contraction, all of which are involved in muscle development and metabolism. This study also measured the meat quality of the breast and leg muscles from the two breeds, which demonstrated superior overall meat quality in Chinese Dagu Chicken compared to the AA+ broiler. Conclusions Our findings show that the meat quality of dual-purpose breeds (Chinese Dagu chicken) is higher than meat-type (AA+ broiler), which may be related to the DEGs regulating muscle development and metabolism. Our findings also provide transcriptomic insights for a comparative analysis of molecular mechanisms underlying muscle development between the two breeds, and have practical implications for the improvement of chicken breeding practices.

Abstract Background Calmodulin-binding transcription activator (CAMTA) is comprised of a group of transcription factors and plays an important role in the Ca2+ signaling pathway, mediating various molecular responses via interactions with other transcription factors and binding to the promoter region of specific genes. Mung beans (Vigna radiata) are one of the most commonly consumed commodities in Asia. To date, CAMTA proteins have not been characterized in this important crop plant. Results Eight paralogous Vr CAMTA genes were identified and found to be distributed on five of the 11 chromosomes. The proteins possessed CG-1 DNA-binding domains with bipartite NLS signals, ankyrin domains, CaM-binding IQ motifs, and CaM-binding domain (CaMBD). The 2 kb upstream regions of VrCAMTA genes contained sequence motifs of abscisic acid-responsive elements (ABRE) and ethylene-responsive elements (ERE), and binding sites for transcription factors of the bZIP and bHLH domains. Analysis of RNA-seq data from a public repository revealed ubiquitous expression of the VrCAMTA genes, as VrCAMTA1 was expressed at the highest level in seedling leaves, whereas VrCAMTA8 was expressed at the lowest level, which agreed with the RT-qPCR analysis performed on the first true leaves. On day four after leaf emergence, all VrCAMTA genes were upregulated, with VrCAMTA1 exhibiting the highest degree of upregulation. In darkness on day 4, upregulation was not observed in most VrCAMTA genes, except VrCAMTA7, for which a low degree of upregulation was found, whereas no difference was found in VrCAMTA8 expression between light and dark conditions. Treatment with calcium ionophores enhanced VrCAMTA expression under light and/or dark conditions at different times after leaf emergence, suggesting that calcium signaling is involved in the light-induced upregulation of VrCAMTA gene expression. Conclusions The expression dependence of nearly all VrCAMTA genes on light and calcium signaling suggests their possible differential but likely complementary roles during the early stages of mung bean growth and development.

Abstract The use of complete organelle genomes, including chloroplast and mitochondrial genomes, is a powerful molecular method for studying biological evolution and gene transfer. However, in the case of Polygonaceae, an important family with numerous edible, medicinal, and ornamental species, the mitochondrial genomes of only three species have been sequenced and analyzed. In this study, we present the mitochondrial and chloroplast genomes of two important Tibetan medicinal plants, Bistorta viviparum and B. macrophyllum. All the organelle genomes are assembled into a single circular structure and contain a common set of 32 protein-coding genes (PCGs). Some genes such as rps2 and ndhF were found to have high nucleotide polymorphism (Pi) in the chloroplast genomes, while cox1, mttB and rps12 showed pronounced Pi values in the mitochondrial genomes. Furthermore, our analysis revealed that most chloroplast genes and mitochondrial PCGs in Polygonaceae plants are under purifying selection. However, a few genes, including the chloroplast gene psaJ and the mitochondrial genes ccmFc, atp8 and nad4, showed positive selection in certain Polygonaceae plants, as indicated by a Ka/Ks ratio greater than one. Structural variation analysis revealed a wealth of differences between the mitochondrial genomes of five Polygonaceae species, with a particularly notable large-scale inversion observed between Reynoutria japonica and Fallopia aubertii. Furthermore, an analysis of the homologous sequences in the chloroplast and mitochondrial genomes revealed that the rps7 has been transferred from the chloroplast to the mitochondrial genome in all five Polygonaceae species. Finally, ecological niche models were constructed for B. viviparum and B. macrophyllum, indicating that mean annual temperature and altitude are the main climatic factors influencing the distribution of both species. Although the current distribution of B. viviparum is significantly wider than that of B. macrophyllum, projections suggest that the optimal growth ranges of both species will expand in the future, with B. macrophyllum potentially exceeding B. viviparum. This study not only contributes to the plastid genome database for Polygonaceae plants, but also provides theoretical insights into the adaptive evolution of these plants.

Abstract Background Mongolian cattle are local breeds in northern China with excellent adaptability to harsh environmental conditions. Adipose tissues play essential roles in tolerance to cold and disease, but the associated cellular and molecular mechanisms are unclear. Methods Single-nucleus RNA sequencing (snRNA-seq) was performed on the adipose tissues from the subcutaneous (SAT), greater omentum (OAT) and perirenal (PAT) of 3 healthy cattle. The adipogenic trajectory was analyzed, and the functional roles of gene of interest were verified in vitro. Results There were different cell subpopulations in adipose tissues. The lipid-deposition adipocytes identified by the PTGER3 marker exhibited outstanding characteristics in SAT. In PAT and OAT, aldosterone was expressed to provide clues for the differential brown adipocytes. Among the DEGs by comparing OAT with SAT and PAT with OAT, C3 was significantly expressed in most of the cell populations in SAT. G0S2, LIPE, LPIN1, PTGER3 and RGCC took part in the adipogenic trajectory from preadipocyte commitment to mature adipocytes. S100A4 expression affected Ca2+ signaling and the expression of UCP1 ~ 3, FABP4 and PTGER3. Conclusion The cell heterogeneity and genes expressed in adipose tissues of Mongolian cattle not only determine the endocrine and energy storage, but contribute to adapt to cold and disease resistance.

Abstract Whole genome sequencing (WGS) has become a gold standard for diagnosing genomic variation. Peripheral blood is a common sample source for the extraction of nucleic acids for Next-Generation Sequencing (NGS) applications. Here, we present an integrated and fully automated device design that uses new concepts of fluid mechanics, heat-mass transfer, and thermodynamics of enzymatic reactions to extract nucleic acids from the blood and perform DNA library preparation from a pre-filled plate. We demonstrate that the presented device effectively extracts dsDNA with an average of 25.03 µg/mL and 25.91 µg/mL yield from citrate-stabilized human peripheral blood stored in Fresh (4 °C) and Frozen (-20 °C) conditions, respectively. Furthermore, our method automatically extracts nucleic acids and creates a high-quality sequence-ready DNA library from blood stabilized with citrate and EDTA for 8 samples simultaneously in a single run with a total operation time of ~ 7 h. Our results show the required coverage and depth of the genome, highlighting an essential application of this device in processing blood samples for genome sequencing. Graphical Abstract

Abstract Klebsiella pneumoniae is a common pathogen capable of causing a wide range of infections. Antibiotic resistance complicates treatment of these infections significantly. We are comparing resistance levels and genotypes among two collections of K. pneumoniae clinical isolates from Alexandria Main University Hospital (AMUH). We used disc diffusion and Minimum Inhibitory Concentration (MIC) by microbroth dilution to assess resistance levels and performed whole genome sequencing (WGS) to describe multilocus sequence types (MLST) and resistance gene presence. Among a collection of 56 K. pneumoniae clinical isolates (19 from 2019 to 37 from 2021), multidrug resistance (MDR) was 33% and 10%, extended drug resistance (XDR) was 24% and 46% and pan-drug resistance (PDR) was 43% and 43%, respectively. We identified 15 MLST STs including two novel types (ST-6118 and ST-6119 ). ST-101 and ST-383 were common between the two collections; ST-101 was the most common genotype in 2019 (28.6%) and ST-147 was most common in 2021 (25%). Ampicillin/sulbactam, amikacin, cefepime, ceftriaxone and ertapenem MICs were significantly higher in 2021. Prevalence of aph(3’) – Ia, aph(3’)-VI, mphA was significantly higher in 2021. The increasing resistance levels and the persistence of some MDR/XDR genotypes is concerning. Understanding mechanisms of resistance will inform infection control and antimicrobial stewardship plans to prevent evolution and spread of XDR and PDR strains.

Abstract Elsinoë species are phytopathogenic fungi that cause serious scab diseases on economically important plants. The disease symptoms arise from the effects of a group of phytotoxins known as elsinochromes, produced via a type-I polyketide synthase (PKS) biosynthetic pathway. The elsinochrome gene cluster was first annotated in Elsinoë fawcettii where the main type-I PKS gene was characterized as EfPKS1. A later study showed that this gene and the associated cluster had not been correctly annotated, and that EfPKS1 was actually the anchor gene of the melanin biosynthetic pathway. A new type-I PKS gene EfETB1 associated with elsinochrome production was also identified. The aim of this study was to identify all type-I PKS genes in the genomes of seven Elsinoë species with the goal of independently verifying the PKS containing clusters for both melanin and elsinochrome production. A total of six type-I PKS classes were identified, although there was variation between the species in the number and type of classes present. Genes similar to the E. fawcettii EfPKS1 and EfETB1 type-I PKS genes were associated with melanin and elsinochrome production respectively in all species. The complete melanin and elsinochrome PKS containing clusters were subsequently annotated in all the species with high levels of synteny across Elsinoë species. This study provides a genus-level overview of type-I PKS distribution in Elsinoë species, including an additional line of support for the annotation of the melanin and elsinochrome PKS containing clusters in these important plant pathogens.

Abstract Background Female ticks remain attached to their host for multiple days to complete a blood meal. This prolonged feeding period is accompanied by a significant increase in the tick’s size and body weight, paralleled by noteworthy changes to the tick midgut. While the midgut is recognized for its established role in blood storage and processing, its importance extends to playing a crucial role in the acquisition, survival, and proliferation of pathogens. Despite this, our overall understanding of tick midgut biology is limited. Results Our transcriptome analysis identified 15,599 putative DNA coding sequences (CDS), which were classified into 26 functional groups. Dimensional and differential expression analyses revealed four primary transcriptional profiles corresponding to unfed, slow-feeding, transitory (from slow- to rapid-feeding), and rapid-feeding stages. Additionally, comparing the current dataset with previously deposited transcriptome from other tick species allowed the identification of commonly expressed transcripts across different feeding stages. Conclusion Our findings provide a detailed temporal resolution of numerous metabolic pathways in the midgut of A. americanum adult females throughout the feeding process, highlighting the dynamic transcriptional regulation of the tick’s midgut as feeding progresses. Furthermore, we identified conserved transcripts across three different tick species that exhibit similar expression patterns. This knowledge not only enhances our understanding of the physiological processes within the tick midgut but also opens up potential avenues for developing control methods that target multiple tick species.

Abstract Background Vibrio cholerae O1 El Tor, the etiological agent responsible for the last cholera pandemic, has become a well-established model organism for which some genetic tools are available. While CRISPRi technology has been applied to V. cholerae, improvements were necessary to upscale it and enable pooled screening by high-throughput sequencing in this bacterium. Results In this study, we present a genome-wide CRISPR-dCas9 screen specifically optimized for the N16961 El Tor model strain of V. cholerae. This approach is characterized by a tight control of dCas9 expression and activity, as well as a streamlined experimental setup. Our library allows the depletion of 3,674 (98.9%) annotated genes from the V. cholerae genome. To confirm its effectiveness, we screened for genes that are essential during exponential growth in rich medium and identified 369 genes for which guides were significantly depleted from the library (log2FC < -2). Remarkably, 82% of these genes had previously been described as hypothetical essential genes in V. cholerae or in a closely related bacterium, V. natriegens. Conclusion We thus validated the robustness and accuracy of our CRISPRi-based approach for assessing gene fitness in a given condition. Our findings highlight the efficacy of the developed CRISPRi platform as a powerful tool for high-throughput functional genomics studies of V. cholerae.

Abstract Klebsiella variicola is considered an emerging pathogen, which may colonize a variety of hosts, including environmental sources. Klebsiella variicola investigated in this study was obtained from an influent wastewater treatment plant in the North-West Province, South Africa. Whole genome sequencing was conducted to unravel the genetic diversity and antibiotic resistance patterns of K. variicola. Whole genome core SNP phylogeny was employed on publicly available 170 genomes. Furthermore, capsule types and antibiotic resistance genes, particularly beta-lactamase and carbapenems genes were investigated from the compared genomes. A 38 099 bp bacteriophage was uncovered alongside with K. variicola genome. Whole genome sequencing revealed that the extended beta-lactamase bla LEN (75.3%) of the beta-lactamase is dominant among compared K. variicola strains. The identified IncF plasmid AA035 confers resistance genes of metal and heat element subtypes, i.e., silver, copper, and tellurium. The capsule type KL107-D1 is a predominant capsule type present in 88.2% of the compared K. variicola genomes. The phage was determined to be integrase-deficient consisting of a fosB gene associated with fosfomycin resistance and clusters with the Wbeta genus Bacillus phage group. In silico analysis showed that the phage genome interacts with B. cereus as opposed to K. variicola strain T2. The phage has anti-repressor proteins involved in the lysis-lysogeny decision. This phage will enhance our understanding of its impact on bacterial dissemination and how it may affect disease development and antibiotic resistance mechanisms in wastewater treatment plants. This study highlights the need for ongoing genomic epidemiological surveillance of environmental K. variicola isolates.

Abstract The nutritional quality and yield of oats (Avena sativa) are often compromised by plant diseases such as red leaf, powdery mildew, and leaf spot. Sugars Will Eventually be Exported Transporters (SWEETs) are newly identified sugar transporters involved in regulating plant growth and stress responses. However, the roles of SWEET genes in biotic stress responses remain uncharacterized in oats. In this study, 13 AsSWEET genes were identified across nine chromosomes of the oat genome, all of which were predicted to contain seven transmembrane regions. Phylogenetic analysis revealed four clades of AsSWEET proteins, with high homology to SWEET proteins in the Poaceae family. Collinearity analysis demonstrated strong relationships between oat and Zea mays SWEETs. Using subcellular localization prediction tools, AsSWEET proteins were predicted to localize to the plasma membrane. Promoter analysis revealed cis-acting elements associated with light response, growth, and stress regulation. Six AsSWEET proteins were predicted to interact in a network centered on AsSWEET1a and AsSWEET11. Gene expression analysis of two oat varieties, ‘ForagePlus’ and ‘Molasses’, indicated significant expression differences in several AsSWEET genes following infection with powdery mildew or leaf spot, including AsSWEET1a, AsSWEET1b, AsSWEET2b, AsSWEET3a, AsSWEET11, and AsSWEET16. These SWEET genes are potential candidates for disease resistance in oats. This study provides a foundation for understanding the regulatory mechanisms of AsSWEET genes, particularly in response to powdery mildew and leaf spot, and offers insights for enhancing oat molecular breeding.

Abstract Background Atrial fibrillation (AF) stands as a prevalent and detrimental arrhythmic disorder, characterized by intricate pathophysiological mechanisms. The availability of reliable and reproducible AF models is pivotal in unraveling the underlying mechanisms of this complex condition. Unfortunately, the researchers are still confronted with the absence of consistent in vitro AF models, hindering progress in this crucial area of research. Methods Human induced pluripotent stem cells derived atrial myocytes (hiPSC-AMs) were generated based on the GiWi methods and were verified by whole-cell patch clamp, immunofluorescent staining, and flow cytometry. Then hiPSC-AMs were employed to establish the AF model by HS. Whole-cell patch clamp technique and calcium imaging were used to identify the AF model. The stability of 29 reference genes was evaluated using delta-Ct, GeNorm, NormFinder, and BestKeeper algorithms; Results HiPSC-AMs displayed atrial myocyte action potentials and expressed the atrial-specific protein MLC-2 A and NR2F2, about 70% of the cardiomyocytes were MLC-2 A positive. After HS, hiPSC-AMs showed a significant increase in beating frequency, a shortened action potential duration, and increased calcium transient frequency. Of the 29 candidate genes, the top five most stably ranked genes were ABL1, RPL37A, POP4, RPL30, and EIF2B1. After normalization using ABL1, KCNJ2 was significantly upregulated in the AF model; Conclusions: In the hiPSC-AMs AF model established by HS, ABL1 provides greater normalization efficiency than commonly used GAPDH.

Abstract Background MYB transcription factors regulate anthocyanin biosynthesis across numerous plant species. However, comprehensive genome-wide investigations regarding the R2R3-MYB gene family and its involvement in regulating anthocyanin biosynthesis in the red and white fruit color morphs of Fragaria pentaphylla remain scarce. Results A total of 101 FpR2R3-MYB genes were identified from the F. pentaphylla genome and were divided into 34 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were particularly conserved among the FpR2R3-MYB genes, especially members within the same subgroup. The FpR2R3-MYB genes were distributed over seven F. pentaphylla chromosomes. Analysis of gene duplication events revealed five pairs of tandem duplication genes and 16 pairs of segmental duplication genes, suggesting that segmental duplications are the major pattern for expansion of the FpR2R3-MYB gene family expansion in F. pentaphylla. Cis-regulatory elements of the FpR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponse. Based on the analysis of the FpR2R3-MYB gene family and transcriptome sequencing (RNA-seq) data, FpMYB9 was identified as a key transcription factor involved in the regulation of anthocyanin synthesis in F. pentaphylla fruits. The expression of FpMYB9 increases significantly during the ripening stage of red fruits, as confirmed by reverse transcription quantitative real-time PCR. In addition, subcellular localization experiments further confirmed the nuclear presence of FpMYB9, supporting its role as a transcription factor involved in anthocyanin biosynthesis. Conclusion Our results showed that the FpR2R3-MYB genes are highly conserved and play important roles in the anthocyanin biosynthesis in F. pentaphylla fruits. Our results also provide a compelling basis for further understanding of the regulatory mechanism underlying the role of FpMYB9 in anthocyanin formation in F. pentaphylla fruits.

Abstract Semen Ziziphi Spinosae (SZS) is a traditional Chinese herbal medicine widely used to treat insomnia and anxiety in clinical practice. Currently, the demand for SZS is increasing every year, but the production of wild SZS is unstable due to environmental factors. Grafting sour jujube scions onto sour jujube or jujube tree stocks can achieve a high production rate within a short period of time. However, the effects of grafting on the quality of SZS have not been reported. This study investigated the differences between wild-type and grafted SZS from three aspects: phenotype, chemical composition, and molecular mechanism. The findings revealed that the grafted specimens were generally larger in morphology and lighter in color than the wild-type samples. The dimensions of both the grafted specimens were generally larger than those of the wild specimens. The HPLC-ELSD results revealed that the three main chemical components in the grafted SZS, namely, spinosin, jujuboside A, and jujuboside B, had higher contents than their wild-type counterparts. Comprehensive transcriptome sequencing analysis and KEGG annotation revealed that DEG enrichment between grafted and wild-type SZS occurred mainly during stress resistance and rootstock scion healing. There were 23 DEGs that may encode enzymes involved in the biosynthetic pathway of flavonoids and 21 genes encoding terpenoid saponins. Further investigation revealed that the expression of the genes C4H, CHS, CHI, and F3’5’H in the flavonoid biosynthesis pat.hway and HMGR, MVK, MVD, and FPPS in the saponin biosynthesis pathway accounted for the difference in quality between grafted and wild SZS. Furthermore, WGCNA identified 15 core genes related to medicinal ingredients between grafted and wild SZS. These results provide support for further research on the differences in the quality of medicinal ingredients between grafted and wild SZS.

Abstract Background Cold is a significant limiting factor in productivity, particularly in northwestern and eastern China. Calcium-Dependent Protein Kinases (CDPKs), a primary calcium signaling sensor in plants, play an important role in their response to cold. Snow lotus (Sasussured involucrata Kar L) is a plant that thrives in harsh climates and grows in northwest China. However, there were no reports on the transcriptome of OE-SikCDPK1 transgenic tobacco in response to cold. Results When exposed to cold stress, OE-SikCDPK1 plants displayed a cold-tolerant phenotype compared to non-transgenic tobacco. Under cold conditions, relative water content reduced, relative conductivity increased, malondialdehyde levels rose, and cold-responsive gene expression increased. The OE-SikCDPK1 gene and non-transgenic tobacco were employed for research purposes. The transcriptome of leaves was sequenced using the HISAT2 sequencing platform, and the data were used to examine gene function annotation and differentially expressed genes (DEGs). 53,022 DEGs in tobacco leaves under cold treatment were obtained. The GO enrichment results revealed that it was enriched for biological-process, defense response and other processes under cold stress. The KEGG pathway enrichment analysis revealed that the metabolic pathways of significant enrichment of DEGs under cold stress mainly involved MAPK signaling pathway transduction. The transcription factor identification results showed that the transcription factors with the largest number of differential expressions under cold treatment were mainly from WRKY, AP2, MYB, bHLH, NAC and other transcription factor families. Conclusion The cold tolerance mechanism of snow lotus SikCDPK1 was comprehensively analyzed at the transcriptional level for the first time using RNA-seq technology. This study demonstrates that SikCDPK1 can respond to cold by participating in the MAPK signaling pathway and regulating the expression levels of transcription factors, including WRKY, AP2, MYB, bHLH, and NAC. These results offer valuable insights for further exploration of the cold tolerance mechanism associated with SikCDPK1.

Abstract Background Mesenchymal stem cells (MSCs) are multipotent stem cells that are under investigation for use in clinical trials because they are capable of self-renewal and differentiating into different cell types under defined conditions. Nonetheless, the therapeutic effects of MSCs have been constrained by low engraftment rates, cell fusion, and cell survival. Various strategies have been explored to improve the therapeutic efficacy of MSCs, with platelet-derived growth factor (PDGF)-BB emerging as a promising candidate. To enhance our comprehension of the impact of PDGF-BB on the gene expression profile and chromosomal accessibility of MSCs, RNA-sequencing and analysis of chromatin accessibility profiles were conducted on three human primary MSCs in culture, both with and without stimulation by PDGF-BB. Results Integrative analysis of gene expression and chromatin accessibility demonstrated that PDGF-BB treatment modified the chromatin accessibility landscape, marking regions for activation or repression through the AP-1 family transcription factors TEAD, CEBP, and RUNX2. These changes in AP-1 transcription factor expression, in turn, led to cell proliferation and differentiation potential towards osteoblasts, adipocytes, or chondrocytes. The degree of MSC differentiation varies among cells isolated from different donors. The presence of an enrichment of exosome-related genes is also noted among all the differentially expressed genes. Conclusions In conclusion, the observed changes in AP-1 transcription factor expression not only induced cellular proliferation and differentiation, but also revealed variations in the degree of MSC differentiation based on donor-specific differences. Moreover, the enrichment of exosome-related genes among differentially expressed genes suggests a potential significant role for PDGF-BB in facilitating intercellular communication.

Abstract Drought stress poses a critical threat to global crop yields and sustainable agriculture. The GASA genes are recognized for their pivotal role in stress tolerance and plant growth, but little is known about how they function in sunflowers. The investigation aimed to identify and elucidate the role of HaGASA genes in conferring sunflowers with drought tolerance. Twenty-seven different HaGASA gene family members were found in this study that were inconsistently located across eleven sunflower chromosomes. Phylogeny analysis revealed that the sunflower HaGASA genes were divided into five subgroups by comparing GASA genes with those from Arabidopsis, peanut, and soybean, with members within each subgroup displaying similar conserved motifs and gene structures. In-silico evaluation of cis-regulatory elements indicated the existence of specific elements associated with stress-responsiveness being the most abundant, followed by hormone, light, and growth-responsive elements. Transcriptomic data from the NCBI database was utilized to assess the HaGASA genes expression profile in different sunflower varieties under drought conditions. The HaGASA genes expression across ten sunflower genotypes under drought stress, revealed 14 differentially expressed HaGASA genes, implying their active role in the plant’s stress response. The expression in different organs revealed that HaGASA2, HaGASA11, HaGASA17, HaGASA19, HaGASA21 and HaGASA26 displayed maximum expression in the stem. Our findings implicate HaGASA genes in mediating sunflower growth maintenance and adaptation to abiotic stress, particularly drought. The findings, taken together, provided a basic understanding of the structure and potential functions of HaGASA genes, setting the framework for further functional investigations into their roles in drought stress mitigation and crop improvement strategies.

Abstract Background Cashmere, named as “soft gold”, derives from the secondary hair follicles (SHFs) of cashmere goat which is vital to Northwest China’s economy. The cytodifferentiation stage (E120), mirroring the complete hair follicle (HF) structure of adult goats and marking a critical phase in SHF development. Therefore, this study aims to enhance the understanding of SHF development and its impact on fiber quality, informing breeding strategies. Results From the scRNA-seq data analysis, the intricate processes and transcriptional dynamics of inner layer cell differentiation of HFs were unveiled in this study. we identified nine cell populations during cytodifferentiation and key structures such as the hair shaft and inner root sheath. And we discovered three main inner layer lineages and seven subpopulations, clarifying their roles in specialization and signaling. Pseudotime mapping analysis showed cell evolution from early stage to mature stages marked by unique gene expressions, and the intermediate stage on the differentiation of each lineage was revealed. The identification and spatial localization of specific transcription factors, such as GATA3, LEF1 and PRDM1, as well as keratin genes highlight regulatory pathways involved in HF development, which was further validated by immunofluorescence. These findings suggested the potential strategies to improve fiber quality, and the discovery of diverse cell types and their developmental molecular mechanisms, particularly in this species-specific context, offered a nuanced view of the regulatory mechanisms driving HF development in cashmere goats. Conclusion Overall, these findings provide a systematic molecular atlas of skin, defining three major branches and cell states of inner layer cells of HF, and determining how the branch-specific transcription factors, keratins, and signals coordinate HF morphogenesis during cytodifferentiation stage. This research not only advances skin tissue research in goats but also holds broader implications for the understanding of HF regeneration and development across various species.

Abstract Background Cryptolaemus montrouzieri (Coccinellidae) is widely utilized as biological control agents in modern agriculture. A comprehensive understanding of its food preference can help guide mass rearing and safety management during field application of pest control. Although some studies have paid attentions to the impacts of prey shift on C. montrouzieri, little is known regarding the role of post-transcriptional regulations in its acclimation to unnatural preys. Results We performed a genome-wide investigation on alternative splicing dynamics in C. montrouzieri in response to the predation transition from natural prey to unnatural ones. When feeding on undesired diets, 402–764 genes were differentially alternative spliced in C. montrouzieri. It is noteworthy that the majority of these genes (> 87%) were not differentially expressed, and these differentially spliced genes regulated distinct biological processes from differentially expressed genes, such as organ development and morphogenesis, locomotory behavior, and homeostasis processes. These suggested the functionally nonredendant role of alternative splicing in modulating physiological and metabolic responses of C. montrouzieri to the shift to undesired preys. In addition, the individuals feeding on aphids were subject to a lower level of changes in splicing than other alternative diets, which might be because of the similar chemical and microbial compositions. Our study further suggested a putative coupling of alternative splicing and nonsense-mediated decay (AS-NMD), which may play an important role in fine-tuning the protein repertoire of C. montrouzieri, and promoting its acclimation to predation changes. Conclusion These findings highlight the key role of alternative splicing in modulating the acclimation of ladybirds to prey shift and provide new genetic clues for the future application of ladybirds in biocontrol.

Abstract Background Multiple sequence alignment (MSA) has proven extremely useful in computational biology, especially in inferring evolutionary relationships via phylogenetic analysis and providing insight into protein structure and function. An alternative to the standard MSA model is partial order alignment (POA), in which aligned sequences are represented as paths in a graph rather than rows in a matrix. While the POA model has proven useful in several applications (e.g. sequencing reads assembly and pangenome structure exploration), we lack efficient visualization tools that could highlight its advantages. Results We propose Sequence Flow – a web application designed to address the above problem. Sequence Flow presents the POA as a Sankey diagram, a kind of graph visualisation typically used for graphs representing flowcharts. Sequence Flow enables interactive alignment exploration, including fragment selection, highlighting a selected group of sequences, modification of the position of graph nodes, structure simplification etc. After adjustment, the visualization can be saved as a high-quality graphic file. Thanks to the use of SanKEY.js – a JavaScript library for creating Sankey diagrams, designed specifically to visualize POAs, Sequence Flow provides satisfactory performance even with large alignments. Conclusions We provide Sankey diagram-based POA visualization tools for both end users (Sequence Flow) and bioinformatic software developers (SanKEY.js). Sequence Flow webservice is available at https://sequenceflow.mimuw.edu.pl/ . The source code for SanKEY.js is available at https://github.com/Krzysiekzd/SanKEY.js and for Sequence Flow at https://github.com/Krzysiekzd/SequenceFlow .

Abstract Background Certain structural variants (SVs) including large-scale genetic copy number variants, as well as copy number-neutral inversions and translocations may not all be resolved by chromosome karyotype studies. The identification of genetic risk factors for Parkinson’s disease (PD) has been primarily focused on the gene-disruptive single nucleotide variants. In contrast, larger SVs, which may significantly influence human phenotypes, have been largely underexplored. Optical genomic mapping (OGM) represents a novel approach that offers greater sensitivity and resolution for detecting SVs. In this study, we used induced pluripotent stem cell (iPSC) lines of patients with PD-linked SNCA and PRKN variants as a proof of concept to (i) show the detection of pathogenic SVs in PD with OGM and (ii) provide a comprehensive screening of genetic abnormalities in iPSCs. Results OGM detected SNCA gene triplication and duplication in patient-derived iPSC lines, which were not identified by long-read sequencing. Additionally, various exon deletions were confirmed by OGM in the PRKN gene of iPSCs, of which exon 3–5 and exon 2 deletions were unable to phase with conventional multiplex-ligation-dependent probe amplification. In terms of chromosomal abnormalities in iPSCs, no gene fusions, no aneuploidy but two balanced inter-chromosomal translocations were detected in one line that were absent in the parental fibroblasts and not identified by routine single nucleotide variant karyotyping. Conclusions In summary, OGM can detect pathogenic SVs in PD-linked genes as well as reveal genomic abnormalities for iPSCs that were not identified by other techniques, which is supportive for OGM’s future use in gene discovery and iPSC line screening.

Abstract Background In the domain of forensic science, the application of kinship identification and mixture deconvolution techniques are of critical importance, providing robust scientific evidence for the resolution of complex cases. Microhaplotypes, as the emerging class of genetic markers, have been widely studied in forensics due to their high polymorphisms and excellent stability. Results and discussion In this research, a novel and high-efficient panel integrating 33 microhaplotype loci along with a sex-determining locus was developed by the next generation sequencing technology. In addition, we also assessed its forensic utility and delved into its capacity for kinship analysis and mixture deconvolution. The average effective number of alleles (Ae) of the 33 microhaplotype loci in the Guizhou Han population was 6.06, and the Ae values of 30 loci were greater than 5. The cumulative power of discrimination and cumulative power of exclusion values of the novel panel in the Guizhou Han population were 1-5.6 × 10− 43 and 1-1.6 × 10− 15, respectively. In the simulated kinship analysis, the panel could effectively distinguish between parent-child, full-sibling, half-sibling, grandfather-grandson, aunt-nephew and unrelated individuals, but uncertainty rates clearly increased when distinguishing between first cousins and unrelated individuals. For the mixtures, the novel panel had demonstrated excellent performance in estimating the number of contributors of mixtures with 1 to 5 contributors in combination with the machine learning methods. Conclusions In summary, we have developed a small and high-efficient panel for forensic genetics, which could provide novel insights into forensic complex kinships testing and mixture deconvolution.

Abstract Background The excavation and utilization of endophytic fungi from medicinal plants is highly important for the development of new drugs. The endophytic fungus Aspergillus tamarii TPD11, which was isolated and obtained by the authors in the previous stage, can produce a variety of polyphyllins with important potential applications in hemostasis, inflammation and antitumor activities; however, the genomic information of TPD11 is still unknown. Results In this study, we sequenced and assembled the whole genome of the endophytic fungus A. tamarii TPD11, resolved the genome evolutionary relationships of 24 Aspergillus strains, and phylogenetic analysis of the genomes of 16 strains revealed the evolutionary differences between Aspergillus and Penicillium and the mechanisms of genome expansion and contraction. CAZy annotation analysis revealed that TPD11 obtains nutrients mainly by ingesting starch from the host plant. TPD11 has a biosynthesis-related gene cluster for the synthesis of squalestatin S1, and the silencing of this biosynthesis-related gene cluster might increase the content of polyphyllin. Annotation of 11 UDP-glycosyltransferase genes helps to further reveal the biosynthetic pathway of polyphyllin. In addition, secondary metabolism gene cluster and CAZy analyses confirmed the potential probiotic, insecticidal and antimicrobial activities of TPD11 on host plants. Conclusions This study reveals the intrinsic mechanism by which endophytic fungi increase the content of polyphyllin, which provides a basis for the synthetic synthesis of the natural product polyphyllin.

Abstract Background The R2R3-MYB transcription factors in plants participate in various physiological and biochemical processes and responds to various external stimuli. Prunus sibirica (known as Siberian apricot) is a drupe tree species that produces extremely high nutritional value kernels. However, it is susceptiblility to frost damage during the flowering period, results in a marked reduction in kernel yield. Results In this study, the MYB gene family of P. sibirica (PsMYB) was systematically analyzed, and 116 R2R3-MYB genes that were distributed unevenly over eight chromosomes were ultimately screened. Phylogenetic analysis divided these 116 genes into 30 subgroups. We discovered that 37 PsMYBs had cold stress-responsive promoters, and six PsMYBs were annotated to be associated with cold response. Intraspecific homology analysis identified segmental duplication as the primary gene amplification mechanism, and homology analysis of the PsMYB genes with those of five other species revealed phylogenetic relationships with Rosaceae species. Protein interaction studies revealed collaborative regulation of the PsMYB proteins with Arabidopsis protein, and transcriptome analysis identified PsMYB genes that were highly expressed at low temperatures. Additionally, the expression levels of 22 PsMYBs in different tissue parts of P. sibirica and under different low-temperature stress conditions were evaluated using quantitative real-time PCR, with the results verifying that PsMYBs are specifically expressed in different plant parts and may be involved in the growth and development of P. sibirica species. Genes upregulated after exposure to low-temperature stress and likely involved in cold response were identified. Conclusion This study lays a foundation for understanding the molecular biology of PsMYBs in P. sibirica and provides a theoretical basis for the future study of transgenic lines with cold resistance during the flowering period of this tree.

Abstract Background Austroasiatic (AA)-speaking populations in northern Thailand are of significant interest due to their status as indigenous descendants and their location at the crossroads of AA prehistoric distribution across Southern China, the Indian Subcontinent, and Mainland Southeast Asia. However, the complexity of ethnic identification can result in inaccuracies regarding the origin and migration history of these populations. To address this, we have conducted a genome-wide SNP analysis of 89 individuals from two Lavue and three Lwa-endonym populations. We then combined our outcomes with previously published data to elucidate the genetic diversity and clustering of AA groups in northern Thailand. Results Our findings align with existing linguistic classifications, revealing different genetic compositions among the three branches of the Mon-Khmer subfamily within the AA family: Monic, Khmuic, and Palaungic. Although the term “Lua” ethnicity is confusingly used to identify ethnic groups belonging to both Khmuic and Palaungic branches, our genomic data indicate that the Khmuic-speaking Lua living on the eastern side of the region are relatively distant from the Palaungic-speaking Lavue and Lwa populations living on the western side. The Lavue populations, primarily inhabiting mountainous areas, exhibit a genetic makeup unique to the AA family, with a close genetic relationship to the Karenic subgroup of the Sino-Tibetan language family. Conversely, the Lwa and Blang populations, residing in lowland river valleys, display genetic signatures resulting from admixture with Tai-Kadai-speaking ethnic groups. Conclusion Utilizing genome-wide SNP markers, our findings indicate genetic heterogeneity among the Lua, Lavue, and Lwa ethnic groups. The intricate interplay of genetics, cultural heritage, and historical influences has shaped these ethnic communities. Our study underscores the importance of accurate ethnic classifications, emphasizing the use of self-identified endonyms, names created and used by the ethnic groups themselves. This approach respects the AA communities in northern Thailand and acknowledges their significant contributions to advancing our understanding of genetic anthropology.

Abstract Background African cattle represent a unique resource of genetic diversity in response to adaptation to numerous environmental challenges. Characterising the genetic landscape of indigenous African cattle and identifying genomic regions and genes of functional importance can contribute to targeted breeding and tackle the loss of genetic diversity. However, pinpointing the adaptive variant and determining underlying functional mechanisms of adaptation remains challenging. Results In this study, we use selection signatures from whole-genome sequence data of eight indigenous African cattle breeds in combination with gene expression and quantitative trait loci (QTL) databases to characterise genomic targets of artificial selection and environmental adaptation and to identify the underlying functional candidate genes. In general, the trait-association analyses of selection signatures suggest the innate and adaptive immune system and production traits as important selection targets. For example, a large genomic region, with selection signatures identified for all breeds except N’Dama, was located on BTA27, including multiple defensin DEFB coding-genes. Out of 22 analysed tissues, genes under putative selection were significantly enriched for those overexpressed in adipose tissue, blood, lung, testis and uterus. Our results further suggest that cis-eQTL are themselves selection targets; for most tissues, we found a positive correlation between allele frequency differences and cis-eQTL effect size, suggesting that positive selection acts directly on regulatory variants. Conclusions By combining selection signatures with information on gene expression and QTL, we were able to reveal compelling candidate selection targets that did not stand out from selection signature results alone (e.g. GIMAP8 for tick resistance and NDUFS3 for heat adaptation). Insights from this study will help to inform breeding and maintain diversity of locally adapted, and hence important, breeds.

Abstract Background Archaea and Bacteria are distinct domains of life that are adapted to a variety of ecological niches. Several genome-based methods have been developed for their accurate classification, yet many aspects of the specific genomic features that determine these differences are not fully understood. In this study, we used publicly available whole-genome sequences from bacteria ( $$N=2546$$ N = 2546 ) and archaea ( $$N=109$$ N = 109 ). From these, a set of genomic features (nucleotide frequencies and proportions, coding sequences (CDS), non-coding, ribosomal and transfer RNA genes (ncRNA, rRNA, tRNA), Chargaff’s, topological entropy and Shannon’s entropy scores) was extracted and used as input data to develop machine learning models for the classification of archaea and bacteria. Results The classification accuracy ranged from 0.993 (Random Forest) to 0.998 (Neural Networks). Over the four models, only 11 examples were misclassified, especially those belonging to the minority class (Archaea). From variable importance, tRNA topological and Shannon’s entropy, nucleotide frequencies in tRNA, rRNA and ncRNA, CDS, tRNA and rRNA Chargaff’s scores have emerged as the top discriminating factors. In particular, tRNA entropy (both topological and Shannon’s) was the most important genomic feature for classification, pointing at the complex interactions between the genetic code, tRNAs and the translational machinery. Conclusions tRNA, rRNA and ncRNA genes emerged as the key genomic elements that underpin the classification of archaea and bacteria. In particular, higher nucleotide diversity was found in tRNA from bacteria compared to archaea. The analysis of the few classification errors reflects the complex phylogenetic relationships between bacteria, archaea and eukaryotes.

Abstract Background Microsatellites are increasingly realized to have biological significance in human genome and health in past decades, the assembled complete reference sequence of human genome T2T-CHM13 brought great help for a comprehensive study of short tandem repeats in the human genome. Results Microsatellites density landscapes of all 24 chromosomes were built here for the first complete reference sequence of human genome T2T-CHM13. These landscapes showed that short tandem repeats (STRs) are prone to aggregate characteristically to form a large number of STRs density peaks. We classified 8,823 High Microsatellites Density Peaks (HMDPs), 35,257 Middle Microsatellites Density Peaks (MMDPs) and 199, 649 Low Microsatellites Density Peaks (LMDPs) on the 24 chromosomes; and also classified the motif types of every microsatellites density peak. These STRs density aggregation peaks are mainly composing of a single motif, and AT is the most dominant motif, followed by AATGG and CCATT motifs. And 514 genomic regions were characterized by microsatellite density feature in the full T2T-CHM13 genome. Conclusions These landscape maps exhibited that microsatellites aggregate in many genomic positions to form a large number of microsatellite density peaks with composing of mainly single motif type in the complete reference genome, indicating that the local microsatellites density varies enormously along the every chromosome of T2T-CHM13.

Abstract The proteasome plays a crucial role in cellular homeostasis by degrading misfolded, damaged, or unnecessary proteins. Understanding the regulatory mechanisms of proteasome activity is vital, particularly the interaction with activators containing the hydrophobic-tyrosine-any amino acid (HbYX) motif. Here, we present ProEnd, a comprehensive database designed to identify and catalog HbYX motif-containing proteins across the tree of life. Using a simple bioinformatics pipeline, we analyzed approximately 73 million proteins from 22,000 reference proteomes in the UniProt/SwissProt database. Our findings reveal the widespread presence of HbYX motifs in diverse organisms, highlighting their evolutionary conservation and functional significance. Notably, we observed an interesting prevalence of these motifs in viral proteomes, suggesting strategic interactions with the host proteasome. As validation two novel HbYX proteins found in this database were experimentally tested by pulldowns, confirming that they directly interact with the proteasome, with one of them directly activating it. ProEnd’s extensive dataset and user-friendly interface enable researchers to explore the potential proteasomal regulator landscape, generating new hypotheses to advance proteasome biology. This resource is set to facilitate the discovery of novel therapeutic targets, enhancing our approach to treating diseases such as neurodegenerative disorders and cancer.

Abstract Background The Caridea, affiliated with Malacostraca, Decapoda, and Pleocyemata, constitute one of the most significant shrimp groups. They are widely distributed across diverse aquatic habitats worldwide, enriching their evolutionary history. In recent years, considerable attention has been focused on the classification and systematic evolution of Caridea, yet controversies still exist regarding the phylogenetic relationships among families. Methods Here, the complete mitochondrial genome (mitogenome) sequences of five caridean species, namely Heterocarpus sibogae, Procletes levicarina, Macrobrachium sp., Latreutes anoplonyx, and Atya gabonensis, were determined using second-generation high-throughput sequencing technology. The basic structural characteristics, nucleotide composition, amino acid content, and codon usage bias of their mitogenomes were analyzed. Selection pressure values of protein-coding genes (PCGs) in species within the families Pandalidae, Palaemonidae, and Atyidae were also computed. Phylogenetic trees based on the nucleotide and amino acid sequences of 13 PCGs from 103 caridean species were constructed, and divergence times for various families within Caridea were estimated. Results The mitogenome of these five caridean species vary in length from 15,782 to 16,420 base pairs, encoding a total of 37 or 38 genes, including 13 PCGs, 2 rRNA genes, and 22 or 23 tRNA genes. Specifically, L. anoplonyx encodes an additional tRNA gene, bringing its total gene count to 38. The base composition of the mitogenomes of these five species exhibited a higher proportion of adenine-thymine (AT) bases. Six start codons and four stop codons were identified across the five species. Analysis of amino acid content and codon usage revealed variations among the five species. Analysis of selective pressure in Pandalidae, Palaemonidae, and Atyidae showed that the Ka/Ks values of PCGs in all three families were less than 1, indicating that purifying selection is influencing on their evolution. Phylogenetic analysis revealed that each family within Caridea is monophyletic. The results of gene rearrangement and phylogenetic analysis demonstrated correlations between these two aspects. Divergence time estimation, supported by fossil records, indicated that the divergence of Caridea species occurred in the Triassic period of the Mesozoic era, with subsequent differentiation into two major lineages during the Jurassic period. Conclusions This study explored the fundamental characteristics and phylogenetic relationships of mitogenomes within the infraorder Caridea, providing valuable insights into their classification, interspecific evolutionary patterns, and the evolutionary status of various Caridea families. The findings provide essential references for identifying shrimp species and detecting significant gene rearrangements within the Caridea infraorder.

Abstract Salinity adaptation is an important issue in aquaculture. Understanding the immediate-early response to salinity stress helps in comprehending this process. In vitro experiments using cell lines can explain cell-independent reactions without the involvement of hormones in vivo. In this study, salinity stress experiments were conducted using cell line derived from the gills of Gymnocypris przewalskii (GPG cell line) to isolate immediate-early response-related genes and miRNAs using transcriptomics, followed by bioinformatics analysis. The results showed that intracellular free Ca2+ appeared to be a key factor in cell sensing and initiating downstream cell signaling in response to external salinity. Additionally, cell apoptosis was the most common feature of salinity stress, with multiple signaling pathways involved in salinity-induced cell apoptosis. Furthermore, MiRNAs played a crucial role in the rapid response to salinity stress by selectively inhibiting the expression of specific genes. Additionally, for the first time in the G. przewalskii genome, Tf2 and TY3 families of transposons were found to have responsive roles to the external salinity stress. This study contributes to a better understanding of osmotic sensing in G. przewalskii and provides theoretical assistance for improving salinity adaptation in aquaculture fish species.

Abstract Background Assembly of metagenomic samples can provide essential information about the mobility potential and taxonomic origin of antibiotic resistance genes (ARGs) and inform interventions to prevent further spread of resistant bacteria. However, similar to other conserved regions, such as ribosomal RNA genes and mobile genetic elements, almost identical ARGs typically occur in multiple genomic contexts across different species, representing a considerable challenge for the assembly process. Usually, this results in many fragmented contigs of unclear origin, complicating the risk assessment of ARG detections. To systematically investigate the impact of this issue on detection, quantification and contextualization of ARGs, we evaluated the performance of different assembly approaches, including genomic-, metagenomic- and transcriptomic-specialized assemblers. We quantified recovery and accuracy rates of each tool for ARGs both from in silico spiked metagenomic samples as well as real samples sequenced using both long- and short-read sequencing technologies. Results The results revealed that none of the investigated tools can accurately capture genomic contexts present in samples of high complexity. The transcriptomic assembler Trinity showed a better performance in terms of reconstructing longer and fewer contigs matching unique genomic contexts, which can be beneficial for deciphering the taxonomic origin of ARGs. The currently commonly used metagenomic assembly tools metaSPAdes and MEGAHIT were able to identify the ARG repertoire but failed to fully recover the diversity of genomic contexts present in a sample. On top of that, in a complex scenario MEGAHIT produced very short contigs, which can lead to considerable underestimation of the resistome in a given sample. Conclusions Our study shows that metaSPAdes and Trinity would be the preferable tools in terms of accuracy to recover correct genomic contexts around ARGs in metagenomic samples characterized by uneven coverages. Overall, the inability of assemblers to reconstruct long ARG-containing contigs has impacts on ARG quantification, suggesting that directly mapping reads to an ARG database should be performed as a complementary strategy to get accurate ARG abundance and diversity measures.

Abstract This research aimed to investigate the effects of dietary fructooligosaccharides (FOS) on attenuating the Aeromonas hydrophila (A. hydrophila)-induced oxidative stress and apoptosis in blunt snout bream Megalobrama amblycephala. Fish were divided into three groups as follows: C1 (Control), T1 (A. hydrophila), and T2 (A. hydrophila + 4 g/kg FOS). The results showed that the activities of antioxidant enzymes increased, the liver morphology had disorderly arrangement, and extensive cell necrosis occurred because of A. hydrophila-infection. While the dietary FOS improved the above-mentioned liver damage. Additionaly, FOS elevated mRNA levels of pro-apoptotic molecules, including caspase-8 and 9, and down-regulated mRNA levels of the anti-apoptotic molecule Bcl-2, which is triggered by A. hydrophila-infection. The transcriptome analysis showed that the oxidative stress-related DEGs pathways were activated in intestine of blunt snout bream by A. hydrophila-infection. The FOS-added group led to the enrichment of more pathways to health. Further WGCNA co-expression network analysis showed that the screened single genes were clustered into 49 modules. The two modules with the highest association to the five traits (10 hub genes) were chosen to build the network by combining the physiological and biochemical characteristic. In summary, this research offers a foundation for the exploring of A. hydrophila-restoration genes in dietary FOS, and also lays a theoretical foundation for aquaculture in the future.

Abstract Prediction of chromatin interactions from DNA sequence has been a significant research challenge in the last couple of years. Several solutions have been proposed, most of which are based on encoder-decoder architecture, where 1D sequence is convoluted, encoded into the latent representation, and then decoded using 2D convolutions into the Hi-C pairwise chromatin spatial proximity matrix. Those methods, while obtaining high correlation scores and improved metrics, produce Hi-C matrices that are artificial - they are blurred due to the deep learning model architecture. In our study, we propose the HiCDiffusion, sequence-only model that addresses this problem. We first train the encoder-decoder neural network and then use it as a component of the diffusion model - where we guide the diffusion using a latent representation of the sequence, as well as the final output from the encoder-decoder. That way, we obtain the high-resolution Hi-C matrices that not only better resemble the experimental results - improving the Fréchet inception distance by an average of 11 times, with the highest improvement of 56 times - but also obtain similar classic metrics to current state-of-the-art encoder-decoder architectures used for the task.

Abstract Background The Plant-specific AP2/ERF gene family encodes proteins involved in various biological and physiological processes. Although the genome of Coptis chinensis Franch, a plant producing benzylisoquinoline alkaloids (BIAs), has been sequenced at the chromosome level, studies on the AP2/ERF gene family in C. chinensis are lacking. Thus, a genome-wide identification of AP2/ERF gene family in C. chinensis was conducted to explore its role in BIAs biosynthesis. Results A total of 96 CcAP2/ERF genes were identified and categorized into five subfamilies, including 43 ERFs, 32 DREBs, 17 AP2s, 3 RAVs, and 1 Soloist, based on their structural domains. These CcAP2/ERF genes were unevenly distributed across nine chromosomes. Analysis of gene duplication events identified 17 CcAP2/ERF gene pairs in the genome, with 7 involved in tandem duplication events and 10 involved in segmental duplicate events, indicating that both types of duplications contributed to the expansion of the AP2/ERF gene family. The Ka/Ks ratio analysis suggested that the CcAP2/ERF gene family underwent strong purifying selection. Two phytohormones, methyl jasmonate and abscisic acid, were identified as potential key inducers of BIAs biosynthesis due to the cis-acting element prediction. Analysis of the spatial transcriptomic data revealed that 28 differentially expressed AP2/ERF genes had the highest or relatively higher expression levels in the rhizome, 17 of which positively correlated with the tissue-specific accumulation of BIAs. Further real-time PCR verification and protein-protein interaction analysis indicated that DREB1B might be one of the central regulators in a highly complex BIAs biosynthesis network. Conclusion These findings provide significant insight into the function of AP2/ERF genes in C. chinensis, particularly in the regulatory network of BIAs biosynthesis in C. chinensis. This study also identifies candidate genes for metabolic engineering to increase BIAs content in C. chinensis.

Abstract Background Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed. Results We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools. Conclusions This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches.

Abstract Background Transcription factors (TFs) bind to different parts of the genome in different types of cells, but it is usually assumed that the inherent DNA-binding preferences of a TF are invariant to cell type. Yet, there are several known examples of TFs that switch their DNA-binding preferences in different cell types, and yet more examples of other mechanisms, such as steric hindrance or cooperative binding, that may result in a “DNA signature” of differential binding. Results To survey this phenomenon systematically, we developed a deep learning method we call SigTFB (Signatures of TF Binding) to detect and quantify cell-type specificity in a TF’s known genomic binding sites. We used ENCODE ChIP-seq data to conduct a wide scale investigation of 169 distinct TFs in up to 14 distinct cell types. SigTFB detected statistically significant DNA binding signatures in approximately two-thirds of TFs, far more than might have been expected from the relatively sparse evidence in prior literature. We found that the presence or absence of a cell-type specific DNA binding signature is distinct from, and indeed largely uncorrelated to, the degree of overlap between ChIP-seq peaks in different cell types, and tended to arise by two mechanisms: using established motifs in different frequencies, and by selective inclusion of motifs for distint TFs. Conclusions While recent results have highlighted cell state features such as chromatin accessibility and gene expression in predicting TF binding, our results emphasize that, for some TFs, the DNA sequences of the binding sites contain substantial cell-type specific motifs.

Abstract Background Mutton quality is closely related to genetic variants and gene expression alterations during growth and development, resulting in differences in nutritional values, flavor, and odor. Results We first evaluated and compared the composition of crude protein, crude fat, cholesterol, amino acid (AA), and fatty acid (FA) in the longissimus dorsi muscle of Guizhou black goats (GZB, n = 5) and Yunshang black goats (YBG, n = 6). The contents of cholesterol and FA related to odor in GZB were significantly lower than that in YBG, while the concentrations of umami amino acids and intramuscular fat were significantly higher in GZB. Furthermore, structural variants (SVs) in the genomes of GZB (n = 30) and YBG (n = 11) were explored. It was found that some regions in Chr 10/12/18 were densely involved with a large number of SVs in the genomes of GZB and YBG. By setting F ST ≥ 0.25, we got 837 stratified SVs, of which 25 SVs (involved in 12 genes, e.g., CORO1A, CLIC6, PCSK2, and TMEM9) were limited in GZB. Functional enrichment analysis of 14 protein-coding genes (e.g., ENPEP, LIPC, ABCA5, and SLC6A15) revealed multiple terms and pathways related with metabolisms of AA, FA, and cholesterol. The SVs (n = 10) obtained by the whole genome resequencing were confirmed in percentages of 36.67 to 86.67% (n = 96) by PCR method. The SVa and SVd polymorphisms indicated a moderate negative correlation with HMGCS1 activity (n = 17). Conclusion This study is the first to comprehensively reveal potential SVs related to mutton nutritional values, flavor, and odor based on genomic compare between two black goat breeds with closely genetic relationship. The SVs generated in this study provide a data resource for deeper studies to understand the genomic characteristics and possible evolutionary outcomes with better nutritional values, flavor and extremely light odor.

Abstract Background The blunt snout bream (Megalobrama amblycephala) is an important economic freshwater fish in China with tender flesh and high nutritional value. With the cultivation of superior new varieties and the expansion of breeding scale, it becomes imperative to employ sex-control technology to cultivate monosexual populations of M. amblycephala, thereby preventing the deterioration of desirable traits. The development of specific markers capable of accurately identifying the sex of M. amblycephala would facilitate the determination of the genetic sex of the breeding population before gonad maturation, thereby expediting the processes of sex-controlled breeding of M. amblycephala. Results A whole-genome re-sequencing was performed for 116 females and 141 males M. amblycephala collected from nine populations. Seven candidate male-specific sequences were identified through comparative analysis of male and female genomes, which were further compared with the sequencing data of 257 individuals, and finally three male-specific sequences were generated. These three sequences were further validated by PCR amplification in 32 males and 32 females to confirm their potential as male-specific molecular markers for M. amblycephala. One of these markers showed potential applicability in M. pellegrini as well, enabling males to be identified using this specific molecular marker. Conclusions The study provides a high-efficiency and cost-effective approach for the genetic sex identification in two species of Megalobrama. The developed markers in this study have great potential in facilitating sex-controlled breeding of M. amblycephala and M. pellegrini, while also contributing valuable insights into the underlying mechanisms of fish sex determination.

Abstract Pheromone communication is a key mechanism by which the reproductive division of labor is maintained within insect communities. Understanding how pheromones evolved to regulate social behavior requires knowledge of the molecular regulation of their production. However, even in cases where pheromones were identified, our understanding of their biosynthesis and molecular regulation remains limited. Bumble bees provide a unique system to explore pheromone biosynthesis since workers produce ester sterility signals in their Dufour’s gland that differ from gyne-specific esters and are not produced by queens. These esters are hypothesized to be produced in the exocrine gland where they are stored, and indeed queens, gynes and workers differ significantly in the expression of Dufour’s gland genes coding to enzymes involved in the biosynthesis of esters. However, a previous transcriptome analysis revealed no gene expression differences in the Dufour’s gland of workers despite differences in both ester production and ovarian activation, suggesting that ester production may be regulated lower down. Proteomics of the Dufour’s gland of queens, gynes, and workers recovered over 2400 proteins and broadly matched the previous RNAseq data. However, more than 100 differentially expressed proteins were found between the worker groups, including key enzymes in fatty acid biosynthesis, indicating that the regulation of reproductive signal biosynthesis in workers is done post-transcription. Overall, our data provide evidence that pheromone biosynthesis in the Dufour’s gland is caste specific, that gynes and workers are likely using different enzymes to make their respective wax esters, and that the regulation on pheromone production in queens, gynes and workers is likely done at different regulatory levels, with workers signals being subjected to regulation at the protein level.

Abstract Non-typhoidal Salmonella (NTS) is one of the top causes of diarrhea worldwide. Ceftriaxone is commonly recommended as the initial treatment option for Salmonella infections due to its antibacterial effectiveness. The objective of this study was to investigate the molecular epidemiological characteristics of NTS and to compare the phenotypic and genotypic profiles of antimicrobial resistance in multidrug-resistant Salmonella strains by sequencing 329 NTS strains collected from a county-level hospital between 2018 and 2021. Multi-locus sequence typing (MLST), antimicrobial resistance genes and plasmid types were identified by BacWGSTdb 2.0 webserver. Phylogenetic analysis of all NTS strains was carried out using Snippy and Gubbins software. The transferability of ceftriaxone resistant plasmids was confirmed through plasmid conjugation assays, and verified by S1-PFGE-Southern blot assays. The predominant serotypes among all NTS strains were Typhimurium (161/329), Enteritidis (49/329) and London (45/329). The most common sequence type observed was ST34 (86/329), followed by ST19 (72/329) and ST11 (47/329). The antimicrobial resistance of Salmonella to a wide range of antimicrobials showed an overall increase. Out of these 37 (11.24%) ceftriaxone-resistant strains, with the majority of them (33/37) being bla CTX−M. The predominant plasmid types identified were IncHI2 (14/21) and IncI1 (6/21), ranging in size from 70 kb to 360 kb. The conjugation efficiency was calculated with the high conjugation efficiency of 1.1 × 10− 5 to 9.3 × 10− 2. The strains varied widely, ranging from 3 to 45,024 single nucleotide polymorphisms (SNPs). There are close linkages observed among the predominant lineage, with an average of 78 SNPs between each pair of ST34 strains. The findings contribute to our understanding of the transmission and resistance mechanisms of multidrug-resistant Salmonella, thereby facilitating the development of effective control strategies.