Evidence for the utility of quantum computing before fault tolerance.

Quantum computing promises to offer substantial speed-ups over its classical counterpart for certain problems. However, the greatest impediment to realizing its full potential is noise that is inherent to these systems. The widely accepted solution to this challenge is the implementation of fault-tolerant quantum circuits, which is out of reach for current processors. Here we report experiments on a noisy 127-qubit processor and demonstrate the measurement of accurate expectation values for circuit volumes at a scale beyond brute-force classical computation. We argue that this represents evidence for the utility of quantum computing in a pre-fault-tolerant era. These experimental results are enabled by advances in the coherence and calibration of a superconducting processor at this scale and the ability to characterize1 and controllably manipulate noise across such a large device. We establish the accuracy of the measured expectation values by comparing them with the output of exactly verifiable circuits. In the regime of strong entanglement, the quantum computer provides correct results for which leading classical approximations such as pure-state-based 1D (matrix product states, MPS) and 2D (isometric tensor network states, isoTNS) tensor network methods2,3 break down. These experiments demonstrate a foundational tool for the realization of near-term quantum applications4,5.

Establishment of a secondary infection laboratory model of Echinococcus shiquicus metacestode using BALB/c mice and Mongolian jirds (Meriones unguiculatus).

Echinococcus shiquicus is peculiar to the Qinghai­Tibet plateau of China. Research on this parasite has mainly focused on epidemiological surveys and life cycle studies. So far, limited laboratory studies have been reported. Here, experimental infection of E. shiquicus metacestode in BALB/c mice and Mongolian jirds (Meriones unguiculatus) was carried out to establish alternative laboratory animal models. Intraperitoneal inoculation of metacestode material containing protoscoleces (PSCs) obtained from infected plateau pikas were conducted on BALB/c mice. Furthermore, metacestode material without PSCs deriving from infected BALB/c mice was intraperitoneally inoculated to Mongolian jirds. Experimental animals were dissected for macroscopic and histopathological examination. The growth of cysts in BALB/c mice was infiltrative, and they invaded the murine entire body. Most of the metacestode cysts were multicystic, but a few were unilocular. The cysts contained sterile vesicles, which had no PSCs. The metacestode materials were able to successfully infect new mice. In the jirds model, E. shiquicus cysts were typically formed freely in the peritoneal cavity; the majority of these cysts were free while a small portion adhered loosely to nearby organs. The proportion of fertile cysts was high, and contained many PSCs. The PSCs produced in Mongolian jirds also successfully infected new ones, which confirms that jirds can serve as an alternative experimental intermediate host. In conclusion, a laboratory animal infection was successfully established for E. shiquicus using BALB/c mice and Mongolian jirds. These results provide new models for the in-depth study of Echinococcus metacestode survival strategy, host interactions and immune escape mechanism.

Update on the genetic diversity and population structure of Echinococcus granulosus in Gansu Province, Tibet Autonomous Region, and Xinjiang Uygur Autonomous Region, Western China, inferred from mitochondrial cox1, nad1, and nad5 sequences.

The identification of additional Echinococcus granulosus sensu lato (s.l.) complex species/genotypes in recent years raises the possibility that there might be more variation among this species in China than is currently understood. The aim of this study was to explore intra- and inter-species variation and population structure of Echinococcus species isolated from sheep in three areas of Western China. Of the isolates, 317, 322, and 326 were successfully amplified and sequenced for cox1, nad1, and nad5 genes, respectively. BLAST analysis revealed that the majority of the isolates were E. granulosus s.s., and using the cox1, nad1, and nad5 genes, respectively, 17, 14, and 11 isolates corresponded to Elodea canadensis (genotype G6/G7). In the three study areas, G1 genotypes were the most prevalent. There were 233 mutation sites along with 129 parsimony informative sites. A transition/transversion ratio of 7.5, 8, and 3.25, respectively, for cox1, nad1, and nad5 genes was obtained. Every mitochondrial gene had intraspecific variations, which were represented in a star-like network with a major haplotype with observable mutations from other distant and minor haplotypes. The Tajima's D value was significantly negative in all populations, indicating a substantial divergence from neutrality and supporting the demographic expansion of E. granulosus s.s. in the study areas. The phylogeny inferred by the maximum likelihood (ML) method using nucleotide sequences of cox1-nad1-nad5 further confirmed their identity. The nodes assigned to the G1, G3, and G6 clades as well as the reference sequences utilized had maximal posterior probability values (1.00). In conclusion, our study confirms the existence of a significant major haplotype of E. granulosus s.s. where G1 is the predominant genotype causing of CE in both livestock and humans in China.

Identification of Na+/K+-ATPase Inhibitor Bufalin as a Novel Pseudorabies Virus Infection Inhibitor In Vitro and In Vivo.

Pseudorabies virus (PRV), an alpha herpesvirus, induces significant economic losses to the swine industry and infects multiple kinds of animals. Therefore, it is of great importance to explore anti-PRV compounds. In this study, to explore the anti-PRV compounds, a library of natural compounds was screened through a cell-based ELISA assay, and it was discovered that bufalin, a Na+/K+-ATPase inhibitor, had a robust inhibitory effect on PRV replication. A time-of-addition experiment and temperature-shift assay showed that bufalin significantly inhibited the entry stage of PRV. NaCl- or KCl-treatment showed that NaCl could enhance the inhibitory effect of bufalin on PRV replication, whereas there was no significant effect under the treatment of KCl. Meanwhile, it was also found that bufalin possessed antiviral activity against other alpha herpesviruses, including human herpes simplex virus type 1 (HSV-1) and chicken Marek's disease virus (MDV). Finally, it was found that bufalin could decrease the viral load in multiple tissues, and reduce the morbidity and mortality in PRV-challenged BALB/c mice. Overall, our findings demonstrated that bufalin has the potential to be developed as an anti-PRV compound.

dCas9-BE3 and dCas12a-BE3 Systems Mediated Base Editing in Kiwifruit Canker Causal Agent Pseudomonas syringae pv. actinidiae.

Pseudomonas syringae pv. actinidiae (Psa) causes bacterial canker of kiwifruit with heavy economic losses. However, little is known about the pathogenic genes of Psa. CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas-mediated genome editing technology has dramatically facilitated the characterization of gene function in various organisms. However, CRISPR genome editing could not be efficiently employed in Psa due to lacking homologous recombination repair. The base editor (BE) system, which depends on CRISPR/Cas, directly induces single nucleoside C to T without homology recombination repair. Here, we used dCas9-BE3 and dCas12a-BE3 systems to create substitutions of C to T and to convert CAG/CAA/CGA codons to stop codons (TAG/TAA/TGA) in Psa. The dCas9-BE3 system-induced single C-to-T conversion frequency of 3 to 10 base positions ranged from 0% to 100%, with a mean of 77%. The dCas12a-BE3 system-induced single C-to-T conversion frequency of 8 to 14 base positions in the spacer region ranged from 0% to 100%, with a mean of 76%. In addition, a relatively saturated Psa gene knockout system covering more than 95% of genes was developed based on dCas9-BE3 and dCas12a-BE3, which could knock out two or three genes at the same time in the Psa genome. We also found that hopF2 and hopAO2 were involved in the Psa virulence of kiwifruit. The HopF2 effector can potentially interact with proteins such as RIN, MKK5, and BAK1, while the HopAO2 effector can potentially interact with the EFR protein to reduce the host's immune response. In conclusion, for the first time, we established a PSA.AH.01 gene knockout library that may promote research on elucidating the gene function and pathogenesis of Psa.

Simple yet effective analysis of waveguide mode symmetry: generalized eigenvalue approach based on Maxwell's equations.

Particular waveguide structures and refractive index distribution can lead to specified degeneracy of eigenmodes. To obtain an accurate understanding of this phenomenon, we propose a simple yet effective approach, i.e., generalized eigenvalue approach based on Maxwell's equations, for the analysis of waveguide mode symmetry. In this method, Maxwell's equations are reformulated into generalized eigenvalue problems. The waveguide eigenmodes are completely determined by the generalized eigenvalue problem given by two matrices (M, N), where M is 6 × 6 waveguide Hamiltonian and N is a constant singular matrix. Close examination shows that N usually commute with the corresponding matrix of a certain symmetry operation, thus the waveguide eigenmode symmetry is essentially determined by M, in contrast to the tedious and complex procedure given in the previous work [Opt. Express25, 29822 (2017)10.1364/OE.25.029822]. Based on this new approach, we discuss several symmetry operations and the corresponding symmetries including chiral, parity-time reversal, rotation symmetry, wherein the constraints of symmetry requirements on material parameters are derived in a much simpler way. In several waveguides with balanced gain and loss, anisotropy, and geometrical symmetry, the analysis of waveguide mode symmetry based on our simple yet effective approach is consistent with previous results, and shows perfect agreement with full-wave simulations.

Pathogenicity comparison between QX-type and Mass-type infectious bronchitis virus to different segments of the oviducts in laying phase.

BACKGROUND: The QX-type infectious bronchitis virus (IBV) has become the predominant genotype worldwide in recent years and has caused serious economic losses to the chicken industry. The most significant feature of QX IBV is that its infection in the early growing stage can cause abnormal oviduct development, resulting in a high proportion of 'false layers' in poultry flocks of laying hens and breeders. However, few studies have evaluated whether infections of QX-type IBV in laying stages can also cause severe pathological changes in the oviduct. METHODS: In this study, 300-day-old specific-pathogen-free chickens were infected either with the QX-type strain QXL or Massachusetts (Mass)-type strain M41 to compare their pathogenicity on different segments of the oviduct. RESULTS: Both the QXL and M41 strains successfully replicated in all segments of the oviduct; however, the QXL strain was more highly distributed in mucosal layer and caused severe lesions in the lamina propria, including interstitial dilation, inflammatory cell infiltration, and distinct expansion of tubular glands. Moreover, the QXL strain induced high expression of proinflammatory cytokines and cytotoxic molecules in the majority of segments in the oviduct. Further research found that the QXL strain may affected the formation of shell membranes and eggshells by inhibiting the expression of type I collagen and CaBP-D28k. CONCLUSIONS: Our results indicate that the QX-type IBV is more pathogenic than Mass-type IBV to oviduct in laying phase. Collectively, these findings provide detailed information on the pathological changes in different segments of the oviduct in laying phase, which could offer a better understanding about the pathogenicity of IBV.

Pathogenicity, colonization, and innate immune response to Pasteurella multocida in rabbits.

BACKGROUND: Pasteurella multocida (P. multocida) infection can cause a series of diseases in different animals and cause huge economic losses to the breeding industry. P. multocida is considered to be one of the most significant pathogens in rabbits. In order to elucidate the pathogenic mechanism and innate immune response of P. multocida, an infection experiment was carried out in this study. RESULTS: Our results showed that the clinical symptoms of rabbits were severe dyspnoea and serous nasal fluid. During the course of the disease, the deaths peaked at 2 days post infection (dpi) and mortality rate was 60%. The pathological changes of the lung, trachea, and thymus were observed. In particular, consolidation and abscesses appeared in lung. Histopathologic changes in rabbits showed edema, hemorrhage, and neutrophil infiltration in the lung. P. multocida can rapidly replicate in a variety of tissues, and the colonization in most of the tested tissues reached the maximum at 2 dpi and then decreased at 3 dpi. The number of P. multocida in lung and thymus remained high level at 3 dpi. Toll-like receptors 2 and 4 signaling pathways were activated after P. multocida infection. The expression of Il1ß, Il6, Il8, and Tnf-α was significantly increased. The expression of most proinflammatory cytokines peaked at 2 dpi and decreased at 3 dpi, and the expression trend of cytokines was consistent with the colonization of P. multocida in rabbit tissues. CONCLUSIONS: The P. multocida can rapidly replicate in various tissues of rabbit and cause bacteremia after infection. TLRs signaling pathways were activated after P. multocida infection, significantly inducing the expression of proinflammatory cytokines, which is might the main cause of respiratory inflammation and septicemia.

Antibiotic resistance of Mycoplasma Synoviae strains isolated in China from 2016 to 2019.

BACKGROUND: In the past decade, Mycoplasma synoviae (M. synoviae) infection has become widely prevalent in China, has caused serious economic losses and has become one of the most important diseases in the chicken industry. Medication is a general approach for the control of M. synoviae infection, but antibiotics are sometimes ineffective in clinical practice. To investigate the sensitivity of M. synoviae to antimicrobials commonly used in the treatment of M. synoviae infection, the antibiotic susceptibility of 32 M. synoviae strains isolated from China from 2016 to 2019 were determined using the minimum inhibitory concentration (MIC) method. RESULTS: All isolates had low MIC values for the combination of lincomycin and spectinomycin, pleuromutilin, and macrolides. However, the M. synoviae isolates displayed variance in MICs for doxycycline hydrochloride with a range of 0.25 to 8 µg/mL, and oxytetracycline hydrochloride with a range of 0.5 to 8 µg/mL. Three and one M. synoviae isolates showed intermediate MIC values to doxycycline hydrochloride and oxytetracycline hydrochloride, respectively. High MIC values for enrofloxacin were detected in all isolates with MICs ranging from 4 to 32 µg/mL. Furthermore, comparison of the parC QRDR identified a mutation at nucleotide position 254 (C254T) resulting in a Thr 85 Ile amino acid change in all M. synoviae isolates and the reference strain ATCC 25204 being resistant to enrofloxacin. Moreover, mutations at Glu 804 Gly and Thr 686 Ala of gyrA QRDR were identified in all M. synoviae isolates and ATCC 25204. The mutation in the QRDR of the parE gene resulted in amino acid changes at positions 197 (Pro to Ser) in 27/32 M. synoviae isolates. CONCLUSION: Three nonsynonymous mutations in gyrA and parE were first identified to be related to enrofloxacin resistance. Our results showed that M. synoviae resistance to enrofloxacin is widespread.

Echinococcus shiquicus in Qinghai-Tibet plateau: population structure and confirmation of additional endemic areas.

Echinococcus shiquicus is currently limited to the Qinghai­Tibet plateau, a large mountainous region in China. Although the zoonotic potential remains unknown, progress is being made on the distribution and intermediate host range. In this study, we report E. shiquicus within Gansu and Qinghai provinces in regions located not only around the central areas but also the southeast edge of the plateau and describe their genetic relationship with previous isolates from the plateau. From 1879 plateau pikas examined, 2.39% (95% CI 1.79­3.18) were infected with E. shiquicus. The highest prevalence of 10.26% (4.06­23.58) was recorded in Makehe town, Qinghai province. Overall the prevalence was marginally higher in Qinghai (2.5%, CI 1.82­3.43) than in Gansu (2%, CI 1.02­3.89). The cox1 and nad1 genes demonstrated high and low haplotype and nucleotide diversities, respectively. The median-joining network constructed by the cox1­nad1 gene sequences demonstrated a star-like configuration with a median vector (unsampled haplotype) occupying the centre of the network. No peculiar distinction or common haplotype was observed in isolates originating from the different provinces. The presence of E. shiquicus in regions of the southeast and northeast edges of the Qinghai­Tibet plateau and high genetic variation warrants more investigation into the haplotype distribution and genetic polymorphism by exploring more informative DNA regions of the mitochondrial genome to provide epidemiologically useful insight into the population structure of E. shiquicus across the plateau and its axis.

Monte Carlo Renormalization Group for Classical Lattice Models with Quenched Disorder.

We extend to quenched-disordered systems the variational scheme for real-space renormalization group calculations that we recently introduced for homogeneous spin Hamiltonians. When disorder is present our approach gives access to the flow of the renormalized Hamiltonian distribution, from which one can compute the critical exponents if the correlations of the renormalized couplings retain finite range. Key to the variational approach is the bias potential found by minimizing a convex functional in statistical mechanics. This potential reduces dramatically the Monte Carlo relaxation time in large disordered systems. We demonstrate the method with applications to the two-dimensional dilute Ising model, the random transverse field quantum Ising chain, and the random field Ising in two- and three-dimensional lattices.

Evaluation of the reproductive system development and egg-laying performance of hens infected with TW I-type infectious bronchitis virus.

The prevalence of TW I-type infectious bronchitis virus (IBV) has been increasing rapidly, and it has become the second most common genotype of IBV in China threatening the poultry industry. In this study, 1-day-old specific-pathogen-free (SPF) chickens infected with TW I-type IBV were continuously observed for 200 days. TW I-type IBV affected the respiratory, urinary, and female reproductive systems, resulting in a mortality rate of 10% as well as a decrease in egg quantity and an increase in inferior eggs. During the monitoring period, serious lesions occurred in the female reproductive system, such as yolk peritonitis, a shortened oviduct, and cysts of different sizes with effusion in the degenerated right oviduct. The infective viruses persisted in vivo for a long time, and due to the stress of laying, virus shedding was detected again after the onset of egg production. Our findings suggest that TW I-type IBV is deadly to chickens and could cause permanent damage to the oviduct, resulting in the poor laying performance of female survivors and decreasing the breeding value and welfare of the infected flock.

Quantum momentum distribution and quantum entanglement in the deep tunneling regime.

In this paper, we consider the momentum operator of a quantum particle directed along the displacement of two of its neighbors. A modified open-path path integral molecular dynamics is presented to sample the distribution of this directional momentum distribution, where we derive and use a new estimator for this distribution. Variationally enhanced sampling is used to obtain this distribution for an example molecule, malonaldehyde, in the very low temperature regime where deep tunneling happens. We find no secondary feature in the directional momentum distribution and that its absence is due to quantum entanglement through a further study of the reduced density matrix.

A multiplex PCR assay for the simultaneous detection of Taenia hydatigena, T. multiceps, T. pisiformis, and Dipylidium caninum infections.

BACKGROUND: Taenia hydatigena, T. multiceps, T. pisiformis, and Dipylidium caninum are four common large and medium-sized tapeworms parasitizing the small intestine of dogs and other canids. These parasites cause serious impact on the health and development of livestock. However, there are, so far, no commercially available molecular diagnostic kits capable of simultaneously detecting all four parasites in dogs. The aim of the study was therefore to develop a multiplex PCR assay that will accurately detect all four cestode infections in one reaction. METHODS: Specific primers for a multiplex PCR were designed based on corresponding mitochondrial genome sequences, and its detection limit was assessed by serial dilutions of the genomic DNAs of tapeworms examined. Furthermore, field samples of dog feces were tested using the developed assay. RESULTS: A multiplex polymerase chain reaction (PCR) assay was developed based on mitochondrial DNA (mtDNA) that accurately and simultaneously identify four cestode species in one reaction using specific fragment sizes of 592, 385, 283, and 190 bp for T. hydatigena, T. multiceps, T. pisiformis, and D. caninum, respectively. The lowest DNA concentration detected was 1 ng for T. hydatigena, T. multiceps and T. pisiformis, and 0.1 ng for D. caninum in a 25 µl reaction system. This assay offers high potential for the rapid detection of these four tapeworms in host feces simultaneously. CONCLUSIONS: This study provides an efficient tool for the simultaneous detection of T. hydatigena, T. multiceps, T. pisiformis, and D. caninum. The assay will be potentially useful in epidemiological studies, diagnosis, and treatment of these four cestodes infections during prevention and control program.

Autophagy induced by avian reovirus enhances viral replication in chickens at the early stage of infection.

BACKGROUND: Avian reovirus (ARV) is an important pathogen that can cause serious disease in poultry. Though several in vitro studies revealed some molecular mechanisms that are responsible for ARV-induced autophagy, it is still largely unknown how ARV manipulates autophagy to promote its own propagation. RESULTS: In this study, we demonstrated that ARV infection triggered autophagy in chicken tissues, evident from the enhancement of LC3-I/-II conversion and the appearance of abundant autophagosomes. Moreover, viral replication and the expression of IL-1ß were coupled with the process of ARV-induced autophagy in the early stage of infection. Furthermore, regulation of autophagy affected the accumulation of LC3-II, the production of ARV and the expression of IL-1ß. CONCLUSIONS: Altogether, our data suggest that ARV induces autophagy, which benefits its replication and dissemination in chicken tissues at the early infection stage.

Transcriptome analysis of avian reovirus-mediated changes in gene expression of normal chicken fibroblast DF-1 cells.

BACKGROUND: Avian reovirus (ARV) is an important poultry pathogen that can cause immunosuppression. In this study, RNA-Seq technology was applied to investigate the transcriptome-wide changes of DF-1 cells upon ARV infection at the middle stage. RESULTS: Total RNA of ARV-infected or mock-infected samples at 10 and 18 h post infection (hpi) was extracted to build RNA-Seq datasets. Analysis of the sequencing data revealed that the expressions of numerous genes were altered, and a panel of differentially expressed genes were confirmed with RT-qPCR. At 10 hpi, 104 genes were down-regulated and 64 were up-regulated, while the expressions of 47 genes were increased and only one was down-regulated, which may play a role in retinoic acid biosynthesis, at 18 hpi in the ARV-infected cells. The similar profiles of up-regulated genes between the two groups of infected cells suggest that ARV infection activated a prolonged antiviral response of host cells. Alternative splicing analysis found no significantly changed events altered by ARV infection. CONCLUSIONS: Overall, the differential expression profile presented in this study can be used to expand our understanding of the comprehensive interactions between ARV and the host cells, and may be helpful for us to reveal the pathogenic mechanism on the molecular level.

Variational Approach to Monte Carlo Renormalization Group.

We present a Monte Carlo method for computing the renormalized coupling constants and the critical exponents within renormalization theory. The scheme, which derives from a variational principle, overcomes critical slowing down, by means of a bias potential that renders the coarse grained variables uncorrelated. The two-dimensional Ising model is used to illustrate the method.

The papain-like protease of avian infectious bronchitis virus has deubiquitinating activity.

Coronavirus papain-like proteases (PLPs) can act as proteases that process virus-encoded large replicase polyproteins and also as deubiquitinating (DUB) enzymes. Like the PLPs of other coronaviruses (CoVs), the avian infectious bronchitis virus (IBV) PLP catalyzes proteolysis of Gly-Gly dipeptide bonds to release mature cleavage products. However, the other functions of the IBV PLP are not well understood. In this study, we found that IBV exhibits strong global DUB activity with significant reductions of the levels of ubiquitin (Ub)-, K48-, and K63-conjugated proteins. The DUB activity exhibited a clear time dependence, with stronger DUB activity in the early stage of viral infection. Furthermore, the IBV replicase-encoded PLP, including the downstream transmembrane (TM) domain, is a DUB enzyme and dramatically reduced the level of Ub-conjugated proteins, while processing both K48- and K63-linked polyubiquitin chains. By contrast, PLP did not cause any reduction of haemagglutinin (HA)-Ub-conjugated proteins. In addition, mutations of the catalytic residues of PLP-TM, Cys1274Ser and His1437Lys, reduced DUB activity against Ub-, K48- and K63- conjugated proteins, indicating that the DUB activity of the PLP-TM wild-type protein is not completely dependent on its catalytic activity. Overall, these results demonstrate that the IBV-encoded PLP-TM functions as a DUB enzyme and suggest that IBV may interfere with the activation of host antiviral signaling pathway by degrading polyubiquitin-associated proteins.

Avian infectious bronchitis virus disrupts the melanoma differentiation associated gene 5 (MDA5) signaling pathway by cleavage of the adaptor protein MAVS.

BACKGROUND: Melanoma differentiation associated gene 5 (MDA5) and retinoic acid-inducible gene-I (RIG-I) selectively sense cytoplasmic viral RNA to induce an antiviral immune response. Infectious bronchitis virus (IBV) is one of the most important infectious agents in chickens, and in chicken cells, it can be recognized by MDA5 to activate interferon production. RIG-I is considered to be absent in chickens. However, the absence of RIG-I in chickens raises the question of whether this protein influences the antiviral immune response against IBV infection. RESULTS: Here, we showed that chicken cells transfected with domestic goose RIG-I (dgRIG-I) exhibited increased IFN-ß activity after IBV infection. We also found that IBV can cleave MAVS, an adaptor protein downstream of RIG-I and MDA5 that acts as a platform for antiviral innate immunity at an early stage of infection. CONCLUSIONS: Although chicken MDA5 (chMDA5) is functionally active during IBV infection, the absence of RIG-I may increase the susceptibility of chickens to IBV infection, and IBV may disrupt the activation of the host antiviral response through the cleavage of MAVS.

The Distribution of 18 Enterotoxin and Enterotoxin-Like Genes in Staphylococcus aureus Strains from Different Sources in East China.

The distribution of 18 staphylococcal enterotoxin (SE) or SE-like (SEl) genes in Staphylococcus aureus strains from different sources in east China was investigated. Among all 496 S. aureus strains, 291 strains carried one or more SE genes. The more frequently occurred genes were sea, seb, seg, selk, sell, selm, selo, and seq; the less frequent occurred genes were sec, selj, and ser. The classic SE genes and the enterotoxin gene cluster (egc) (seg, sei, selm, seln, selo, and/or selu) accounted for 25.67% and 61.68% of all detected genes, respectively. There were three gene clusters (egc, sea-sek-seq, and sed-sej-ser), of which the egc cluster was the important one that could generate novel complexes, and the sea-sek-seq cluster was a close relative to the hospital-acquired methicillin-resistant S. aureus. The SE gene distributions were different among strains of different sources and formed diverse toxin gene profiles. The human- and foodborne-origin strains harbored classic and novel SE and SEl genes, whereas animal-origin strains harbored egc and other novel SE and SEl genes mainly. The foodborne- and human-origin strains were the main dangerous factors of classic staphylococcal foodborne poisoning, whereas the strains (especially from animals) that carried egc and other novel genes mainly should be new potential dangerous factors for food safety.