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.

A study of genetic variants associated with skin traits in the Vietnamese population.

BACKGROUND: Most skin-related traits have been studied in Caucasian genetic backgrounds. A comprehensive study on skin-associated genetic effects on underrepresented populations such as Vietnam is needed to fill the gaps in the field. OBJECTIVES: We aimed to develop a computational pipeline to predict the effect of genetic factors on skin traits using public data (GWAS catalogs and whole-genome sequencing (WGS) data from the 1000 Genomes Project-1KGP) and in-house Vietnamese data (WGS and genotyping by SNP array). Also, we compared the genetic predispositions of 25 skin-related traits of Vietnamese population to others to acquire population-specific insights regarding skin health. METHODS: Vietnamese cohorts of whole-genome sequencing (WGS) of 1008 healthy individuals for the reference and 96 genotyping samples (which do not have any skin cutaneous issues) by Infinium Asian Screening Array-24 v1.0 BeadChip were employed to predict skin-associated genetic variants of 25 skin-related and micronutrient requirement traits in population analysis and correlation analysis. Simultaneously, we compared the landscape of cutaneous issues of Vietnamese people with other populations by assessing their genetic profiles. RESULTS: The skin-related genetic profile of Vietnamese cohorts was similar at most to East Asian cohorts (JPT: Fst = 0.036, CHB: Fst = 0.031, CHS: Fst = 0.027, CDX: Fst = 0.025) in the population study. In addition, we identified pairs of skin traits at high risk of frequent co-occurrence (such as skin aging and wrinkles (r = 0.45, p = 1.50e-5) or collagen degradation and moisturizing (r = 0.35, p = 1.1e-3)). CONCLUSION: This is the first investigation in Vietnam to explore genetic variants of facial skin. These findings could improve inadequate skin-related genetic diversity in the currently published database.

LmTag: functional-enrichment and imputation-aware tag SNP selection for population-specific genotyping arrays.

Despite the rapid development of sequencing technology, single-nucleotide polymorphism (SNP) arrays are still the most cost-effective genotyping solutions for large-scale genomic research and applications. Recent years have witnessed the rapid development of numerous genotyping platforms of different sizes and designs, but population-specific platforms are still lacking, especially for those in developing countries. SNP arrays designed for these countries should be cost-effective (small size), yet incorporate key information needed to associate genotypes with traits. A key design principle for most current platforms is to improve genome-wide imputation so that more SNPs not included in the array (imputed SNPs) can be predicted. However, current tag SNP selection methods mostly focus on imputation accuracy and coverage, but not the functional content of the array. It is those functional SNPs that are most likely associated with traits. Here, we propose LmTag, a novel method for tag SNP selection that not only improves imputation performance but also prioritizes highly functional SNP markers. We apply LmTag on a wide range of populations using both public and in-house whole-genome sequencing databases. Our results show that LmTag improved both functional marker prioritization and genome-wide imputation accuracy compared to existing methods. This novel approach could contribute to the next generation genotyping arrays that provide excellent imputation capability as well as facilitate array-based functional genetic studies. Such arrays are particularly suitable for under-represented populations in developing countries or non-model species, where little genomics data are available while investment in genome sequencing or high-density SNP arrays is limited. $\textrm{LmTag}$ is available at: https://github.com/datngu/LmTag.

Assessing polygenic risk score models for applications in populations with under-represented genomics data: an example of Vietnam.

Most polygenic risk score (PRS)models have been based on data from populations of European origins (accounting for the majority of the large genomics datasets, e.g. >78% in the UK Biobank and >85% in the GTEx project). Although several large-scale Asian biobanks were initiated (e.g. Japanese, Korean, Han Chinese biobanks), most other Asian countries have little or near-zero genomics data. To implement PRS models for under-represented populations, we explored transfer learning approaches, assuming that information from existing large datasets can compensate for the small sample size that can be feasibly obtained in developing countries, like Vietnam. Here, we benchmark 13 common PRS methods in meta-population strategy (combining individual genotype data from multiple populations) and multi-population strategy (combining summary statistics from multiple populations). Our results highlight the complementarity of different populations and the choice of methods should depend on the target population. Based on these results, we discussed a set of guidelines to help users select the best method for their datasets. We developed a robust and comprehensive software to allow for benchmarking comparisons between methods and proposed a computational framework for improving PRS performance in a dataset with a small sample size. This work is expected to inform the development of genomics applications in under-represented populations. PRSUP framework is available at: https://github.com/BiomedicalMachineLearning/VGP.

Investigation of target sequencing of SARS-CoV-2 and immunogenic GWAS profiling in host cells of COVID-19 in Vietnam.

BACKGROUND: A global pandemic has been declared for coronavirus disease 2019 (COVID-19), which has serious impacts on human health and healthcare systems in the affected areas, including Vietnam. None of the previous studies have a framework to provide summary statistics of the virus variants and assess the severity associated with virus proteins and host cells in COVID-19 patients in Vietnam. METHOD: In this paper, we comprehensively investigated SARS-CoV-2 variants and immune responses in COVID-19 patients. We provided summary statistics of target sequences of SARS-CoV-2 in Vietnam and other countries for data scientists to use in downstream analysis for therapeutic targets. For host cells, we proposed a predictive model of the severity of COVID-19 based on public datasets of hospitalization status in Vietnam, incorporating a polygenic risk score. This score uses immunogenic SNP biomarkers as indicators of COVID-19 severity. RESULT: We identified that the Delta variant of SARS-CoV-2 is most prevalent in southern areas of Vietnam and it is different from other areas in the world using various data sources. Our predictive models of COVID-19 severity had high accuracy (Random Forest AUC = 0.81, Elastic Net AUC = 0.7, and SVM AUC = 0.69) and showed that the use of polygenic risk scores increased the models' predictive capabilities. CONCLUSION: We provided a comprehensive analysis for COVID-19 severity in Vietnam. This investigation is not only helpful for COVID-19 treatment in therapeutic target studies, but also could influence further research on the disease progression and personalized clinical outcomes.

Leveraging known genomic variants to improve detection of variants, especially close-by Indels.

Motivation: The detection of genomic variants has great significance in genomics, bioinformatics, biomedical research and its applications. However, despite a lot of effort, Indels and structural variants are still under-characterized compared to SNPs. Current approaches based on next-generation sequencing data usually require large numbers of reads (high coverage) to be able to detect such types of variants accurately. However Indels, especially those close to each other, are still hard to detect accurately. Results: We introduce a novel approach that leverages known variant information, e.g. provided by dbSNP, dbVar, ExAC or the 1000 Genomes Project, to improve sensitivity of detecting variants, especially close-by Indels. In our approach, the standard reference genome and the known variants are combined to build a meta-reference, which is expected to be probabilistically closer to the subject genomes than the standard reference. An alignment algorithm, which can take into account known variant information, is developed to accurately align reads to the meta-reference. This strategy resulted in accurate alignment and variant calling even with low coverage data. We showed that compared to popular methods such as GATK and SAMtools, our method significantly improves the sensitivity of detecting variants, especially Indels that are close to each other. In particular, our method was able to call these close-by Indels at a 15-20% higher sensitivity than other methods at low coverage, and still get 1-5% higher sensitivity at high coverage, at competitive precision. These results were validated using simulated data with variant profiles extracted from the 1000 Genomes Project data, and real data from the Illumina Platinum Genomes Project and ExAC database. Our finding suggests that by incorporating known variant information in an appropriate manner, sensitive variant calling is possible at a low cost. Availability and implementation: Implementation can be found in our public code repository https://github.com/namsyvo/IVC. Supplementary information: Supplementary data are available at Bioinformatics online.

How genome complexity can explain the difficulty of aligning reads to genomes.

BACKGROUND: Although it is frequently observed that aligning short reads to genomes becomes harder if they contain complex repeat patterns, there has not been much effort to quantify the relationship between complexity of genomes and difficulty of short-read alignment. Existing measures of sequence complexity seem unsuitable for the understanding and quantification of this relationship. RESULTS: We investigated several measures of complexity and found that length-sensitive measures of complexity had the highest correlation to accuracy of alignment. In particular, the rate of distinct substrings of length k, where k is similar to the read length, correlated very highly to alignment performance in terms of precision and recall. We showed how to compute this measure efficiently in linear time, making it useful in practice to estimate quickly the difficulty of alignment for new genomes without having to align reads to them first. We showed how the length-sensitive measures could provide additional information for choosing aligners that would align consistently accurately on new genomes. CONCLUSIONS: We formally established a connection between genome complexity and the accuracy of short-read aligners. The relationship between genome complexity and alignment accuracy provides additional useful information for selecting suitable aligners for new genomes. Further, this work suggests that the complexity of genomes sometimes should be thought of in terms of specific computational problems, such as the alignment of short reads to genomes.

Exploiting dependencies of pairwise comparison outcomes to predict patterns of gene response.

BACKGROUND: The analysis of gene expression has played an important role in medical and bioinformatics research. Although it is known that a large number of samples is needed to determine the patterns of gene expression accurately, practical designs of gene expression studies occasionally have insufficient numbers of samples, making it difficult to ascertain true response patterns of variantly expressed genes. RESULTS: We describe an approach to cope with the challenge of predicting true orders of gene response to treatments. We show that true patterns of gene response must be orderable sets. In experiments with few samples, we modify the conventional pairwise comparison tests and increase the significance level α intelligently to deduce orderable patterns, which are most likely true orders of gene response. Additionally, motivated by the fact that a gene can be involved in multiple biological functions, our method further resamples experimental replicates and predicts multiple response patterns for each gene. CONCLUSIONS: This method can be useful in designing cost-effective experiments with small sample sizes. Patterns of highly-variantly expressed genes can be predicted by varying α intelligently. Furthermore, clusters are labeled meaningfully with patterns that describe precisely how genes in such clusters respond to treatments.

RandAL: a randomized approach to aligning DNA sequences to reference genomes.

BACKGROUND: The alignment of short reads generated by next-generation sequencers to genomes is an important problem in many biomedical and bioinformatics applications. Although many proposed methods work very well on narrow ranges of read lengths, they tend to suffer in performance and alignment quality for reads outside of these ranges. RESULTS: We introduce RandAL, a novel method that aligns DNA sequences to reference genomes. Our approach utilizes two FM indices to facilitate efficient bidirectional searching, a pruning heuristic to speed up the computing of edit distances, and most importantly, a randomized strategy that enables effective estimation of key parameters. Extensive comparisons showed that RandAL outperformed popular aligners in most instances and was unique in its consistent and accurate performance over a wide range of read lengths and error rates. The software package is publicly available at https://github.com/namsyvo/RandAL. CONCLUSIONS: RandAL promises to align effectively and accurately short reads that come from a variety of technologies with different read lengths and rates of sequencing error.

Whole-Genome Sequencing Reveals Temporal Trends in Antibiotic Resistance Genes in Escherichia coli Causing Pediatric Urinary Tract Infections in Central Vietnam.

(1) Background: Pediatric urinary tract infections (UTIs) pose significant challenges due to drug-resistant Escherichia coli (E. coli) strains. This study utilizes whole-genome sequencing to analyze temporal trends in antibiotic resistance genes (ARGs) in clinical E. coli isolates from pediatric UTI cases in central Vietnam. (2) Methods: We conducted whole-genome sequencing on 71 E. coli isolates collected from pediatric UTI patients between 2018 and 2020. ARGs were identified, and their prevalence over time was analyzed. Statistical tests were used to correlate ARG presence with antibiotic resistance. (3) Results: Of the 47 E. coli isolates with complete data, 40 distinct ARGs were identified, with a median of 10 resistance genes per isolate. A significant increase in the total number of ARGs per isolate was observed over time, from an average of 8.88 before June 2019 to 11.63 after. Notably, the prevalence of the aadA2 gene (aminoglycoside resistance) rose from 0% to 26.7%, and that of the blaNDM-5 gene (beta-lactam and carbapenem resistance) increased from 0% to 23.3%. Key correlations include blaEC with cephalosporin resistance, blaNDM-5 with carbapenem resistance, and sul2 with sulfamethoxazole/trimethoprim resistance. (4) Conclusions: Whole-genome sequencing reveals complex and evolving antibiotic resistance patterns in pediatric E. coli UTIs in central Vietnam, with a marked increase in ARG prevalence over time. Continuous surveillance and targeted treatments are essential to address these trends. Understanding genetic foundations is crucial for effective intervention strategies.

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.

The therapeutic landscape for COVID-19 and post-COVID-19 medications from genetic profiling of the Vietnamese population and a predictive model of drug-drug interaction for comorbid COVID-19 patients.

Despite the raised awareness of the role of pharmacogenomic (PGx) in personalized medicines for COVID-19, data for COVID-19 drugs is extremely scarce and not even a publication on this topic for post-COVID-19 medications to date. In the current study, we investigated the genetic variations associated with COVID-19 and post-COVID-19 therapies by using whole genome sequencing data of the 1000 Vietnamese Genomes Project (1KVG) in comparison with other populations retrieved from the 1000 Genomes Project Phase 3 (1KGP3) and the Genome Aggregation Database (gnomAD). Moreover, we also evaluated the risk of drug interactions in comorbid COVID-19 and post-COVID-19 patients based on pharmacogenomic profiles of drugs using a computational approach. For COVID-19 therapies, variants related to the response of two causal treatment agents (tolicizumab and ritonavir) and antithrombotic drugs are common in the Vietnamese cohort. Regarding post-COVID-19, drugs for mental manipulations possess the highest number of clinical annotated variants carried by Vietnamese individuals. Among the superpopulations, East Asian populations shared the most similar genetic structure with the Vietnamese population, whereas the African population showed the most difference. Comorbid patients are at an increased drug-drug interaction (DDI) risk when suffering from COVID-19 and after recovering as well due to a large number of potential DDIs which have been identified. Our results presented the population-specific understanding of the pharmacogenomic aspect of COVID-19 and post-COVID-19 therapy to optimize therapeutic outcomes and promote personalized medicine strategy. We also partly clarified the higher risk in COVID-19 patients with underlying conditions by assessing the potential drug interactions.

A comprehensive evaluation of polygenic score and genotype imputation performances of human SNP arrays in diverse populations.

Regardless of the overwhelming use of next-generation sequencing technologies, microarray-based genotyping combined with the imputation of untyped variants remains a cost-effective means to interrogate genetic variations across the human genome. This technology is widely used in genome-wide association studies (GWAS) at bio-bank scales, and more recently, in polygenic score (PGS) analysis to predict and stratify disease risk. Over the last decade, human genotyping arrays have undergone a tremendous growth in both number and content making a comprehensive evaluation of their performances became more important. Here, we performed a comprehensive performance assessment for 23 available human genotyping arrays in 6 ancestry groups using diverse public and in-house datasets. The analyses focus on performance estimation of derived imputation (in terms of accuracy and coverage) and PGS (in terms of concordance to PGS estimated from whole-genome sequencing data) in three different traits and diseases. We found that the arrays with a higher number of SNPs are not necessarily the ones with higher imputation performance, but the arrays that are well-optimized for the targeted population could provide very good imputation performance. In addition, PGS estimated by imputed SNP array data is highly correlated to PGS estimated by whole-genome sequencing data in most cases. When optimal arrays are used, the correlations of PGS between two types of data are higher than 0.97, but interestingly, arrays with high density can result in lower PGS performance. Our results suggest the importance of properly selecting a suitable genotyping array for PGS applications. Finally, we developed a web tool that provides interactive analyses of tag SNP contents and imputation performance based on population and genomic regions of interest. This study would act as a practical guide for researchers to design their genotyping arrays-based studies. The tool is available at: https://genome.vinbigdata.org/tools/saa/ .

BRCA1/2 Mutations in Vietnamese Patients with Hereditary Breast and Ovarian Cancer Syndrome.

(1) Background: Individuals with BRCA1/2 gene mutations are at increased risk of breast and ovarian cancer. The prevalence of BRCA1/2 mutations varies by race and ethnicity, and the prevalence and the risks associated with most BRCA1/2 mutations has not been unknown in the Vietnamese population. We herein screen the entire BRCA1 and BRCA2 genes for breast and ovarian cancer patients with a family history of breast cancer and ovarian cancer, thereby, suggesting a risk score associated with carrier status and history for aiding personalized treatment; (2) Methods: Between December 2017 and December 2019, Vietnamese patients who had a pathological diagnosis of breast and epithelial ovarian cancer were followed up, prospectively, after treatment from two large institutions in Vietnam. Blood samples from 33 Vietnamese patients with hereditary breast and ovarian cancers (HBOC) syndrome were collected and analyzed using Next Generation Sequencing; (3) Results: Eleven types of mutations in both BRCA1 (in nine patients) and BRCA2 (in three patients) were detected, two of which (BRCA1:p.Tyr1666Ter and BRCA2:p.Ser1341Ter) have not been previously documented in the literature. Seven out of 19 patient's relatives had BRCA1/2 gene mutations. All selected patients were counselled about the likelihood of cancer rising and prophylactic screening and procedures. The study established a risk score associated with the cohorts based on carrier status and family history; (4) Conclusions: Our findings suggested the implications for the planning of a screening programme for BRCA1 and BRCA2 genes testing in breast and ovarian cancer patients and genetic screening in their relatives. BRCA1/2 mutation carriers without cancer should have early and regular cancer screening, and prophylactic measures. This study could be beneficial for a diverse group in a large population-specific cohort, related to HBOC Syndrome.

Review on Databases and Bioinformatic Approaches on Pharmacogenomics of Adverse Drug Reactions.

Pharmacogenomics has been used effectively in studying adverse drug reactions by determining the person-specific genetic factors associated with individual response to a drug. Current approaches have revealed the significant importance of sequencing technologies and sequence analysis strategies for interpreting the contribution of genetic variation in developing adverse reactions. Advance in next generation sequencing and platform brings new opportunities in validating the genetic candidates in certain reactions, and could be used to develop the preemptive tests to predict the outcome of the variation in a personal response to a drug. With the highly accumulated available data recently, the in silico approach with data analysis and modeling plays as other important alternatives which significantly support the final decisions in the transformation from research to clinical applications such as diagnosis and treatments for various types of adverse responses.