Pasa: leveraging population pangenome graph to scaffold prokaryote genome assemblies.

Whole genome sequencing has increasingly become the essential method for studying the genetic mechanisms of antimicrobial resistance and for surveillance of drug-resistant bacterial pathogens. The majority of bacterial genomes sequenced to date have been sequenced with Illumina sequencing technology, owing to its high-throughput, excellent sequence accuracy, and low cost. However, because of the short-read nature of the technology, these assemblies are fragmented into large numbers of contigs, hindering the obtaining of full information of the genome. We develop Pasa, a graph-based algorithm that utilizes the pangenome graph and the assembly graph information to improve scaffolding quality. By leveraging the population information of the bacteria species, Pasa is able to utilize the linkage information of the gene families of the species to resolve the contig graph of the assembly. We show that our method outperforms the current state of the arts in terms of accuracy, and at the same time, is computationally efficient to be applied to a large number of existing draft assemblies.

AMRViz enables seamless genomics analysis and visualization of antimicrobial resistance.

We have developed AMRViz, a toolkit for analyzing, visualizing, and managing bacterial genomics samples. The toolkit is bundled with the current best practice analysis pipeline allowing researchers to perform comprehensive analysis of a collection of samples directly from raw sequencing data with a single command line. The analysis results in a report showing the genome structure, genome annotations, antibiotic resistance and virulence profile for each sample. The pan-genome of all samples of the collection is analyzed to identify core- and accessory-genes. Phylogenies of the whole genome as well as all gene clusters are also generated. The toolkit provides a web-based visualization dashboard allowing researchers to interactively examine various aspects of the analysis results. Availability: AMRViz is implemented in Python and NodeJS, and is publicly available under open source MIT license at https://github.com/amromics/amrviz .

AMRomics: a scalable workflow to analyze large microbial genome collections.

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

The interplay of membrane fluidity, acyl chain order and area per lipid on the partitioning of two antidepressants paroxetine and sertraline.

Selective serotonin reuptake inhibitors (SSRIs) are among the popular drugs for treating depression and mental disorders. Membrane fluidity has previously been considered as the main factor in modulating the membrane partitioning of SSRIs, while other biophysical properties, such as the acyl chain order and area per lipid, were often neglected. Varying the lipid membrane composition and temperature can significantly modify the physical phase and, in turn, affect its fluidity, acyl chain order and area per lipid. Here, we investigate the role of membrane fluidity, acyl chain order and area per lipid in the partitioning of two SSRIs, paroxetine (PAX) and sertraline (SER). The model membranes were either POPC : SM (1 : 1 mol ratio) or POPC : SM : Chol (1 : 1 : 1 mol ratio) and studied in the temperature range of 25-45 °C. The order parameters and area per lipid in the two lipid mixtures were calculated using molecular dynamics simulations. The membrane partitioning of PAX and SER was determined via second derivative spectrophotometry. In a lower temperature range (25-32 °C), membrane fluidity favors the SSRI partitioning into Lo/Ld POPC:SM:Chol. In a higher temperature range (37-45 °C), the interplay between membrane fluidity, acyl chain order and area per lipid favors drug partitioning into Ld POPC:SM. The findings offer indication for the inconsistent distribution of SSRIs in tissues as well as the possible interaction of SSRIs with lipid domains and membrane-bound proteins.

Evaluation of the expression levels of BRAFV600E mRNA in primary tumors of thyroid cancer using an ultrasensitive mutation assay.

BACKGROUND: The BRAFV600E gene encodes for the mutant BRAFV600E protein, which triggers downstream oncogenic signaling in thyroid cancer. Since most currently available methods have focused on detecting BRAFV600E mutations in tumor DNA, there is limited information about the level of BRAFV600E mRNA in primary tumors of thyroid cancer, and the diagnostic relevance of these RNA mutations is not known. METHODS: Sixty-two patients with thyroid cancer and non-malignant thyroid disease were included in the study. Armed with an ultrasensitive technique for mRNA-based mutation analysis based on a two step RT-qPCR method, we analysed the expression levels of the mutated BRAFV600E mRNA in formalin-fixed paraffin-embedded samples of thyroid tissues. Sanger sequencing for detection of BRAFV600E DNA was performed in parallel for comparison and normalization of BRAFV600E mRNA expression levels. RESULTS: The mRNA-based mutation detection assay enables detection of the BRAFV600E mRNA transcripts in a 10,000-fold excess of wildtype BRAF counterparts. While BRAFV600E mutations could be detected by Sanger sequencing in 13 out of 32 malignant thyroid cancer FFPE tissue samples, the mRNA-based assay detected mutations in additionally 5 cases, improving the detection rate from 40.6 to 56.3%. Furthermore, we observed a surprisingly large, 3-log variability, in the expression level of the BRAFV600E mRNA in FFPE samples of thyroid cancer tissue. CONCLUSIONS: The expression levels of BRAFV600E mRNA was characterized in the primary tumors of thyroid cancer using an ultrasensitive mRNA-based mutation assay. Our data inspires further studies on the prognostic and diagnostic relevance of the BRAFV600E mRNA levels as a molecular biomarker for the diagnosis and monitoring of various genetic and malignant diseases.

Extendable blocking probe in reverse transcription for analysis of RNA variants with superior selectivity.

Here we provide the first strategy to use a competitive Extendable Blocking Probe (ExBP) for allele-specific priming with superior selectivity at the stage of reverse transcription. In order to analyze highly similar RNA variants, a reverse-transcriptase primer whose sequence matches a specific variant selectively primes only that variant, whereas mismatch priming to the alternative variant is suppressed by virtue of hybridization and subsequent extension of the perfectly matched ExBP on that alternative variant template to form a cDNA-RNA hybrid. This hybrid will render the alternative RNA template unavailable for mismatch priming initiated by the specific primer in a hot-start protocol of reverse transcription when the temperature decreases to a level where such mismatch priming could occur. The ExBP-based reverse transcription assay detected BRAF and KRAS mutations in at least 1000-fold excess of wild-type RNA and detection was linear over a 4-log dynamic range. This novel strategy not only reveals the presence or absence of rare mutations with an exceptionally high selectivity, but also provides a convenient tool for accurate determination of RNA variants in different settings, such as quantification of allele-specific expression.

Innovative Semi-Nested Realtime PCR Assay with Extendable Blocking Probe for Enhanced Analysis of SEPT9 Methylation in Colorectal Cancer.

(1) Background: The detection of methylated SEPT9 (mSEPT9) in plasma is a promising approach to non-invasive colorectal cancer (CRC) screening. Traditional approaches have limitations in sensitivity and cost-effectiveness, particularly in resource-limited settings. (2) Methods: We developed a semi-nested realtime PCR assay utilizing extendable blocking probes (ExBP) to enhance the detection of low-level mSEPT9 based on DNA melting. This assay allows for the discrimination of mSEPT9 in the presence of high concentrations of non-methylated SEPT9 (up to 100,000 times higher). (3) Results: The assay demonstrated a sensitivity of 73.91% and specificity of 80%, showcasing its ability to detect very low levels of methylated DNA effectively. The innovative use of ExBP without costly modified probes simplifies the assay setup and reduces the overall costs, enhancing its applicability in diverse clinical settings. (4) Conclusions: This novel assay significantly improves the detection of mSEPT9, offering a potential advance in CRC screening and monitoring. Its cost-efficiency and high sensitivity make it particularly suitable for the early detection and management of CRC, especially in settings with limited resources. Future studies are encouraged to validate this assay in larger populations to establish its clinical benefits and practical utility.

Novel Semi-Nested Real-Time PCR Assay Leveraging Extendable Blocking Probes for Improved SHOX2 Methylation Analysis in Lung Cancer.

Lung cancer is the leading cause of cancer deaths globally, necessitating effective early detection methods. Traditional diagnostics like low-dose computed tomography (LDCT) often yield high false positive rates. SHOX2 gene methylation has emerged as a promising biomarker. This study aimed to develop and validate a novel semi-nested real-time PCR assay enhancing sensitivity and specificity for detecting SHOX2 methylation using extendable blocking probes (ExBPs). The assay integrates a semi-nested PCR approach with ExBPs, enhancing the detection of low-abundance methylated SHOX2 DNA amidst unmethylated sequences. It was tested on spiked samples with varied methylation levels and on clinical samples from lung cancer patients and individuals with benign lung conditions. The assay detected methylated SHOX2 DNA down to 0.01%. Clinical evaluations confirmed its ability to effectively differentiate between lung cancer patients and those with benign conditions, demonstrating enhanced sensitivity and specificity. The use of ExBPs minimized non-target sequence amplification, crucial for reducing false positives. The novel semi-nested real-time PCR assay offers a cost-effective, highly sensitive, and specific method for detecting SHOX2 methylation, enhancing early lung cancer detection and monitoring, particularly valuable in resource-limited settings.

Efficient inference of large prokaryotic pangenomes with PanTA.

Pangenome inference is an indispensable step in bacterial genomics, yet its scalability poses a challenge due to the rapid growth of genomic collections. This paper presents PanTA, a software package designed for constructing pangenomes of large bacterial datasets, showing unprecedented efficiency levels multiple times higher than existing tools. PanTA introduces a novel mechanism to construct the pangenome progressively without rebuilding the accumulated collection from scratch. The progressive mode is shown to consume orders of magnitude less computational resources than existing solutions in managing growing datasets. The software is open source and is publicly available at https://github.com/amromics/panta and at 10.6084/m9.figshare.23724705 .

PanKA: Leveraging population pangenome to predict antibiotic resistance.

Machine learning has the potential to be a powerful tool in the fight against antimicrobial resistance (AMR), a critical global health issue. Machine learning can identify resistance mechanisms from DNA sequence data without prior knowledge. The first step in building a machine learning model is a feature extraction from sequencing data. Traditional methods like single nucleotide polymorphism (SNP) calling and k-mer counting yield numerous, often redundant features, complicating prediction and analysis. In this paper, we propose PanKA, a method using the pangenome to extract a concise set of relevant features for predicting AMR. PanKA not only enables fast model training and prediction but also improves accuracy. Applied to the Escherichia coli and Klebsiella pneumoniae bacterial species, our model is more accurate than conventional and state-of-the-art methods in predicting AMR.

Fibroblast growth factor-21 (FGF21) regulates low-density lipoprotein receptor (LDLR) levels in cells via the E3-ubiquitin ligase Mylip/Idol and the Canopy2 (Cnpy2)/Mylip-interacting saposin-like protein (Msap).

The LDLR is a critical factor in the regulation of blood cholesterol levels that are altered in different human diseases. The level of LDLR in the cell is regulated by both transcriptional and post-transcriptional events. The E3 ubiquitin ligase, myosin regulatory light chain-interacting protein (Mylip)/inducible degrader of the LDL-R (Idol) was shown to induce degradation of LDLR via protein ubiquitination. We have here studied novel factors and mechanisms that may regulate Mylip/Idol in human hepatocyte cells and in mouse macrophages. We observed that FGF21 that is present in serum in different conditions reduced Mylip/Idol at the RNA and protein level, and increased LDLR levels and stability in the cells. FGF21 also enhanced expression of Canopy2 (Cnpy2)/MIR-interacting Saposin-like protein (Msap) that is known to interact with Mylip/Idol. Overexpression of Cnpy2/Msap increased LDLRs, and knockdown experiments showed that Cnpy2/Msap is crucial for the FGF21 effect on LDLRs. Experiments using DiI-labeled LDL particles showed that FGF21 increased lipoprotein uptake and the effect of FGF21 was additive to that of statins. Our results are consistent with an important role of FGF21 and Cnpy2/Msap in the regulation of LDLRs in cultured cells, which warrants further studies using human samples.

Selective detection of HBV pre-genomic RNA in chronic hepatitis B patients using a novel RT-PCR assay.

In chronic hepatitis B (CHB) patients, quantification of HBV pgRNA in plasma has the potential to provide information on disease prognosis and liver injury or histopathology. However, current methods for detecting HBV pgRNA present technical difficulties due to the co-existence of HBV DNA in plasma samples. We have successfully established a novel one-step RT-PCR assay that allows selective quantification of HBV pgRNA. Two cohorts of participants were recruited for assay validation, including treatment-naïve patients with CHB and HBeAg-positive CHB patients who were treated with Tenofovir and monitored for 6 months to assess the predictive value of baseline HBV RNA for HBeAg seroclearance. Statistical analysis was performed using MedCalc version 20.019 software. The novel selective one-step RT-PCR assay for detecting HBV pgRNA was validated with a limit of detection of 100 copies/mL. The assay was able to selectively measure HBV pgRNA even in the presence of excess HBV rcDNA. In treatment-naïve CHB patients, HBV pgRNA levels were significantly lower than HBV DNA concentration. Serum HBV DNA levels and HBeAg status were positively associated with HBV pgRNA. Baseline serum HBV pgRNA levels were found to be strong predictors of HBeAg seroclearance after 6 months of Tenofovir treatment. The study presents a novel RT-PCR assay that allows accurate measurement of plasma HBV pgRNA in chronic hepatitis B patients, even in the presence of excess HBV DNA. The assay is highly selective and represents a significant advancement with potential for further breakthroughs in understanding the clinical significance of HBV pgRNA.

Multiple heat pulses during PCR extension enabling amplification of GC-rich sequences and reducing amplification bias.

PCR amplification over GC-rich and/or long repetitive sequences is challenging because of thermo-stable structures resulting from incomplete denaturation, reannealing, and self-annealing of target sequences. These structures block the DNA polymerase during the extension step, leading to formation of incomplete extension products and favoring amplification of nonspecific products rather than specific ones. We have introduced multiple heat pulses in the extension step of a PCR cycling protocol to temporarily destabilize such blocking structures, in order to enhance DNA polymerase extension over GC-rich sequences. With this novel type of protocol, we were able to amplify all expansions of CGG repeats in five Fragile X cell lines, as well as extremely GC-rich nonrepetitive segments of the GNAQ and GP1BB genes. The longest Fragile X expansion contained 940 CGG repeats, corresponding to about 2.8 kilo bases of 100% GC content. For the GNAQ and GP1BB genes, different length PCR products in the range of 700 bases to 2 kilobases could be amplified without addition of cosolvents. As this technique improves the balance of amplification efficiencies between GC-rich target sequences of different length, we were able to amplify all of the allelic expansions even in the presence of the unexpanded allele.

Formation of lipid raft nanodomains in homogeneous ternary lipid mixture of POPC/DPSM/cholesterol: Theoretical insights.

Lipid rafts, in biological membranes, are cholesterol-rich nanodomains that regulate many protein activities and cellular processes. Understanding the formation of the lipid-raft nanodomains helps us elucidate many complex interactions in the cell. In this study, the formation of lipid-raft nanodomains in a ternary palmitoyl-oleoyl-phosphatidylcholine/stearoyl-sphingomyelin/cholesterol (POPC/DPSM/Chol) lipid mixture, the most realistic surrogate model for biological membranes, has been successfully observed for the first time in-silico using microsecond timescale molecular dynamics simulations. The model reveals the formation of cholesterol-induced nanodomains with raft-like characteristics and their underlying mechanism: the cholesterol molecules segregate themselves into cholesterol nanodomains and then enrich the cholesterol-rich domain with sphingomyelin molecules to form a lipid raft thanks to the weak bonding of cholesterol with sphingomyelin. Besides, it is found that the increase in cholesterol concentration enhances the biophysical properties (e.g., bilayer thickness, area per lipid headgroup, and order parameter) of the lipid raft nanodomains. Such findings suggest that the POPC/DPSM/Chol bilayer is a suitable model to fundamentally extend the nanodomain evolution to investigate their lifetime and kinetics as well as to study protein-lipid interaction, protein-protein interaction, and selection of therapeutic molecules in the presence of lipid rafts.

Clinical Epidemiology Characteristics and Antibiotic Resistance Associated with Urinary Tract Infections Caused by E. coli.

Introduction: In individuals with urinary tract infections, Escherichia coli (E. coli) is an ubiquitous causative agent and antibiotic resistance is on the rise throughout the world. Therefore, early diagnosis and appropriate choice of antimicrobials are essential. The purpose of our study is to describe some of the clinical and epidemiological characteristics and the laboratory test results of children treated in our hospital for urinary tract infections caused by E. coli. Methods: The study included 128 patients from 2 months to 15 years of age with urinary tract infections caused by E. coli and treated at the Haiphong Children's Hospital during the periods of 2011-2013 and 2018-2020. Results: During the two study periods, 57 and 71 cases, respectively, were included. The most common clinical symptom was fever in 40 and 46 cases, respectively. The proportion of E. coli's resistance to ampicillin increased from 85.3% in 2011-2013 to 97.1% in 2018-2020. In 2011-2013, 70.5% of E. coli isolates were resistant to cotrimoxazole, which increased to 81.4% during 2018-2020. During both periods, E. coli was highly sensitive to amikacin, at 87% and 95.5%, respectively. In 2018-2020, carbapenems (meropenem and imipenem) and piperacillin were also effective against E. coli. Conclusion: Our study revealed that high fever was the most prevalent clinical characteristic in urinary tract infections caused by E. coli in children and E. coli was mostly resistant to ampicillin, nalidixic acid, and cotrimoxazole but was highly sensitive to ciprofloxacin, amikacin, piperacillin, meropenem, and imipenem.

Semi-nested RT-PCR enables sensitive and high-throughput detection of SARS-CoV-2 based on melting analysis.

BACKGROUND: Asymptomatic transmission was found to be the Achilles' heel of the symptom-based screening strategy, necessitating the implementation of mass testing to efficiently contain the transmission of COVID-19 pandemic. However, the global shortage of molecular reagents and the low throughput of available realtime PCR facilities were major limiting factors. METHODS: A novel semi-nested and heptaplex (7-plex) RT-PCR assay with melting analysis for detection of SARS-CoV-2 RNA has been established for either individual testing or 96-sample pooled testing. The complex melting spectrum collected from the heptaplex RT-PCR amplicons was interpreted with the support of an artificial intelligence algorithm for the detection of SARS-CoV-2 RNA. The analytical and clinical performance of the semi-nested RT-PCR assay was evaluated using RNAs synthesized in-vitro and those isolated from nasopharyngeal samples. RESULTS: The LOD of the assay for individual testing was estimated to be 7.2 copies/reaction. Clinical performance evaluation indicated a sensitivity of 100% (95% CI: 97.83-100) and a specificity of 99.87% (95% CI: 99.55-99.98). More importantly, the assay supports a breakthrough sample pooling method, which makes possible parallel screening of up to 96 samples in one real-time PCR well without loss of sensitivity. As a result, up to 8,820 individual pre-amplified samples could be screened for SARS-CoV-2 within each 96-well plate of realtime PCR using the pooled testing procedure. CONCLUSION: The novel semi-nested RT-PCR assay provides a solution for highly multiplex (7-plex) detection of SARS-CoV-2 and enables 96-sample pooled detection for increase of testing capacity. .

Establishment of Recombinase Polymerase Amplification assay for rapid and sensitive detection of Orientia tsutsugamushi in Southeast Asia.

Scrub typhus, caused by Orientia tsutsugamushi, is a common fever in parts of Southern and Southeast Asia. As delayed diagnosis of scrub typhus leads to inappropriate treatment and high mortality rates, of up to 70%, sensitive and rapid detection of O. tsutsugamushi is required for timely and appropriate treatment. Molecular assays, such as PCR and real-time PCR, have been shown to be more sensitive than conventional immunoassay, however, they are only available in centralized laboratories. In contrast to PCR assays, Recombinase Polymerase Amplification (RPA) is conducted under a constant temperature ranging from 24°C to 45°C. Therefore, this technology is very promising for nucleic acid testing in the field, and in resource-limited areas. An RPA assay for the detection of O. tsutsugamushi based on the target gene encoding for the 47 kDa outer membrane protein has been reported, but the primer and probe sequences of this assay are suboptimal for detection of the majority of recently published sequences of O. tsutsugamushi isolates from Southeast Asia. We have established a real-time RPA assay with primer and probe sequences that are optimized for most Southeast Asia's isolates of O. tsutsugamushi. As a result, the new RPA assay showed better performance than the previous assay in detecting O. tsutsugamushi in clinical samples of scrub typhus cases found in Vietnam. The specificity of RPA assay was also evaluated using genomic DNA from microorganisms commonly encountered in the differential diagnosis of scrub typhus, and blood samples from healthy controls and O. tsutsugamushi negative confirmed cases.

MAPK inhibitors induce serine peptidase inhibitor Kazal type 1 (SPINK1) secretion in BRAF V600E-mutant colorectal adenocarcinoma.

The mitogen-activated protein kinase (MAPK) pathway plays a central role in colorectal cancers (CRC). In particular, BRAF V600E-mutant tumors, which represent around 10% of CRCs, are refractory to current therapies. Overexpression and secretion of serine peptidase inhibitor Kazal type 1 (SPINK1) are observed in around 50% of CRCs, and its serum level can be used as a biomarker for poor prognosis. Utilizing a recently developed extendable blocking probe assay, we analyzed the BRAF mutation status in a CRC patient cohort (N = 571) using tissue-derived RNA as the starting material. From the same RNA samples, we measured the relative SPINK1 expression levels using a quantitative real-time PCR method. Expression of mutant BRAF V600E correlated with poor prognosis, as did low expression of SPINK1 mRNA. Further, BRAF V600E correlated negatively with SPINK1 levels. In order to investigate the effect of MAPK pathway-targeted therapies on SPINK1 secretion, we conducted in vitro studies using both wild-type and V600E CRC cell lines. BRAF inhibitor vemurafenib, and subsequent MAPK pathway inhibitors trametinib and SCH772984, significantly increased SPINK1 secretion in V600E CRC cell lines Colo205 and HT-29 with a concomitant decrease in trypsin-1 and -2 secretion. Notably, no SPINK1 increase or trypsin-1 decrease was observed in BRAF wild-type CRC cell line Caco-2 in response to MAPK pathway inhibitors. In further mechanistic studies, we observed that only trametinib was able to diminish completely both MEK and ERK phosphorylation in the V600E CRC cells. Furthermore, the key regulator of integrated stress response, activating transcription factor 4 (ATF-4), was downregulated both at mRNA and at protein level in response to trametinib treatment. In conclusion, these data suggest that sustained inhibition of not only MAPK pathway activation, but also ATF-4 and trypsin, might be beneficial in the therapy of BRAF V600E-mutant CRC and that SPINK1 levels may serve as an indicator of therapy response.

No Tissue Expression of KRAS or BRAF Mutations in 61 Adult Patients Treated for Esophageal Atresia in Early Childhood.

BACKGROUND: Previous studies have reported an association among esophageal atresia (EA), Barrett's esophagus, and esophageal adenocarcinoma later in life. OBJECTIVE: The objective of the article is to evaluate KRAS and BRAF mutations as potential genetic markers for early detection of malignant transformation, we used an ultrasensitive technique to detect tissue expression of KRAS and BRAF mutations in endoscopic biopsies from 61 adult patients under follow-up after treatment for EA. MATERIALS AND METHODS: RNA was extracted from 112 fresh-frozen endoscopic tissue biopsies from 61 adult patients treated for EA in early childhood. RNA was reverse transcribed using the extendable blocking probe reverse transcription method. KRAS codons 12 and 13, as well as BRAF mutations were detected by quantitative polymerase chain reaction. RESULTS: No mutations of KRAS codon 12, KRAS codon 13, or BRAF were found in 112 endoscopic biopsy samples from 61 patients. CONCLUSION: Despite the presence of histological findings indicating long-standing gastroesophageal reflux in 25%, as well as symptomatic gastroesophageal reflux in more than 40%, there was no detectable tissue expression of KRAS or BRAF mutations in this cohort of patients.