Analysis Of Whole Genome Sequencing Data Research Articles

A rare genetic variant in APEX1 is associated with familial amyotrophic lateral sclerosis with slow progression.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by loss of motor neurons and progressive muscle weakness. We aimed to identify the pathogenic genetic variants in familial ALS (fALS) pedigrees and to elucidate their impact on the disease phenotype. Through the analysis of whole-genome sequencing data of 34 fALS probands that screened negative for mutations in the most common ALS-causing genes, we identified a rare missense variant in APEX1 (NM_001641.4: c.22G > A, p.Gly8Arg) associated with ALS in one pedigree. Fluorescence microscopy images using green fluorescent protein (GFP)-fusion proteins suggested that this amino acid substitution could cause an impairment in nuclear localization of the protein. We described the clinical characteristics of this cohort analyzed and found that patients carrying this variant exhibit lower motor neuron onset and prolonged survival. The relation between APEX1 and ALS occurrence has been elusive despite evidence of a neuroprotective role for the gene. This study provides evidence linking an APEX1 variant with fALS and information on the distinct clinical manifestation. This study contributes to the understanding of the genetic basis of ALS, as well as a potential mechanism leading to loss of neurons, highlighting possible opportunities of targeted treatment harnessing the DNA repair process or ameliorating the oxidative stress.

Genome-scale analysis of Mycobacterium avium complex isolates from Portugal reveals extensive genetic diversity

Opportunist infections caused by nontuberculous mycobacteria (NTM) have emerged as a significant public health problem. Among these, species of the Mycobacterium avium complex (MAC) are the main responsible for the increase in the number of human disease cases. In order to address the current needs in the detection and surveillance of MAC disease cases, we evaluated different species classification methodologies (BLASTn-based marker-gene approach, Kraken v2, rMLST and MLST databases) and their congruence with a core-SNP phylogenetic approach, based on whole genome sequencing (WGS) data. For this purpose, we used a collection of 142 MAC isolates from Portuguese patients diagnosed between 2014 and 2022. The marker-gene approach (based on the rpoB, hsp65 and groEL genes), showed the best results, allowing the identification of the 142 MAC isolates to the species/subspecies level (M. avium subsp. hominissuis, M. intracellulare, M. intracellulare subsp. chimaera, M. intracellulare subsp. yongonense, M. marseillence and M. colombiense). Additionally, we performed drug susceptibility testing that confirmed clarithromycin efficacy as a first-line treatment for MAC disease, as 93 % of the Portuguese isolates were susceptible. Using a core-SNP approach we also performed an in-depth phylogenetic analysis within each identified species group, and despite the high genetic diversity within the MAC species, we were able to clearly distinguish all the species/subspecies and identify genetic clusters with epidemiological potential.We highlight not only the need for the standardization of an appropriate genotyping approach for species identification and management of MAC disease, but also a more robust large-scale WGS data analysis, in a One Health perspective, in order to identify potential routes of transmission.

Mycobacterium tuberculosis evolution from monoresistance to pre-extensive drug resistance during a prolonged household outbreak

Whole genome sequencing (WGS) is sensitive tool for the analysis of tuberculosis transmission and drug-resistance. We used WGS to analyze the Mycobacterium tuberculosis evolution from isoniazid monoresistance to MDR/preXDR during a prolonged household outbreak. The outbreak started with a isoniazid resistant strain (katG S315T mutation) and evolve in two cases to pre-XDR phenotype (with mutations in katG, rpoB, embB, pncA and gyrA genes). Based on WGS data and epidemiological interview we proposed a possible chain of transmission an evolution of the strains.Similar intra-patient and inter-patient acquisition of variability was observed, making difficult to distinguish reinfection or reactivation. Analysis of WGS data together with epidemiological clinical history are discussed in order to distinguish between prolonged infections or transition from latency to reactivation. Classical interview and clinical history taking should be consider to fully understanding WGS data. With a still low incidence of TB cases, Uruguay could use universal WGS of all isolates to reduce time of diagnosis, detect outbreaks and perform public actions to reduce TB incidence.

Whole-Genome Sequencing of Newly Emerged Fungal Pathogen Aspergillus Lentulus and Its Azole Resistance Gene Prediction.

Aims: Aspergillus lentulus is an important newly recorded species in the A. fumigatus complex and its resistance to azole drugs and the high mortality rate of infected individuals have emerged as problems. Comprehensive understanding of the A. lentulus is limited due to lack of genome-wide fine mapping data. The aim of this study was to investigate the A. lentulus signature at the molecular level, analyze the genome-wide profile of this strain, and predict its possible genes that execute azole resistance. Methods: In this study, a whole-genome sequencing of a clinically isolated A. lentulus strain (named A. lentulus PWCAL1) was studied by PacBio Sequel sequencing platform. Azole resistance genes were predicted based on whole-genome sequencing data analysis, gene function annotation, comparative genomic analysis, and BLASTP alignment using the Mycology Antifungal Resistance Database to comprehensively understanding the genome-wide features, pathogenicity, and resistance mechanisms of A. lentulus. Results: The results of whole-genome sequencing demonstrated that the total length of A. lentulus PWCAL1 genome was 31255105 bp, the GC content was 49.24%, and 6883 coding genes were predicted. A total of 4565, 1824, and 6405 genes were annotated in the Gene Ontology, Clusters of Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes databases, respectively. In the Pathogen Host Interactions Database and the Database of Fungal Virulence Factors, 949 and 259 interacting virulence factors were identified, respectively, with the main virulence factor-mutant virulence phenotype, being enriched in reduced virulence. Comparative genomic analysis showed that there were 5456 consensus core genes in this strain and four closely related strains of A. fumigatus complex, which were mainly involved in human diseases, metabolism, organismal systems, etc. Among the three aligned A. lentulus strains, the number of unique genes of this bacterium was the highest with a number of 171, and these genes were mainly associated with carbohydrate metabolism and cell growth and death. Resistance gene prediction demonstrated that the A5653 gene of this bacterium had F46Y/N248T double point mutations on the CYP51A gene, but no tandem repeat mutations in the promoter region were detected. Furthermore, 12 genes belonging to the fungal multidrug resistance ATP-binding cassette (ABC) transporters were identified based on the complete genome sequence and phylogenetic analysis of A. lentulus, which belonged to the ALDp subfamily, the PDR subfamily (AtrB, CDR1, and CDR2), and the MDR subfamily (MDR1), respectively, and there were four genes that are annotated to the major facilitator superfamily multidrug transporter. Further phylogenetic tree classification of the ABC transporter subfamilies predicted in the nine selected A. fumigatus complex strains showed that these putative ABC proteins were divided into two main clusters, which belonged to the PDR (CDR1, CDR2, AtrB, and AtrF) and MDR subfamilies (MDR1, MDR2, and MDR3). The distribution of these ABC proteins varies among different species of the A. fumigatus complex. Conclusions: The main result obtained from this study for the whole genome of A. lentulus provide new insights into better understanding the biological characteristics, pathogenicity, and resistance mechanisms of this bacterium. In this study, two resistance mechanisms, which include CYP51A gene mutation and multidrug-resistant ABC transporter, were predicted in a single isolate. Based on the predicted CYP51A-F46Y/N248T site mutation combination, we speculate that the CYP51A gene of A. lentulus may be partially responsible for azole resistance. Based on the predicted ABC transporter family genes, we hypothesize that resistance to multiple azoles in A. lentulus is mediated, at least in part, by these ABC transporters with resistance.

Whole-genome Sequencing Association Analysis of Quantitative Platelet Traits in A Large Cohort of β-thalassemia.

Platelet acts as a crucial monitoring indicator for hypercoagulability and thrombosis and a key target for drug regulation. Genotype-phenotype association studies have confirmed that platelet traits are quantitatively regulated by multiple genes. However, there is currently a lack of genetic studies on the heterogeneity of platelet traits in β-thalassemia under hypercoagulable state. Here, we studied the phenotypic heterogeneity of platelet count (PLT) and mean platelet volume (MPV) in 1020 β-thalassemia patients. We further performed a functionally informed whole genome sequencing association analysis of common variants and rare variants (RVs) for PLT and MPV in 916 patients through integrative analysis of whole-genome sequencing data and functional annotation data. Extreme phenotypic heterogeneity of platelet traits was observed in β-thalassemia patients. Additionally, the common variant based gene-level analysis identified the novel gene of RNF144B associated with MPV. The RV analysis identified several novel associations in both coding and noncoding genome, including missense RVs of PPP2R5C associated with PLT and missense RVs of TSSK1B associated with MPV. In conclusion, we performed a comprehensive and systematic whole genome scan of platelet traits in the β-thalassemia cohort, demonstrating the specificity of genetic regulation of platelet traits in the context of β-thalassemia, providing potential targets for intervention.

Easy analysis of bacterial whole-genome sequencing data for clinical microbiologists using open-source Galaxy platform: Characterization of ESBL-producing Enterobacterales from bloodstream infections

ObjectivesClinical microbiologists require easy-to-use open access tools with graphical interfaces to perform bacterial whole-genome sequencing (WGS) in routine practice. This study aimed to build a bioinformatics pipeline on the open-source Galaxy platform, facilitating comprehensive and reproducible analysis of bacterial WGS data in a few steps. We then used it to characterize our local epidemiology of ESBL-producing Enterobacterales isolated from patients with bacteremia. MethodsWe built a bioinformatics pipeline consisting of the following sequential tools: Fastp (input data trimming); FastQC (read quality control); SPAdes (genome assembly); Quast (quality control of genome assembly); Prokka (gene annotation); Staramr (ResFinder database) and ABRicate (CARD database) for antimicrobial resistance (AMR) gene screening and molecular strain typing. Paired-end short read WGS data from all ESBL-producing Enterobacterales strains isolated from patients with bacteremia over one year were analysed. ResultsThe Galaxy platform does not require command line tools. The bioinformatics pipeline was constructed within one hour. It only required uploading fastq files and facilitated systematization of de novo assembly of genomes, MLST typing, and AMR gene screening in one step. Among the 66 ESBL-producing strains analysed, the two most frequent ESBL genes were blaCTX−M-15 (62.1%) and blaCTX−M-27 (13.6%). ConclusionsThe open-access Galaxy platform provides a graphical interface and easy-to-use tools suitable for routine use in clinical microbiology laboratories without bioinformatics specialists. We believe that this platform will facilitate fast and low-cost analysis of bacterial WGS data, especially in resource-limited settings.

A plasmid-mediated type III secretion system associated with invasiveness and diarrheagenicity of Providencia rustigianii.

We have recently described a clinical isolate of Providencia rustigianii strain JH-1 carrying the genes for cytolethal distending toxin (CDT) in a conjugative plasmid. A cdtB mutant of strain JH-1, which lost CDT activity, was still found to retain invasiveness and diarrheagenicity. The strain was subjected to phenotypic and genetic analyses including whole genome sequencing (WGS) to explore the genetic determinants of the observed invasiveness and diarrheagenic properties. Analysis and annotation of WGS data revealed the presence of two distinct type III secretion systems (T3SS) in strain JH-1, one of which was located on the chromosome designated as cT3SS (3,992,833 bp) and the other on a mega-plasmid designated as pT3SS (168,819 bp). Comparative genomic analysis revealed that cT3SS is generally conserved in Providencia spp. but pT3SS was limited to a subset of Providencia spp., carrying cdt genes. Strain JH-1 was found to invade HeLa cells and induce fluid accumulation with characteristic pathological lesions in rabbit ileal loops. Remarkably, these phenomena were associated with the pT3SS but not cT3SS. The plasmid could be transferred by conjugation from strain JH-1 to other strains of P. rustigianii, Providencia rettgeri, and Escherichia coli with concomitant transfer of these virulence properties. This is the first report of a functional and mobile T3SS in P. rustigianii and its association with invasiveness and diarrheagenicity of this bacterium. These data suggest that P. rustigianii and other CDT-producing Providencia strains might carry T3SS and exert their diarrheagenic effect by exploiting the T3SS nano-machinery.IMPORTANCEThe precise mechanism of virulence of Providencia rustigianii is unclear, although some strains produce cytolethal distending toxin as a putative virulence factor. We have detected the presence of a type III secretion system (T3SS) for the first time on a plasmid in a P. rustigianii strain. Plasmid-mediated T3SS seems to be directly involved in virulence of P. rustigianii and may serve as a means of horizontal transfer of T3SS genes. Our results may have implication in understanding the mechanism of emergence of new pathogenic strains of P. rustigianii.

Association mapping with a diverse population of Puccinia graminis f. sp. tritici identified avirulence loci interacting with the barley Rpg1 stem rust resistance gene

BackgroundWheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is an important disease of barley and wheat. A diverse sexual Pgt population from the Pacific Northwest (PNW) region of the US contains a high proportion of individuals with virulence on the barley stem rust resistance (R) gene, Rpg1. However, the evolutionary mechanisms of this virulence on Rpg1 are mysterious considering that Rpg1 had not been deployed in the region and the gene had remained remarkably durable in the Midwestern US and prairie provinces of Canada.Methods and resultsTo identify AvrRpg1 effectors, genome wide association studies (GWAS) were performed using 113 Pgt isolates collected from the PNW (n = 89 isolates) and Midwest (n = 24 isolates) regions of the US. Disease phenotype data were generated on two barley lines Morex and the Golden Promise transgenic (H228.2c) that carry the Rpg1 gene. Genotype data was generated by whole genome sequencing (WGS) of 96 isolates (PNW = 89 isolates and Midwest = 7 isolates) and RNA sequencing (RNAseq) data from 17 Midwestern isolates. Utilizing ~1.2 million SNPs generated from WGS and phenotype data (n = 96 isolates) on the transgenic line H228.2c, 53 marker trait associations (MTAs) were identified. Utilizing ~140 K common SNPs generated from combined analysis of WGS and RNAseq data, two significant MTAs were identified using the cv Morex phenotyping data. The 55 MTAs defined two distinct avirulence loci, on supercontig 2.30 and supercontig 2.11 of the Pgt reference genome of Pgt isolate CRL 75-36-700-3. The major avirulence locus designated AvrRpg1A was identified with the GWAS using both barley lines and was delimited to a 35 kb interval on supercontig 2.30 containing four candidate genes (PGTG_10878, PGTG_10884, PGTG_10885, and PGTG_10886). The minor avirulence locus designated AvrRpg1B identified with cv Morex contained a single candidate gene (PGTG_05433). AvrRpg1A haplotype analysis provided strong evidence that a dominant avirulence gene underlies the locus.ConclusionsThe association analysis identified strong candidate AvrRpg1 genes. Further analysis to validate the AvrRpg1 genes will fill knowledge gaps in our understanding of rust effector biology and the evolution and mechanism/s of Pgt virulence on Rpg1.

Genome analysis of Salmonella enterica serovar enteritis strains isolated from poultry and humans in Burkina Faso.

Whole-genome sequencing (WGS) was used to characterize four Salmonella enterica Enteritidis isolates from poultry (n=2) and human (n=2) from Ouagadougou, Burkina Faso. Antimicrobial resistance genes, chromosomal mutations, and mobile genetic elements were identified by analysis of WGS data using sequence homology.

RSV-GenoScan: An automated pipeline for whole-genome human respiratory syncytial virus (RSV) sequence analysis

BackgroundAdvances in high-throughput sequencing (HTS) technologies and reductions in sequencing costs have revolutionised the study of genomics and molecular biology by making whole-genome sequencing (WGS) accessible to many laboratories. However, the analysis of WGS data requires significant computational effort, which is the major drawback in implementing WGS as a routine laboratory technique. ObjectiveAutomated pipelines have been developed to overcome this issue, but they do not exist for all organisms. This is the case for human respiratory syncytial virus (RSV), which is a leading cause of lower respiratory tract infections in infants, the elderly, and immunocompromised adults. ResultsWe present RSV-GenoScan, a fast and easy-to-use pipeline for WGS analysis of RSV generated by HTS on Illumina or Nanopore platforms. RSV-GenoScan automates the WGS analysis steps directly from the raw sequence data. The pipeline filters the sequence data, maps the reads to the RSV reference genomes, generates a consensus sequence, identifies the RSV subgroup, and lists amino acid mutations, insertions and deletions in the F and G viral genes. This enables the rapid identification of mutations in these coding genes that are known to confer resistance to monoclonal antibodies. AvailabilityRSV-GenoScan is freely available at https://github.com/AlexandreD-bio/RSV-GenoScan.

Abstract 3887: The origins and clinical impact of extrachromosomal DNA across 39 cancers

Abstract Extrachromosomal DNA (ecDNA) is emerging as a major contributor to cancer drug resistance and poor outcome in many cancer types. To understand the underlying genetic and mutational processes associated with ecDNA formation and presence, we conducted an analysis of whole-genome sequencing data from 14,778 cancer patients representing 39 tumor types from the Genomics England cohort. Our findings revealed focal amplifications driven by ecDNA in 17.1% of tumor samples, exhibiting varying frequencies across tissues. Notably, liposarcomas (55.4%) and HER2+ breast cancers (46.4%) displayed higher rates compared to lower rates found in pancreatic adenocarcinomas (2%). In addition to oncogenic amplifications such as HER2 derived from ecDNA during cancer evolution, ecDNAs that did not contain oncogenes were enriched for immunomodulatory genes (e.g., CXCL6, LRRC32, and HLA-G). Immunomodulatory ecDNAs exhibited reduced T cell infiltration. In the context of intrinsic and environmental mutational processes, deamination clock-like, tobacco smoking, and APOBEC cytidine deaminase signatures positively correlated with ecDNA presence. Conversely, mismatch repair deficiency (MMRd) and POLD1/POLE deficiency signatures exhibited a negative correlation. Mutational processes associated with tobacco exposure, UV light, and spontaneous deamination primarily occurred before ecDNA formation, while signatures of homologous recombination repair deficiency tended to manifest post-ecDNA formation.Notably, ecDNA was enriched at metastatic sites and significantly associated with shorter overall survival in multivariable analysis (HR 1.31, 95% CI 1.17-1.47), underscoring its clinical relevance. This study enhances our understanding of the mutational processes and tissue contexts that shape the ecDNA landscape in human cancer and highlights its implications for clinical outcomes. Citation Format: Chris Bailey, Oriol Pich, Kerstin Thol, Jens Luebeck, Georgia Stavrou, Thomas B. Watkins, Natasha E. Weiser, Wei-Ting Lu, Jeanette Kittel, Robert Bentham, Natasha Sharma, Roberto Salgado, Dave Moore, Maria Litovchenko, King L. Hung, Alex J. Cornish, Bhargavi Dameracharla, Serena Nik-Zainal, Vineet Bafna, Howard Chang, Nicholas McGranahan, Nnennaya Kanu, Genomics England Consortium, Paul S. Mischel, Mariam Jamal-Hanjani, Charles Swanton. The origins and clinical impact of extrachromosomal DNA across 39 cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 3887.

Phenotypic and genotypic characterization of Neisseria gonorrhoeae isolates among individuals at high risk for sexually transmitted diseases in Zurich, Switzerland.

While ceftriaxone resistance remains scarce in Switzerland, global Neisseria gonorrhoeae (NG) antimicrobial resistance poses an urgent threat. This study describes clinical characteristics in MSM (men who have sex with men) diagnosed with NG infection and analyses NG resistance by phenotypic and genotypic means. Data of MSM enrolled in three clinical cohorts with a positive polymerase chain reaction test (PCR) for NG were analysed between January 2019 and December 2021 and linked with antibiotic susceptibility testing. Bacterial isolates were subjected to whole genome sequencing (WGS). Of 142 participants, 141 (99%) were MSM and 118 (84%) living with HIV. Participants were treated with ceftriaxone (N = 79), azithromycin (N = 2), or a combination of both (N = 61). No clinical or microbiological failures were observed. From 182 positive PCR samples taken, 23 were available for detailed analysis. Based on minimal inhibitory concentrations (MICs), all isolates were susceptible to ceftriaxone, gentamicin, cefixime, cefpodoxime, ertapenem, zoliflodacin, and spectinomycin. Resistance to azithromycin, tetracyclines and ciprofloxacin was observed in 10 (43%), 23 (100%) and 11 (48%) of the cases, respectively. Analysis of WGS data revealed combinations of resistance determinants that matched with the corresponding phenotypic resistance pattern of each isolate. Among the MSM diagnosed with NG mainly acquired in Switzerland, ceftriaxone MICs were low for a subset of bacterial isolates studied and no treatment failures were observed. For azithromycin, high occurrences of in vitro resistance were found. Gentamicin, cefixime, cefpodoxime, ertapenem, spectinomycin, and zoliflodacin displayed excellent in vitro activity against the 23 isolates underscoring their potential as alternative agents to ceftriaxone.

Genomic comparison between Mycobacterium bovis and Mycobacterium microti and in silico analysis of peptide-based biomarkers for serodiagnosis.

In recent years, there has been an increase in the number of reported cases of Mycobacterium microti infection in various animals, which can interfere with the ante-mortem diagnosis of animal tuberculosis caused by Mycobacterium bovis. In this study, whole genome sequencing (WGS) was used to search for protein-coding genes to distinguish M. microti from M. bovis. In addition, the population structure of the available M. microti genomic WGS datasets is described, including three novel Belgian isolates from infections in alpacas. Candidate genes were identified by examining the presence of the regions of difference and by a pan-genome analysis of the available WGS data. A total of 80 genes showed presence-absence variation between the two species, including genes encoding Proline-Glutamate (PE), Proline-Proline-Glutamate (PPE), and Polymorphic GC-Rich Sequence (PE-PGRS) proteins involved in virulence and host interaction. Filtering based on predicted subcellular localization, sequence homology and predicted antigenicity resulted in 28 proteins out of 80 that were predicted to be potential antigens. As synthetic peptides are less costly and variable than recombinant proteins, an in silico approach was performed to identify linear and discontinuous B-cell epitopes in the selected proteins. From the 28 proteins, 157 B-cell epitope-based peptides were identified that discriminated between M. bovis and M. microti species. Although confirmation by in vitro testing is still required, these candidate synthetic peptides containing B-cell epitopes could potentially be used in serological tests to differentiate cases of M. bovis from M. microti infection, thus reducing misdiagnosis in animal tuberculosis surveillance.

Streptococcus pneumoniae serotypes 8 and 22F causing invasive pneumococcal disease in the Czech Republic in 2014-2020: whole genome sequencing (WGS) analysis.

An analysis is presented of whole genome data of Streptococcus pneumoniae serotypes 8 and 22F isolated in the Czech Republic from invasive pneumococcal disease (IPD) in 2014-2020. New multivalent pneumococcal conjugate vaccines (PCVs) are effective against these serotypes. Recently, serotypes 8 and 22F have been among the leading causes of IPD in the Czech Republic.S. pneumoniae isolates from the Czech Republic were compared with those of the same serotypes recovered in other countries in the same period and available in the international database PubMLST. Isolates from IPD of serotypes 8 (22 isolates) and 22F (21 isolates) recovered in the Czech Republic in 2014-2020 were subjected to whole genome sequencing (WGS). The genomes were analysed and compared using the international database PubMLST. Most of the studied Czech serotype 8 isolates belong to two main subpopulations. The first subpopulation, dominated by ST-53 isolates, is part of a highly abundant group of genetically close European and non-European isolates that are clearly separated on the phylogenetic network. The second subpopulation of Czech serotype 8 isolates (dominated by ST-404) is more genetically variable and forms a separate lineage on the global phylogenetic network, with no other European isolates. Czech isolates of serotype 22F are a homogeneous population with a clear predominance of ST-433, which belongs to a genetically close European population. The analysis of WGS data of IPD isolates of serotypes 8 and 22F provided a detailed insight into the genetic relationships between the Czech populations of these serotypes. It also allowed comparison of the Czech populations with the matched populations from other European and non-European countries. The obtained results add to the body of knowledge about the spread of genetic lineages causing IPD in the Czech Republic in the post-vaccination period and provide a basis for considering whether the use of the new multivalent PCVs in the Czech Republic would be beneficial.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds

BackgroundThe genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome.ResultsWe identified several thousand variants that had significantly steep clines (‘SCV’) across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle.ConclusionsTo our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Whole-Genome Sequencing-Based Profiling of Antimicrobial Resistance Genes and Core-Genome Multilocus Sequence Typing of Campylobacter jejuni from Different Sources in Lithuania.

Campylobacter jejuni is known as one of the main causative agents of gastroenteritis in humans worldwide, and the rise of antimicrobial resistance (AMR) in Campylobacter is a growing public health challenge of special concern. Whole-genome sequencing (WGS) was used to characterize genetic determinants of AMR in 53 C. jejuni isolates from dairy cattle, broiler products, wild birds, and humans in Lithuania. The WGS-based study revealed 26 C. jejuni AMR markers that conferred resistance to various antimicrobials. Genetic markers associated with resistance to beta-lactamases, tetracycline, and aminoglycosides were found in 79.3%, 28.3%, and 9.4% of C. jejuni isolates, respectively. Additionally, genetic markers associated with multidrug resistance (MDR) were found in 90.6% of C. jejuni isolates. The WGS data analysis revealed that a common mutation in the quinolone resistance-determining region (QRDR) was R285K (854G > A) at 86.8%, followed by A312T (934G > A) at 83% and T86I (257C > T) at 71.7%. The phenotypic resistance analysis performed with the agar dilution method revealed that ciprofloxacin (CIP) (90.6%), ceftriaxone (CRO) (67.9%), and tetracycline (TET) (45.3%) were the predominant AMR patterns. MDR was detected in 41.5% (22/53) of the isolates tested. Fifty-seven virulence genes were identified in all C. jejuni isolates; most of these genes were associated with motility (n = 28) and chemotaxis (n = 10). Additionally, all C. jejuni isolates harbored virulence genes related to adhesion, invasion, LOS, LPS, CPS, transportation, and CDT. In total, 16 sequence types (STs) and 11 clonal complexes (CC) were identified based on core-genome MLST (cgMLST) analysis. The data analysis revealed distinct diversity depending on phenotypic and genotypic antimicrobial resistance of C. jejuni.

Genetic determinants underlying the progressive phenotype of β-lactam/β-lactamase inhibitor resistance in Escherichia coli.

Currently, whole-genome sequencing (WGS) data have not shown strong concordance with Escherichia coli susceptibility profiles to the commonly used β-lactam/β-lactamase inhibitor (BL/BLI) combinations: ampicillin-sulbactam (SAM), amoxicillin-clavulanate (AMC), and piperacillin-tazobactam (TZP). Progressive resistance to these BL/BLIs in the absence of cephalosporin resistance, also known as extended-spectrum resistance to BL/BLI (ESRI), has been suggested to primarily result from increased copy numbers of blaTEM variants, which is not routinely assessed in WGS data. We sought to determine whether addition of gene amplification could improve genotype-phenotype associations through WGS analysis of 147 E. coli bacteremia isolates with increasing categories of BL/BLI non-susceptibility ranging from ampicillin (AMP) susceptibility to being fully resistant to all three BL/BLIs. Consistent with a key role of blaTEM in ESRI, 112/134 strains (84%) with at least ampicillin non-susceptibility encoded blaTEM. Evidence of blaTEM amplification (i.e., blaTEM gene copy number estimates > 2×) was present in 40/112 (36%) strains. There were positive correlations between blaTEM copy numbers with minimum inhibitory concentrations of AMC and TZP (P < 0.05) but not for SAM (P = 0.09). The diversity of β-lactam resistance mechanisms, including non-ceftriaxone hydrolyzing blaCTX-M variants, blaOXA-1, and ampC and blaTEM strong promoter mutations, was greater in AMC- and TZP-non-susceptible strains but rarely observed within SAM- and AMP-non-susceptible isolates. Our study indicates that comprehensive analysis of WGS data, including β-lactamase-encoding gene amplification, can help categorize E. coli with AMC or TZP non-susceptibility but that discerning the transition from SAM susceptibility to SAM non-susceptibility using genetic data requires further refinement. IMPORTANCE The increased feasibility of whole-genome sequencing has generated significant interest in using such molecular diagnostic approaches to characterize difficult-to-treat, antimicrobial-resistant (AMR) infections. Nevertheless, there are current limitations in the accurate prediction of AMR phenotypes based on existing AMR gene database approaches, which primarily correlate a phenotype with the presence/absence of a single AMR gene. Our study utilized a large cohort of cephalosporin-susceptible Escherichia coli bacteremia samples to determine how increasing the dosage of narrow-spectrum β-lactamase-encoding genes in conjunction with other diverse β-lactam/β-lactamase inhibitor (BL/BLI) genetic determinants contributes to progressively more severe BL/BLI phenotypes. We were able to characterize the complexity of the genetic mechanisms underlying progressive BL/BLI resistance including the critical role of β-lactamase encoding gene amplification. For the diverse array of AMR phenotypes with complex mechanisms involving multiple genomic factors, our study provides an example of how composite risk scores may improve understanding of AMR genotype/phenotype correlations.

ABCA4 c.6480-35A>G, a novel branchpoint variant associated with Stargardt disease.

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 (ABCA4) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone-rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants. Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro. The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48. Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699_768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans. For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases. Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis.

Use of genome sequencing to hunt for cryptic second-hit variants: analysis of 31 cases recruited to the 100 000 Genomes Project

BackgroundCurrent clinical testing methods used to uncover the genetic basis of rare disease have inherent limitations, which can lead to causative pathogenic variants being missed. Within the rare disease arm...

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.