Short-sequence DNA repeats in prokaryotic genomes

Mechanisms and implications of phenotypic switching in bacterial pathogens

Bacteria encounter various stressful conditions within a variety of dynamic environments, which they must overcome for survival. One way they achieve this is by developing phenotypic heterogeneity to introduce diversity within their population. Such distinct subpopulations can arise through endogenous fluctuations in regulatory components, wherein bacteria can express diverse phenotypes and switch between them, sometimes in a heritable and reversible manner. This switching may also lead to antigenic variation, enabling pathogenic bacteria to evade the host immune response. Therefore, phenotypic heterogeneity plays a significant role in microbial pathogenesis, immune evasion, antibiotic resistance, host niche tissue establishment, and environmental persistence. This heterogeneity can result from stochastic and responsive switches, as well as various genetic and epigenetic mechanisms. The development of phenotypic heterogeneity may create clonal populations that differ in their level of virulence, contribute to the formation of biofilms, and allow for antibiotic persistence within select morphological variants. This review delves into the current understanding of the molecular switching mechanisms underlying phenotypic heterogeneity, highlighting their roles in establishing infections caused by select bacterial pathogens.

Synthetic strategy for the production of electrically polarized polyvinylidene fluoride-trifluoroethylene-co-polymer osseo-functionalized with hydroxyapatite scaffold

The physiological mechanism of bone tissue regeneration is intricately organized and involves several cell types, intracellular, and extracellular molecular signaling networks. To overcome the drawbacks of autografts and allografts, a number of synthetically produced scaffolds have been manufactured by integrating ceramics, polymers, and their hybrid-composites. Considering the fact that natural bone is composed primarily of collagen and hydroxyapatite, ceramic-polymer composite materials seem to be the most viable alternative to bone implants. Here, in this experimental study, copolymer PVDF-TrFE has been amalgamated with HA ceramics to produce composite scaffolds as bone implants. In order to fabricate PVDF-TrFE-HA (polyvinylidene fluoride-trifluoroethylene-hydroxyapatite) composite scaffolds, solvent casting-particulate leaching technique was devised. Two scaffold specimens were produced, with different PVDF-TrFE and HA molar ratios (70:30 and 50:50), and then electrically polarized to observe the subsequent polarization impact on the tissue growth and the suppression of bacterial cell proliferation. Both the specimens underwent characterization to analyze their biocompatibility and bactericidal activities. The bacterial culture of Pseudomonas aeruginosa (P. aeruginosa) and Staphylococcus aureus (S. aureus) bacteria on the composites was studied to understand the antibacterial characteristics. Moreover, MG63 cells cultured on these as-formed composites provided information about osteogenesis. Improved osteogenesis and antibacterial efficacy were observed on both the composites. However, the composite with 70 wt% PVDF-TrFE and 30 wt% HA showed a higher bactericidal effect as well as osteogenesis. It was found that PVDF-TrFE-HA-based biomaterials have the potential for bone tissue engineering applications.

We are doing it wrong: Putting homology before phylogeny in cyanobacterial taxonomy

The rapid expansion of whole genome sequencing in bacterial taxonomy has revealed deep evolutionary relationships and speciation signals, but assembly methods often miss true nucleotide diversity in the ribosomal operons. Though it lacks sufficient phylogenetic signal at the species level, the 16S ribosomal RNA gene is still much used in bacterial taxonomy. In cyanobacterial taxonomy, comparisons of 16S-23S Internal Transcribed Spacer (ITS) regions are used to bridge this information gap. Although ITS rRNA region analyses are routinely being used to identify species, researchers often do not identify orthologous operons, which leads to improper comparisons. No method for delineating orthologous operon copies from paralogous ones has been established. A new method for recognizing orthologous ribosomal operons by quantifying the conserved paired nucleotides in a helical domain of the ITS, has been developed. The D1' Index quantifies differences in the ratio of pyrimidines to purines in paired nucleotide sequences of this helix. Comparing 111 operon sequences from 89 strains of Brasilonema, four orthologous operon types were identified. Plotting D1' Index values against the length of helices produced clear separation of orthologs. Most orthologous operons in this study were observed both with and without tRNA genes present. We hypothesize that genomic rearrangement, not gene duplication, is responsible for the variation among orthologs. This new method will allow cyanobacterial taxonomists to utilize ITS rRNA region data more correctly, preventing erroneous taxonomic hypotheses. Moreover, this work could assist genomicists in identifying and preserving evident sequence variability in ribosomal operons, which is an important proxy for evolution in prokaryotes.

Comparative analysis of multiple consensus genomes of the same strain of Marek's disease virus reveals intrastrain variation

Current strategies to understand the molecular basis of Marek's disease virus (MDV) virulence primarily consist of cataloging divergent nucleotides between strains with different phenotypes. However, most comparative genomic studies of MDV rely on previously published consensus genomes despite the confirmed existence of MDV strains as mixed viral populations. To assess the reliability of interstrain genomic comparisons relying on published consensus genomes of MDV, we obtained two additional consensus genomes of vaccine strain CVI988 (Rispens) and two additional consensus genomes of the very virulent strain Md5 by sequencing viral stocks and cultured field isolates. In conjunction with the published genomes of CVI988 and Md5, this allowed us to perform three-way comparisons between multiple consensus genomes of the same strain. We found that consensus genomes of CVI988 can vary in as many as 236 positions involving 13 open reading frames (ORFs). By contrast, we found that Md5 genomes varied only in 11 positions involving a single ORF. Notably, we were able to identify 3 single-nucleotide polymorphisms (SNPs) in the unique long region and 16 SNPs in the unique short (US) region of CVI988GenBank.BAC that were not present in either CVI988Pirbright.lab or CVI988USDA.PA.field. Recombination analyses of field strains previously described as natural recombinants of CVI988 yielded no evidence of crossover events in the US region when either CVI988Pirbright.lab or CVI988USDA.PA.field were used to represent CVI988 instead of CVI988GenBank.BAC. We were also able to confirm that both CVI988 and Md5 populations were mixed, exhibiting a total of 29 and 27 high-confidence minor variant positions, respectively. However, we did not find any evidence of minor variants in the positions corresponding to the 19 SNPs in the unique regions of CVI988GenBank.BAC. Taken together, our findings suggest that continued reliance on the same published consensus genome of CVI988 may have led to an overestimation of genomic divergence between CVI988 and virulent strains and that multiple consensus genomes per strain may be necessary to ensure the accuracy of interstrain genomic comparisons.

Belt and braces: Two escape ways to maintain the cassette reservoir of large chromosomal integrons

Integrons are adaptive devices that capture, stockpile, shuffle and express gene cassettes thereby sampling combinatorial phenotypic diversity. Some integrons called sedentary chromosomal integrons (SCIs) can be massive structures containing hundreds of cassettes. Since most of these cassettes are non-expressed, it is not clear how they remain stable over long evolutionary timescales. Recently, it was found that the experimental inversion of the SCI of Vibrio cholerae led to a dramatic increase of the cassette excision rate associated with a fitness defect. Here, we question the evolutionary sustainability of this apparently counter selected genetic context. Through experimental evolution, we find that the integrase is rapidly inactivated and that the inverted SCI can recover its original orientation by homologous recombination between two insertion sequences (ISs) present in the array. These two outcomes of SCI inversion restore the normal growth and prevent the loss of cassettes, enabling SCIs to retain their roles as reservoirs of functions. These results illustrate a nice interplay between gene orientation, genome rearrangement, bacterial fitness and demonstrate how integrons can benefit from their embedded ISs.

Pathogenomes of Shiga Toxin Positive and Negative Escherichia coli O157:H7 Strains TT12A and TT12B: Comprehensive Phylogenomic Analysis Using Closed Genomes

Shiga toxin-producing Escherichia coli are zoonotic pathogens that cause food-borne human disease. Among these, the O157:H7 serotype has evolved from an enteropathogenic O55:H7 ancestor through the displacement of the somatic gene cluster and recurrent toxigenic conversion by Shiga toxin-converting bacteriophages. However, atypical strains that lack the Shiga toxin, the characteristic virulence hallmark, are circulating in this lineage. For this study, we analyzed the pathogenome and virulence inventories of the stx+ strain, TT12A, isolated from a patient with hemorrhagic colitis, and its respective co-isolated stx- strain, TT12B. Sequencing the genomes to closure proved critical to the cataloguing of subtle strain differentiating sequence and structural polymorphisms at a high-level of phylogenetic accuracy and resolution. Phylogenomic profiling revealed SNP and MLST profiles similar to the near clonal outbreak isolates. Their prophage inventories, however, were notably different. The attenuated atypical non-shigatoxigenic status of TT12B is explained by the absence of both the ΦStx1a- and ΦStx2a-prophages carried by TT12A, and we also recorded further alterations in the non-Stx prophage complement. Phenotypic characterization indicated that culture growth was directly impacted by the strains' distinct lytic phage complement. Altogether, our phylogenomic and phenotypic analyses show that these intimately related isogenic strains are on divergent Stx(+/stx-) evolutionary paths.

CRISPR genotyping methods: Tracing the evolution from spoligotyping to machine learning

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems provide prokaryotes with adaptive immunity defenses against foreign genetic invaders. The identification of CRISPR-Cas function is among the most impactful discoveries of recent decades that have shaped the development of genome editing in various organisms paving the way for a plethora of promising applications in biotechnology and health. Even before the discovery of CRISPR-Cas biological role, the particular structure of CRISPR loci has been explored for epidemiological genotyping of bacterial pathogens. CRISPR-Cas loci are arranged in CRISPR arrays of mostly identical direct repeats intercalated with invader-derived spacers and an operon of cas genes encoding the Cas protein components. Each small CRISPR RNA (crRNA) encoded within the CRISPR array constitutes a key functional unit of this RNA-based CRISPR-Cas defense system guiding the Cas effector proteins toward the foreign nucleic acids for their destruction. The information acquired from prior invader encounters and stored within CRISPR arrays turns out to be extremely valuable in tracing the microevolution and epidemiology of major bacterial pathogens. We review here the history of CRISPR-based typing strategies highlighting the first PCR-based methods that have set the stage for recent developments of high-throughput sequencing and machine learning-based approaches. A great amount of whole genome sequencing and metagenomic data accumulated in recent years opens up new avenues for combining experimental and computational approaches of high-resolution CRISPR-based typing.

LUSTR: a new customizable tool for calling genome-wide germline and somatic short tandem repeat variants

BACKGROUND

Short tandem repeats (STRs) are widely distributed across the human genome and are associated with numerous neurological disorders. However, the extent that STRs contribute to disease is likely under-estimated because of the challenges calling these variants in short read next generation sequencing data. Several computational tools have been developed for STR variant calling, but none fully address all of the complexities associated with this variant class.

RESULTS

Here we introduce LUSTR which is designed to address some of the challenges associated with STR variant calling by enabling more flexibility in defining STR loci, allowing for customizable modules to tailor analyses, and expanding the capability to call somatic and multiallelic STR variants. LUSTR is a user-friendly and easily customizable tool for targeted or unbiased genome-wide STR variant screening that can use either predefined or novel genome builds. Using both simulated and real data sets, we demonstrated that LUSTR accurately infers germline and somatic STR expansions in individuals with and without diseases.

CONCLUSIONS

LUSTR offers a powerful and user-friendly approach that allows for the identification of STR variants and can facilitate more comprehensive studies evaluating the role of pathogenic STR variants across human diseases.

Comparative Analysis of Multiple Consensus Genomes of the Same Strain of Marek's Disease Virus Reveals Intrastrain Variation

Current strategies to understand the molecular basis of Marek's disease virus (MDV) virulence primarily consist of cataloguing divergent nucleotides between strains with different phenotypes. However, each MDV strain is typically represented by a single consensus genome despite the confirmed existence of mixed viral populations. To assess the reliability of single-consensus interstrain genomic comparisons, we obtained two additional consensus genomes of vaccine strain CVI988 (Rispens) and two additional consensus genomes of the very virulent strain Md5 by sequencing viral stocks and cultured field isolates. In conjunction with the published genomes of CVI988 and Md5, this allowed us to perform 3-way comparisons between consensus genomes of the same strain. We found that consensus genomes of CVI988 can vary in as many as 236 positions involving 13 open reading frames (ORFs). In contrast, we found that Md5 genomes varied only in 11 positions involving a single ORF. Phylogenomic analyses showed all three Md5 consensus genomes clustered closely together, while also showing that CVI988 GenBank.BAC diverged from CVI988 Pirbright.lab and CVI988 USDA.PA.field . Comparison of CVI988 consensus genomes revealed 19 SNPs in the unique regions of CVI988 GenBank.BAC that were not present in either CVI988 Pirbright.lab or CVI988 USDA.PA.field . Finally, we evaluated the genomic heterogeneity of CVI988 and Md5 populations by identifying positions with >2% read support for alternative alleles in two ultra-deeply sequenced samples. We were able to confirm that both populations of CVI988 and Md5 were mixed, exhibiting a total of 29 and 27 high-confidence minor variant positions, respectively. We did not find any evidence of minor variants in the positions corresponding to the 19 SNPs in the unique regions of CVI988 GenBank.BAC . Taken together, our findings confirm that consensus genomes of the same strain of MDV can vary and suggest that multiple consensus genomes per strain are needed in order to maximize the accuracy of interstrain genomic comparisons.

Breaking the Ice: A Review of Phages in Polar Ecosystems

Bacteriophages, or phages, are viruses that infect and replicate within bacterial hosts, playing a significant role in regulating microbial populations and ecosystem dynamics. However, phages from extreme environments such as polar regions remain relatively understudied due to challenges such as restricted ecosystem access and low biomass. Understanding the diversity, structure, and functions of polar phages is crucial for advancing our knowledge of the microbial ecology and biogeochemistry of these environments. In this review, we will explore the current state of knowledge on phages from the Arctic and Antarctic, focusing on insights gained from -omic studies, phage isolation, and virus-like particle abundance data. Metagenomic studies of polar environments have revealed a high diversity of phages with unique genetic characteristics, providing insights into their evolutionary and ecological roles. Phage isolation studies have identified novel phage-host interactions and contributed to the discovery of new phage species. Virus-like particle abundance and lysis rate data, on the other hand, have highlighted the importance of phages in regulating bacterial populations and nutrient cycling in polar environments. Overall, this review aims to provide a comprehensive overview of the current state of knowledge about polar phages, and by synthesizing these different sources of information, we can better understand the diversity, dynamics, and functions of polar phages in the context of ongoing climate change, which will help to predict how polar ecosystems and residing phages may respond to future environmental perturbations.

Validating the Utility of Multilocus Variable Number Tandem-repeat Analysis (MLVA) as a Subtyping Strategy to Monitor Listeria monocytogenes In-built Food Processing Environments

Listeria monocytogenes is a serious human pathogen and an enduring challenge to control for the ready-to-eat food processing industry. Cost-effective tools that can be deployed by commercial or in-house laboratories to rapidly investigate and resolve contamination events in the built food processing environment are of value to the food industry. Multilocus variable number tandem-repeat analysis (MLVA) is a molecular subtyping method, which along with other same-generation methods such as pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST) is being superseded in disease tracking and outbreak investigations by whole-genome sequencing (WGS). In this paper, it is demonstrated that MLVA can continue to play a valuable role as a valid, fast, simple, and cost-effective method to identify and track Listeria monocytogenes subtypes in factory environments, with the method being highly congruent with MLST. Although MLVA does not have the discriminatory power of WGS to identify truly persistent clones, with careful interpretation of results alongside isolate metadata, it remains a powerful tool in situations and locations where WGS may not be readily available to food business operators.

Evolution of genetic architecture and gene regulation in biphenyl/PCB-degrading bacteria

A variety of bacteria in the environment can utilize xenobiotic compounds as a source of carbon and energy. The bacterial strains degrading xenobiotics are suitable models to investigate the adaptation and evolutionary processes of bacteria because they appear to have emerged relatively soon after the release of these compounds into the natural environment. Analyses of bacterial genome sequences indicate that horizontal gene transfer (HGT) is the most important contributor to the bacterial evolution of genetic architecture. Further, host bacteria that can use energy effectively by controlling the expression of organized gene clusters involved in xenobiotic degradation will have a survival advantage in harsh xenobiotic-rich environments. In this review, we summarize the current understanding of evolutionary mechanisms operative in bacteria, with a focus on biphenyl/PCB-degrading bacteria. We then discuss metagenomic approaches that are useful for such investigation.

Clonal isolates of Treponema pallidum subsp. pallidum Nichols provide evidence for the occurrence of microevolution during experimental rabbit infection and in vitro culture

The recent development of a system for long-term in vitro culture of the syphilis spirochete, Treponema pallidum subsp. pallidum, has introduced the possibility of detailed genetic analysis of this bacterium. In this study, the in vitro culture system was used to isolate and characterize clonal populations of T. pallidum subsp. pallidum Nichols, the most widely studied strain. In limiting dilutions experiments, it was possible to establish cultures with inocula as low as 0.5 T. pallidum per well despite the long generation time (~35 to 40 hours) of this organism. Six Nichols strain clones isolated by limiting dilution were characterized in detail. All clones exhibited indistinguishable morphology and motility, highly similar in vitro multiplication rates, and comparable infectivity in the rabbit model (ID50 ≤ 100 bacteria). Genomic sequencing revealed sequence heterogeneity in the form of insertions or deletions at 5 sites, single nucleotide variations at 20 sites, and polynucleotide (polyG/C) tract length differences at 22 locations. Genomic sequences of the uncloned Nichols strain preparations propagated in rabbits or in vitro cultures exhibited substantial heterogeneity at these locations, indicating coexistence of many varied 'clonotypes' within these populations. Nearly all genetic variations were specific for the Nichols strain and were not detected in the >280 T. pallidum genomic sequences that are currently available. We hypothesize that these Nichols strain-specific sequence variations arose independently either during human infection or within the 110 years since the strain's initial isolation, and thus represent examples of microevolution and divergence.

Molecular Epidemiology of Ralstonia pseudosolanacearum Phylotype I Strains in the Southwest Indian Ocean Region and Their Relatedness to African Strains

The increasing requirement for developing tools enabling fine strain traceability responsible for epidemics is tightly linked with the need to understand factors shaping pathogen populations and their environmental interactions. Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is one of the most important plant diseases in tropical and subtropical regions. Sadly, little, outdated, or no information on its epidemiology is reported in the literature, although alarming outbreaks are regularly reported as disasters. A large set of phylotype I isolates (n = 2,608) was retrieved from diseased plants in fields across the Southwest Indian Ocean (SWIO) and Africa. This collection enabled further assessment of the epidemiological discriminating power of the previously published RS1-MLVA14 scheme. Thirteen markers were validated and characterized as not equally informative. Most had little infra-sequevar polymorphism, and their performance depended on the sequevar. Strong correlation was found with a previous multilocus sequence typing scheme. However, 2 to 3% of sequevars were not correctly assigned through endoglucanase gene sequence. Discriminant analysis of principal components (DAPC) revealed four groups with strong phylogenetic relatedness to sequevars 31, 33, and 18. Phylotype I-31 isolates were highly prevalent in the SWIO and Africa, but their dissemination pathways remain unclear. Tanzania and Mauritius showed the greatest diversity of RSSC strains, as the four DAPC groups were retrieved. Mauritius was the sole territory harboring a vast phylogenetic diversity and all DAPC groups. More research is still needed to understand the high prevalence of phylotype I-31 at such a large geographic scale.

A reversible mutation in a genomic hotspot saves bacterial swarms from extinction

Microbial adaptation to changing environmental conditions is frequently mediated by hypermutable sequences. Here we demonstrate that such a hypermutable hotspot within a gene encoding a flagellar unit of Paenibacillus glucanolyticus generated spontaneous non-swarming mutants with increased stress resistance. These mutants, which survived conditions that eliminated wild-type cultures, could be carried by their swarming siblings when the colony spread, consequently increasing their numbers at the spreading edge. Of interest, the hypermutable nature of the aforementioned sequence enabled the non-swarming mutants to serve as "seeds" for a new generation of wild-type cells through reversion of the mutation. Using a mathematical model, we examined the survival dynamics of P. glucanolyticus colonies under fluctuating environments. Our experimental and theoretical results suggest that the non-swarming, stress-resistant mutants can save the colony from extinction. Notably, we identified this hypermutable sequence in flagellar genes of additional Paenibacillus species, suggesting that this phenomenon could be wide-spread and ecologically important.

A minisatellite-based MLVA for deciphering the global epidemiology of the bacterial cassava pathogen Xanthomonas phaseoli pv. manihotis

Cassava Bacterial Blight (CBB) is a destructive disease widely distributed in the different areas where this crop is grown. Populations studies have been performed at local and national scales revealing a geographical genetic structure with temporal variations. A global epidemiology analysis of its causal agent Xanthomonas phaseoli pv. manihotis (Xpm) is needed to better understand the expansion of the disease for improving the monitoring of CBB. We targeted new tandem repeat (TR) loci with large repeat units, i.e. minisatellites, that we multiplexed in a scheme of Multi-Locus Variable number of TR Analysis (MLVA-8). This genotyping scheme separated 31 multilocus haplotypes in three clusters of single-locus variants and a singleton within a worldwide collection of 93 Xpm strains isolated over a period of fifty years. The major MLVA-8 cluster 1 grouped strains originating from all countries, except the unique Chinese strain. On the contrary, all the Xpm strains genotyped using the previously developed MLVA-14 microsatellite scheme were separated as unique haplotypes. We further propose an MLVA-12 scheme which takes advantage of combining TR loci with different mutation rates: the eight minisatellites and four faster evolving microsatellite markers, for global epidemiological surveillance. This MLVA-12 scheme identified 78 haplotypes and separated most of the strains in groups of double-locus variants (DLV) supporting some phylogenetic relationships. DLV groups were subdivided into closely related clusters of strains most often sharing the same geographical origin and isolated over a short period, supporting epidemiological relationships. The main MLVA-12 DLV group#1 was composed by strains from South America and all the African strains. The MLVA-12 scheme combining both minisatellite and microsatellite loci with different discriminatory power is expected to increase the accuracy of the phylogenetic signal and to minimize the homoplasy effects. Further investigation of the global epidemiology of Xpm will be helpful for a better control of CBB worldwide.

In vitro effects of diazabicyclooctane β-lactamase inhibitors relebactam and nacubactam against three subspecies of Mycobacterium abscessus complex

BACKGROUND

Mycobacterium abscessus complex (MABC) pulmonary disease is notoriously difficult to treat due to intrinsic resistance to many common antibiotics. MABC is β-lactam-resistant as it produces class A β-lactamases, such as bla_Mab, which are inhibited by diazabicyclooctane (DBO) β-lactamase inhibitors.

OBJECTIVES

To investigate the microbiological effects of the combination of β-lactam and DBO β-lactamase inhibitors (relebactam and nacubactam) against MABC and determine if the effects are associated with the MABC subspecies and colony morphotype.

METHODS

The antimicrobial susceptibility of three type strains and 20 clinical isolates of MABC to the combination of seven β-lactams with relebactam or nacubactam was evaluated using broth microdilution checkerboard assays. For these strains, expression levels of bla_Mab were assessed using quantitative real-time polymerase chain reaction and genotypic diversity was evaluated using 18-locus variable number tandem repeat assay.

RESULTS

Relebactam and nacubactam lowered the minimum inhibitory concentrations of β-lactams, particularly imipenem, meropenem, and tebipenem, against MABC. There was no difference in efficacy of combination treatment between three subspecies, but rough morphotypes tended to be less susceptible than smooth morphotypes. There were no differences in bla_Mab expression levels and genotypic diversity between the morphotypes.

CONCLUSIONS

The combination of β-lactam with relebactam or nacubactam improved the efficacy of β-lactams against all MABC subspecies, but higher concentrations of β-lactams were needed for rough morphotypes.

From genome structure to function: insights into structural variation in microbiology

Structural variation in bacterial genomes is an important evolutionary driver. Genomic rearrangements, such as inversions, duplications, and insertions, can regulate gene expression and promote niche adaptation. Importantly, many of these variations are reversible and preprogrammed to generate heterogeneity. While many tools have been developed to detect structural variation in eukaryotic genomes, variation in bacterial genomes and metagenomes remains understudied. However, recent advances in genome sequencing technology and the development of new bioinformatic pipelines hold promise in further understanding microbial genomics.

Development of a Novel Diagnostic Tool for Cercospora Species Based on BOX-PCR System

The genus Cercospora contains many devastating plant pathogens linked to leaf spot diseases afflicting various plants. Identification of Cercospora species based on morphology or host plant association has proven unreliable due to simple morphology and wide host range in many cases; hence, multi-gene DNA sequence data are essential for accurate species identification. Considering the complexity and cost involved in application of multi-locus DNA phylogenetic approaches for species delineation in Cercospora; rapid and cost-effective methods are urgently needed for species recognition. In this study, we applied rep-PCR (repetitive-sequence based polymerase chain reaction) fingerprinting methods referred to as BOX-PCR to differentiate species of Cercospora. Cluster analysis of the banding patterns of 52 Cercospora strains indicated the ability of BOX-PCR technique using BOXA1R primer to generate species-specific DNA fingerprints from all the tested strains. Since this technique was able to discriminate between all the 20 examined Cercospora species during this study, which corresponded well to the species identified based on multi-gene DNA sequence data, our findings revealed the efficiency of BOX-PCR system as a suitable complementary method for molecular identification of the genus Cercospora at species level.

Assays for Identification and Differentiation of Brucella Species: A Review

Brucellosis is one of the most important and widespread bacterial zoonoses worldwide. Cases are reported annually across the range of known infectious species of the genus Brucella. Globally, Brucella melitensis, primarily hosted by domestic sheep and goats, affects large proportions of livestock herds, and frequently spills over into humans. While some species, such as Brucella abortus, are well controlled in livestock in areas of North America, the Greater Yellowstone Ecosystem supports the species in native wild ungulates with occasional spillover to livestock. Elsewhere in North America, other Brucella species still infect domestic dogs and feral swine, with some associated human cases. Brucella spp. patterns vary across space globally with B. abortus and B. melitensis the most important for livestock control. A myriad of other species within the genus infect a wide range of marine mammals, wildlife, rodents, and even frogs. Infection in humans from these others varies with geography and bacterial species. Control in humans is primarily achieved through livestock vaccination and culling and requires accurate and rapid species confirmation; vaccination is Brucella spp.-specific and typically targets single livestock species for distribution. Traditional bacteriology methods are slow (some media can take up to 21 days for bacterial growth) and often lack the specificity of molecular techniques. Here, we summarize the molecular techniques for confirming and identifying specific Brucella species and provide recommendations for selecting the appropriate methods based on need, sensitivity, and laboratory capabilities/technology. As vaccination/culling approaches are costly and logistically challenging, proper diagnostics and species identification are critical tools for targeting surveillance and control.

Multiple-Locus Variable-Number Tandem Repeat Analysis Genotypes of Listeria monocytogenes Isolated from Farms, Abattoirs, and Retail in Gauteng Province, South Africa

ABSTRACT

The use of multiple-locus variable-number analysis (MLVA) of tandem repeats (TRs) for subtyping Listeria monocytogenes has proven to be reliable and fast. This study determined the MLVA genotypes of 60 isolates of L. monocytogenes recovered from cattle farms, abattoirs, and retail outlets in Gauteng province, South Africa. The distribution of the 60 L. monocytogenes isolates analyzed by type of sample was as follows: raw beef (28, 46.7%), ready-to-eat beef products (9, 15.0%), beef carcass swabs (9, 15.0%), cattle environment (6, 10.0%), and cattle feces (8, 13.3%). The serogroups of the isolates were determined using PCR and the MLVA genotypes based on six selected loci. The frequency of the 60 serogroups detected was as follows: 1/2a-3a (IIa) (27, 45.0%); 4b-4d-4e (1Vb) (24, 40.0%); 1/2c-3c (IIc) (8, 13.3%); and 1/2b-3b (IIb) (1, 1.7%). MLVA successfully clustered genetically related isolates and differentiated nonrelated isolates, irrespective of their sources, sample types, and serogroups, as demonstrated by 16 MLVA pattern types detected. For serogroup 4b-4d-4e (IVb), there was no variation in TRs LM-TR2, LM-TR4, and LM-TR6, which each contained only one allele (02, 00, and 93, respectively). However, across the sources and sample types of isolates, there was variation in serogroup 4b-4d-4e (IVb): LM-TR1 contained 00, 03, and 05; LM-TR3 contained 14, 20, and 22; and LM-TR5 contained 14, 21, and 25. Similar patterns of variation in the TRs were detected in the other serogroups (1/2a-3a, 1/2b-3b, and 1/2c-3c). BioNumeric data analysis identified at least five types in Gauteng province. MLVA epidemiologically clustered the related isolates and differentiated unrelated isolates.

HIGHLIGHTS

Identification and validation of SSR markers for Xanthomonas axonopodis pv. punicae an incitant of bacterial blight of pomegranate

This study reports genome wide characterization and development of first set of microsatellite markers through in silico analysis of eight sequenced Xanthomonas axonopodis pv. punicae strains available in the public database. SSR survey resulted in identification of ~ 4638 perfect SSRs, with mean marker frequency 901 SSRs/Mb and densitiy of 11,006 bp/Mb aross the eight genomes. Frequency distribution graphs revealed hexa-nucleotide repeats were more prominent fowllowed by tri-, tetra-, di- and penta-nucleotides in the analysed genomes. We desinged 2927 SSR primers that are specific to the strain LMG 859 and ePCR confirmed on seven other Xap genomes. This resulted in identification of 542 informative SSRs that are producing single amplicons, from which 66 primers were successfully validated through wet lab experiments on eight Xap isolates of pomegranate. Furthermore, utility of these SSRs were demostrated by analysing molecular diversity among 22 Xap isolates using 20 Xap_SSR primers. SSRs revealed moderate genetic diversity among Xap isolates (61%) and grouped 11 isolates that are repersenting six different states into one cluster. This proved the earlier evidence of wider spread of ST3 type Xap acoss India using Multi locus Sequence Typing (MLST) technique. In summary, Xap_SSR will serve as powerful genomics tools that would helps in monitoring of population dynamics, taxonomy, epidomology and quarantine aspects in bacterial blight pathogen through development of microsatellite based Multilocus Variable number of Tandem repeat analysis (MLVA) in future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03209-z.

Population Diversity of 'Candidatus Liberibacter asiaticus' and Diaphorina citri in Sichuan: A Case Study for Huanglongbing Monitoring and Interception

Citrus huanglongbing (HLB) is present in 10 provinces in China and is associated with 'Candidatus Liberibacter asiaticus' (CLas), which is transmitted by the Asian citrus psyllid (Diaphorina citri, ACP). To date, HLB and ACP have expanded to Yibin city of Sichuan Province, posing an imminent threat to the citrus belt of the upper and middle reaches of the Yangtze River, an important late-maturing citrus-producing area in China. To understand the epidemiological route of CLas and ACP in newly invaded regions of Sichuan and thereby better establish an HLB interception zone ranging from Leibo to Yibin, we evaluated the molecular variability of 19 CLas draft genomes from citrus or dodder (Cuscuta campestris). They include three type-specific prophage loci, three variable number tandem repeat loci, a miniature inverted-repeat transposable element, and population diversity of 44 ACP mitochondrial genomes. The results indicated that CLas isolates in the newly invaded area (Pingshan) were more diverse than those in the HLB endemic areas (Leibo and Ningnan). Phylogenetic analysis based on mitochondrial genomes demonstrated that ACPs in Leibo, Pingshan, and Xuzhou (rural areas) represent a new mitochondrial group (MG4), distinguished by the three unique single-nucleotide polymorphisms in cox1, nad4, and cytb. However, the ACPs sampled from the urban areas of Cuiping and Xuzhou belonged to the southeastern China group (MG2-1). Altogether, our study revealed multiple sources of ACP and CLas in the HLB interception zone and proposed their transmission route. This study contributes to the formulation of precise HLB prevention and control strategies in the HLB interception zone in Sichuan and could be useful for HLB management efforts in other regions.

Using unique ORFan genes as strain-specific identifiers for Escherichia coli

BACKGROUND

Bacterial identification at the strain level is a much-needed, but arduous and challenging task. This study aimed to develop a method for identifying and differentiating individual strains among multiple strains of the same bacterial species. The set used for testing the method consisted of 17 Escherichia coli strains picked from a collection of strains isolated in Germany, Spain, the United Kingdom and Vietnam from humans, cattle, swine, wild boars, and chickens. We targeted unique or rare ORFan genes to address the problem of selective and specific strain identification. These ORFan genes, exclusive to each strain, served as templates for developing strain-specific primers.

RESULTS

Most of the experimental strains (14 out of 17) possessed unique ORFan genes that were used to develop strain-specific primers. The remaining three strains were identified by combining a PCR for a rare gene with a selection step for isolating the experimental strains. Multiplex PCR allowed the successful identification of the strains both in vitro in spiked faecal material in addition to in vivo after experimental infections of pigs and recovery of bacteria from faecal material. In addition, primers for qPCR were also developed and quantitative readout from faecal samples after experimental infection was also possible.

CONCLUSIONS

The method described in this manuscript using strain-specific unique genes to identify single strains in a mixture of strains proved itself efficient and reliable in detecting and following individual strains both in vitro and in vivo, representing a fast and inexpensive alternative to more costly methods.

CRISPRs in the human genome are differentially expressed between malignant and normal adjacent to tumor tissue

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) have been identified in bacteria, archaea and mitochondria of plants, but not in eukaryotes. Here, we report the discovery of 12,572 putative CRISPRs randomly distributed across the human chromosomes, which we termed hCRISPRs. By using available transcriptome datasets, we demonstrate that hCRISPRs are distinctively expressed as small non-coding RNAs (sncRNAs) in cell lines and human tissues. Moreover, expression patterns thereof enabled us to distinguish normal from malignant tissues. In prostate cancer, we confirmed the differential hCRISPR expression between normal adjacent and malignant primary prostate tissue by RT-qPCR and demonstrate that the SHERLOCK and DETECTR dipstick tools are suitable to detect these sncRNAs. We anticipate that the discovery of CRISPRs in the human genome can be further exploited for diagnostic purposes in cancer and other medical conditions, which certainly will lead to the development of point-of-care tests based on the differential expression of the hCRISPRs.

Discovery by metagenomics of a functional tandem repeat sequence that controls gene expression in bacteria

The ability to degrade exogenous compounds is acquired by adaptive processes of microorganisms when they are exposed to compounds that are foreign to their existing enzyme systems. Previously, we reported that simultaneous point mutations and mobile genetic elements cause the evolution and optimization of the degradation systems for aromatic compounds. In the present study, we propose another element with this role—tandem repeats. The novel metagenomic tandem repeat (MTR) sequence T(G/A)ACATG(A/C)T was identified in the 5′-untranslated regions of catechol 2,3-dioxygenase (C23O)-encoding genes by metagenomic analysis. Recombinant Escherichia coli carrying a C23O gene with various numbers of MTRs exhibited increased C23O protein expression and enzyme activity compared with cells expressing the C23O gene without MTRs. Real-time reverse transcription-PCR showed that changes in the numbers of MTRs affected the levels of detectable C23O mRNA in the E. coli host. Furthermore, the mRNAs transcribed from C23O genes containing various numbers of MTRs had longer half-lives than those transcribed from a C23O gene without MTRs. Thus, MTRs would affect the translation efficiency of the gene expression system. MTRs may change the expression levels of their downstream genes for adaptation to a fluctuating environment.

Sequencing Independent Molecular Typing of Staphylococcus aureus Isolates: Approach for Infection Control and Clonal Characterization

Staphylococcus aureus is a major bacterial human pathogen that causes a wide variety of clinical manifestations. The main aim of the presented study was to determine and optimize a novel sequencing independent approach that enables molecular typing of S. aureus isolates and elucidates the transmission of emergent clones between patients. In total, 987 S. aureus isolates including both methicillin-resistant S. aureus (MRSA) and methicillin-sensitive S. aureus (MSSA) isolates were used to evaluate the novel typing approach combining high-resolution melting (HRM) analysis of multilocus sequence typing (MLST) genes (mini-MLST) and spa gene (spa-HRM). The novel approach's discriminatory ability was evaluated by whole-genome sequencing (WGS). The clonal relatedness of tested isolates was set by the BURP and BURST approach using spa and MLST data, respectively. Mini-MLST classified the S. aureus isolates into 38 clusters, followed by spa-HRM classifying the isolates into 101 clusters. The WGS proved HRM-based methods to effectively differentiate between related S. aureus isolates. Visualizing evolutionary relationships among different spa-types provided by the BURP algorithm showed comparable results to MLST/mini-MLST clonal clusters. We proved that the combination of mini-MLST and spa-HRM is rapid, reproducible, and cost-efficient. In addition to high discriminatory ability, the correlation between spa evolutionary relationships and mini-MLST clustering allows the variability in population structure to be monitored. IMPORTANCE Rapid and cost-effective molecular typing tools for Staphylococcus aureus epidemiological applications such as transmission tracking, source attribution and outbreak investigations are highly desirable. High-resolution melting based methods are effective alternative to those based on sequencing. Their good reproducibility and easy performance allow prospective typing of large set of isolates while reaching great discriminatory power. In this study, we established a new epidemiological approach to S. aureus typing. This scheme has the potential to greatly improve epidemiological investigations of S. aureus.

Epigenetics and Helicobacter pylori

Epigenetics regulates gene expression, cell type development during differentiation, and the cell response to environmental stimuli. To survive, bacteria need to evade the host immune response. Bacteria, including Helicobacter pylori (Hp), reach this target epigenetically, altering the chromatin of the host cells, in addition to several more approaches, such as DNA mutation and recombination. This review shows that Hp prevalently silences the genes of the human gastric mucosa by DNA methylation. Epigenetics includes different mechanisms. However, DNA methylation persists after DNA replication and therefore is frequently associated with the inheritance of repressed genes. Chromatin modification can be transmitted to daughter cells leading to heritable changes in gene expression. Aberrant epigenetic alteration of the gastric mucosa DNA remains the principal cause of gastric cancer. Numerous methylated genes have been found in cancer as well as in precancerous lesions of Hp-infected patients. These methylated genes inactivate tumor-suppressor genes. It is time for us to complain about our genetic and epigenetic makeups for our diseases.

Genotyping of Chlamydia abortus using multiple loci variable number of tandem repeats analysis technique

Background

The correlation between various factors (geographical region, clinical incidence, and host type) and the genomic heterogeneity has been shown in several bacterial strains including Chlamydia abortus.

Methods

The aim of this study was to survey the predominant types of C. abortus strains isolated from ruminants in Iran by the multiple loci variable number of tandem repeats (VNTR) analysis (MLVA) method. C. abortus infection was evaluated in a total of 117 aborted fetuses by real-time PCR. The isolation was done via the inoculation of the positive samples in chicken embryo and the L929 cell line. Genotyping was carried out by MLVA typing technique.

Results

Forty samples (34.2%) were detected with C. abortus infection; however, chlamydial infection in ruminants of Charmahal/Bakhtiari (3 bovines and 35 sheep) was higher than that of Khuzestan (2 sheep). All MLVA types (MT1-MT8) were detected in the collected samples from Charmahal/Bakhtiari but only 2 types (MT1 and MT3) were reported in samples from Khuzestan. The main MT type was MT1 (32% of aborted fetuses). Although in this study only 9 cow samples were investigated, they possessed similar clusters to those obtained from sheep (MT1 and MT6).

Variation of type in sheep samples (MT1 to MT8) was more than that of bovine samples (MT1, and MT6).

Conclusion

By this research revealed that C.abortus was responsible for a significant percentage of ruminant abortion in two studied regions. The main MT type was MT1 (32% of aborted fetuses) and also 7 different genotypes were involved in infections. So it is concluded that diversity in C.abortus genotyping is high in two regions.

Salmonella Typhi from Northwest Pakistan: Molecular Strain Typing and Drug Resistance Signature

Objective: To determine the molecular strain typing and drug resistance pattern of Salmonella enterica serovar Typhi prevalent in Northwest Pakistan. Methodology: A total of 2,138 blood samples of suspected typhoid patients from Northwest Pakistan were collected followed by identification of Salmonella Typhi through biochemical, serological, and species-specific fliC-d gene amplification. These isolates were typed by variable-number tandem repeat (VNTR) profiling and investigated for drug resistance. Results: The overall prevalence of Salmonella Typhi was found to be 8.8% (n = 189). Thirty different VNTR strain types of Salmonella Typhi were detected and the most prevalent strain types were T1 and T4, whereas T27 was less prevalent strain. Among the 189 isolates 175 (92.5%) isolates were multidrug resistant, whereas 12 (5.8%) isolates were extensively drug resistant. Resistance to imipenem in Salmonella Typhi was not observed. Most of the isolates have genes encoding for resistance to fluoroquinolones, including gyrA (n = 164), gyrB (n = 160), parC (n = 164), parE (n = 160), ac(6')-ib-cr (n = 163), qnrS (n = 15), and qnrB (n = 3). Similarly, chloramphinicol (cat; n = 147), azithromycin (msrA; n = 3), and co-trimoxazole (dfrA7; n = 145) resistance genes were detected among Salmonella Typhi isolates. Conclusion: In this study, T1 and T4 type Salmonella Typhi strains were predominantly prevalent in Northwest Pakistan. Antibiotic resistance among Salmonella Typhi isolates were observed. Findings of the study would be helpful to devise an appropriate antibiotic policy to control the emergence of drug-resistant Salmonella Typhi in Pakistan.

Multilocus Variable-Number Tandem-Repeat Analysis of Enterohemorrhagic Escherichia coli Serogroups O157, O26, and O111 Based on a De Novo Look-Up Table Constructed by Regression Analysis

Multilocus variable-number tandem-repeat analysis (MLVA) is a widely accepted molecular typing tool for enterohemorrhagic Escherichia coli (EHEC). However, ensuring the accuracy of MLVA data among multiple laboratories remains difficult. We developed a method of constructing adjusted look-up tables, which are necessary for MLVA profiling, at each laboratory using a regression analysis based on electrophoresis data from 24 in-house reference strains. On performing MLVA against 51 EHEC O157 isolates, the repeat numbers of 46 isolates were determined accurately using the look-up table with a 99% prediction interval, an outcome superior to that when using a 95% prediction interval. For the remaining five isolates, although the electrophoresis size fell outside the look-up table, we were able to predict the repeat number accurately by extrapolation or the nearest values of the look-up table. Our approach provides more accurate results than a nonadjusted conventional look-up table for calibrating MLVA profiles.

STRs: Ancient Architectures of the Genome beyond the Sequence

Short tandem repeats (STRs) are commonly defined as short runs of repetitive nucleotides, consisting of tandemly repeating 2-6- bp motif units, which are ubiquitously distributed throughout genomes. Functional STRs are polymorphic in the population, and their variations influence gene expression, which subsequently may result in pathogenic phenotypes. To understand STR phenotypic effects and their functional roles, we describe four different mutational mechanisms including the unequal crossing-over model, gene conversion, retrotransposition mechanism and replication slippage. Due to the multi-allelic nature, small length, abundance, high variability, codominant inheritance, nearly neutral evolution, extensive genome coverage and simple assaying of STRs, these markers are widely used in various types of biological research, including population genetics studies, genome mapping, molecular epidemiology, paternity analysis and gene flow studies. In this review, we focus on the current knowledge regarding STR genomic distribution, function, mutation and applications.

Correlating Genotyping Data of Coxiella burnetii with Genomic Groups

Coxiella burnetii is a zoonotic pathogen that resides in wild and domesticated animals across the globe and causes a febrile illness, Q fever, in humans. Several distinct genetic lineages or genomic groups have been shown to exist, with evidence for different virulence potential of these lineages. Multispacer Sequence Typing (MST) and Multiple-Locus Variable number tandem repeat Analysis (MLVA) are being used to genotype strains. However, it is unclear how these typing schemes correlate with each other or with the classification into different genomic groups. Here, we created extensive databases for published MLVA and MST genotypes of C. burnetii and analysed the associated metadata, revealing associations between animal host and human disease type. We established a new classification scheme that assigns both MST and MLVA genotypes to a genomic group and which revealed additional sub-lineages in two genomic groups. Finally, we report a novel, rapid genomotyping method for assigning an isolate into a genomic group based on the Cox51 spacer sequence. We conclude that by pooling and streamlining existing datasets, associations between genotype and clinical outcome or host source were identified, which in combination with our novel genomotyping method, should enable an estimation of the disease potential of new C. burnetii isolates.

Molecular Markers Reveal Epidemiological Patterns and Evolutionary Histories of the Human Pathogenic Cryptococcus

The human pathogenic Cryptococcus species are the main agents of fungal meningitis in humans and the causes of other diseases collectively called cryptococcosis. There are at least eight evolutionary divergent lineages among these agents, with different lineages showing different geographic and/or ecological distributions. In this review, we describe the main strain typing methods that have been used to analyze the human pathogenic Cryptococcus and discuss how molecular markers derived from the various strain typing methods have impacted our understanding of not only cryptococcal epidemiology but also its evolutionary histories. These methods include serotyping, multilocus enzyme electrophoresis, electrophoretic karyotyping, random amplified polymorphic DNA, restriction fragment length polymorphism, PCR-fingerprinting, amplified fragment length polymorphism, multilocus microsatellite typing, single locus and multilocus sequence typing, matrix-assisted laser desorption/ionization time of flight mass spectrometry, and whole genome sequencing. The major findings and the advantages and disadvantages of each method are discussed. Together, while controversies remain, these strain typing methods have helped reveal (i) the broad phylogenetic pattern among these agents, (ii) the centers of origins for several lineages and their dispersal patterns, (iii) the distributions of genetic variation among geographic regions and ecological niches, (iv) recent hybridization among several lineages, and (v) specific mutations during infections within individual patients. However, significant challenges remain. Multilocus sequence typing and whole genome sequencing are emerging as the gold standards for continued strain typing and epidemiological investigations of cryptococcosis.

Detection and application of novel SSR markers from transcriptome data for Astronium fraxinifolium Schott, a threatened Brazilian tree species

Astronium fraxinifolium is an endangered tree species from Brazil. Due to its significance in environmental reforestation, as well as the continued exploitation of its wood, it is necessary to develop management programs that support the conservation of the species. Simple sequence repeats (SSR) or microsatellite markers are widely used in population genetic studies across a range of diverse organisms. In this study, we present the first SSR markers developed for A. fraxinifolium as well as their frequency and distribution based on transcriptome data. From transcriptome data, we identified more than 100 thousand sequences presenting microsatellites, with a predominant distribution of trinucleotide repeats. From the initial screening, we selected 20 microsatellite loci which were validated and evaluated for genetic indices in two natural populations. All loci were polymorphic, ranging from four to 11 alleles per locus. The observed and expected heterozygosities ranged from 0 to 1.0 and from 0.533 to 1.0, respectively, while the genetic differentiation (G_ST = 0.363) was greater within than between populations. The developed SSR loci from RNA-Seq data provides a foundation for future studies on genetic diversity and population structure, mating system, and gene flow for A. fraxinifolium populations and related species, aiming at conservation and management.

Trends in Molecular Diagnosis and Diversity Studies for Phytosanitary Regulated Xanthomonas

Bacteria in the genus Xanthomonas infect a wide range of crops and wild plants, with most species responsible for plant diseases that have a global economic and environmental impact on the seed, plant, and food trade. Infections by Xanthomonas spp. cause a wide variety of non-specific symptoms, making their identification difficult. The coexistence of phylogenetically close strains, but drastically different in their phenotype, poses an added challenge to diagnosis. Data on future climate change scenarios predict an increase in the severity of epidemics and a geographical expansion of pathogens, increasing pressure on plant health services. In this context, the effectiveness of integrated disease management strategies strongly depends on the availability of rapid, sensitive, and specific diagnostic methods. The accumulation of genomic information in recent years has facilitated the identification of new DNA markers, a cornerstone for the development of more sensitive and specific methods. Nevertheless, the challenges that the taxonomic complexity of this genus represents in terms of diagnosis together with the fact that within the same bacterial species, groups of strains may interact with distinct host species demonstrate that there is still a long way to go. In this review, we describe and discuss the current molecular-based methods for the diagnosis and detection of regulated Xanthomonas, taxonomic and diversity studies in Xanthomonas and genomic approaches for molecular diagnosis.

Molecular Diversity and Pathogenicity of Ralstonia solanacearum Species Complex Associated With Bacterial Wilt of Potato in Rwanda

Bacterial wilt (BW), caused by Ralstonia solanacearum species complex (RSSC), leads to substantial potato yield losses in Rwanda. Studies were conducted to (i) determine the molecular diversity of RSSC strains associated with BW of potato, (ii) generate an RSSC distribution map for epidemiological inferences, and (iii) test the pathogenicity of predominant RSSC phylotypes on six commercial potato cultivars. In surveys conducted in 2018 and 2019, tubers from wilting potato plants were collected for pathogen isolation. DNA was extracted from 95 presumptive RSSC strain colonies. The pathogen was phylotyped by multiplex PCR and typed at sequevar level. Phylotype II sequevar 1 strains were then haplotyped using multilocus tandem repeat sequence typing (TRST) schemes. Pathogenicity of one phylotype II strain and two phylotype III strains were tested on cultivars Kinigi, Kirundo, Victoria, Kazeneza, Twihaze, and Cruza. Two RSSC phylotypes were identified, phylotype II (95.79%, n = 91) and phylotype III (4.21%, n = 4). This is the first report of phylotype III strains from Rwanda. Phylotype II strains were identified as sequevar 1 and distributed across potato growing regions in the country. The TRST scheme identified 14 TRST haplotypes within the phylotype II sequevar 1 strains with moderate diversity index (HGDI = 0.55). Mapping of TRST haplotypes revealed that a single TRST '8-5-12-7-5' haplotype plays an important epidemiological role in BW of potato in Rwanda. None of the cultivars had complete resistance to the tested phylotypes; the level of susceptibility varied among cultivars. Cultivar Cruza, which is less susceptible to phylotype II and III strains, is recommended when planting potatoes in the fields with history of BW.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

BACKGROUND

Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes.

RESULTS

We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats.

CONCLUSION

This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

Genome (in)stability at tandem repeats

Repeat sequences account for over half of the human genome and represent a significant source of variation that underlies physiological and pathological states. Yet, their study has been hindered due to limitations in short-reads sequencing technology and difficulties in assembly. A important category of repetitive DNA in the human genome is comprised of tandem repeats (TRs), where repetitive units are arranged in a head-to-tail pattern. Compared to other regions of the genome, TRs carry between 10 and 10,000 fold higher mutation rate. There are several mutagenic mechanisms that can give rise to this propensity toward instability, but their precise contribution remains speculative. Given the high degree of homology between these sequences and their arrangement in tandem, once damaged, TRs have an intrinsic propensity to undergo aberrant recombination with non-allelic exchange and generate harmful rearrangements that may undermine the stability of the entire genome. The dynamic mutagenesis at TRs has been found to underlie individual polymorphism associated with neurodegenerative and neuromuscular disorders, as well as complex genetic diseases like cancer and diabetes. Here, we review our current understanding of the surveillance and repair mechanisms operating within these regions, and we describe how alterations in these protective processes can readily trigger mutational signatures found at TRs, ultimately resulting in the pathological correlation between TRs instability and human diseases. Finally, we provide a viewpoint to counter the detrimental effects that TRs pose in light of their selection and conservation, as important drivers of human evolution.

A positive perspective on DNA methylation: regulatory functions of DNA methylation outside of host defense in Gram-positive bacteria

The presence of post-replicative DNA methylation is pervasive among both prokaryotic and eukaryotic organisms. In bacteria, the study of DNA methylation has largely been in the context of restriction-modification systems, where DNA methylation serves to safeguard the chromosome against restriction endonuclease cleavage intended for invading DNA. There has been a growing recognition that the methyltransferase component of restriction-modification systems can also regulate gene expression, with important contributions to virulence factor gene expression in bacterial pathogens. Outside of restriction-modification systems, DNA methylation from orphan methyltransferases, which lack cognate restriction endonucleases, has been shown to regulate important processes, including DNA replication, DNA mismatch repair, and the regulation of gene expression. The majority of research and review articles have been focused on DNA methylation in the context of Gram-negative bacteria, with emphasis toward Escherichia coli, Caulobacter crescentus, and related Proteobacteria. Here we summarize the epigenetic functions of DNA methylation outside of host defense in Gram-positive bacteria, with a focus on the regulatory effects of both phase variable methyltransferases and DNA methyltransferases from traditional restriction-modification systems.

Unique Features of Tandem Repeats in Bacteria

DNA tandem repeats, or satellites, are well described in eukaryotic species, but little is known about their prevalence across prokaryotes. Here, we performed the most complete characterization to date of satellites in bacteria. We identified 121,638 satellites from 12,233 fully sequenced and assembled bacterial genomes with a very uneven distribution. We also determined the families of satellites which have a related sequence. There are 85 genomes that are particularly satellite rich and contain several families of satellites of yet unknown function. Interestingly, we only found two main types of noncoding satellites, depending on their repeat sizes, 22/44 or 52 nucleotides (nt). An intriguing feature is the constant size of the repeats in the genomes of different species, whereas their sequences show no conservation. Individual species also have several families of satellites with the same repeat length and different sequences. This result is in marked contrast with previous findings in eukaryotes, where noncoding satellites of many sizes are found in any species investigated. We describe in greater detail these noncoding satellites in the spirochete Leptospira interrogans and in several bacilli. These satellites undoubtedly play a specific role in the species which have acquired them. We discuss the possibility that they represent binding sites for transcription factors not previously described or that they are involved in the stabilization of the nucleoid through interaction with proteins.IMPORTANCE We found an enigmatic group of noncoding satellites in 85 bacterial genomes with a constant repeat size but variable sequence. This pattern of DNA organization is unique and had not been previously described in bacteria. These findings strongly suggest that satellite size in some bacteria is under strong selective constraints and thus that satellites are very likely to play a fundamental role. We also provide a list and properties of all satellites in 12,233 genomes, which may be used for further genomic analysis.

A Simple Blocking PCR-Based Method for the Synthesis of High-Copy dsDNA Tandem Repeats

DNA tandem repeats are frequently found in eukaryotic genomes. High-copy DNA repeats can serve as building blocks of complex DNA structures, but the in vitro synthesis of DNA repeats has been challenging due to complicated procedures and the high cost. Here, a new, simple method is developed using the strategy of blocking polymerase chain reaction for highly efficient DNA repeat expansion (BPRE). With BPRE, dsDNA fragments composed of more than 40 copies of the repeat sequence can be quickly produced, while the cost is reduced by at least 90%. As a typical application, reannealing of the dsDNA repeats generates elastic hydrogels, which shows a high capacity for doxycycline absorption and prolonged release.

Site-Specific Recombination - How Simple DNA Inversions Produce Complex Phenotypic Heterogeneity in Bacterial Populations

Many bacterial species generate phenotypically heterogeneous subpopulations as a strategy for ensuring the survival of the population as a whole - an environmental stress that eradicates one subpopulation may leave other phenotypic groups unharmed, allowing the lineage to continue. Phase variation, a process that functions as an ON/OFF switch for gene expression, is one way that bacteria achieve phenotypic heterogeneity. Phase variation occurs stochastically and reversibly, and in the presence of a selective pressure the advantageous phenotype(s) predominates in the population. Phase variation can occur through multiple genetic and epigenetic mechanisms. This review focuses on conservative site-specific recombination that generates reversible DNA inversions as a genetic mechanism mediating phase variation. Recent studies have sparked a renewed interest in phase variation mediated through DNA inversion, revealing a high level of complexity beyond simple ON/OFF switching, including unusual modes of gene regulation, and highlighting an underappreciation of the use of these mechanisms by bacteria.

Structures and stability of simple DNA repeats from bacteria

DNA is a fundamentally important molecule for all cellular organisms due to its biological role as the store of hereditary, genetic information. On the one hand, genomic DNA is very stable, both in chemical and biological contexts, and this assists its genetic functions. On the other hand, it is also a dynamic molecule, and constant changes in its structure and sequence drive many biological processes, including adaptation and evolution of organisms. DNA genomes contain significant amounts of repetitive sequences, which have divergent functions in the complex processes that involve DNA, including replication, recombination, repair, and transcription. Through their involvement in these processes, repetitive DNA sequences influence the genetic instability and evolution of DNA molecules and they are located non-randomly in all genomes. Mechanisms that influence such genetic instability have been studied in many organisms, including within human genomes where they are linked to various human diseases. Here, we review our understanding of short, simple DNA repeats across a diverse range of bacteria, comparing the prevalence of repetitive DNA sequences in different genomes. We describe the range of DNA structures that have been observed in such repeats, focusing on their propensity to form local, non-B-DNA structures. Finally, we discuss the biological significance of such unusual DNA structures and relate this to studies where the impacts of DNA metabolism on genetic stability are linked to human diseases. Overall, we show that simple DNA repeats in bacteria serve as excellent and tractable experimental models for biochemical studies of their cellular functions and influences.

Genetic heterogeneity of the Spy1336/R28-Spy1337 virulence axis in Streptococcus pyogenes and effect on gene transcript levels and pathogenesis

Streptococcus pyogenes is a strict human pathogen responsible for more than 700 million infections annually worldwide. Strains of serotype M28 S. pyogenes are typically among the five more abundant types causing invasive infections and pharyngitis in adults and children. Type M28 strains also have an unusual propensity to cause puerperal sepsis and neonatal disease. We recently discovered that a one-nucleotide indel in an intergenic homopolymeric tract located between genes Spy1336/R28 and Spy1337 altered virulence in a mouse model of infection. In the present study, we analyzed size variation in this homopolymeric tract and determined the extent of heterogeneity in the number of tandemly-repeated 79-amino acid domains in the coding region of Spy1336/R28 in large samples of strains recovered from humans with invasive infections. Both repeat sequence elements are highly polymorphic in natural populations of M28 strains. Variation in the homopolymeric tract results in (i) changes in transcript levels of Spy1336/R28 and Spy1337 in vitro, (ii) differences in virulence in a mouse model of necrotizing myositis, and (iii) global transcriptome changes as shown by RNAseq analysis of isogenic mutant strains. Variation in the number of tandem repeats in the coding sequence of Spy1336/R28 is responsible for size variation of R28 protein in natural populations. Isogenic mutant strains in which genes encoding R28 or transcriptional regulator Spy1337 are inactivated are significantly less virulent in a nonhuman primate model of necrotizing myositis. Our findings provide impetus for additional studies addressing the role of R28 and Spy1337 variation in pathogen-host interactions.

Ribosomal protein L4 of Lactobacillus rhamnosus LRB alters resistance to macrolides and other antibiotics

Lactobacillus rhamnosus is an important lactic acid bacterium that is predominantly used as a probiotic supplement. This bacterium secretes immunomodulatory and antibacterial peptides that are necessary for the probiotic trait. This organism also occupies diverse ecological niches, such as gastrointestinal tracts and the oral cavity. Several studies have shown that L. rhamnosus is prone to spontaneous genome rearrangement irrespective of the ecological origins. We previously characterized an oral isolate of L. rhamnosus, LRB, which is genetically closely related to the widely used probiotic strain L. rhamnosus LGG. In this study, we isolated a nontargeted mutant that was particularly sensitive to acid stress. Using next generation sequencing, we further mapped the putative mutations in the genome and found that the mutant had acquired a deletion of 75 base pairs in the rplD gene that encodes the large ribosomal subunit L4. The mutant had a growth defect at 37°C and at ambient temperature. Further antibiotic sensitivity analyses indicated that the mutant is relatively more resistant to erythromycin and chloramphenicol; two antibiotics that target the 50S subunit. In contrast, the mutant was more sensitive to tetracycline, which targets the 30S subunit. Thus, it appears that nontargeted mutations could significantly alter the antibiotic resistance profile of L. rhamnosus. Our study raises concern that probiotic use of L. rhamnosus should be carefully monitored to avoid unintended consequences.

Evaluation of genetic diversity among aquatic and fecal isolates of Escherichia coli using multilocus variable number of tandem repeat analysis

In developing countries like India, fecal pollution of surface waters is a major threat to public and environmental health. The aim of the study was to assess serological, phylogenetic and molecular diversity among aquatic Escherichia coli isolates from Yamuna river and their comparison with the animal fecal isolates. A total of 97 E. coli isolates from Yamuna river and domesticated animals were characterized by multilocus variable number tandem repeat analysis (MLVA) using four VNTR loci. The pathogenicity of these strains by serological and phylogenetic analysis was also determined. E. coli strains were differentiated into 53 distinct MLVA types with high discriminatory power, Simpson's index of 0.95 (95% CI 0.923-0.978). Cluster analysis and population modeling using minimum spanning tree suggested a possible epidemiological linkage among aquatic and fecal isolates. The study also reported the presence of highly diverse and pathogenic serotypes belonging to STEC and EPEC strains, particularly O157 and high prevalence of pathogenic phylogroups (phylogroup, B2 and D). The presence of such a high molecular heterogeneity among aquatic and fecal E. coli isolates emphasizes upon the need to develop proper fecal pollution abatement strategies for Indian natural bodies.

Genetic Variability Based on Tandem Repeat Numbers in a Genomic Locus of 'Candidatus Liberibacter asiaticus' Prevalent in North East India

The genetic variability of 'Candidatus Liberibacter asiaticus' (CLas) population associated with huanglongbing (HLB) disease of citrus in North Eastern (NE) region of India, a geographically locked region, and home for the diversity of many citrus species was analyzed on the basis of tandem repeat numbers (TRN) in variable CLIBASIA_01645 genomic loci. Fifty-five CLas strains sampled from different groves of NE Hill (NEH) region of India were in single amplicon group, but there was remarkable genetic variability in TRNs. The TRN in HLB-associated CLas strains varied from 0-21 and two novel repeat motifs were also identified. Among the NE population of CLas, TRN5 and TRN9 were most frequent (total frequency of 36.36%) followed by TRN4 (14.55%) and TRN6, TNR7 with a frequency of 12.73% each. Class II type CLas genotypes (5 < TRN ≤ 10) had highest prevalence (frequency of 60.00%) in the samples characterized in present study. Class I (TRN ≤ 5) genotypes were second highest prevalent (29.09%) in the NEH region. Further analysis of genetic diversity parameters using Nei's measure (H value) indicated wide genetic diversity in the CLas strains of NE India (H value of 0.58-0.86). Manipur CLas strains had highest genetic variability (0.86) as compared to Eastern, Southern and Central India. The R₁₀ values (TRN ≤ 10/TRN > 10) of NE CLas population was 10.43 (73/7), higher from other regions of India. Present study conclusively reported the occurrence of high genetic variability in TRN of CLas population in North East Indian citrus groves which have evolved to adapt to the specific ecological niche.