Whole Genome Sequencing Research Articles

New clinical treatment plan for acute lymphoblastic leukemia after whole genome sequencing: An update of the development of new technologies for genome editing

Acute lymphoblastic leukemia (ALL) is now being reclassified into newly identified subtypes due to genetic abnormalities, Philadelphia chromosome‐like B‐lineage ALL is one of the new high‐risk subtypes characterized by genetic alterations that activate various signaling pathways, including those involving cytokine receptors, tyrosine kinases, and epigenetic modifiers. Philadelphia chromosome‐like ALL is essentially heterogeneous; however, deletion mutations in the IKZF1 gene encoding the transcription factor IKAROS underlie many cases as a key factor inducing aggressive phenotypes and poor treatment responses. WGS (whole genome sequencing) studies of ALL patients and ethnically matched controls also identified inherited genetic variations in lymphoid neoplasm‐related genes, which are likely to increase ALL susceptibility. On the contrary, Genome editing technologies offers new opportunities for tackling diseases, and genome-editing tools such as CRISPR-Cas9 are an important means of advancing functional studies of ALL through the incorporation, elimination, and modification of somatic mutations. These findings are directly relevant to clinical hematology, and further studies on this aspect could contribute to accurate diagnosis, effective monitoring of residual disease, and patient‐oriented therapies.

Identification of a novel first-generation HIV-1 circulating recombinant form (CRF152_DG) among people living with HIV in Karachi, Pakistan.

In Pakistan, the genetic diversity of human immunodeficiency virus type 1 (HIV-1) is becoming increasingly complex, compared to the early years of the epidemic that started after the detection of the first cases of HIV-1 in 1987 in Karachi. Based on the available molecular studies, two dominant HIV-1 clades, sub-subtype A1 and CRF02_AG, have been found to co-circulate with other clades, namely B, C, D, G, CRF01_AE, CRF35_A1D, and CRF56_cpx, in various urban areas of Pakistan. Several novel recombinant forms have also been detected. This first report of CRF152_DG highlights the complex nature of the HIV epidemic in Pakistan and emphasizes the importance of continual molecular surveillance (ideally based on whole-genome sequences) of HIV.

Lavender Harbors More Viruses than Previously Thought: First Report of Raspberry Ringspot Virus and Phlox Virus M in Lavandula × intermedia.

Plants of the genus Lavandula are thought to be rarely infected by viruses. To date, only alfalfa mosaic virus, cucumber mosaic virus, tobacco mosaic virus, and tomato spotted wilt virus have been reported in this host. In this study, we identified for the first time raspberry ringspot virus (RpRSV) and phlox virus M (PhlVM) in lavender using herbaceous indexing, enzyme-linked immunosorbent assay, and high-throughput sequencing. Nearly complete genome sequences for both viruses were determined. Phylogenetic and serological characterizations suggest that the obtained RpRSV isolate is a raspberry strain. A preliminary survey of 166 samples indicated RpRSV was spread only in the lavender cultivar 'Grosso', while PhlVM was detected in multiple lavender cultivars. Although RpRSV raspberry strain may have spread throughout Auckland and nearby areas in New Zealand, it is very likely restricted to the genus Lavandula or even to the cultivar 'Grosso' due to the absence or limited occurrence of the nematode vector. Interestingly, all infected lavender plants, regardless of their infection status (by RpRSV, PhlVM, or both) were asymptomatic. RpRSV is an important virus that infects horticultural crops including grapevine, cherry, berry fruits, and rose. It remains on the list of regulated pests in New Zealand. RpRSV testing is mandatory for imported Fragaria, Prunus, Ribes, Rosa, Rubus, and Vitis nursery stock and seeds for sowing, while this is not required for Lavandula importation. Our study revealed that lavender could play a role not only as a reservoir but also as an uncontrolled import pathway of viruses that pose a threat to New Zealand's primary industries.

Population structure and history of Mycobacterium bovis European 3 clonal complex reveal transmission across ecological corridors of unrecognized importance in Portugal.

Efforts to strengthen surveillance and control of animal tuberculosis (TB) are ongoing worlwide. Here, we developed an eco-phylodynamic framework based on discrete phylogenetic approaches informed by M. bovis whole-genome sequence data representing a multi-host transmission system at the livestock-wildlife interface, within a rich ecological landscape in Portugal, to understand transmission processes and translate this knowledge into disease management benefits. We find evidence for the co-circulation of several M. bovis clades, with frequent transmission of the Eu3 lineage among cattle and wildlife populations. Most transition events between different ecological settings took place toward host, climate and land use gradients, underscoring animal TB expansion and a potential corridor of unrecognized importance for M. bovis maintenance. Results stress that animal TB is an established wildlife disease without ecological barriers, showing that control measures in place are insufficient to prevent long-distance transmission and spillover across multi-host communities, demanding new interventions targeting livestock-wildlife interactions.

Improving nitrofurantoin resistance prediction in Escherichia coli from whole-genome sequence by integrating NfsA/B enzyme assays.

Nitrofurantoin resistance in Escherichia coli is primarily caused by mutations damaging two enzymes, NfsA and NfsB. Studies based on small isolate collections with defined nitrofurantoin MICs have found significant random genetic drift in nfsA and nfsB, making it extremely difficult to predict nitrofurantoin resistance from whole-genome sequence (WGS) where both genes are not obviously disrupted by nonsense or frameshift mutations or insertional inactivation. Here, we report a WGS survey of 200 oqxAB-negative E. coli from community urine samples, of which 34 were nitrofurantoin resistant. We characterized individual non-synonymous mutations seen in nfsA and nfsB among this collection using complementation cloning and NfsA/B enzyme assays in cell extracts. We definitively identified R203C, H11Y, W212R, A112E, and A112T in NfsA and R121C, Q142H, F84S, P163H, W46R, K57E, and V191G in NfsB as amino acid substitutions that reduce enzyme activity sufficiently to cause resistance. In contrast, E58D, I117T, K141E, L157F, A172S, G187D, and A188V in NfsA and G66D, M75I, V93A, and A174E in NfsB are functionally silent in this context. We identified that 9/166 (5.4%) nitrofurantoin-susceptible isolates were "pre-resistant," defined as having loss of function mutations in nfsA or nfsB. Finally, using NfsA/B enzyme assays and proteomics, we demonstrated that 9/34 (26.5%) ribE wild-type nitrofurantoin-resistant isolates also carried functionally wild-type nfsB or nfsB/nfsA. In these cases, NfsA/B activity was reduced through downregulated gene expression. Our biological understanding of nitrofurantoin resistance is greatly improved by this analysis but is still insufficient to allow its reliable prediction from WGS data.

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

BackgroundTuberculosis 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.MethodsEpidemiologically 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.ResultsWe 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.ConclusionThe 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.

Insights into the biocontrol and plant growth promotion functions of Bacillus altitudinis strain KRS010 against Verticillium dahliae

BackgroundVerticillium wilt, caused by the fungus Verticillium dahliae, is a soil-borne vascular fungal disease, which has caused great losses to cotton yield and quality worldwide. The strain KRS010 was isolated from the seed of Verticillium wilt-resistant Gossypium hirsutum cultivar “Zhongzhimian No. 2.”ResultsThe strain KRS010 has a broad-spectrum antifungal activity to various pathogenic fungi as Verticillium dahliae, Botrytis cinerea, Fusarium spp., Colletotrichum spp., and Magnaporthe oryzae, of which the inhibition rate of V. dahliae mycelial growth was 73.97% and 84.39% respectively through confrontation test and volatile organic compounds (VOCs) treatments. The strain was identified as Bacillus altitudinis by phylogenetic analysis based on complete genome sequences, and the strain physio-biochemical characteristics were detected, including growth-promoting ability and active enzymes. Moreover, the control efficiency of KRS010 against Verticillium wilt of cotton was 93.59%. After treatment with KRS010 culture, the biomass of V. dahliae was reduced. The biomass of V. dahliae in the control group (Vd991 alone) was 30.76-folds higher than that in the treatment group (KRS010+Vd991). From a molecular biological aspect, KRS010 could trigger plant immunity by inducing systemic resistance (ISR) activated by salicylic acid (SA) and jasmonic acid (JA) signaling pathways. Its extracellular metabolites and VOCs inhibited the melanin biosynthesis of V. dahliae. In addition, KRS010 had been characterized as the ability to promote plant growth.ConclusionsThis study indicated that B. altitudinis KRS010 is a beneficial microbe with a potential for controlling Verticillium wilt of cotton, as well as promoting plant growth.

Genome analyses of colistin-resistant high-risk blaNDM-5 producing Klebsiella pneumoniae ST147 and Pseudomonas aeruginosa ST235 and ST357 in clinical settings

BackgroundColistin is a last-resort antibiotic used in extreme cases of multi-drug resistant (MDR) Gram-negative bacterial infections. Colistin resistance has increased in recent years and often goes undetected due to the inefficiency of predominantly used standard antibiotic susceptibility tests (AST). To address this challenge, we aimed to detect the prevalence of colistin resistance strains through both Vitek®2 and broth micro-dilution. We investigated 1748 blood, tracheal aspirate, and pleural fluid samples from the Intensive Care Unit (ICU), Neonatal Intensive Care Unit (NICU), and Tuberculosis and Respiratory Disease centre (TBRD) in an India hospital. Whole-genome sequencing (WGS) of extremely drug-resitant (XDR) and pan-drug resistant (PDR) strains revealed the resistance mechanisms through the Resistance Gene Identifier (RGI.v6.0.0) and Snippy.v4.6.0. Abricate.v1.0.1, PlasmidFinder.v2.1, MobileElementFinder.v1.0.3 etc. detected virulence factors, and mobile genetic elements associated to uncover the pathogenecity and the role of horizontal gene transfer (HGT).ResultsThis study reveals compelling insights into colistin resistance among global high-risk clinical isolates: Klebsiella pneumoniae ST147 (16/20), Pseudomonas aeruginosa ST235 (3/20), and ST357 (1/20). Vitek®2 found 6 colistin-resistant strains (minimum inhibitory concentrations, MIC = 4 μg/mL), while broth microdilution identified 48 (MIC = 32–128 μg/mL), adhering to CLSI guidelines. Despite the absence of mobile colistin resistance (mcr) genes, mechanisms underlying colistin resistance included mgrB deletion, phosphoethanolamine transferases arnT, eptB, ompA, and mutations in pmrB (T246A, R256G) and eptA (V50L, A135P, I138V, C27F) in K. pneumoniae. P. aeruginosa harbored phosphoethanolamine transferases basS/pmrb, basR, arnA, cprR, cprS, alongside pmrB (G362S), and parS (H398R) mutations. Both strains carried diverse clinically relevant antimicrobial resistance genes (ARGs), including plasmid-mediated blaNDM-5 (K. pneumoniae ST147) and chromosomally mediated blaNDM-1 (P. aeruginosa ST357).ConclusionThe global surge in MDR, XDR and PDR bacteria necessitates last-resort antibiotics such as colistin. However, escalating resistance, particularly to colistin, presents a critical challenge. Inefficient colistin resistance detection methods, including Vitek2, alongside limited surveillance resources, accentuate the need for improved strategies. Whole-genome sequencing revealed alarming colistin resistance among K. pneumoniae and P. aeruginosa in an Indian hospital. The identification of XDR and PDR strains underscores urgency for enhanced surveillance and infection control. SNP analysis elucidated resistance mechanisms, highlighting the complexity of combatting resistance.

Genome-wide association analysis uncovers rice blast resistance alleles of Ptr and Pia

A critical step to maximize the usefulness of genome-wide association studies (GWAS) in plant breeding is the identification and validation of candidate genes underlying genetic associations. This is of particular importance in disease resistance breeding where allelic variants of resistance genes often confer resistance to distinct populations, or races, of a pathogen. Here, we perform a genome-wide association analysis of rice blast resistance in 500 genetically diverse rice accessions. To facilitate candidate gene identification, we produce de-novo genome assemblies of ten rice accessions with various rice blast resistance associations. These genome assemblies facilitate the identification and functional validation of novel alleles of the rice blast resistance genes Ptr and Pia. We uncover an allelic series for the unusual Ptr rice blast resistance gene, and additional alleles of the Pia resistance genes RGA4 and RGA5. By linking these associations to three thousand rice genomes we provide a useful tool to inform future rice blast breeding efforts. Our work shows that GWAS in combination with whole-genome sequencing is a powerful tool for gene cloning and to facilitate selection of specific resistance alleles for plant breeding.

Complete genome sequence of Bacillus subtilis NA05 (=NBRC 116153), a phosphate-solubilizing bacterium isolated from crop field soil.

We report the complete genome sequence of the phosphate-solubilizing bacterium Bacillus subtilis NA05 (=NBRC 116153), consisting of a circular chromosome of ~3.8 M bp and two circular plasmids. The data presented here provide further insight into the genetic and functional potential of B. subtilis and the mechanism of phosphate solubilization.

Complete genome sequence of the linezolid-resistant clinical Enterococcus faecalis N23-3408 linked to a livestock lineage in Switzerland.

Here, we report the complete genome of human clinical linezolid-resistant Enterococcus faecalis N23-3408. N23-3408 harbored a 59.5 kb plasmid with antimicrobial resistance genes cat, erm(B), fexA, optrA, tet(L), and tet(M). Closely related E. faecalis harboring this plasmid was previously obtained from livestock animals and pet food in Switzerland.

Complete genome sequence of potential plant growth-promoting Bacillus altitudinis strain AIMST-CREST03 isolated from paddy field bulk soil.

We present the complete genome of a potential plant growth-promoting bacteria Bacillus altitudinis AIMST-CREST03 isolated from a high-yielding paddy plot. The genome is 3,669,202 bp in size with a GC content of 41%. Annotation predicted 3,327 coding sequences, including several genes required for plant growth promotion.

Whole-genome sequence of six Rhizobium laguerreae strains.

Rhizobium laguerreae is regarded as a promising candidate for biofertilization of legume plants worldwide through its high efficiency in symbiosis. In this paper, we report high-quality sequences of six R. laguerreae strains with total genome completeness from 93.5% to 97.5%.

Fragmentomics features of ovarian cancer.

Ovarian cancer (OC) is a major cause of cancer mortality in women worldwide. Due to the occult onset of OC, its nonspecific clinical symptoms in the early phase, and a lack of effective early diagnostic tools, most OC patients are diagnosed at an advanced stage. In this study, shallow whole-genome sequencing was utilized to characterize fragmentomics features of circulating tumor DNA (ctDNA) in OC patients. By applying a machine learning model, multiclass fragmentomics data achieved a mean area under the curve (AUC) of 0.97 (95% CI 0.962-0.976) for diagnosing OC. OC scores derived from this model strongly correlated with the disease stage. Further comparative analysis of OC scores illustrated that the fragmentomics-based technology provided additional clinical benefits over the traditional serum biomarkers cancer antigen 125 (CA125) and the Risk of Ovarian Malignancy Algorithm (ROMA) index. In conclusion, fragmentomics features in ctDNA are potential biomarkers for the accurate diagnosis of OC.

Research progress on genetics in cardioembolic stroke.

Cardioembolic stroke, characterized by severe illness, poor prognosis, and high recurrence rate, is one of the important causes of ischemic stroke. In the field of genetic research, numerous genes associated with cardioembolic stroke have been identified, and their potential in predicting disease risk and evaluating risk factors has been progressively explored. Here, we provide an overview of the latest advancements in genetics for cardioembolic stroke, including genome-wide association studies, copy number variation studies, whole-genome sequencing studies. Furthermore, we also summarize the application of genetic datasets in polygenic risk score and Mendelian randomization. The aim of this overview is to provide insights and references from multiple perspectives for future investigations on the genetic information for cardioembolic stroke.

Whole-genome sequence of the strictly anaerobic bacterial strain SANA belonging to the family Gottschalkiaceae, isolated from a xenic culture of an anaerobic protist.

An anaerobic bacterial strain SANA was isolated from a xenic culture of an anaerobic heterolobosean protist which was obtained from a saline lake in Japan. Its draft genome comprises 1 circular chromosome (3,490,293 bp), harboring 3,275 predicted protein-coding and 73 tRNA-encoding genes and 8 rRNA operons.

Detection and evaluation of susceptibility to antibiotics in non-hydrogen sulfide-producing antibiotic-resistant soil microbe: Pseudomonas guariconensis.

Antimicrobial resistance in bacteria is a global threat that can make antibacterial treatments ineffective. One well-known method of antibiotic resistance and a common defensive mechanism in many harmful bacteria is the synthesis of endogenous hydrogen sulfide (H2S) in bacteria. In this study, soil bacteria were screened using the lead acetate agar test and the triple sugar iron test to determine that they were non-endogenous H2S producers. This was further validated by full genome analysis of the identified organism against the gene sequences of H2S-producing genes. Antibacterial resistance of the bacteria was phenotypically analyzed using the Kirby-Bauer disk diffusion method. Then, the effect of exogenous H2S on the antibiotic-resistant bacteria was checked in sodium sulfide, leading to antibiotic re-sensitization.

Genetic characterization of bovine coronavirus strain isolated in Inner Mongolia of China

BackgroundBovine coronavirus (BCoV) is implicated in severe diarrhea in calves and contributes to the bovine respiratory disease complex; it shares a close relationship with human coronavirus. Similar to other coronaviruses, remarkable variability was found in the genome and biology of the BCoV. In 2022, samples of feces were collected from a cattle farm. A virus was isolated from 7-day-old newborn calves. In this study, we present the genetic characteristics of a new BCoV isolate. The complete genomic, spike protein, and nucleocapsid protein gene sequences of the BCoV strain, along with those of other coronaviruses, were obtained from the GenBank database. Genetic analysis was conducted using MEGA7.0 and the Neighbor-Joining (NJ) method. The reference strains’ related genes were retrieved from GenBank for comparison and analysis using DNAMAN.ResultsThe phylogenetic tree and whole genome consistency analysis showed that it belonged to the GIIb subgroup, which is epidemic in Asia and America, and was quite similar to the Chinese strains in the same cluster. Significantly, the S gene was highly consistent with QH1 (MH810151.1) isolated from yak. This suggests that the strain may have originated from interspecies transmission involving mutations of wild strains. The N gene was conserved and showed high sequence identity with the epidemic strains in China and the USA.ConclusionsGenetic characterization suggests that the isolated strain could be a new mutant from a wild-type lineage, which is in the same cluster as most Chinese epidemic strains but on a new branch.

Isolation of a CTX-M-55 (ESBL)-Producing Escherichia coli Strain of the Global ST6448 Clone from a Captive Orangutan in the USA.

The purpose of this study was to determine if orangutans (Pongo spp.) living in captivity at a zoo in Wisconsin were colonized with antimicrobial-resistant bacteria and, if found, to identify underlying genetic mechanisms contributing to their resistant phenotypes. We hypothesize that since antimicrobial-resistant bacteria are so prevalent within humans, the animals could also be carriers of such strains given the daily contact between the animals and the zoo staff that care for them. To test this theory, fecal samples from two orangutans were examined for resistant bacteria by inoculation on HardyCHROM™ ESBL and HardyCHROM™ CRE agars. Isolates were identified using MALDI-TOF mass spectrometry and antimicrobial susceptibility testing was performed using a Microscan autoSCAN-4 System. An isolate was selected for additional characterization, including whole genome sequencing (WGS). Using the Type (Strain) Genome Server (TYGS) the bacterium was identified as Escherichia coli. The sequence type identified was (ST/phylogenetic group/β-lactamase): ST6448/B1/CTX-M-55.

Negligible effects of read trimming on the accuracy of germline short variant calling in the human genome

Background Next generation sequencing (NGS) has become a standard tool in the molecular diagnostics of Mendelian disease, and the precision of such diagnostics is greatly affected by the accuracy of variant calling from sequencing data. Recently, we have comprehensively evaluated the performance of multiple variant calling pipelines. However, no systematic analysis of the effects of read trimming on variant discovery with modern variant calling software has yet been performed. Methods In this work, we systematically evaluated the effects of adapters on the performance of 8 variant calling and filtering methods using 14 standard reference Genome-in-a-Bottle (GIAB) samples. Variant calls were compared to the ground truth variant sets, and the effect of adapter trimming with different tools was assessed using major performance metrics (precision, recall, and F1 score). Results We show that adapter trimming has no effect on the accuracy of the best-performing variant callers (e.g., DeepVariant) on whole-genome sequencing (WGS) data. For whole-exome sequencing (WES) datasets subtle improvement of accuracy was observed in some of the samples. In high-coverage WES data (~200x mean coverage), adapter removal allowed for discovery of 2-4 additional true positive variants in only two out of seven datasets tested. Moreover, this effect was not dependent on the median insert size and proportion of adapter sequences in reads. Surprisingly, the effect of trimming on variant calling was reversed when moderate coverage (~80-100x) WES data was used. Finally, we show that some of the recently developed machine learning-based variant callers demonstrate greater dependence on the presence of adapters in reads. Conclusions Taken together, our results indicate that adapter removal is unnecessary when calling germline variants, but suggest that preprocessing methods should be carefully chosen when developing and using machine learning-based variant analysis methods.