bsgenova: an accurate, robust, and fast genotype caller for bisulfite-sequencing data.

BACKGROUND: Bisulfite sequencing (BS-Seq) is a fundamental technique for characterizing DNA methylation profiles. Genotype calling from bisulfite-converted BS-Seq data allows allele-specific methylation analysis and the concurrent exploration of genetic and epigenetic profiles. Despite various methods have been proposed, single nucleotide polymorphisms (SNPs) calling from BS-Seq data, particularly for SNPs on chromosome X and in the presence of contaminative data, poses ongoing challenges. RESULTS: We introduce bsgenova, a novel SNP caller tailored for bisulfite sequencing data, employing a Bayesian multinomial model. The performance of bsgenova is assessed by comparing SNPs called from real-world BS-Seq data with those from corresponding whole-genome sequencing (WGS) data across three human cell lines. bsgenova is both sensitive and precise, especially for chromosome X, compared with three existing methods. Moreover, in the presence of low-quality reads, bsgenova outperforms other methods notably. In addition, bsgenova is meticulously implemented, leveraging matrix imputation and multi-process parallelization. Compared to existing methods, bsgenova stands out for its speed and efficiency in memory and disk usage. Furthermore, bsgenova integrates bsextractor, a methylation extractor, enhancing its flexibility and expanding its utility. CONCLUSIONS: We introduce bsgenova for SNP calling from bisulfite-sequencing data. The source code is available at https://github.com/hippo-yf/bsgenova under license GPL-3.0.

Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Comparative Genomics and Characterization of Shigella flexneri Isolated from Urban Wastewater.

Shigella species are a group of highly transmissible Gram-negative pathogens. Increasing reports of infection with extensively drug-resistant varieties of this stomach bug has convinced the World Health Organization to prioritize Shigella for novel therapeutic interventions. We herein coupled the whole-genome sequencing of a natural isolate of Shigella flexneri with a pangenome ana-lysis to characterize pathogen genomics within this species, which will provide us with an insight into its existing genomic diversity and highlight the root causes behind the emergence of quick vaccine escape variants. The isolated novel strain of S. flexneri contained ~4,500 protein-coding genes, 57 of which imparted resistance to antibiotics. A comparative pan-genomic ana-lysis revealed genomic variability of ~64%, the shared conservation of core genes in central metabolic processes, and the enrichment of unique/accessory genes in virulence and defense mechanisms that contributed to much of the observed antimicrobial resistance (AMR). A pathway ana-lysis of the core genome mapped 22 genes to 2 antimicrobial resistance pathways, with the bulk coding for multidrug efflux pumps and two component regulatory systems that are considered to work synergistically towards the development of resistance phenotypes. The prospective evolvability of Shigella species as witnessed by the marked difference in genomic content, the strain-specific essentiality of unique/accessory genes, and the inclusion of a potent resistance mechanism within the core genome, strengthens the possibility of novel serotypes emerging in the near future and emphasizes the importance of tracking down genomic diversity in drug/vaccine design and AMR governance.

European Autism GEnomics Registry (EAGER): protocol for a multicentre cohort study and registry.

INTRODUCTION: Autism is a common neurodevelopmental condition with a complex genetic aetiology that includes contributions from monogenic and polygenic factors. Many autistic people have unmet healthcare needs that could be served by genomics-informed research and clinical trials. The primary aim of the European Autism GEnomics Registry (EAGER) is to establish a registry of participants with a diagnosis of autism or an associated rare genetic condition who have undergone whole-genome sequencing. The registry can facilitate recruitment for future clinical trials and research studies, based on genetic, clinical and phenotypic profiles, as well as participant preferences. The secondary aim of EAGER is to investigate the association between mental and physical health characteristics and participants' genetic profiles. METHODS AND ANALYSIS: EAGER is a European multisite cohort study and registry and is part of the AIMS-2-TRIALS consortium. EAGER was developed with input from the AIMS-2-TRIALS Autism Representatives and representatives from the rare genetic conditions community. 1500 participants with a diagnosis of autism or an associated rare genetic condition will be recruited at 13 sites across 8 countries. Participants will be given a blood or saliva sample for whole-genome sequencing and answer a series of online questionnaires. Participants may also consent to the study to access pre-existing clinical data. Participants will be added to the EAGER registry and data will be shared externally through established AIMS-2-TRIALS mechanisms. ETHICS AND DISSEMINATION: To date, EAGER has received full ethical approval for 11 out of the 13 sites in the UK (REC 23/SC/0022), Germany (S-375/2023), Portugal (CE-085/2023), Spain (HCB/2023/0038, PIC-164-22), Sweden (Dnr 2023-06737-01), Ireland (230907) and Italy (CET_62/2023, CEL-IRCCS OASI/24-01-2024/EM01, EM 2024-13/1032 EAGER). Findings will be disseminated via scientific publications and conferences but also beyond to participants and the wider community (eg, the AIMS-2-TRIALS website, stakeholder meetings, newsletters).

Whole genome sequencing analysis on antibiotic-resistant Escherichia coli isolated from pig farms in Banten Province, Indonesia.

IMPORTANCE: The emergence and rapid increase in the incidence of multidrug-resistant (MDR) bacteria in pig farms has become a serious concern and reduced the choice of effective antibiotics. OBJECTIVE: This study analyzed the phylogenetics and diversity of antibiotic resistance genes (ARGs) and molecularly identified the source of ARGs in antibiotic-resistant Escherichia coli isolated from pig farms in Banten Province, Indonesia. METHODS: Forty-four antibiotic-resistant E. coli isolates from fecal samples from 44 pig farms in Banten Province, Indonesia, were used as samples. The samples were categorized into 14 clusters. Sequencing was performed using the Oxford Nanopore Technologies MinION platform, with barcoding before sequencing with Nanopore Rapid sequencing gDNA-barcoding (SQK-RBK110.96) according to manufacturing procedures. ARG detection was conducted using ResFinder, and the plasmid replicon was determined using PlasmidFinder. RESULTS: Three phylogenetic leaves of E. coli were identified in the pig farming cluster in Banten Province. The E. coli isolates exhibited potential resistance to nine classes of antibiotics. Fifty-one ARGs were identified across all isolates, with each cluster carrying a minimum of 10 ARGs. The ant(3'')-Ia and qnrS1 genes were present in all isolates. ARGs in the E. coli pig farming cluster originated mainly from plasmids, accounting for an average of 89.4%. CONCLUSIONS AND RELEVANCE: The elevated potential for MDR events, coupled with the dominance of ARGs originating from plasmids, increases the risk of ARG spread among bacterial populations in animals, humans, and the environment.

Morphological, Metabolomic and Genomic Evidences on Drought Stress Protective Functioning of the Endophyte Bacillus safensis Ni7.

The metabolomic and genomic characterization of an endophytic Bacillus safensis Ni7 was carried out in this study. This strain has previously been isolated from the xerophytic plant Nerium indicum L. and reported to enhance the drought tolerance in Capsicum annuum L. seedlings. The effects of drought stress on the morphology, biofilm production, and metabolite production of B. safensis Ni7 are analyzed in the current study. From the results obtained, the organism was found to have multiple strategies such as aggregation and clumping, robust biofilm production, and increased production of surfactin homologues under the drought induced condition when compared to non-stressed condition. Further the whole genome sequencing (WGS) based analysis has demonstrated B. safensis Ni7 to have a genome size of 3,671,999 bp, N50 value of 3,527,239, and a mean G+C content of 41.58%. Interestingly the organism was observed to have the presence of various stress-responsive genes (13, 20U, 16U,160, 39, 17M, 18, 26, and ctc) and genes responsible for surfactin production (srfAA, srfAB, srfAC, and srfAD), biofilm production (epsD, epsE, epsF, epsG, epsH, epsI, epsK, epsL, epsM, epsN, and pel), chemotaxis (cheB_1, cheB_2, cheB_3, cheW_1, cheW_2 cheR, cheD, cheC, cheA, cheY, cheV, and cheB_4), flagella synthesis (flgG_1, flgG_2, flgG_3, flgC, and flgB) as supportive to the drought tolerance. Besides these, the genes responsible for plant growth promotion (PGP), including the genes for nitrogen (nasA, nasB, nasC, nasD, and nasE) and sulfur assimilation (cysL_1&L_2, cysI) and genes for phosphate solubilization (phoA, phoP_1& phoP_2, and phoR) could also be predicted. Along with the same, the genes for catalase, superoxide dismutase, protein homeostasis, cellular fitness, osmoprotectants production, and protein folding could also be predicted from its WGS data. Further pan-genome analysis with plant associated B. safensis strains available in the public databases revealed B. safensis Ni7 to have the presence of a total of 5391 gene clusters. Among these, 3207 genes were identified as core genes, 954 as shell genes and 1230 as cloud genes. This variation in gene content could be taken as an indication of evolution of strains of Bacillus safensis as per specific conditions and hence in the case of B. safensis Ni7 its role in habitat adaptation of plant is well expected. This diversity in endophytic bacterial genes may attribute its role to support the plant system to cope up with stress conditions. Overall, the study provides genomic evidence on Bacillus safensis Ni7 as a stress alleviating microbial partner in plants.

Whole-genome sequencing reveals genomic diversity and selection signatures in Xia'nan cattle.

BACKGROUND: The crossbreeding of specialized beef cattle breeds with Chinese indigenous cattle is a common method of genetic improvement. Xia'nan cattle, a crossbreed of Charolais and Nanyang cattle, is China's first specialized beef cattle breed with independent intellectual property rights. After more than two decades of selective breeding, Xia'nan cattle exhibit a robust physique, good environmental adaptability, good tolerance to coarse feed, and high meat production rates. This study analyzed the population genetic structure, genetic diversity, and genomic variations of Xia'nan cattle using whole-genome sequencing data from 30 Xia'nan cattle and 178 published cattle genomic data. RESULT: The ancestry estimating composition analysis showed that the ancestry proportions for Xia'nan cattle were mainly Charolais with a small amount of Nanyang cattle. Through the genetic diversity studies (nucleotide diversity and linkage disequilibrium decay), we found that the genomic diversity of Xia'nan cattle is higher than that of specialized beef cattle breeds in Europe but lower than that of Chinese native cattle. Then, we used four methods to detect genome candidate regions influencing the excellent traits of Xia'nan cattle. Among the detected results, 42 genes (θπ and CLR) and 131 genes (FST and XP-EHH) were detected by two different detection strategies. In addition, we found a region in BTA8 with strong selection signals. Finally, we conducted functional annotation on the detected genes and found that these genes may influence body development (NR6A1), meat quality traits (MCCC1), growth traits (WSCD1, TMEM68, MFN1, NCKAP5), and immunity (IL11RA, CNTFR, CCL27, SLAMF1, SLAMF7, NAA35, and GOLM1). CONCLUSION: We elucidated the genomic features and population structure of Xia'nan cattle and detected some selection signals in genomic regions potentially associated with crucial economic traits in Xia'nan cattle. This research provided a basis for further breeding improvements in Xia'nan cattle and served as a reference for genetic enhancements in other crossbreed cattle.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes.

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multilocus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Whole-genome analysis of human group A rotaviruses in 1980s Japan and evolutionary assessment of global Wa-like strains across half a century.

Historically, the Wa-like strains of human group A rotavirus (RVA) have been major causes of gastroenteritis. However, since the 2010s, the circulation of non-Wa-like strains has been increasingly reported, indicating a shift in the molecular epidemiology of RVA. Although understanding RVA evolution requires the analysis of both current and historical strains, comprehensive pre-1980's sequencing data are scarce globally. We determined the whole-genome sequences of representative strains from six RVA gastroenteritis outbreaks observed at an infant home in Sapporo, Japan, between 1981 and 1989. These outbreaks were mainly caused by G1 or G3 Wa-like strains, resembling strains from the United States in the 1970s-1980s and from Malawi in the 1990s. Phylogenetic analysis of these infant home strains, together with Wa-like strains collected worldwide from the 1970s to 2020, revealed a notable trend: pre-2010 strains diverged into multiple lineages in many genomic segments, whereas post-2010 strains tended to converge into a single lineage. However, Bayesian skyline plot indicated near-constant effective population sizes from the 1970s to 2020, and selection pressure analysis identified positive selection only at amino acid 75 of NSP2. These results suggest that evidence supporting the influence of rotavirus vaccines, introduced globally since 2006, on Wa-like RVA molecular evolution is lacking at present, and phylogenetic analysis may simply reflect natural fluctuations in RVA molecular evolution. Evaluating the long-term impact of RV vaccines on the molecular evolution of RVA requires sustained surveillance.

Exploring the potential of Oxford Nanopore Technologies sequencing for Mycobacterium tuberculosis sequencing: An assessment of R10 flowcells and V14 chemistry.

Oxford Nanopore Technologies (ONT) sequencing is a promising technology. We assessed the performance of the new ONT R10 flowcells and V14 rapid sequencing chemistry for Mtb whole genome sequencing of Mycobacterium tuberculosis (Mtb) DNA extracted from clinical primary liquid cultures (CPLCs). Using the recommended protocols for MinION Mk1C, R10.4.1 MinION flowcells, and the ONT Rapid Sequencing Kit V14 on six CPLC samples, we obtained a pooled library yield of 10.9 ng/µl, generated 1.94 Gb of sequenced bases and 214k reads after 48h in a first sequencing run. Only half (49%) of all generated reads met the Phred Quality score threshold (>8). To assess if the low data output and sequence quality were due to impurities present in DNA extracted directly from CPLCs, we added a pre-library preparation bead-clean-up step and included purified DNA obtained from an Mtb subculture as a control sample in a second sequencing run. The library yield for DNA extracted from four CPLCs and one Mtb subculture (control) was similar (10.0 ng/µl), 2.38 Gb of bases and 822k reads were produced. The quality was slightly better with 66% of the produced reads having a Phred Quality >8. A third run of DNA from six CPLCs with bead clean-up pre-processing produced a low library yield (±1 Gb of bases, 166k reads) of low quality (51% of reads with a Phred Quality score >8). A median depth of coverage above 10× was only achieved for five of 17 (29%) sequenced libraries. Compared to Illumina WGS of the same samples, accurate lineage predictions and full drug resistance profiles from the generated ONT data could not be determined by TBProfiler. Further optimization of the V14 ONT rapid sequencing chemistry and library preparation protocol is needed for clinical Mtb WGS applications.

Klenkia sesuvii sp. nov., isolated from leaves of halophyte Sesuvium portulacastrum.

During a study on the diversity of culturable actinobacteria from coastal halophytes in Thailand, strain LSe6-5T was isolated from leaves of sea purslane (Sesuvium portulacastrum L.), and a polyphasic approach was employed to determine its taxonomic position. The 16S rRNA gene sequences analysis indicated that the strain was most closely related to Klenkia brasiliensis Tu 6233T (99.2 %), Klenkia marina YIM M13156T (99.1 %), and Klenkia terrae PB261T (98.7 %). The genome of strain LSe6-5T was estimated to be 4.33 Mbp in size, with DNA G+C contents of 74.3%. A phylogenomic tree based on whole-genome sequences revealed that strain LSe6-5T formed a clade with Klenkia marina DSM 45722T, indicating their close relationship. However, the average nucleotide identity (ANI)-blast, ANI-MUMmer, and dDDH values between strain LSe6-5T with K. marina DSM 45722T (87.1, 88.9, and 33.0 %) were below the thresholds of 95-96 % ANI and 70 % dDDH for identifying a novel species. Furthermore, strain LSe6-5T showed morphological and chemotaxonomic characteristics of the genus Klenkia. Cells were motile, rod-shaped, and Gram-stain-positive. Optimal growth of strain LSe6-5T occurred at 28 °C, pH 7.0, and 0-3 % NaCl. The whole-cell hydrolysates contained meso-diaminopimelic acid as the diagnostic diamino acid, with galactose, glucose, mannose, and ribose as whole-cell sugars. The predominant menaquinones were MK-9(H4) and MK-9(H0). The polar lipid profile was composed of diphosphatidylglycerol, hydroxyphosphatidylethanolamine, phosphatidylinositol, glycophosphatidylinositol, an unidentified phospholipid, and an unidentified lipid. Major cellular fatty acids were iso-C15 : 0, iso-C16 : 0, and iso-C17 : 0. From the distinct phylogenetic position and combination of genotypic and phenotypic characteristics, it is supported that strain LSe6-5T represents a novel species of the genus Klenkia, for which the name Klenkia sesuvii sp. nov. is proposed. The type strain is strain LSe6-5T (=TBRC 16417T= NBRC 115929T).

Prevalence of ST1049-KL5 carbapenem-resistant Klebsiella pneumoniae with a blaKPC-2 and blaNDM-1 co-carrying hypertransmissible IncM1 plasmid.

Infection caused by KPC and NDM carbapenemases co-producing Klebsiella pneumoniae (KPC_NDM_CRKP) poses serious public health concerns. Here, we elucidate the prevalence of a hypertransmissible lncM1 plasmid, pKPC_NDM, co-carrying blaKPC-2 and blaNDM-1 genes in sequence type 1049 K_locus 5 (ST1049-KL5) KPC_NDM_CRKP isolates. Genetic and clonal relatedness analyses using pulsed-field gel electrophoresis, single nucleotide polymorphism analysis and core genome multilocus sequence typing suggested clonal dissemination of ST1049-KL5 KPC_NDM_CRKP strains in our hospital. Whole genome sequencing identified an identical 76,517 bp- blaKPC-2 and blaNDM-1 genes co-carrying IncM1 plasmid pKPC_NDM and a pLVPK-like hypervirulent plasmid in all ST1049-KL5 KPC_NDM_CRKP isolates. pKPC_NDM shared 100% identity with a previously sequenced plasmid CRKP35_unnamed4, demonstrating high transferability in conjugation assay, with conjugation frequencies reaching 10-4 and 10-5 in Escherichia coli and K. pneumoniae recipients, respectively. It also maintained favorable stability and flexible compatibility, with retention rates exceeding 80% after 10 days of continuous passage, and could be compatible with pre-existing blaKPC- or blaNDM-carrying plasmids in recipient strains. This study summarizes the characteristics of KPC_NDM_CRKP outbreaks and highlights the importance of ongoing surveillance and infection control strategies to address the challenges posed by ST1049 K. pneumoniae strains.

Genome-wide association study of BNT162b2 vaccine-related myocarditis identifies potential predisposing functional areas in Hong Kong adolescents.

Vaccine-related myocarditis associated with the BNT162b2 vaccine is a rare complication, with a higher risk observed in male adolescents. However, the contribution of genetic factors to this condition remains uncertain. In this study, we conducted a comprehensive genetic association analysis in a cohort of 43 Hong Kong Chinese adolescents who were diagnosed with myocarditis shortly after receiving the BNT162b2 mRNA COVID-19 vaccine. A comparison of whole-genome sequencing data was performed between the confirmed myocarditis cases and a control group of 481 healthy individuals. To narrow down potential genomic regions of interest, we employed a novel clustering approach called ClusterAnalyzer, which prioritised 2,182 genomic regions overlapping with 1,499 genes for further investigation. Our pathway analysis revealed significant enrichment of these genes in functions related to cardiac conduction, ion channel activity, plasma membrane adhesion, and axonogenesis. These findings suggest a potential genetic predisposition in these specific functional areas that may contribute to the observed side effect of the vaccine. Nevertheless, further validation through larger-scale studies is imperative to confirm these findings. Given the increasing prominence of mRNA vaccines as a promising strategy for disease prevention and treatment, understanding the genetic factors associated with vaccine-related myocarditis assumes paramount importance. Our study provides valuable insights that significantly advance our understanding in this regard and serve as a valuable foundation for future research endeavours in this field.

Genome sequencing of captive white tigers from Bangladesh.

OBJECTIVES: The Bengal tiger Panthera tigris tigris, is an emblematic animal for Bangladesh. Despite being the apex predator in the wild, their number is decreasing due to anthropogenic activities such as hunting, urbanization, expansion of agriculture and deforestation. By contrast, captive tigers are flourishing due to practical conservation efforts. Breeding within the small captive population can produce inbreeding depression and genetic bottlenecks, which may limit the success of conservation efforts. Despite past decades of research, a comprehensive database on genetic variation in the captive and wild Bengal tigers in Bangladesh still needs to be included. Therefore, this research aimed to investigate the White Bengal tiger genome to create a resource for future studies to understand variation underlying important functional traits. DATA DESCRIPTION: Blood samples from Chattogram Zoo were collected for three white Bengal tigers. Genomic DNA for all collected samples were extracted using a commercial DNA extraction kit. Whole genome sequencing was performed using a DNBseq platform. We generated 77 Gb of whole-genome sequencing (WGS) data for three white Bengal tigers (Average 11X coverage/sample). The data we generated will establish a paradigm for tiger research in Bangladesh by providing a genomic resource for future functional studies on the Bengal white tiger.

Whole genome sequences of 70 indigenous Ethiopian cattle.

Indigenous animal genetic resources play a crucial role in preserving global genetic diversity and supporting the livelihoods of millions of people. In Ethiopia, the majority of the cattle population consists of indigenous breeds. Understanding the genetic architecture of these cattle breeds is essential for effective management and conservation efforts. In this study, we sequenced DNA samples from 70 animals from seven indigenous cattle breeds, generating about two terabytes of pair-end reads with an average coverage of 14X. The sequencing data were pre-processed and mapped to the cattle reference genome (ARS-UCD1.2) with an alignment rate of 99.2%. Finally, the variant calling process produced approximately 35 million high-quality SNPs. These data provide a deeper understanding of the genetic landscape, facilitate the identification of causal mutations, and enable the exploration of evolutionary patterns to assist cattle improvement and sustainable utilization, particularly in the face of unpredictable climate changes.

Emergence of Eurasian Avian-Like Swine Influenza A (H1N1) virus in a child in Shandong Province, China.

BACKGROUND: Influenza A virus infections can occur in multiple species. Eurasian avian-like swine influenza A (H1N1) viruses (EAS-H1N1) are predominant in swine and occasionally infect humans. A Eurasian avian-like swine influenza A (H1N1) virus was isolated from a boy who was suffering from fever; this strain was designated A/Shandong-binzhou/01/2021 (H1N1). The aims of this study were to investigate the characteristics of this virus and to draw attention to the need for surveillance of influenza virus infection in swine and humans. METHODS: Throat-swab specimens were collected and subjected to real-time fluorescent quantitative polymerase chain reaction (RT‒PCR). Positive clinical specimens were inoculated onto Madin-Darby canine kidney (MDCK) cells to isolate the virus, which was confirmed by a haemagglutination assay. Then, whole-genome sequencing was carried out using an Illumina MiSeq platform, and phylogenetic analysis was performed with MEGA X software. RESULTS: RT‒PCR revealed that the throat-swab specimens were positive for EAS-H1N1, and the virus was subsequently successfully isolated from MDCK cells; this strain was named A/Shandong-binzhou/01/2021 (H1N1). Whole-genome sequencing and phylogenetic analysis revealed that A/Shandong-binzhou/01/2021 (H1N1) is a novel triple-reassortant EAS-H1N1 lineage that contains gene segments from EAS-H1N1 (HA and NA), triple-reassortant swine influenza H1N2 virus (NS) and A(H1N1) pdm09 viruses (PB2, PB1, PA, NP and MP). CONCLUSIONS: The isolation and analysis of the A/Shandong-binzhou/01/2021 (H1N1) virus provide further evidence that EAS-H1N1 poses a threat to human health, and greater attention should be given to the surveillance of influenza virus infections in swine and humans.

Whole genome sequencing of HER2-positive metastatic extramammary Paget's disease: a case report.

BACKGROUND: Extramammary Paget's disease (EMPD) is a rare cancer that occurs within the epithelium of the skin, arising predominantly in areas with high apocrine gland concentration such as the vulva, scrotum, penis and perianal regions. Here, we aim to integrate clinicopathological data with genomic analysis of aggressive, rapidly-progressing de novo metastatic EMPD responding to HER2-directed treatment in combination with other agents, to attain a more comprehensive understanding of the disease landscape. METHODS: Immunohistochemical staining on the scrotal wall tumor and bone marrow metastasis demonstrated HER2 overexpression. Whole genome sequencing of the tumor and matched blood was performed. RESULTS: Notable copy number gains (log2FC > 0.9) on chromosomes 7 and 8 were detected (n = 81), with 92.6% of these unique genes specifically located on chromosome 8. Prominent cancer-associated genes include ZNF703, HOOK3, DDHD2, LSM1, NSD3, ADAM9, BRF2, KAT6A and FGFR1. Interestingly, ERBB2 gene did not exhibit high copy number gain (log2FC = 0.4) although 90% of tumor cells stained HER2-positive. Enrichment in pathways associated with transforming growth factor-beta (TGFß) (FDR = 0.0376, Enrichment Ratio = 8.12) and fibroblast growth factor receptor (FGFR1) signaling (FDR = 0.0082, Enrichment Ratio = 2.3) was detected. Amplicon structure analysis revealed that this was a simple-linear amplification event. CONCLUSION: Whole genome sequencing revealed the underlying copy number variation landscape in HER2-positive metastatic EMPD. The presence of alternative signalling pathways and genetic variants suggests potential interactions with HER2 signalling, which possibly contributed to the HER2 overexpression and observed response to HER2-directed therapy combined with other agents in a comprehensive treatment regimen.

Unraveling tuberculosis patient cluster transmission chains: integrating WGS-based network with clinical and epidemiological insights.

Background: Tuberculosis remains a global health threat, and the World Health Organization reports a limited reduction in disease incidence rates, including both new and relapse cases. Therefore, studies targeting tuberculosis transmission chains and recurrent episodes are crucial for developing the most effective control measures. Herein, multiple tuberculosis clusters were retrospectively investigated by integrating patients' epidemiological and clinical information with median-joining networks recreated based on whole genome sequencing (WGS) data of Mycobacterium tuberculosis isolates. Methods: Epidemiologically linked tuberculosis patient clusters were identified during the source case investigation for pediatric tuberculosis patients. Only M. tuberculosis isolate DNA samples with previously determined spoligotypes identical within clusters were subjected to WGS and further median-joining network recreation. Relevant clinical and epidemiological data were obtained from patient medical records. Results: We investigated 18 clusters comprising 100 active tuberculosis patients 29 of whom were children at the time of diagnosis; nine patients experienced recurrent episodes. M. tuberculosis isolates of studied clusters belonged to Lineages 2 (sub-lineage 2.2.1) and 4 (sub-lineages 4.3.3, 4.1.2.1, 4.8, and 4.2.1), while sub-lineage 4.3.3 (LAM) was the most abundant. Isolates of six clusters were drug-resistant. Within clusters, the maximum genetic distance between closely related isolates was only 5-11 single nucleotide variants (SNVs). Recreated median-joining networks, integrated with patients' diagnoses, specimen collection dates, sputum smear microscopy, and epidemiological investigation results indicated transmission directions within clusters and long periods of latent infection. It also facilitated the identification of potential infection sources for pediatric patients and recurrent active tuberculosis episodes refuting the reactivation possibility despite the small genetic distance of ≤5 SNVs between isolates. However, unidentified active tuberculosis cases within the cluster, the variable mycobacterial mutation rate in dormant and active states, and low M. tuberculosis genetic variability inferred precise transmission chain delineation. In some cases, heterozygous SNVs with an allelic frequency of 10-73% proved valuable in identifying direct transmission events. Conclusion: The complex approach of integrating tuberculosis cluster WGS-data-based median-joining networks with relevant epidemiological and clinical data proved valuable in delineating epidemiologically linked patient transmission chains and deciphering causes of recurrent tuberculosis episodes within clusters.

Antimicrobial resistance and genome characteristics of Salmonella enteritidis from Huzhou, China.

Salmonella enteritidis is a main pathogen responsible for sporadic outbreaks of gastroenteritis, and therefore is an important public health problem. This study investigated the drug resistance and genomic characteristics of S. enteritidis isolated from clinical and food sources in Huzhou, Zhejiang Province, China, from February 1, 2021, to December 30, 2023. In total, 43 S. enteritidis strains isolated during the study period were subjected to virulence gene, drug resistance gene, genetic correlation, antibiotic resistance, and multilocus sequence typing analyses. All 43 isolates were identified as ST11, and contained 108 virulence-related genes. Drug sensitivity analysis of the 43 isolates showed resistance rates of 100% to nalidixic acid and 90.70% to ampicillin and ampicillin/sulbactam. Multidrug resistance is a serious issue, with 81.40% of strains resistant to three or more antibacterial drugs. Genome sequencing indicated that S. enteritidis possessed 23 drug resistance genes, of which 14 were common to all 43 isolates. Phylogenetic analysis based on core genome single-nucleotide polymorphisms divided the 43 S. enteritidis strains into three clusters, with the 10 samples from an outbreak forming an independent branch located in cluster 3.

Genomic surveillance of SARS-CoV-2 evolution by a centralised pipeline and weekly focused sequencing, Austria, January 2021 to March 2023.

BackgroundThe COVID-19 pandemic was largely driven by genetic mutations of SARS-CoV-2, leading in some instances to enhanced infectiousness of the virus or its capacity to evade the host immune system. To closely monitor SARS-CoV-2 evolution and resulting variants at genomic-level, an innovative pipeline termed SARSeq was developed in Austria.AimWe discuss technical aspects of the SARSeq pipeline, describe its performance and present noteworthy results it enabled during the pandemic in Austria.MethodsThe SARSeq pipeline was set up as a collaboration between private and public clinical diagnostic laboratories, a public health agency, and an academic institution. Representative SARS-CoV-2 positive specimens from each of the nine Austrian provinces were obtained from SARS-CoV-2 testing laboratories and processed centrally in an academic setting for S-gene sequencing and analysis.ResultsSARS-CoV-2 sequences from up to 2,880 cases weekly resulted in 222,784 characterised case samples in January 2021-March 2023. Consequently, Austria delivered the fourth densest genomic surveillance worldwide in a very resource-efficient manner. While most SARS-CoV-2 variants during the study showed comparable kinetic behaviour in all of Austria, some, like Beta, had a more focused spread. This highlighted multifaceted aspects of local population-level acquired immunity. The nationwide surveillance system enabled reliable nowcasting. Measured early growth kinetics of variants were predictive of later incidence peaks.ConclusionWith low automation, labour, and cost requirements, SARSeq is adaptable to monitor other pathogens and advantageous even for resource-limited countries. This multiplexed genomic surveillance system has potential as a rapid response tool for future emerging threats.