An overview of next generation sequencing strategies and genomics tools used for tuberculosis research.

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research.

Distinct oral DNA viral signatures in rheumatoid arthritis: a Pilot study.

Background: Despite evidence linking viruses and oral microbiome to rheumatoid arthritis (RA), limited whole genome sequencing research has been conducted on the oral virome (a viral component of the microbiome) of untreated RA patients. This pilot research seeks to address this knowledge gap by comparing the oral virome of untreated rheumatoid arthritis patients (RAs) and healthy individuals (HCs). Method: Whole genome DNA sequence of saliva samples from 45 participants including 21 RAs and 24 age and gender matched HCs was obtained from the BioProject: PRJEB6997. Metaphlan3 pipeline and LEfSe analysis were used for the viral signature detection. Wilcoxon pairwise test and ROC analysis were used to validate and predict signatures. Results: RA exhibits higher alpha diversity compared to HCs. Callitrichine gammaherpesvirus 3, Human gammaherpesvirus 4 (EBV), Murid betaherpesvirus 8, and Suid alphaherpesvirus 1 were enriched in RAs, while Aotine betaherpesvirus 1 from the Cytomegalovirus genus was enriched in HCs. In addition, Saccharomyces cerevisiae killer virus M1 (ScV-M1) was found to be enriched in RAs, whereas bacteriophage Hk97virus (Siphoviridae) and Cd119virus (Myoviridae) were enriched in HCs. Conclusion: This study identifies significant DNA oral viral signatures at species level as potential biomarkers for the early detection and diagnosis of rheumatoid arthritis.

Whole Genome-Sequencing and Phylogenetic Analysis of Bacillus anthracis from 2019-2023 in Ningxia Hui Autonomous Region, China.

BACKGROUND: Since 2019, the incidence of anthrax in the Ningxia Hui Autonomous Region has increased significantly compared with previous years, so in this situation the anthrax in the Ningxia region not only had a detrimental impact on public health, but also inflicted significant economic repercussions. Therefore, we conducted a molecular epidemiological study of 20 strains from 2019-2023 isolates. This study investigated the origin of Bacillus anthracis and its genetic diversity. METHODS: We conducted canonical single-nucleotide polymorphisms (CanSNPs) typing and whole genome sequencing based on the extracted nucleic acid of Bacillus anthracis. Based on the whole genome drafts, we studied the genomic characteristics of 20 isolates. Meanwhile, we performed phylogenetic studies based on genome-wide core single-nucleotide polymorphisms (SNPs) using MEGA's Maximum Likelihood (ML) method and core-genome-based multilocus sequence typing (cgMLST) of the core genomes of these strains using BioNumerics' minimum spanning tree (MST) model. RESULTS: The 20 isolates were categorized into sub-lineages A.Br.001/002, and comparative genomic analyses of these strains with other isolates from other parts of the world showed that the strains from Ningxia were correlated with isolates from Europe, Indonesia, Georgia (USA), and Beijing (China). For the 20 isolates in Ningxia, the genetic relationship of the isolates isolated from the same year or region was relatively close. CONCLUSION: The A.Br.001/002 subgroup was the dominant endemic strain in Ningxia. The genetic relationship and phylogenesis between isolates from Ningxia and strains from Europe and Indonesia suggest that anthrax spread around the globe through ancient trade routes.

Whole-genome informed circulating tumor DNA analysis by multiplex digital PCR for disease monitoring in B-cell lymphomas: a proof-of-concept study.

Introduction: Analyzing liquid biopsies for tumor-specific aberrations can facilitate detection of measurable residual disease (MRD) during treatment and at follow-up. In this study, we assessed the clinical potential of using whole-genome sequencing (WGS) of lymphomas at diagnosis to identify patient-specific structural (SVs) and single nucleotide variants (SNVs) to enable longitudinal, multi-targeted droplet digital PCR analysis (ddPCR) of cell-free DNA (cfDNA). Methods: In 9 patients with B-cell lymphoma (diffuse large B-cell lymphoma and follicular lymphoma), comprehensive genomic profiling at diagnosis was performed by 30X WGS of paired tumor and normal specimens. Patient-specific multiplex ddPCR (m-ddPCR) assays were designed for simultaneous detection of multiple SNVs, indels and/or SVs, with a detection sensitivity of 0.0025% for SV assays and 0.02% for SNVs/indel assays. M-ddPCR was applied to analyze cfDNA isolated from serially collected plasma at clinically critical timepoints during primary and/or relapse treatment and at follow-up. Results: A total of 164 SNVs/indels were identified by WGS including 30 variants known to be functionally relevant in lymphoma pathogenesis. The most frequently mutated genes included KMT2D, PIM1, SOCS1 and BCL2. WGS analysis further identified recurrent SVs including t(14;18)(q32;q21) (IGH::BCL2), and t(6;14)(p25;q32) (IGH::IRF4). Plasma analysis at diagnosis showed positive circulating tumor DNA (ctDNA) levels in 88% of patients and the ctDNA burden correlated with baseline clinical parameters (LDH and sedimentation rate, p-value <0.01). While clearance of ctDNA levels after primary treatment cycle 1 was observed in 3/6 patients, all patients analyzed at final evaluation of primary treatment showed negative ctDNA, hence correlating with PET-CT imaging. One patient with positive ctDNA at interim also displayed detectable ctDNA (average variant allele frequency (VAF) 6.9%) in the follow-up plasma sample collected 2 years after final evaluation of primary treatment and 25 weeks before clinical manifestation of relapse. Conclusion: In summary, we demonstrate that multi-targeted cfDNA analysis, using a combination of SNVs/indels and SVs candidates identified by WGS analysis, provides a sensitive tool for MRD monitoring and can detect lymphoma relapse earlier than clinical manifestation.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton.

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Mining the genome of Bacillus velezensis FS26 for probiotic markers and secondary metabolites with antimicrobial properties against aquaculture pathogens.

Bacillus velezensis FS26 is a bacterium from the genus Bacillus that has been proven as a potential probiotic in aquaculture with a good antagonistic effect on Aeromonas spp. and Vibrio spp. Whole-genome sequencing (WGS) allows a comprehensive and in-depth analysis at the molecular level, and it is becoming an increasingly significant technique in aquaculture research. Although numerous probiotic genomes have been sequenced and investigated recently, there are minimal data on in silico analysis of B. velezensis as a probiotic bacterium isolated from aquaculture sources. Thus, this study aims to analyse the general genome characteristics and probiotic markers from the B. velezensis FS26 genome with secondary metabolites predicted against aquaculture pathogens. The B. velezensis FS26 genome (GenBank Accession: JAOPEO000000000) assembly proved to be of high quality, with eight contigs containing 3,926,371 bp and an average G + C content of 46.5%. According to antiSMASH analysis, five clusters of secondary metabolites from the B. velezensis FS26 genome showed 100% similarity. These clusters include Cluster 2 (bacilysin), Cluster 6 (bacillibactin), Cluster 7 (fengycin), Cluster 8 (bacillaene), and Cluster 9 (macrolactin H), which signify promising antibacterial, antifungal, and anticyanobacterial agents against pathogens in aquaculture. The probiotic markers of B. velezensis FS26 genome for adhesion capability in the hosts' intestine, as well as the acid and bile salt-tolerant genes, were also detected through the Prokaryotic Genome Annotation System (Prokka) annotation pipeline. These results are in agreement with our previous in vitro data, suggesting that the in silico investigation facilitates establishing B. velezensis FS26 as a beneficial probiotic for use in aquaculture.

Unveiling the Human Brain Virome in Brodmann Area 46: Novel Insights Into Dysbiosis and Its Association With Schizophrenia.

Research suggests a potential role of the oral-neuro and gut-brain axes in schizophrenia, involving non-brain microbiomes such as salivary and gut microbiomes. However, the blood-brain barrier effectively prevents microorganism entry. Additionally, despite approximately 8% of the human genome consisting of retroviruses and the established link between viral infections and schizophrenia, the presence of a resident virome (a viral component of the microbiome) in the brain and its association with mental disorders remain unexplored. METHODS: Whole-genome sequencing raw data from postmortem Brodmann Area 46 (BA46) tissue from 49 individuals (20 healthy controls [HCs], 29 with schizophrenia [SCZs]) obtained from the NCBI SRA database from BioProject: PRJNA422380.Virome profiles were retrieved using Metaphlan3, and viral signatures were identified using linear discriminant analysis effect size (LEfSe). Mann-Whitney tests and receiver operating characteristic curve validated the viral signatures. RESULTS: In BA46, 30 distinct species representing 9 phyla, 10 classes, 10 orders, 13 families, and 19 genera were identified. HCs exhibited greater alpha diversity, and there were significant differences in beta diversity between the groups. LEfSe analysis highlighted distinct viral levels, including Escherichia virus Lambda, Escherichia virus phiV10, Human endogenous retrovirus K, Taterapox virus, Alcelaphine gammaherpesvirus 1, and Bovine gammaherpesvirus 4 in HCs, while Glypta fumiferanae ichnovirus and unknown virus showed higher levels in schizophrenia. CONCLUSION: This is the first study to identify a human brain virome associated with schizophrenia in BA46. Brain virome dysbiosis may be associated with mental illness, and viral signatures may serve as biomarkers for the early detection of schizophrenia.

Favored single nucleotide variants identified using whole genome Re-sequencing of Austrian and Chinese cattle breeds.

Cattle have been essential for the development of human civilization since their first domestication few thousand years ago. Since then, they have spread across vast geographic areas following human activities. Throughout generations, the cattle genome has been shaped with detectable signals induced by various evolutionary processes, such as natural and human selection processes and demographic events. Identifying such signals, called selection signatures, is one of the primary goals of population genetics. Previous studies used various selection signature methods and normalized the outputs score using specific windows, in kbp or based on the number of SNPs, to identify the candidate regions. The recent method of iSAFE claimed for high accuracy in pinpointing the candidate SNPs. In this study, we analyzed whole-genome resequencing (WGS) data of ten individuals from Austrian Fleckvieh (Bos taurus) and fifty individuals from 14 Chinese indigenous breeds (Bos taurus, Bos taurus indicus, and admixed). Individual WGS reads were aligned to the cattle reference genome of ARS. UCD1.2 and subsequently undergone single nucleotide variants (SNVs) calling pipeline using GATK. Using these SNVs, we examined the population structure using principal component and admixture analysis. Then we refined selection signature candidates using the iSAFE program and compared it with the classical iHS approach. Additionally, we run Fst population differentiation from these two cattle groups. We found gradual changes of taurine in north China to admixed and indicine to the south. Based on the population structure and the number of individuals, we grouped samples to Fleckvieh, three Chinese taurines (Kazakh, Mongolian, Yanbian), admixed individuals (CHBI_Med), indicine individuals (CHBI_Low), and a combination of admixed and indicine (CHBI) for performing iSAFE and iHS tests. There were more significant SNVs identified using iSAFE than the iHS for the candidate of positive selection and more detectable signals in taurine than in indicine individuals. However, combining admixed and indicine individuals decreased the iSAFE signals. From both within-population tests, significant SNVs are linked to the olfactory receptors, production, reproduction, and temperament traits in taurine cattle, while heat and parasites tolerance in the admixed individuals. Fst test suggests similar patterns of population differentiation between Fleckvieh and three Chinese taurine breeds against CHBI. Nevertheless, there are genes shared only among the Chinese taurine, such as PAX5, affecting coat color, which might drive the differences between these yellowish coated breeds, and those in the greater Far East region.

Clarification of relationship between single-nucleotide polymorphism panels of Shiga toxin-producing Escherichia coli O157:H7/H- strains.

Eighty strains of enterohemorrhagic Escherichia coli O157:H7/H- were analyzed by three single-nucleotide polymorphism (SNP) panels using whole-genome sequencing data. The partial concordance of SNP types among the different SNP panels was observed on minimum spanning trees reconstructed with SNP data. As for lineage I/II strains, some of the clade 7 strains belonged to one unique SNP type as determined by three panels, suggesting that clade 7 should be divided into at least two genotypes, namely, the unique type and the rest. In addition, clade 8 contained two unique genotypes, which was consistent with the previous prediction. Similarly, for lineage II, clade 12 should be divided into three genotype strains. In contrast, many strains of several clades belonging to lineage I were clustered into the same node on each minimum spanning tree upon testing with the three SNP panels. Previous studies reported that lineage I diverged more recently than lineages I/II and II. Such low diversity in lineage I in this study may have arisen because this lineage has not accumulated SNPs because of its relatively recent divergence. Based on the concordance observed in this study, some of the previously published O157 genotype distribution data were successfully interpreted to clarify the clade distribution, which was well supported by previous literature.

Clonal serotype 1c multidrug-resistant Shigella flexneri detected in multiple institutions by sentinel-site sequencing.

Shigella flexneri is a major diarrhoeal pathogen, and the emergence of multidrug-resistant S. flexneri is of public health concern. We report the detection of a clonal cluster of multidrug-resistant serotype 1c (7a) S. flexneri in Singapore in April 2022. Long-read whole-genome sequence analysis found five S. flexneri isolates to be clonal and harboring the extended-spectrum ß-lactamases bla CTX-M-15 and bla TEM-1. The isolates were phenotypically resistant to ceftriaxone and had intermediate susceptibility to ciprofloxacin. The S. flexneri clonal cluster was first detected in a tertiary hospital diagnostic laboratory (sentinel-site), to which the S. flexneri isolates were sent from other hospitals for routine serogrouping. Long-read whole-genome sequence analysis was performed in the sentinel-site near real-time in view of the unusually high number of S. flexneri isolates received within a short time frame. This study demonstrates that near real-time sentinel-site sequence-based surveillance of convenience samples can detect possible clonal outbreak clusters and may provide alerts useful for public health mitigations at the earliest possible opportunity.

Whole-Genome Sequencing and Genomic Variant Analysis of Kazakh Individuals.

Kazakhstan, the ninth-largest country in the world, is located along the Great Silk Road and connects Europe with Asia. Historically, its territory has been inhabited by nomadic tribes, and modern-day Kazakhstan is a multiethnic country with a dominant Kazakh population. We sequenced and analyzed the genomes of five ethnic Kazakhs at high coverage using the Illumina HiSeq2000 next-generation sequencing platform. The five Kazakhs yielded a total number of base pairs ranging from 87,308,581,400 to 107,526,741,301. On average, 99.06% were properly mapped. Based on the Het/Hom and Ti/Tv ratios, the quality of the genomic data ranged from 1.35 to 1.49 and from 2.07 to 2.08, respectively. Genetic variants were identified and annotated. Functional analysis of the genetic variants identified several variants that were associated with higher risks of metabolic and neurogenerative diseases. The present study showed high levels of genetic admixture of Kazakhs that were comparable to those of other Central Asians. These whole-genome sequence data of healthy Kazakhs could contribute significantly to biomedical studies of common diseases as their findings could allow better insight into the genotype-phenotype relations at the population level.

First outbreak of multidrug-resistant tuberculosis (MDR-TB) in Denmark involving six Danish-born cases.

BACKGROUND: Denmark is a low-incidence country for tuberculosis (TB) and multidrug-resistant (MDR) TB at 5 and 0.05 cases per 100,000 population, respectively. Until 2018, the transmission of MDR-TB was nonexistent except for a few pairwise related family cases. In this study, we describe the first MDR-TB outbreak in Denmark. METHODS: On the basis of genotyping of all Mycobacterium tuberculosis (Mtb) culture-positive cases in Denmark spanning 3 decades, 6 molecular- and epidemiologically linked Danish-born cases were identified as the first cluster of an MDR-TB in Denmark. The primary case was diagnosed posthumously in 2010 followed by 5 epidemiologically linked cases from 2018 to 2019. RESULTS AND CONCLUSION: Through a combination of routine Mtb genotyping and clinical epidemiological surveillance data, we identified the first Danish MDR-TB outbreak spanning 10 years and were able to disclose the specific transmission pathways in detail, which helped guide the outbreak investigations. The occurrence of an MDR-TB outbreak in a resource-rich low TB incidence setting such as Denmark highlights the importance of a collaborative control system combining classic contact tracing; timely identification of drug-resistant TB through rapid diagnostics; and a close collaboration between clinicians and classical- and molecular epidemiologists for the benefit of TB control.

Telomere length analysis in amyotrophic lateral sclerosis using large-scale whole genome sequence data.

Background: Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by the loss of upper and lower motor neurons, leading to progressive weakness of voluntary muscles, with death following from neuromuscular respiratory failure, typically within 3 to 5 years. There is a strong genetic contribution to ALS risk. In 10% or more, a family history of ALS or frontotemporal dementia is obtained, and the Mendelian genes responsible for ALS in such families have now been identified in about 50% of cases. Only about 14% of apparently sporadic ALS is explained by known genetic variation, suggesting that other forms of genetic variation are important. Telomeres maintain DNA integrity during cellular replication, differ between sexes, and shorten naturally with age. Sex and age are risk factors for ALS and we therefore investigated telomere length in ALS. Methods: Samples were from Project MinE, an international ALS whole genome sequencing consortium that includes phenotype data. For validation we used donated brain samples from motor cortex from people with ALS and controls. Ancestry and relatedness were evaluated by principal components analysis and relationship matrices of DNA microarray data. Whole genome sequence data were from Illumina HiSeq platforms and aligned using the Isaac pipeline. TelSeq was used to quantify telomere length using whole genome sequence data. We tested the association of telomere length with ALS and ALS survival using Cox regression. Results: There were 6,580 whole genome sequences, reducing to 6,195 samples (4,315 from people with ALS and 1,880 controls) after quality control, and 159 brain samples (106 ALS, 53 controls). Accounting for age and sex, there was a 20% (95% CI 14%, 25%) increase of telomere length in people with ALS compared to controls (p = 1.1 × 10-12), validated in the brain samples (p = 0.03). Those with shorter telomeres had a 10% increase in median survival (p = 5.0×10-7). Although there was no difference in telomere length between sporadic ALS and familial ALS (p=0.64), telomere length in 334 people with ALS due to expanded C9orf72 repeats was shorter than in those without expanded C9orf72 repeats (p = 5.0×10-4). Discussion: Although telomeres shorten with age, longer telomeres are a risk factor for ALS and worsen prognosis. Longer telomeres are associated with ALS.

Systematic Evaluation of Whole-Genome Sequencing Based Prediction of Antimicrobial Resistance in Campylobacter jejuni and C. coli.

The increasing prevalence of antimicrobial resistance (AMR) in Campylobacter spp. is a global concern. This study evaluated the use of whole-genome sequencing (WGS) to predict AMR in Campylobacter jejuni and C. coli. A panel of 271 isolates recovered from Canadian poultry was used to compare AMR genotype to antimicrobial susceptibility testing (AST) results (azithromycin, ciprofloxacin, erythromycin, gentamicin, tetracycline, florfenicol, nalidixic acid, telithromycin, and clindamycin). The presence of antibiotic resistance genes (ARGs) was determined for each isolate using five computational approaches to evaluate the effect of: ARG screening software, input data (i.e., raw reads, draft genome assemblies), genome coverage and genome assembly software. Overall, concordance between the genotype and phenotype was influenced by the computational pipelines, level of genome coverage and the type of ARG but not by input data. For example, three of the pipelines showed a 99% agreement between detection of a tet(O) gene and tetracycline resistance, whereas agreement between the detection of tet(O) and TET resistance was 98 and 93% for two pipelines. Overall, higher levels of genome coverage were needed to reliably detect some ARGs; for example, at 15X coverage a tet(O) gene was detected in >70% of the genomes, compared to <60% of the genomes for bla(OXA). No genes associated with florfenicol or gentamicin resistance were found in the set of strains included in this study, consistent with AST results. Macrolide and fluoroquinolone resistance was associated 100% with mutations in the 23S rRNA (A2075G) and gyrA (T86I) genes, respectively. A lower association between a A2075G 23S rRNA gene mutation and resistance to clindamycin and telithromycin (92.8 and 78.6%, respectively) was found. While WGS is an effective approach to predicting AMR in Campylobacter, this study demonstrated the impact that computational pipelines, genome coverage and the genes can have on the reliable identification of an AMR genotype.

Complexity of Genomic Epidemiology of Carbapenem-Resistant Klebsiella pneumoniae Isolates in Colombia Urges the Reinforcement of Whole Genome Sequencing-Based Surveillance Programs.

BACKGROUND: Carbapenem-resistant Klebsiella pneumoniae (CRKP) is an emerging public health problem. This study explores the specifics of CRKP epidemiology in Colombia based on whole genome sequencing (WGS) of the National Reference Laboratory at Instituto Nacional de Salud (INS)'s 2013-2017 sample collection. METHODS: A total of 425 CRKP isolates from 21 departments were analyzed by HiSeq-X10®Illumina high-throughput sequencing. Bioinformatic analysis was performed, primarily using the pipelines developed collaboratively by the National Institute for Health Research Global Health Research Unit (GHRU) on Genomic Surveillance of Antimicrobial Resistance (AMR), and AGROSAVIA. RESULTS: Of the 425 CRKP isolates, 91.5% were carbapenemase-producing strains. The data support a recent expansion and the endemicity of CRKP in Colombia with the circulation of 7 high-risk clones, the most frequent being CG258 (48.39% of isolates). We identified genes encoding carbapenemases blaKPC-3, blaKPC-2, blaNDM-1, blaNDM-9, blaVIM-2, blaVIM-4, and blaVIM-24, and various mobile genetic elements (MGE). The virulence of CRKP isolates was low, but colibactin (clb3) was present in 25.2% of isolates, and a hypervirulent CRKP clone (CG380) was reported for the first time in Colombia. ST258, ST512, and ST4851 were characterized by low levels of diversity in the core genome (ANI > 99.9%). CONCLUSIONS: The study outlines complex CRKP epidemiology in Colombia. CG258 expanded clonally and carries specific carbapenemases in specific MGEs, while the other high-risk clones (CG147, CG307, and CG152) present a more diverse complement of carbapenemases. The specifics of the Colombian situation stress the importance of WGS-based surveillance to monitor evolutionary trends of sequence types (STs), MGE, and resistance and virulence genes.

Copy Number Variation Identification on 3,800 Alzheimer's Disease Whole Genome Sequencing Data from the Alzheimer's Disease Sequencing Project.

Alzheimer's Disease (AD) is a progressive neurologic disease and the most common form of dementia. While the causes of AD are not completely understood, genetics plays a key role in the etiology of AD, and thus finding genetic factors holds the potential to uncover novel AD mechanisms. For this study, we focus on copy number variation (CNV) detection and burden analysis. Leveraging whole-genome sequence (WGS) data released by Alzheimer's Disease Sequencing Project (ADSP), we developed a scalable bioinformatics pipeline to identify CNVs. This pipeline was applied to 1,737 AD cases and 2,063 cognitively normal controls. As a result, we observed 237,306 and 42,767 deletions and duplications, respectively, with an average of 2,255 deletions and 1,820 duplications per subject. The burden tests show that Non-Hispanic-White cases on average have 16 more duplications than controls do (p-value 2e-6), and Hispanic cases have larger deletions than controls do (p-value 6.8e-5).

Whole genome global insight of antibiotic resistance gene repertoire and virulome of high - risk multidrug-resistant Uropathogenic Escherichiacoli.

Elucidation of genetic determinants via whole genome sequence (WGS) analyses can help understand the high risk multidrug-resistant (MDR) Uropathogenic Escherichia coli (UPEC) associated with urinary tract infections (UTI) and its evasion strategies from treatment. We investigated the WGS of 30 UPEC strains from UTI samples across the world (2016-2019) and found 25 UPEC strains carrying 2-23 antibiotic resistance genes (ARGs) scattered across 1-3 plasmids per strain. Different ARGs (blaTEM, blaCTXM, blaNDM, blaOXA, blaCMY) encoding extended-spectrum beta-lactamases (TEM, CTXM, CMY) and carbapenemases (NDM, OXA) were found in 24/30, ARGs encoding aminoglycoside modifying enzymes (AAC, APH, AAD) variants in 23/30, trimethoprim ARGs (dfrA17, dfrA12, dfrA5, dfrB4 variants) encoding dihydrofolate reductase in 19/30 and sulfonamide ARGs (sul1, sul2, sul3) encoding dihydropteroate synthase and macrolide ARGs (mph1) encoding macrolide 2' phosphotransferase in 15/30 UPEC strains. Collectively the ARGs were distributed in different combinations in 40 plasmids across UPEC strains with 20 plasmids displaying co-occurrence of multiple ARGs conferring resistance to beta lactam, aminoglycoside, sulfonamide, trimethoprim and macrolide antibiotics. These resistance plasmids belonged to seven incompatibility groups (IncF, IncI, IncC, IncH, IncN, IncB and Col), with IncFI and IncFII being the predominant resistance plasmids. Additionally, we observed co-occurrence of specific mutation pattern in quinolone resistance determining region (QRDR) viz., DNA gyrase (gyrA: S83L, D87N), and topoisomerase IV (parC: S80I, E84V; parE: I529L) in 18/30 strains. The strains also harbored diverse virulence genes, such as fimH, gad, iss, iha, ireA, iroN, cnf1 and san. Multilocus sequence typing (MLST) reconfirmed ST131(n = 10) as the predominant global high-risk clonal strain causing UTI. In summary, our findings contribute to better understand the plasmid mediated ARGs and its encoded enzymes that may contribute in antibiotic inactivation/modification or alteration in the antibiotic target site in high risk MDR hypervirulent UPEC strains causing UTI. The study reinforces the need to characterize and design appropriate inhibitors to counterattack different enzymes and devise strategies to curtail resistance plasmid.

Translation of continuous artificial selection on phenotype into genotype during rice breeding programs.

Understanding genetic diversity among local populations is a primary goal of modern crop breeding programs. Here, we demonstrated the genetic relationships of rice varieties in Hokkaido, Japan, one of the northern limits of rice cultivation around the world. Furthermore, artificial selection during rice breeding programs has been characterized using genome sequences. We utilized 8,565 single nucleotide polymorphisms and insertion/deletion markers distributed across the genome in genotype-by-sequencing for genetic diversity analyses. Phylogenetics, genetic population structure, and principal component analysis showed that a total of 110 varieties were classified into four distinct clusters according to different populations geographically and historically. Furthermore, the genome sequences of 19 rice varieties along with historic representations in Hokkaido, nucleotide diversity and FST values in each cluster revealed that artificial selection of elite phenotypes focused on chromosomal regions. These results clearly demonstrated the history of the selections on agronomic traits as genome sequences among current rice varieties from Hokkaido.

Genetic Analysis of Tryptophan Metabolism Genes in Sporadic Amyotrophic Lateral Sclerosis.

The essential amino acid tryptophan (TRP) is the initiating metabolite of the kynurenine pathway (KP), which can be upregulated by inflammatory conditions in cells. Neuroinflammation-triggered activation of the KP and excessive production of the KP metabolite quinolinic acid are common features of multiple neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). In addition to its role in the KP, genes involved in TRP metabolism, including its incorporation into proteins, and synthesis of the neurotransmitter serotonin, have also been genetically and functionally linked to these diseases. ALS is a late onset neurodegenerative disease that is classified as familial or sporadic, depending on the presence or absence of a family history of the disease. Heritability estimates support a genetic basis for all ALS, including the sporadic form of the disease. However, the genetic basis of sporadic ALS (SALS) is complex, with the presence of multiple gene variants acting to increase disease susceptibility and is further complicated by interaction with potential environmental factors. We aimed to determine the genetic contribution of 18 genes involved in TRP metabolism, including protein synthesis, serotonin synthesis and the KP, by interrogating whole-genome sequencing data from 614 Australian sporadic ALS cases. Five genes in the KP (AFMID, CCBL1, GOT2, KYNU, HAAO) were found to have either novel protein-altering variants, and/or a burden of rare protein-altering variants in SALS cases compared to controls. Four genes involved in TRP metabolism for protein synthesis (WARS) and serotonin synthesis (TPH1, TPH2, MAOA) were also found to carry novel variants and/or gene burden. These variants may represent ALS risk factors that act to alter the KP and lead to neuroinflammation. These findings provide further evidence for the role of TRP metabolism, the KP and neuroinflammation in ALS disease pathobiology.