Whole Genome Sequencing Research Articles

Co-occurrence of triple carbapenemase genes, bla VIM-2, bla NDM-1, and bla OXA-48 in Enterobacter hormaechei clinical isolates –first report from Croatia

Two Enterobacter hormaechei isolates harbouring three carbapenemase genes each, were isolated from two patients from different ICUs at University Hospital Centre Zagreb, Croatia, which is to our knowledge, the first report of triple carbapenemase (bla VIM-2, bla NDM-1, and bla OXA-48 ) co-existence in E. hormachei strains and also among Enterobacterales members in Croatia. Antimicrobial susceptibility testing showed susceptibility only to colistin and amikacin. The production of carbapenemases was phenotypically tested by immunochromatographic assay and confirmed by PCR. Detailed analysis by Whole Genome Sequencing (WGS) of short reads by Illumina and long reads by Oxford Nanopore Technologies (ONT) was additionally performed and showed that both isolates belonged to ST200. They were separated by 98 Single Nucleotide Polymorphisms (SNPs) having variations in the number of bla VIM-2 genes on the chromosome, the number of bla NDM-1 genes on the plasmid, non-identical bla NDM-1 plasmids, different plasmid content in general, and only one isolate carried a 94 kb prophage.

Unveiling novel Neocosmospora species from Thai mangroves as potent biocontrol agents against Colletotrichum species.

Neocosmospora species are saprobes, endophytes, and pathogens belonging to the family Nectriaceae. This study aims to investigate the taxonomy, biosynthetic potential and application of three newly isolated Neocosmospora species from mangrove habitats in the southern part of Thailand using phylogeny, bioactivity screening, genome sequencing and bioinformatics analysis. Detailed descriptions, illustrations, and a multi-locus phylogenetic tree with large subunit ribosomal DNA (LSU), internal transcribed spacer (ITS), translation elongation factor 1-alpha (ef1-α), and RNA polymerase II second largest subunit (RPB2) regions showing the placement of three fungal strains, MFLUCC 17-0253, MFLUCC 17-0257, and MFLUCC 17-0259 clustered within the Neocosmospora clade with strong statistical support. Fungal crude extracts of the new species Neocosmospora mangrovei MFLUCC 17-0253 exhibited strong antifungal activity to control Colletotrichum truncatum CG-0064, while N. ferruginea MFLUCC 17- 0259 exhibited only moderate antifungal activity toward C. acutatum CC-0036. Thus, N. mangrovei MFLUCC 17-0253 was sequenced by Oxford nanopore technology. The bioinformatics analysis revealed that 49.17 Mb genome of this fungus harbors 41 potential biosynthetic gene clusters. Two fungal isolates of Neocosmospora and a new species of N. mangrovei were reported in this study. These fungal strains showed activity against pathogenic fungi causing anthracnose in chili. In addition, full genome sequencing and bioinformatics analysis of N. mangrovei MFLUCC 17-0253 was obtained.

Lactic Acid Bacteria and Bacillus subtilis as Potential Protective Cultures for Biopreservation in the Food Industry

The use of bacteria and/or their compounds is an alternative to the use of positive-list additives that the food industry is using as a tool to meet consumer demands for more natural, long-shelf-life, and healthy products, in short, to offer clean label foods. The aim of this study is to investigate the suitability of cell-free supernatants (CFSs) from Qualified Presumption of Safety strains as bioprotective cultures. Out of an initial screening panel of about 200 isolates, strains Pediococcus acidilactici CNTA 1059, Lactiplantibacillus plantarum CNTA 600, Levilactobacillus brevis CNTA 1374, and Bacillus subtilis CNTA 517 demonstrated strong antimicrobial activity against, especially, Gram-positive bacteria. The CFSs of these four strains showed minimum inhibitory concentration values between 0.15% and 5% against Listeria monocytogenes and Lentilactobacillus parabuchneri. None of the four selected strains exhibited acquired resistance to target antibiotics, and the non-toxigenicity of all the CFSs was demonstrated. In the case of the three lactic acid bacteria, the presence of bacteriocin-like inhibitory substances was confirmed following the decline in antimicrobial activity due to treatment with proteases. Regarding B. subtilis, biosynthetic gene clusters for different bacteriocin-like substances, including protease-resistant lipoproteins, were found via whole-genome sequencing. In addition, all of the CFSs exhibited stable antimicrobial activity at a wide range of temperatures (70–121 °C) used for the pasteurization and sterilization of food products, with a loss of antimicrobial activity ranging from 3% to 28%. These results point to the possibility that CFSs from these strains could be used in the food industry as a biocontrol tool to develop new products.

A new Finnish flavor of feline coat coloration, "salmiak," is associated with a 95-kb deletion downstream of the KIT gene.

Cats with a distinctive white hair pattern of unknown molecular cause have been discovered in the Finnish domestic cat population. Based on the unique appearance of these cats, we have named this phenotype salmiak ("salty licorice"). The use of a commercially available panel test to genotype four salmiak-colored cats revealed the absence of all known variants associated with white-haired phenotypic loci: full White (W), Spotting (Ws) and the Birman white Gloves associated (wg) allele of the KIT proto-oncogene (KIT) gene. Whole-genome sequencing on two salmiak-colored cats was conducted to search for candidate causal variants in the KIT gene. Despite a lack of coding variants, visual inspection of the short read alignments revealed a large ~95 kb deletion located ~65 kb downstream of the KIT gene in the salmiak cats. Additional PCR genotyping of 180 domestic cats and three salmiak-colored cats confirmed the homozygous derived variant genotype fully concordant with the salmiak phenotype. We suggest the newly identified variant be designated as wsal for "w salmiak".

Detection of macropodid alphaherpesvirus 2 infection and lesions in sudden death of a captive Virginia opossum and a water opossum

Macropodid alphaherpesvirus 2 (MaAHV2) is best described in macropods and has been implicated in outbreaks among captive marsupial populations in Australia. Natural disease caused by herpesviruses has not been reported previously in opossum species, to our knowledge. One Virginia opossum ( Didelphis virginiana) and 1 water opossum ( Chironectes minimus) were submitted for postmortem examination from a zoo that housed 6 opossums, all of which died within several weeks. Red kangaroos ( Macropus rufus) and red-necked wallabies ( Macropus rufogriseus) were also present at the facility. Liver samples from both opossums were submitted for transmission electron microscopy and whole-genome sequencing. Microscopically, both opossums had multifocal necrosis in the liver and lung, with intranuclear inclusion bodies within hepatocytes and pneumocytes. Another significant finding in the Virginia opossum was sepsis, with isolation of Streptococcus didelphis from various organs. Ultrastructural analysis of formalin-fixed liver tissue identified herpesviral replication complexes in both opossums; negative-stain electron microscopy of unfixed liver tissue repeatedly yielded a negative result. The herpesvirus had >99% nucleotide identity with MaAHV2. These 2 cases indicate that both opossum species are susceptible to MaAHV2 infection, and the outbreak has implications for mixed-species facilities that house macropods.

Characterization of the Complete Mitochondrial Genome of Salassa thespis (Lepidoptera: Saturniidae) and Comparison with Other Lepidoptera Species

There are many species of Lepidoptera, but few complete mitochondrial genomes of Lepidoptera have been included in databases. Here, the complete mitochondrial genome sequence of Salassa thespis was isolated and characterized. It was 15,302 bp in length and contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. Among the 13 PCGs, the initiation codon of cytochrome c oxidase subunit 1 (cox1) was CGA, and the rest were ATN. The cox1 and cox2 genes had an incomplete stop codon T, while the rest terminated with TAA. Codon usage analysis showed that Phe, Ile, Leu and Asn were the most frequent amino acids, while Trp was the least. Like other Lepidopterans, some conserved motifs were found in the A + T-rich region, including a 17 bp poly-T guided by ATAGA, the AT-rich area and a poly-A element. Bayesian inference and maximum likelihood phylogenetic tree analysis based on 13 PCGs of S. thespis confirmed that it belonged to the Saturniidae family and showed the following relationship: (S. thespis + (Rhodinia fugax + Samia canningi)). The enrichment of mitochondrial DNA provides reference information for the study of the evolution and diversity of Lepidoptera insects.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing.

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Unveiling the emergence of SARS-CoV-2 JN.1 sub-variant: Insights from the first cases at Charles Nicolle Hospital, Tunisia.

The JN.1 sub-variant is a new variant of the SARS-CoV-2 Omicron strain, derived from the BA.2.86 sub-variant. It was first detected in late 2023 and has quickly spread to many countries, becoming the most prevalent variant in some regions. JN.1 exhibits a unique mutation (L455S) in the spike protein compared to the BA.2.86 lineage, which may affect its transmissibility and immune evasion capabilities. JN.1 has been designated as a "variant of interest" by the World Health Organization due to its rapidly increasing spread and is being closely monitored for its impact on the COVID-19 pandemic. This study describes the emergence of SARS-CoV-2 JN.1 sub-variant in Tunisia, and reports its mutation profiles.Nasopharyngeal samples collected over a four-month period (October 2023 to January 2024) were subjected to RNA extraction and real-time RT-PCR confirmation of SARS-CoV-2 infection. The whole-genome sequencing was performed by an iSeq 100 sequencer and COVIDSeq kit reagents (Illumina, USA).Mutation analysis, using the NextClade platform and GISAID database, revealed the presence of JN.1 in 15 out of 80 positive cases (18.75%) during the study period.The emergence of JN.1 highlights the ongoing evolution of SARS-CoV-2 and the need for continued surveillance and research to better understand the characteristics and impact of emerging variants.

Emergence of multidrug-resistant Staphylococcus haemolyticus in Neonatal Intensive Care Unit in Southern France, a Genomic Study

An emergence of multidrug-resistant (MDR) Staphylococcus haemolyticus has been observed in the neonatal intensive care unit (NICU) of Nîmes University Hospital in southern France. A case-control analysis was conducted on 96 neonates, to identify risk factors associated with S. haemolyticus infection, focusing on clinical outcomes. Forty-eight MDR S. haemolyticus strains, isolated from neonates between October 2019 and July 2022, were investigated using routine in vitro procedures and whole-genome sequencing. Additionally, five S. haemolyticus isolates from adult patients were sequenced to identify clusters circulating within the hospital environment. The incidence of neonatal S. haemolyticus was significantly associated with low birth weight, lower gestational age, and central catheter use (p<0.001). Sepsis was the most frequent clinical manifestation in this series (20/46, 43.5%) and was associated with five deaths. Based on whole-genome analysis, three S. haemolyticus genotypes were predicted: ST1 (6/53, 11%), ST25 (3/53, 5.7%), and ST29 (44/53, 83%), which included the subcluster II-A, predominantly emerging in the neonatal department. All strains were profiled in silico to be resistant to methicillin, erythromycin, aminoglycosides, and fluoroquinolones, consistent with in vitro antibiotic susceptibility tests. Moreover, in silico prediction of biofilm formation and virulence-encoding genes supported the association of ST29 with severe clinical outcomes, while the persistence in the NICU could be explained by the presence of antiseptic and heavy metal resistance-encoding genes. The clonality of S. haemolyticus ST29 subcluster II-A isolates confirms healthcare transmission causing severe infections. Based on these results, reinforced hygiene measures are necessary to eradicate the nosocomial transmission of MDR strains.

A fine-scale genetic map of the Japanese population.

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.

Retrospective Genotyping of Enteroviruses Using a Diagnostic Nanopore Sequencing Workflow

Enteroviruses are among the most common viruses pathogenic to humans. They are associated with various forms of disease, ranging from mild respiratory illness to severe neurological diseases. In recent years, an increasing number of isolated cases of children developing meningitis or encephalitis as a result of enterovirus infection have been reported, as well as discrete enterovirus D68 outbreaks in North America in 2014 and 2016. We developed an assay to rapidly genotype enteroviruses by sequencing a region within the VP1 gene using nanopore Flongles. We retrospectively analyzed enterovirus-/rhinovirus-positive clinical samples from the Zurich, Switzerland area mainly collected during two seasons in 2019/2020 and 2021/2022. Respiratory, cerebrospinal fluid, and stool samples were analyzed. Whole-genome sequencing was performed on samples with ambiguous genotyping results and enterovirus D68-positive samples. Out of 255 isolates, a total of 95 different genotypes were found. A difference in the prevalence of enterovirus and rhinovirus infections was observed for both sample type and age group. In particular, children aged 0–4 years showed a higher frequency of enterovirus infections. Comparing the respiratory seasons, a higher prevalence was found, especially for enterovirus A and rhinovirus A after the SARS-CoV-2 pandemic. The enterovirus genotyping workflow provides a rapid diagnostic tool for individual analysis and continuous enterovirus surveillance.

Identification and Safety Assessment of Enterococcus casseliflavus KB1733 Isolated from Traditional Japanese Pickle Based on Whole-Genome Sequencing Analysis and Preclinical Toxicity Studies

The present study involves the precise identification and safety evaluation of Enterococcus casseliflavus KB1733, previously identified using 16S rRNA analysis, through whole-genome sequencing, phenotypic analysis, and preclinical toxicity studies. Analyses based on the genome sequencing data confirm the identity of KB1733 as E. casseliflavus and show that the genes related to vancomycin resistance are only present on the chromosome, while no virulence factor genes are present on the chromosome or plasmid. Phenotypic analyses of antibiotic resistance and hemolytic activity also indicated no safety concerns. A bacterial reverse mutation test showed there was no increase in revertant colonies of heat-killed KB1733. An acute toxicity test employing heat-killed KB1733 at a dose of 2000 mg/kg body weight in rats resulted in no deaths and no weight gain or other abnormalities in the general condition of the animals, with renal depression foci and renal cysts only occurring at the same frequency as in the control. Taking the background data into consideration, the effects on the kidneys observed in the current study were not caused by KB1733. Our findings suggest that KB1733 is non-pathogenic to humans/animals, although further studies involving repeated oral toxicity tests and/or clinical tests are required.

Complete genome sequences of three Porphyromonas gingivalis strains, isolated from patient with esophageal squamous cell carcinoma and healthy individual in China.

Here, we report the complete genome sequences of three Porphyromonas gingivalis, one from patient with esophageal cancer (LyEC01), and the other two from periodontally healthy individuals (LyG-1 and LyG-2) in 2021 and 2023. The genome sizes of LyEC01, LyG-1, and LyG-2 were 2,408,275, 2,411,440, and 2,411,481 bp, respectively.

Identification of novel genes associated with atherosclerosis in Bama miniature pig.

Atherosclerosis is a chronic cardiovascular disease of great concern. However, it is difficult to establish a direct connection between conventional small animal models and clinical practice. The pig's genome, physiology, and anatomy reflect human biology better than other laboratory animals, which is crucial for studying the pathogenesis of atherosclerosis. We used whole-genome sequencing data from nine Bama minipigs to perform a genome-wide linkage analysis, and further used bioinformatic tools to filter and identify underlying candidate genes. Candidate gene function prediction was performed using the online prediction tool STRING 12.0. Immunohistochemistry and immunofluorescence were used to detect the expression of proteins encoded by candidate genes. We mapped differential single nucleotide polymorphisms (SNPs) to genes and obtained a total of 102 differential genes, then we used GO and KEGG pathway enrichment analysis to identify four candidate genes, including SLA-1, SLA-2, SLA-3, and TAP2. nsSNPs cause changes in the primary and tertiary structures of SLA-I and TAP2 proteins, the primary structures of these two proteins have undergone amino acid changes, and the tertiary structures also show slight changes. In addition, immunohistochemistry and immunofluorescence results showed that the expression changes of TAP2 protein in coronary arteries showed a trend of increasing from the middle layer to the inner layer. We have identified SLA-I and TAP2 as potential susceptibility genes of atherosclerosis, highlighting the importance of antigen processing and immune response in atherogenesis.

Whole-genome sequencing of Ganoderma boninense, the causal agent of basal stem rot disease in oil palm, via combined short- and long-read sequencing

The hemibiotrophic Basidiomycete pathogen Ganoderma boninense (Gb) is the dominant causal agent of oil palm basal stem rot disease. Here, we report a complete chromosomal genome map of Gb using a combination of short-read Illumina and long-read Pacific Biosciences (PacBio) sequencing platforms combined with chromatin conformation capture data from the Chicago and Hi-C platforms. The genome was 55.87 Mb in length and assembled to a high contiguity (N50: 304.34 kb) of 12 chromosomes built from 112 scaffolds, with a total of only 4.34 Mb (~ 7.77%) remaining unplaced. The final assemblies were evaluated for completeness of the genome by using Benchmarking Universal Single Copy Orthologs (BUSCO) v4.1.4, and based on 4464 total BUSCO polyporales group searches, the assemblies yielded 4264 (95.52%) of the conserved orthologs as complete and only a few fragmented BUSCO of 42 (0.94%) as well as a missing BUSCO of 158 (3.53%). Genome annotation predicted a total of 21,074 coding genes, with a GC content ratio of 59.2%. The genome features were analyzed with different databases, which revealed 2471 Gene Ontology/GO (11.72%), 5418 KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthologous/KO (25.71%), 13,913 Cluster of Orthologous Groups of proteins/COG (66.02%), 60 ABC transporter (0.28%), 1049 Carbohydrate-Active Enzymes/CAZy (4.98%), 4005 pathogen–host interactions/PHI (19%), and 515 fungal transcription factor/FTFD (2.44%) genes. The results obtained in this study provide deep insight for further studies in the future.

Phage–Bacterial Interaction Alters Phenotypes Associated with Virulence in Acinetobacter baumannii

Bacteriophages exert strong selection on their bacterial hosts to evolve resistance. At the same time, the fitness costs on bacteria following phage resistance may change their virulence, which may affect the therapeutic outcomes of phage therapy. In this study, we set out to assess the costs of phage resistance on the in vitro virulence of priority 1 nosocomial pathogenic bacterium, Acinetobacter baumannii. By subjecting phage-resistant variant Ev5-WHG of A. baumannii WHG40004 to several in vitro virulence profiles, we found that its resistance to phage is associated with reduced fitness in host microenvironments. Also, the mutant exhibited impaired adhesion and invasion to mammalian cells, as well as increased susceptibility to macrophage phagocytosis. Furthermore, the whole-genome sequencing of the mutant revealed that there exist multiple mutations which may play a role in phage resistance and altered virulence. Altogether, this study demonstrates that resistance to phage can significantly alter phenotypes associated with virulence in Acinetobacter baumannii.

Short beak and dwarfism syndrome among Pekin ducks: First detection, full genome sequencing, and immunohistochemical signals of novel goose parvovirus in tongue tissue.

Novel goose parvovirus (NGPV) is continuously threatening the global duck industry, as it causes short beak and dwarfism syndrome among different duck breeds. In this study, we investigated the viral pathogenesis in the tongue of affected ducks, as a new approach for deeper understanding of the syndrome. Seventy-three, 14- to 60-day-old commercial Pekin ducks were clinically examined. Thirty tissue pools of intestine and tongue (15 per tissue) were submitted for molecular identification. Clinical signs in the examined ducks were suggestive of parvovirus infection. All examined ducks had short beaks. Necrotic, swollen, and congested protruding tongues were recorded in adult ducks (37/73, 51%). Tongue protrusion without any marked congestion or swelling was observed in 20-day-old ducklings (13/73, 18%), and no tongue protrusion was observed in 15-day-old ducklings (23/73, 32%). Microscopically, the protruding tongues of adult ducks showed necrosis of the superficial epithelial layer with vacuolar degeneration. Glossitis was present in the nonprotruding tongues of young ducks, which was characterized by multifocal lymphoplasmacytic aggregates and edema in the propria submucosa. Immunohistochemical examination displayed parvovirus immunolabeling, mainly in the tongue propria submucosa. Based on polymerase chain reaction, goose parvovirus was detected in 9 out of 15 tongue sample pools (60%). Next-generation sequencing confirmed the presence of a variant goose parvovirus that is globally named NGPV and closely related to Chinese NGPV isolates. Novel insights are being gained from the study of NGPV pathogenesis in the tongue based on molecular and immunohistochemical identification.

Functional variants in a TTTG microsatellite on 15q26.1 cause familial nonautoimmune thyroid abnormalities.

Insufficient thyroid hormone production in newborns is referred to as congenital hypothyroidism. Multinodular goiter (MNG), characterized by an enlarged thyroid gland with multiple nodules, is usually seen in adults and is recognized as a separate disorder from congenital hypothyroidism. Here we performed a linkage analysis of a family with both nongoitrous congenital hypothyroidism and MNG and identified a signal at 15q26.1. Follow-up analyses with whole-genome sequencing and genetic screening in congenital hypothyroidism and MNG cohorts showed that changes in a noncoding TTTG microsatellite on 15q26.1 were frequently observed in congenital hypothyroidism (137 in 989) and MNG (3 in 33) compared with controls (3 in 38,722). Characterization of the noncoding variants with epigenomic data and in vitro experiments suggested that the microsatellite is located in a thyroid-specific transcriptional repressor, and its activity is disrupted by the variants. Collectively, we presented genetic evidence linking nongoitrous congenital hypothyroidism and MNG, providing unique insights into thyroid abnormalities.

Complete genome sequence of two chromosomes of Vibrio metoecus strain ZF102 isolated from the abdominal cavity of moribund laboratory zebrafish (Danio rerio).

Vibrio metoecus was isolated from the abdominal cavity of moribund laboratory zebrafish. We report complete genomic sequences of V. metoecus strain ZF102 that has two circular chromosomes of 2,872,299 and 1,170,691 bp and two plasmids of 5,265 and 2,361 bp.

Biography of Vitis genomics: recent advances and prospective

Abstract Grape genome is the basis for grape studies and breeding, which is also important for grape industries. In the last two decades, more than 44 grape genomes have been sequenced. Based on these genomes, researchers have made a substantial progress in the mechanism of biotic and abiotic resistance, berry quality formation, and breeding strategies. In addition, these works provide essential data for future pan-genome analysis. Apart from the de novo assembled genomes, more than six whole-genome sequencing projects have provided datasets comprising almost 5000 accessions. Based on these datasets, researchers explored domestications and origin of grape and clarified the gene flow that occurred during the dispersed history. Moreover, GWAS and other methods were used to identify more than 900 genes related to the resistance, quality, and developmental phases of grape. These findings have benefit grape studies and provide some basis for smart genomic selection breeding. Otherwise, the grape genome has played a great role in grape studies and industries, and the importance of genomics will increase sharply in the future.