Comparative Genomics Approach Research Articles

Development and evaluation of genomics informed real-time PCR assays for the detection and strain typing of Mycobacterium avium subsp. paratuberculosis.

This study aimed to identify specific genomic targets for the detection and strain typing of Map and analyse their sensitivity and specificity, and detect Map directly from faeces. A comparative genomics approach was used to identify specific genomic targets for the detection and strain typing of Map. A Map specific qPCR using the primer pair 7132 that targets a DNA segregation ATPase protein was able to detect all strains of Map and is more sensitive than the current Johne's disease PCR assays with a sensitivity of 0.0002 fg µl-1. A strain specific qPCR using the Atsa primer pair that targets the arylsulfase gene was able to differentiate between Type S and Type C strains of Map and was more sensitive than the IS1311 PCR and REA with a sensitivity of 40 fg µl-1 and was specific for Type S Map. Both assays successfully detected Map directly from faeces. This study developed and validated two genomics informed qPCR assays, 7132B Map and Atsa Type S and found both assays to be highly specific and sensitive for the detection of Map from culture and directly from faeces. This is the first time that a probe-based qPCR has been designed and developed for Map strain typing, which will greatly improve the response time during outbreak investigations.

Genomic Insights into Cyanide Biodegradation in the Pseudomonas Genus.

Molecular studies about cyanide biodegradation have been mainly focused on the hydrolytic pathways catalyzed by the cyanide dihydratase CynD or the nitrilase NitC. In some Pseudomonas strains, the assimilation of cyanide has been linked to NitC, such as the cyanotrophic model strain Pseudomonas pseudoalcaligenes CECT 5344, which has been recently reclassified as Pseudomonas oleovorans CECT 5344. In this work, a phylogenomic approach established a more precise taxonomic position of the strain CECT 5344 within the species P. oleovorans. Furthermore, a pan-genomic analysis of P. oleovorans and other species with cyanotrophic strains, such as P. fluorescens and P. monteilii, allowed for the comparison and identification of the cioAB and mqoAB genes involved in cyanide resistance, and the nitC and cynS genes required for the assimilation of cyanide or cyanate, respectively. While cyanide resistance genes presented a high frequency among the analyzed genomes, genes responsible for cyanide or cyanate assimilation were identified in a considerably lower proportion. According to the results obtained in this work, an in silico approach based on a comparative genomic approach can be considered as an agile strategy for the bioprospection of putative cyanotrophic bacteria and for the identification of new genes putatively involved in cyanide biodegradation.

High-throughput prioritization of target proteins for development of new antileishmanial compounds

Leishmaniasis, a vector-borne disease, is caused by the infection of Leishmania spp., obligate intracellular protozoan parasites. Presently, human vaccines are unavailable, and the primary treatment relies heavily on systemic drugs, often presenting with suboptimal formulations and substantial toxicity, making new drugs a high priority for LMIC countries burdened by the disease, but a low priority in the agenda of most pharmaceutical companies due to unattractive profit margins. New ways to accelerate the discovery of new, or the repositioning of existing drugs, are needed. To address this challenge, our study aimed to identify potential protein targets shared among clinically-relevant Leishmania species. We employed a subtractive proteomics and comparative genomics approach, integrating high-throughput multi-omics data to classify these targets based on different druggability metrics. This effort resulted in the ranking of 6502 ortholog groups of protein targets across 14 pathogenic Leishmania species. Among the top 20 highly ranked groups, metabolic processes known to be attractive drug targets, including the ubiquitination pathway, aminoacyl-tRNA synthetases, and purine synthesis, were rediscovered. Additionally, we unveiled novel promising targets such as the nicotinate phosphoribosyltransferase enzyme and dihydrolipoamide succinyltransferases. These groups exhibited appealing druggability features, including less than 40% sequence identity to the human host proteome, predicted essentiality, structural classification as highly druggable or druggable, and expression levels above the 50th percentile in the amastigote form. The resources presented in this work also represent a comprehensive collection of integrated data regarding trypanosomatid biology.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition

Background Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a “cut-and-paste” style, and retrotransposons which propagate in a “copy-and-paste” fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a new type of MEs, which we call “retro-DNAs”, for their combined characteristics of these two fundamentally different ME classes. Methods A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset. Results From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the “TT/AAAA” insertion site motif, suggesting the use of the L1-based target-primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes. Conclusions Although small in number, the identification of these retro-DNAs reveals a new mechanism for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

National genomic epidemiology investigation revealed the spread of carbapenem-resistant Escherichia coli in healthy populations and the impact on public health

BackgroundCarbapenem-resistant Escherichia coli (CREC) has been considered as WHO priority pathogens, causing a great public health concern globally. While CREC from patients has been thoroughly investigated, the prevalence and underlying risks of CREC in healthy populations have been overlooked. Systematic research on the prevalence of CREC in healthy individuals was conducted here. We aimed to characterize CREC collected from healthy populations in China between 2020 and 2022 and to compare the genomes of CREC isolates isolated from healthy individuals and clinical patients.MethodsWe present a nationwide investigation of CREC isolates among healthy populations in China, employing robust molecular and genomic analyses. Antimicrobial susceptibility testing, whole-genome sequencing, and bioinformatics were utilized to analyze a cohort of CREC isolates (n = 113) obtained from fecal samples of 5 064 healthy individuals. Representative plasmids were extracted for third-generation nanopore sequencing. We previously collected 113 non-duplicate CREC isolates (59 in 2018, 54 in 2020) collected from ICU patients in 15 provinces and municipalities in China, and these clinical isolates were used to compare with the isolates in this study. Furthermore, we employ comparative genomics approaches to elucidate molecular variations and potential correlations between clinical and non-clinical CREC isolates.ResultsA total of 147 CREC isolates were identified from 5 064 samples collected across 11 provinces in China. These isolates were classified into 64 known sequence types (STs), but no dominant STs were observed. In total, seven carbapenemase genes were detected with blaNDM-5 (n = 116) being the most prevalent one. Genetic environments and plasmid backbones of blaNDM were conserved in CREC isolated from healthy individuals. Furthermore, we compared clinical and healthy human-originated CRECs, revealing noteworthy distinctions in 23 resistance genes, including blaNDM-1, blaNDM-5, and blaKPC (χ2 test, p < 0.05). Clinical isolates contained more virulence factors associated with iron uptake, adhesion, and invasion than those obtained from healthy individuals. Notably, CREC isolates generally found healthy people are detected in hospitalized patients.ConclusionsOur findings underscore the significance of healthy populations-derived CRECs as a crucial reservoir of antibiotic resistance genes (ARGs). This highlights the need for ongoing monitoring of CREC isolates in healthy populations to accurately assess the potential risks posed by clinical CREC isolates.

A comparative genomics approach reveals a local genetic signature of Leishmania tropica in Morocco.

In Morocco, cutaneous leishmaniasis (CL) caused by Leishmania (L.) tropica is an important health problem. Despite the high incidence of CL in the country, the genomic heterogeneity of these parasites is still incompletely understood. In this study, we sequenced the genomes of 14 Moroccan isolates of L. tropica collected from confirmed cases of CL to investigate their genomic heterogeneity. Comparative genomics analyses were conducted by applying the recently established Genome Instability Pipeline (GIP), which allowed us to conduct phylogenomic and principal components analyses (PCA), and to assess genomic variations at the levels of the karyotype, gene copy number, single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) variants. Read-depth analyses revealed a mostly disomic karyotype, with the exception of the stable tetrasomy of chromosome 31. In contrast, we identified important gene copy number variations across all isolates, which affect known virulence genes and thus were probably selected in the field. SNP-based cluster analysis of the 14 isolates revealed a core group of 12 strains that formed a tight cluster and shared 45.1 % (87 751) of SNPs, as well as two strains (M3015, Ltr_16) that clustered separately from each other and the core group, suggesting the circulation of genetically highly diverse strains in Morocco. Phylogenetic analysis, which compared our 14 L. tropica isolates against 40 published genomes of L. tropica from a diverse array of locations, confirmed the genetic difference of our Moroccan isolates from all other isolates examined. In conclusion, our results indicate potential regional variations in SNP profiles that may differentiate Moroccan L. tropica from other L. tropica strains circulating in endemic countries in the Middle East. Our report paves the way for future research with a larger number of strains that will allow correlation of diverse phenotypes (resistance to treatments, virulence) and origins (geography, host species, year of isolation) to defined genomic signals such as gene copy number variations or SNP profiles that may represent interesting biomarker candidates.

RSSC-Lineage Multiplex PCR Assay Detects and Differentiates Ralstonia solanacearum, R. pseudosolanacearum, R. syzygii, and the R3bv2 Subgroup

Bacterial wilt strains in the Ralstonia solanacearum species complex (RSSC) pose serious threats to economically important crops worldwide. In 2014, Safni et al. proposed the reclassification of the RSSC into three genomospecies: R. solanacearum (Rsol), R. pseudosolanacearum (Rpseu), and R. syzygii (Rsyz). The revision requires the proper identification of strains for diagnostic and epidemiological studies. In response, we developed the inexpensive and user-friendly RSSC-Lineage Multiplex PCR, which effectively detects plant-pathogenic Ralstonia strains in general and also distinguishes between Rpseu, Rsol, Rsyz, and the high-security Select Agent “race 3 biovar 2” subgroup of Rsol, also known as the phylotype IIB-1 potato brown rot pandemic lineage. Genomes were retrieved from the NCBI GenBank database and screened for unique gene regions using OrthoMCL and other comparative genomic approaches. Specific primers were designed for each genomospecies, Ralstonia in general, and “race 3 biovar 2.” AT-rich flaps were added at the 5ʹ position of each primer to optimize the reaction thermodynamics. The specificity was tested in silico using the NCBI GenBank genome database and an in-house database. The in vitro specificity and accuracy of the tool was validated with 113 representative Ralstonia strains and 24 strains from other genera. The assay is highly specific, generating neither false positives nor false negatives. Primer set detection limits ranged from 10 to 100 pg. The assay also detected and differentiated strains from naturally and artificially inoculated plant hosts. This tool is highly specific, reliable, and economical for culture characterization, diagnostics, surveys, quarantine decisions, and epidemiological studies. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Genomic insights into the evolution, pathogenicity, and extensively drug-resistance of emerging pathogens Kluyvera and Phytobacter.

Kluyvera is a Gram-negative, flagellated, motile bacillus within the Enterobacteriaceae. The case reports of clinical infections shed light on the importance of this organism as an emerging opportunistic pathogen. The genus Phytobacter, which often be misidentified with Kluyvera, is also an important clinically relevant member of the Enterobacteriaceae. However, the identification of Kluyvera and Phytobacter is problematic, and their phylogenetic relationship remains unclear. Here, 81 strains of Kluyvera and 16 strains of Phytobacter were collected. A series of comparative genomics approaches were applied to the phylogenetic relationship reconstruction, virulence related genes profiles description, and antibiotic resistance genes prediction. Using average nucleotide identity (ANI) and in silico DNA-DNA hybridization (isDDH), we offered reliable species designations of 97 strains, in which 40 (41.24%) strains were incorrectly labeled. A new Phytobacter genomospecies-1 were defined. Phytobacter and Kluyvera show great genome plasticity and inclusiveness, which may be related to their diverse ecological niches. An intergenomic distances threshold of 0.15875 was used for taxonomy reassignments at the phylogenomic-group level. Further principal coordinates analysis (PCoA) revealed 11 core genes of Kluyvera (pelX, mdtL, bglC, pcak-1, uhpB, ddpA-2, pdxY, oppD-1, cptA, yidZ, csbX) that could be served as potential identification targets. Meanwhile, the Phytobacter specific virulence genes clbS, csgA-C, fliS, hsiB1_vipA and hsiC1_vipB, were found to differentiate from Kluyvera. We concluded that the evolution rate of Kluyvera was 5.25E-6, approximately three times higher than that of Phytobacter. Additionally, the co-existence of ESBLs and carbapenem resistance genes were present in approximately 40% strains, suggesting the potential development of extensively drug-resistant or even fully drug-resistant strains. This work provided a better understanding of the differences between closely related species Kluyvera and Phytobacter. Their genomes exhibited great genome plasticity and inclusiveness. They not only possess a potential pathogenicity threat, but also a risk of multi-drug resistance. The emerging pathogens Kluyvera and Phytobacter warrant close attention.

Compensatory evolution in NusG improves fitness of drug-resistant M. tuberculosis

Drug-resistant bacteria are emerging as a global threat, despite frequently being less fit than their drug-susceptible ancestors1–8. Here we sought to define the mechanisms that drive or buffer the fitness cost of rifampicin resistance (RifR) in the bacterial pathogen Mycobacterium tuberculosis (Mtb). Rifampicin inhibits RNA polymerase (RNAP) and is a cornerstone of modern short-course tuberculosis therapy9,10. However, RifR Mtb accounts for one-quarter of all deaths due to drug-resistant bacteria11,12. We took a comparative functional genomics approach to define processes that are differentially vulnerable to CRISPR interference (CRISPRi) inhibition in RifR Mtb. Among other hits, we found that the universally conserved transcription factor NusG is crucial for the fitness of RifR Mtb. In contrast to its role in Escherichia coli, Mtb NusG has an essential RNAP pro-pausing function mediated by distinct contacts with RNAP and the DNA13. We find this pro-pausing NusG–RNAP interface to be under positive selection in clinical RifR Mtb isolates. Mutations in the NusG–RNAP interface reduce pro-pausing activity and increase fitness of RifR Mtb. Collectively, these results define excessive RNAP pausing as a molecular mechanism that drives the fitness cost of RifR in Mtb, identify a new mechanism of compensation to overcome this cost, suggest rational approaches to exacerbate the fitness cost, and, more broadly, could inform new therapeutic approaches to develop drug combinations to slow the evolution of RifR in Mtb.

Comparative genomics of N-acetyl-5-methoxytryptamine members in four Prunus species with insights into bud dormancy and abiotic stress responses in Prunus avium.

This study provides novel insights into the evolution, diversification, and functions of melatonin biosynthesis genes in Prunus species, highlighting their potential role in regulating bud dormancy and abiotic stresses. The biosynthesis of melatonin (MEL) in plants is primarily governed by enzymatic reactions involving key enzymes such as serotonin N-acetyltransferase (SNAT), tryptamine 5-hydroxylase (T5H), N-acetylserotonin methyltransferase (ASMT) and tryptophan decarboxylase (TDC). In this study, we analyzed Melatonin genes in four Prunus species such as Prunus avium (Pavi), Prunus pusilliflora (Ppus), Prunus serulata (Pser), and Prunus persica (Pper) based on comparative genomics approach. Among the four Prunus species, a total of 29 TDCs, 998 T5Hs, 16 SNATs, and 115 ASMTs within the genome of four Prunus genomes. A thorough investigation of melatonin-related genes was carried out using systematic biological methods and comparative genomics. Through phylogenetic analysis, orthologous clusters, Go enrichment, syntenic relationship, and gene duplication analysis, we discovered both similarities and variations in Melatonin genes among these Prunus species. Additionally, our study revealed the existence of unique subgroup members in the Melatonin genes of these species, which were distinct from those found in Arabidopsis genes. Furthermore, the transcriptomic expression analysis revealed the potential significance of melatonin genes in bud dormancy regulation and abiotic stresses. Our extensive results offer valuable perspectives on the evolutionary patterns, intricate expansion, and functions of PavMEL genes. Given their promising attributes, PavTDCs, PavT5H, PavNAT, and three PavASMT genes warrant in-depth exploration as prime candidates for manipulating dormancy in sweet cherry. This was done to lay the foundation for future explorations into the structural and functional aspects of these factors in Prunus species. This study offers significant insights into the functions of ASMT, SNAT, T5H, and TDC genes and sheds light on their roles in Prunus avium. Moreover, it established a robust foundation for further exploration functional characterization of melatonin genes in fruit species.

Comparative genomics reveals the adaptation of ammonia-oxidising Thaumarchaeota to arid soils.

Thaumarchaeota are predominant in oligotrophic habitats such as deserts and arid soils, but their adaptations to these arid conditions are not well understood. In this study, we assembled 23 Thaumarchaeota genomes from arid and semi-arid soils collected from the Inner Mongolia Steppe and the Qinghai-Tibet Plateau. Using a comparative genomics approach, integrated with 614 Thaumarchaeota genomes from public databases, we identified the traits and evolutionary forces that contribute to their adaptations to aridity. Our results showed that the newly assembled genomes represent an early diverging group within the lineage of ammonia-oxidising Thaumarchaeota. While the genomic functions previously identified in arid soil lineages were conserved across terrestrial, shallow-ocean and deep-ocean lineages, several traits likely contribute to Thaumarchaeota's adaptation to aridity. These include chlorite dismutase, arsenate reductase, V-type ATPase and genes dealing with oxidative stresses. The acquisition and loss of traits at the last common ancestor of arid soil lineages may have facilitated the specialisation of Thaumarchaeota in arid soils. Additionally, the acquisition of unique adaptive traits, such as a urea transporter, Ca2+ :H+ antiporter, mannosyl-3-phosphoglycerate synthase and phosphatase, DNA end-binding protein Ku and phage shock protein A, further distinguishes arid soil Thaumarchaeota. This study provides evidence for the adaptations of Thaumarchaeota to arid soil, enhancing our understanding of the nitrogen and carbon cycling driven by Thaumarchaeota in drylands.

Cis-regulatory control of mammalian Trps1 gene expression.

TRPS1 serves as the causative gene for tricho-rhino phalangeal syndrome, known for its craniofacial and skeletal abnormalities. The Trps1 gene encodes a protein that represses Wnt signaling through strong interactions with Wnt signaling inhibitors. The identification of genomic cis-acting regulatory sequences governing Trps1 expression is crucial for understanding its role in embryogenesis. Nevertheless, to date, no investigations have been conducted concerning these aspects of Trps1.To identify deeply conserved noncoding elements (CNEs) within the Trps1 locus, we employed a comparative genomics approach, utilizing slowly evolving fish such as coelacanth and spotted gar. These analyses resulted in the identification of eight CNEs in the intronic region of the Trps1 gene. Functional characterization of these CNEs in zebrafish revealed their regulatory potential in various tissues, including pectoral fins, heart, and pharyngeal arches. RNA in-situ hybridization experiments revealed concordance between the reporter expression pattern induced by the identified set of CNEs and the spatial expression pattern of the trps1 gene in zebrafish. Comparative in vivo data from zebrafish and mice for CNE7/hs919 revealed conserved functions of these enhancers. Each of these eight CNEs was further investigated in cell line-based reporter assays, revealing their repressive potential. Taken together, in vivo and in vitro assays suggest a context-dependent dual functionality for the identified set of Trps1-associated CNE enhancers.This functionally characterized set of CNE-enhancers will contribute to a more comprehensive understanding of the developmental roles of Trps1 and can aid in the identification of noncoding DNA variants associated with human diseases.

Clostridium strain FAM25158, a unique endospore-forming bacterium related to Clostridium tyrobutyricum and isolated from Emmental cheese shows low tolerance to salt.

The genus Clostridium is a large and diverse group of species that can cause food spoilage, including late blowing defect (LBD) in cheese. In this study, we investigated the taxonomic status of strain FAM25158 isolated from Emmental cheese with LBD using a polyphasic taxonomic and comparative genomic approach. A 16S rRNA gene sequence phylogeny suggested affiliation to the Clostridium sensu stricto cluster, with Clostridium tyrobutyricum DSM 2637T being the closest related type strain (99.16% sequence similarity). Average Nucleotide Identity (ANI) analysis revealed that strain FAM25158 is at the species threshold with C. tyrobutyricum, with ANI values ranging from 94.70 to 95.26%, while the digital DNA-DNA hybridization values were below the recommended threshold, suggesting that FAM25158 is significantly different from C. tyrobutyricum at the genomic level. Moreover, comparative genomic analysis between FAM25158 and its four closest C. tyrobutyricum relatives revealed a diversity of metabolic pathways, with FAM25158 differing from other C. tyrobutyricum strains by the presence of genes such as scrA, srcB, and scrK, responsible for sucrose utilization, and the absence of many important functional genes associated with cold and osmolality adaptation, which was further supported by phenotypic analyses. Surprisingly, strain FAM25158 exhibited unique physiologic traits, such as an optimal growth temperature of 30°C, in contrast to its closest relatives, C. tyrobutyricum species with an optimal growth temperature of 37°C. Additionally, the growth of FAM25158 was inhibited at NaCl concentrations higher than 0.5%, a remarkable observation considering its origin from cheese. While the results of this study provide novel information on the genetic content of strain FAM25158, the relationship between its genetic content and the observed phenotype remains a topic requiring further investigation.

A comparative genomic approach to decipher the mutations associated with Nipah viral human isolates from southeast Asia.

Multiple outbreaks over two decades and a high mortality rate have emphasized the Nipah virus (NiV) as a priority research area. The study focuses on identifying the mutational landscape in sequences from NiV human isolates from different geographical regions. Thirty-seven NiV genomes of human samples from Malaysia, Bangladesh, and India were subjected to phylogeny and metagenomic analysis to decipher the genome variability using MEGA11 software and the meta-CATS web server. Using the Single-Likelihood Ancestor Counting method, the synonymous and nonsynonymous mutations among NiV genes were identified. Further, the nonsynonymous variations were used to identify mutations in all the NiV proteins. The NiV isolates were categorized into NiV-M, NiV-B, and NiV-I clades based on phylogenetic analysis. Metagenomic analysis revealed 1636 variations in the noncoding and coding regions of the genomes of the three clades of NiV. Further analysis of nonsynonymous mutations showed the phosphoprotein to be highly mutating, whereas the matrix protein was stable. Deciphering the mutation pattern using a comparative genomics approach for human isolates provided valuable insight into the stability of NiV proteins which can be further used for understanding variations in host-pathogen interaction and developing effective therapeutic measures.

A Critical Appraisal of DNA Transfer from Plants to Parasitic Cyst Nematodes.

Plant-parasitic nematodes are one of the most economically important pests of crops. It is widely accepted that horizontal gene transfer-the natural acquisition of foreign genes in parasitic nematodes-contributes to parasitism. However, an apparent paradox has emerged from horizontal gene transfer analyses: On the one hand, distantly related organisms with very dissimilar genetic structures (i.e. bacteria), and only transient interactions with nematodes as far as we know, dominate the list of putative donors, while on the other hand, considerably more closely related organisms (i.e. the host plant), with similar genetic structure (i.e. introns) and documented long-term associations with nematodes, are rare among the list of putative donors. Given that these nematodes ingest cytoplasm from a living plant cell for several weeks, there seems to be a conspicuous absence of plant-derived cases. Here, we used comparative genomic approaches to evaluate possible plant-derived horizontal gene transfer events in plant parasitic nematodes. Our evidence supports a cautionary message for plant-derived horizontal gene transfer cases in the sugar beet cyst nematode, Heterodera schachtii. We propose a 4-step model for horizontal gene transfer from plant to parasite in order to evaluate why the absence of plant-derived horizontal gene transfer cases is observed. We find that the plant genome is mobilized by the nematode during infection, but that uptake of the said "mobilome" is the first major barrier to horizontal gene transfer from host to nematode. These results provide new insight into our understanding of the prevalence/role of nucleic acid exchange in the arms race between plants and plant parasites.

Multidisciplinary Methods for Screening Toxic Proteins from Phages and Their Potential Molecular Targets.

This chapter presents a comprehensive methodology for the identification, characterization, and functional analyses of potentially toxic hypothetical proteins of unknown function (toxHPUFs) in phages. The methods begin with in vivo toxicity verification of toxHPUFs in bacterial hosts, utilizing conventional drop tests and following growth curves. Computational methods for structural and functional predictions of toxHPUFs are outlined, incorporating the use of tools such as Phyre2, HHpred, and AlphaFold2. To ascertain potential targets, a comparative genomic approach is described using bioinformatics toolkits for sequence alignment and functional annotation. Moreover, steps are provided to predict protein-protein interactions and visualizing these using PyMOL. The culmination of these methods equips researchers with an effective pipeline to identify and analyze toxHPUFs and their potential targets, laying the groundwork for future experimental confirmations.

Comparative Genomics of Clostridioides difficile.

Clostridioides difficile, a Gram-positive spore-forming anaerobic bacterium, has rapidly emerged as the leading cause of nosocomial diarrhoea in hospitals. The availability of large numbers of genome sequences, mainly due to the use of next-generation sequencing methods, has undoubtedly shown their immense advantages in the determination of C. difficile population structure. The implementation of fine-scale comparative genomic approaches has paved the way for global transmission and recurrence studies, as well as more targeted studies, such as the PaLoc or CRISPR/Cas systems. In this chapter, we provide an overview of recent and significant findings on C. difficile using comparative genomic studies with implications for epidemiology, infection control and understanding of the evolution of C. difficile.

Speeding genomic island discovery through systematic design of reference database composition.

Genomic islands (GIs) are mobile genetic elements that integrate site-specifically into bacterial chromosomes, bearing genes that affect phenotypes such as pathogenicity and metabolism. GIs typically occur sporadically among related bacterial strains, enabling comparative genomic approaches to GI identification. For a candidate GI in a query genome, the number of reference genomes with a precise deletion of the GI serves as a support value for the GI. Our comparative software for GI identification was slowed by our original use of large reference genome databases (DBs). Here we explore smaller species-focused DBs. With increasing DB size, recovery of our reliable prophage GI calls reached a plateau, while recovery of less reliable GI calls (FPs) increased rapidly as DB sizes exceeded ~500 genomes; i.e., overlarge DBs can increase FP rates. Paradoxically, relative to prophages, FPs were both more frequently supported only by genomes outside the species and more frequently supported only by genomes inside the species; this may be due to their generally lower support values. Setting a DB size limit for our SMAll Ranked Tailored (SMART) DB design speeded runtime ~65-fold. Strictly intra-species DBs would tend to lower yields of prophages for small species (with few genomes available); simulations with large species showed that this could be partially overcome by reaching outside the species to closely related taxa, without an FP burden. Employing such taxonomic outreach in DB design generated redundancy in the DB set; as few as 2984 DBs were needed to cover all 47894 prokaryotic species. Runtime decreased dramatically with SMART DB design, with only minor losses of prophages. We also describe potential utility in other comparative genomics projects.

Cuticular proteins in codling moth (Cydia pomonella) respond to insecticide and temperature stress

The insect cuticle consists of chitin and cuticular proteins (CPs), which stabilize the body shape and provide an effective physical barrier against the external environment. They are also potential target sites for developing environmentally friendly insect management through the utilization of physiology-based methods. The codling moth, Cydia pomonella, is a pest afflicting fruit orchards worldwide. This study used a comparative genomic approach, whole-genome resequencing, and transcriptome data to understand the role that CPs played in the environmental adaptation of the codling moth. A total of 182 putative CPs were identified in the codling moth genome, which were classified into 12 CP families. 119 CPR genes, including 54 RR-1, 60 RR-2, and 5 RR-3 genes were identified and accounted for 65.4% of the total CPs. Eight and seven gene clusters are formed in RR1 and RR2 subfamily and the ancestor-descendant relationship was explained. Five CPAP genes were highly expressed during the egg stage and exposed to high temperature, which indicated their potential role in aiding codling moth eggs in acclimating to varying external heat conditions. Moreover, six CPs belonging to the CPR and CPLCP families were identified in association with insecticide resistance by population resequencing. Their expression levels increased after exposure to insecticides, suggesting they might be involved in codling moth resistance to the insecticides azinphos-methyl or deltamethrin. Our results provide insight into the evolution of codling moth CPs and their association with high temperature adaptation and insecticide resistance, and provide an additional information required for further analysis of CPs in environmental adaptation.

Comparative genomics reveals evolutionary drivers of the dietary shift in Hemiptera.

Hemiptera insects exhibit a close relationship to plants and demonstrate a diverse range of dietary preferences, encompassing phytophagy as the predominant feeding habit while a minority engages in carnivorous or haematophagous behaviour. To counteract the challenges posed by phytophagous insects, plants have developed an array of toxic compounds, causing significant evolutionary selection pressure on these insects. In this study, we employed a comparative genomics approach to analyse the expansion and contraction of gene families specific to phytophagous insect lineages, along with their adaptive evolutionary traits, utilising representative species from the Hemiptera order. Our investigation revealed substantial expansions of gene families within the phytophagous lineages, especially in the Pentatomomorpha branch represented by Oncopeltus fasciatus and Riptortus pedestris. Notably, these expansions of gene families encoding enzymes are potentially involved in hemipteran-plant interactions. Moreover, the adaptive evolutionary analysis of these lineages revealed a higher prevalence of adaptively evolved genes in the Pentatomomorpha branch. The observed branch-specific gene expansions and adaptive evolution likely contribute significantly to the diversification of species within Hemiptera. These results help enhance our understanding of the genomic characteristics of the evolution of different feeding habits in hemipteran insects.