Genomic and transcriptomic landscape to decipher the genetic basis of hyperpigmentation in Lanping black-boned sheep (Ovis aries).

BACKGROUND: Lanping black-boned sheep (LPB) represent a distinctive mammalian species characterized by hyperpigmentation, resulting in black bone and muscle features, in contrast to their conventional counterparts exhibiting red muscle and white bone. The genetic basis underlying LPB hyperpigmentation has remained enigmatic. METHODS: In this study, we conducted whole-genome sequencing of 100 LPB and 50 Lanping normal sheep (LPN), and integrated this data with 421 sequenced datasets from wild and domestic sheep, shedding light on the genetic backdrop and genomic variations associated with LPB. Furthermore, we performed comparative RNA-Seq analysis using liver sample to pinpoint genes implicated in the pigmentation process. We generated a comprehensive dataset comprising 97,944,357 SNPs from 571 sheep, facilitating an in-depth exploration of genetic factors. RESULTS: Population genetic structure analysis revealed that the LPB breed traces its origin back to LPN, having evolved into a distinct breed. The integration of positively selected genes with differentially expressed genes identified two candidates, ERBB4 and ROR1, potentially linked to LPB hyperpigmentation. Comparative analysis of ERBB4 and ROR1 mRNA relative expression levels in liver, spleen, and kidney tissues of LPB, in comparison to Diqing sheep, revealed significant upregulation, except for ERBB4 in the liver. Gene expression heatmaps further underscored marked allelic frequency disparities in different populations. CONCLUSION: Our findings establish the evolutionary lineage of the LPB breed from LPN and underscore the involvement of ERBB4 and ROR1 genes in melanin synthesis. These results enhance our comprehension of the molecular basis of hyperpigmentation and contribute to a more comprehensive depiction of sheep diversity.

Evidence for the occurrence of Acinetobacter faecalis in cattle feces and its emended description.

This study provides an emended description of Acinetobacter faecalis, a species previously described based on a single isolate (YIM 103518T) from elephant feces in China. Our emended description is based on 15 novel isolates conspecific with the A. faecalis type strain, obtained from eight cattle farms in the Czech Republic. The A. faecalis strains have relatively small genomes (≈2.5-2.7 Mbp), with a GC content of 36.3-36.7 mol%. Core genome-based phylogenetic analysis showed that the 15 strains, together with the type strain of A. faecalis, form a distinct and internally coherent phylogroup within the genus. Pairwise genomic ANIb values for the 16 A. faecalis strains were 97.32-99.04 %, while ANIb values between the genomes of the 16 strains and those of the other Acinetobacter spp. were ≤ 86.2 %. Analysis of whole-cell MALDI-TOF mass spectra supported the distinctness and cohesiveness of the taxon. The A. faecalis strains could be differentiated from the other validly named Acinetobacter spp. by the absence of hemolytic activity along with their ability to grow at 37 °C and on L-aspartate, ethanol, and L-glutamate but not at 41 °C or on adipate or 2,3-butanediol. Reduced susceptibility to sulfamethoxazole, trimethoprim and/or streptomycin was shown in eight strains, along with the presence of corresponding antibiotic resistance genes. In conclusion, this study provides a comprehensive description of A. faecalis and demonstrates its occurrence in cattle feces. Though the ecological role of A. faecalis remains unknown, our results show its ability to acquire antibiotic resistance genes, likely as an adaptation to antibiotic selection pressure in livestock farms.

Descriptive epidemiology and phylogenetic analysis of highly pathogenic avian influenza H5N1 clade 2.3.4.4b in British Columbia (B.C.) and the Yukon, Canada, September 2022 to June 2023.

Surveillance data from wildlife and poultry was used to describe the spread of highly pathogenic avian influenza (HPAI) H5N1 clade 2.3.4.4b in British Columbia (B.C.) and the Yukon, Canada from September 2022 - June 2023 compared to the first 'wave' of the outbreak in this region, which occurred April - August 2022, after the initial viral introduction. Although the number of HPAI-positive poultry farms and wildlife samples was greater in 'Wave 2', cases were more tightly clustered in southwestern B.C. and the most commonly affected species differed, likely due to an influx of overwintering waterfowl in the area. Eight HPAI genetic clusters, representing seven genotypes and two inter-continental viral incursions, were detected, with significant variation in the relative abundance of each cluster between the waves. Phylogenetic data suggests multiple spillover events from wild birds to poultry and mammals but could not rule out transmission among farms and among mammals.

EFSA statement on the requirements for whole genome sequence analysis of microorganisms intentionally used in the food chain.

Microorganisms, genetically modified or not, may be used in the food chain either as active agents, biomasses or as production organisms of substances of interest. The placement of such microorganisms or their derived substances/products in the European market may be subject to a premarket authorisation process. The authorisation process requires a risk assessment in order to establish the safety and/or the efficacy of the microorganism(s) when used in the food chain as such, as biomasses or as production strains. This includes a full molecular characterisation of the microorganism(s) under assessment. For certain regulated products, the use of whole genome sequence (WGS) data of the microorganism is established as a requirement for the risk assessment. In this regard, data obtained from WGS analysis can provide information on the unambiguous taxonomic identification of the strains, on the presence of genes of concern (e.g. those encoding virulence factors, resistance to antimicrobials of clinical relevance for humans and animals, production of harmful metabolites or of clinically relevant antimicrobials) and on the characterisation of genetic modification(s) (where relevant). This document provides recommendations to applicants on how to describe and report the results of WGS analyses in the context of an application for market authorisation of a regulated product. Indications are given on how to perform genome sequencing and the quality criteria/thresholds that should be reached, as well as the data and relevant information that need to be reported, if required. This updated document replaces the EFSA 2021 Statement and reflects the current knowledge in technologies and methodologies to be used to generate and analyse WGS data for the risk assessment of microorganisms.

An increased prevalence of carbapenem-resistant hypervirulent Klebsiella pneumoniae associated with the COVID-19 pandemic.

BACKGROUND: Klebsiella pneumoniae (Kp) is a common community-acquired and nosocomial pathogen. Carbapenem-resistant and hypervirulent (CR-hvKp) variants can emerge rapidly within healthcare facilities and impacted by other infectious agents such as COVID-19 virus. METHODS: To understand the impact of COVID-19 virus on the prevalence of CR-hvKp, we accessed Kp genomes with corresponding metadata from GenBank. Sequence types (STs), antimicrobial resistance genes, and virulence genes, and those scores and CR-hvKp were identified. We analyzed population diversity and phylogenetic characteristics of five most common STs, measured the prevalence of CR-hvKp, identified CR-hvKp subtypes, and determined associations between carbapenem resistance gene subtypes with STs and plasmid types. These variables were compared pre- and during the COVID-19 pandemic. FINDINGS: The proportion of CR-hvKp isolates increased within multiple STs in different continents during the COVID-19 pandemic and persistent CR-hvKp subtypes were found in common STs. blaKPC was dominant in CG258, blaKPC-2 was detected in 97 % of the ST11 CR-hvKp, blaNDM subtypes were prominent in ST147 (87.4 %) and ST307 (70.8 %); blaOXA-48 and its subtypes were prevalent in ST15 (80.5 %). The possession of carbapenemase genes was different among subclades from different origins in different periods of time within each ST. IncFIB/IncHI1B hybrid plasmids contained virulence genes and carbapenemase genes and were predominant in ST147 (67.37 %) and ST307 (56.25 %). INTERPRETATION: The prevalence of CR-hvKp increased during the COVID-19 pandemic, which was evident by an increase in local endemic clones. This process was facilitated by the convergence of plasmids containing carbapenemase genes and virulence genes. These findings have implications for the appropriate use of antimicrobials and infection prevention and control during outbreaks of respiratory viruses and pandemic management.

Antimicrobial resistance profiles and genome characteristics of Klebsiella isolated from the faeces of neonates in the neonatal intensive care unit.

Introduction. Klebsiella spp. are important bacteria that colonize the human intestine, especially in preterm infants; they can induce local and systemic disease under specific circumstances, including inflammatory bowel disease, necrotizing enterocolitis and colorectal cancer.Hypothesis. Klebsiella spp. colonized in the intestine of the neonates in the neonatal intensive care unit (NICU) may be associated with disease and antibiotic resistance, which will be hazardous to the children.Aim. Our aim was to know about the prevalence, antimicrobial resistance and genome characteristics of Klebsiella spp. in neonate carriers.Methodology. Genome sequencing and analysis, and antimicrobial susceptibility testing were mainly performed in this study.Results. The isolation rates of Klebsiella spp. strains were 3.7% (16/436) in 2014 and 4.3% (18/420) in 2021. Cases with intestinal-colonized Klebsiella spp. were mainly infants with low birth weights or those with pneumonia or hyperbilirubinemia. According to the core-pan genomic analysis, 34 stains showed gene polymorphism and a sequence type (ST) of an emerging high-risk clone (ST11). Eight strains (23.5%) were found to be resistant to 2 or more antibiotics, and 46 genes/gene families along with nine plasmids were identified that conferred resistance to antibiotics. In particular, the two strains were multidrug-resistant. Strain A1256 that is related to Klebsiella quasipneumoniae subsp. similipneumoniae was uncommon, carrying two plasmids similar to IncFII and IncX3 that included five antibiotic resistance genes.Conclusion. The prevention and control of neonatal Klebsiella spp. colonization in the NICU should be strengthened by paying increased attention to preventing antimicrobial resistance in neonates.

Optimizing Polyhydroxyalkanoate production using a novel Bacillus paranthracis isolate: A response surface methodology approach.

Microorganisms have emerged as promising resources for producing economical and sustainable bioproducts like Polyhydroxyalkanoate (PHA), a biodegradable polymer that can replace synthetic plastics. In this study, we screened a novel isolate, Bacillus paranthracis RSKS-3 strain, to produce PHA from sewage water, identifying it using Whole Genome Sequence. This study represents the first report on optimizing PHA production using B. paranthracis RSKS-3, employing Design Expert 12.0 software. Our findings reveal that four factors (temperature, inoculum size, potassium dihydrogen phosphate, and magnesium sulfate) significantly affect PHA production in the Plackett-Burman design experiment. Through Response Surface Methodology, we optimized PHA production to 0.647 g/L with specific values for potassium dihydrogen phosphate (0.55 %), inoculum size (3 %), magnesium sulfate (0.055 %), and a temperature of 35 °C, in agreement with the predicted value of 0.630 g/L. This optimization resulted in a substantial 13.29-fold increase in PHA production from 0.34 g/L to 4.52 g/L, underscoring the promising role of B. paranthracis RSKS-3 in eco-friendly PHA production and advancing sustainable bioproduct development.

First Data on WGS-Based Typing and Antimicrobial Resistance of Human Salmonella Enteritidis Isolates in Greece.

Salmonella enterica subsp. enterica serotype Enteritidis (S. Enteritidis) is one of the major causes of foodborne infections and is responsible for many national and multi-country foodborne outbreaks worldwide. In Greece, human salmonellosis is a mandatory notifiable disease, with laboratory surveillance being on a voluntary basis. This study aims to provide the first insights into the genetic characteristics and antimicrobial resistance profiles of 47 S. Enteritidis human isolates using whole-genome sequencing (WGS) technology. The S. Enteritidis population was mainly resistant to fluoroquinolones due to gyrA point mutations, whereas one isolate presented a multi-resistant plasmid-mediated phenotype. ST11 was the most frequent sequence type, and phylogenetic analysis through the cgMLST and SNP methods revealed considerable genetic diversity. Regarding virulence factors, 8 out of the 24 known SPIs and C63PI were detected. Due to the observed variability between countries, it is of utmost importance to record the circulating S. Enteritidis strains' structure and genomic epidemiology at the national level. WGS is a valuable tool that is revolutionizing our approach to Salmonella by providing a deeper understanding of these pathogens and their impact on human health.

Variant calling in polyploids for population and quantitative genetics.

Advancements in genome assembly and sequencing technology have made whole genome sequence (WGS) data and reference genomes accessible to study polyploid species. Compared to popular reduced-representation sequencing approaches, the genome-wide coverage and greater marker density provided by WGS data can greatly improve our understanding of polyploid species and polyploid biology. However, biological features that make polyploid species interesting also pose challenges in read mapping, variant identification, and genotype estimation. Accounting for characteristics in variant calling like allelic dosage uncertainty, homology between subgenomes, and variance in chromosome inheritance mode can reduce errors. Here, I discuss the challenges of variant calling in polyploid WGS data and discuss where potential solutions can be integrated into a standard variant calling pipeline.

Isolation and Identification of Lactic Acid Bacteria from Environmental Samples.

Isolation is the first and crucial step to study possible roles of lactic acid bacteria (LAB) in environment, especially in food fermentation. It is also important to use the organisms for further application. LABs are diverse bacterial group and have diverse growth characteristics. Culture condition of LAB is thus varied, and selection of a suitable culture medium is essential for the purposes. Identification is also important step, since certain desirable and undesirable characteristics are shared within species. Identification was classically carried out by phenotypic characteristics but is usually performed for DNA sequence-based approaches. 16S rRNA gene sequencing is generally used for the identification, and sequencing of housekeeping genes is used when needed. In addition, identification based on whole-genome sequence similarities was well established during the last decade. Here, we describe methods for isolation and identification of LAB briefly.

Bioinformatics Analysis of the Genome of Rhodococcus rhodochrous IEGM 1362, an (-)-Isopulegol Biotransformer.

A genome of Rhodococcus rhodochrous IEGM 1362 was sequenced and annotated. This strain can transform monoterpene alcohol (-)-isopulegol with the formation of two novel pharmacologically promising metabolites. Nine genes encoding cytochrome P450, presumably involved in (-)-isopulegol transformation, were found in the genome of R. rhodochrous IEGM 1362. Primers and PCR conditions for their detection were selected. The obtained data can be used for the further investigation of genes encoding enzymes involved in monoterpene biotransformation.

Whole-genome resequencing identifies candidate genes associated with heat adaptation in chickens.

The wide distribution and diverse varieties of chickens make them important models for studying genetic adaptation. The aim of this study was to identify genes that alter heat adaptation in commercial chicken breeds by comparing genetic differences between tropical and cold-resistant chickens. We analyzed whole-genome resequencing data of 186 chickens across various regions in Asia, including the following breeds: Bian chickens (B), Dagu chickens (DG), Beijing-You chickens (BY), and Gallus gallus jabouillei from China; Gallus gallus murghi from India; Vietnam native chickens (VN); Thailand native chickens (TN) and Gallus gallus spadiceus from Thailand; and Indonesia native chickens (IN), Gallus gallus gallus, and Gallus gallus bankiva from Indonesia. In total, 5,454,765 SNPs were identified for further analyses. Population genetic structure analysis revealed that each local chicken breed had undergone independent evolution. Additionally, when K = 5, B, BY, and DG chickens shared a common ancestor and exhibited high levels of inbreeding, suggesting that northern cold-resistant chickens are likely the result of artificial selection. In contrast, the runs of homozygosity (ROH) and the ROH-based genomic inbreeding coefficient (FROH) results for IN, TN, and VN chickens showed low levels of inbreeding. Low population differentiation index values indicated low differentiation levels, suggesting low genetic diversity in tropical chickens, implying increased vulnerability to environmental changes, decreased adaptability, and disease resistance. Whole-genome selection sweep analysis revealed 69 candidate genes, including LGR4, G6PC, and NBR1, between tropical and cold-resistant chickens. The genes were further subjected to GO and KEGG enrichment analyses, revealing that most of the genes were primarily enriched in biological synthesis processes, metabolic processes, central nervous system development, ion transmembrane transport, and the Wnt signaling pathway. Our study identified heat adaptation genes and their functions in chickens that primarily affect chickens in high-temperature environments through metabolic pathways. These heat-resistance genes provide a theoretical basis for improving the heat-adaptation capacity of commercial chicken breeds.

Comprehensive physiological and genomic characterization of a potential probiotic strain, Lactiplantibacillus plantarum ILSF15, isolated from the gut of tribes of Odisha, India.

Characterizing probiotic features of organisms isolated from diverse environments can lead to the discovery of novel strains with promising functional features and health attributes. The present study attempts to characterize a novel probiotic strain isolated from the gut of the tribal population of Odisha, India. Based on 16S rRNA-based phylogeny, the strain was identified as a species of the Lactiplantibacillus genus and was named Lactiplantibacillus plantarum strain ILSF15. The current investigation focuses on elucidating this strain's genetic and physiological properties associated with probiotic attributes such as biosafety risk, host adaptation/survival traits, and beneficial functional features. The novel strain was observed, in vitro, exhibiting features such as acid/bile tolerance, adhesion to the host enteric epithelial cells, cholesterol assimilation, and pathogen exclusion, indicating its ability to survive the harsh environment of the human GIT and resist the growth of harmful microorganisms. Additionally, the L. plantarum ILSF15 strain was found to harbor genes associated with the metabolism and synthesis of various bioactive molecules, including amino acids, carbohydrates, lipids, and vitamins, highlighting the organism's ability to efficiently utilize diverse resources and contribute to the host's nutrition and health. Several genes involved in host adaptation/survival strategies and host-microbe interactions were also identified from the ILSF15 genome. Moreover, L. plantarum strains, in general, were found to have an open pangenome characterized by high genetic diversity and the absence of specific lineages associated with particular habitats, signifying its versatile nature and potential applications in probiotic and functional food industries.

Whole-genome sequence and annotation of Diaporthe australafricana strain CBS 118571.

We report the whole-genome sequence of Diaporthe australafricana Crous & J.M. van Niekerkusing using Oxford Nanopore long-read sequencing and Illumina short-read sequencing. The hybrid genome consists of 11 contigs with a total length of 53.509 Mb, and a GC content of 52.40%.

Whole-genome resequencing of native and imported dairy goat identifies genes associated with productivity and immunity.

Understanding the differences in genetic variation between local Chinese dairy goat breeds and imported breeds can help germplasm innovation and molecular breeding. However, the research is limited in this area. In this study, whole-genome resequencing data from 134 individuals of both local and imported dairy goat breeds were analyzed, and their differences in genomic genetic variation, genetic diversity, and population structure were subsequently identified. We also screened candidate genes associated with important traits of dairy goats such as milk production (STK3, GHR, PRELID3B), reproduction (ATP5E), growth and development (CTSZ, GHR), and immune function (CTSZ, NELFCD). Furthermore, we examined allele frequency distributions for the genes of interest and found significant differences between the two populations. This study provides valuable resources for the study of genetic diversity in dairy goats and lays the foundation for the selective breeding of dairy goats in the future.

An overview of next generation sequencing strategies and genomics tools used for tuberculosis research.

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research.

Complete genome sequences of paired isogenic Burkholderia pseudomallei isolated from a Thai pediatric patient with a urinary tract infection.

Burkholderia pseudomallei is the causative agent of melioidosis, the disease endemic in Southeast Asia and northern Australia. We report complete genome sequences of paired isogenic B. pseudomallei isolated from a 12-year-old Thai male presenting with acute urinary tract infection before (SCBP001) and after (SCBP007) a decrease in susceptibility to ceftazidime.

Transposase-Assisted RNA/DNA Hybrid Co-Tagmentation for Target Meta-Virome of Foodborne Viruses.

Foodborne diseases are major public health problems globally. Metagenomics has emerged as a widely used tool for pathogen screening. In this study, we conducted an updated Tn5 transposase-assisted RNA/DNA hybrid co-tagmentation (TRACE) library construction approach. To address the detection of prevalent known foodborne viruses and the discovery of unknown pathogens, we employed both specific primers and oligo-T primers during reverse transcription. The method was validated using clinical samples confirmed by RT-qPCR and compared with standard RNA-seq library construction methods. The mapping-based approach enabled the retrieval of nearly complete genomes (>95%) for the majority of virus genome segments (86 out of 88, 97.73%), with a mean coverage depth of 21,494.53× (ranging from 77.94× to 55,688.58×). Co-infection phenomena involving prevalent genotypes of Norovirus with Astrovirus and Human betaherpesvirus 6B were observed in two samples. The updated TRACE-seq exhibited superior performance in viral reads percentages compared to standard RNA-seq library preparation methods. This updated method has expanded its target pathogens beyond solely Norovirus to include other prevalent foodborne viruses. The feasibility and potential effectiveness of this approach were then evaluated as an alternative method for surveilling foodborne viruses, thus paving the way for further exploration into whole-genome sequencing of viruses.

Unraveling the Genomic Evolution of Dengue Virus Serotype 1: A Case Study from Yantai, China.

In August 2023, we identified a case of dengue fever in Yantai City, which was imported from Xishuangbanna, China. To investigate its evolutionary history and population dynamics, we utilized the metatranscriptomic method to obtain the virus' whole genome sequence. Together with 367 selected dengue virus whole genome sequences from the NCBI database, we constructed a time-scaled Maximum Clade Credibility (MCC) tree. We found that our sequence exhibited a high homology with a sequence of DENV1 (OR418422.1) uploaded by the Guangzhou Center for Disease Control and Prevention in 2023, with an estimated divergence time around 2019 (95% HPD: 2017-2023), coinciding with the emergence of SARS-CoV-2. The DENV strain obtained in this study belongs to genotype I of DENV1. Its ancestors experienced a global epidemic around 2005 (95% HPD: 2002-2010), and its progeny strains have spread extensively in Southeast Asia and China since around 2007 (95% HPD: 2006-2011). The Bayesian skyline plot indicates that the current population of DENV1 has not been affected by SARS-CoV-2 and is expected to maintain stable transmission. Hence, it is imperative to track and monitor its epidemiological trends and genomic variations to prevent potential large-scale outbreaks in the post-SARS-CoV-2 era.

Relevance of genetic causes and environmental adaptation of Cronobacter spp. isolated from infant and follow-up formula production factories and retailed products in China: A 7-year period of continuous surveillance based on genome-wide analysis.

The possible contamination routes, environmental adaptation, and genetic basis of Cronobacter spp. in infant and follow-up formula production factories and retailed products in mainland China have been determined by laboratory studies and whole-genome comparative analysis in a 7-year nationwide continuous surveillance spanning from 2012 to 2018. The 2-year continuous multicenter surveillance of the production process (conducted in 2013 and 2014) revealed that the source of Cronobacter spp. in the dry-blending process was the raw dry ingredients and manufacturing environment (particularly in the vibro sieve and vacuum cleaner), while in the combined process, the main contamination source was identified as the packing room. It is important to note that, according to the contamination control knowledge obtained from the production process surveillance, the contamination rate of retail powdered infant formula (PIF) and follow-up formula (FUF) products in China decreased significantly from 2016 onward, after improving the hygiene management practices in factories. The prevalence of Cronobacter spp. in retailed PIF and FUF in China in 2018 was dramatically reduced from 1.55 % (61/3925, in 2012) to an average as low as 0.17 % (13/7655 in 2018). Phenotype determination and genomic analysis were performed on a total of 90 Cronobacter spp. isolates obtained from the surveillance. Of the 90 isolates, only two showed resistance to either cefazolin or cefoxitin. The multilocus sequence typing results revealed that C. sakazakii sequence type 1 (ST1), ST37, and C. malonaticus ST7 were the dominant sequence types (STs) collected from the production factories, while C. sakazakii ST1, ST4, ST64, and ST8 were the main STs detected in the retailed PIF and FUF nationwide. One C. sakazakii ST4 isolate (1.1 %, 1/90) had strong biofilm-forming ability and 13 isolates (14.4 %, 13/90) had weak biofilm-forming ability. Genomic analysis revealed that Cronobacter spp. have a relatively stable core-genome and an increasing pan-genome size. Plasmid IncFIB (pCTU3) was prevalent in this genus and some contained 14 antibacterial biocide- and metal-resistance genes (BMRGs) including copper, silver, and arsenic resistant genes. Plasmid IncN_1 was predicted to contain 6 ARGs. This is the first time that a multi-drug resistance IncN_1 type plasmid has been reported in Cronobacter spp. Genomic variations with respect to BMRGs, virulence genes, antimicrobial resistance genes (ARGs), and genes involved in biofilm formation were observed among strains of this genus. There were apparent differences in copies of bcsG and flgJ between the biofilm-forming group and non-biofilm-forming group, indicating that these two genes play key roles in biofilm formation. The findings of this study have improved our understanding of the contamination characteristics and genetic basis of Cronobacter spp. in PIF and FUF and their production environment in China and provide important guidance to reduce contamination with this pathogen during the production of PIF and FUF.