Genetic Diversification Research Articles

Diversification of Bourbon Virus in New York State.

Bourbon virus (BRBV, family Orthomyxoviridae) is a tickborne virus recently detected in the United States (US). BRBV was first identified from a fatal human case in 2014 in Bourbon County, Kansas. Enhanced surveillance in Kansas and Missouri implicated Amblyomma americanum as the primary vector for BRBV. Historically, BRBV was only detected in the lower midwestern US, but since 2020 it has been reported in North Carolina, Virginia, New Jersey, and New York State (NYS). This study aimed to elucidate genetic and phenotypic characteristics of BRBV strains from NYS through whole genome sequencing and the assessment of replication kinetics in mammalian cultures and A. americanum nymphs. Sequence analysis revealed the existence of two divergent BRBV clades circulating in NYS. BRBV NY21-2143 is closely related to the midwestern BRBV strains but has unique substitutions in the glycoprotein. Two other NYS BRBV strains, BRBV NY21-1814 and BRBV NY21-2666, form a distinct clade unique from previously sequenced BRBV strains. Phenotypic diversification was also detected in NYS BRBV strains compared to each other and midwestern BRBV strains, with BRBV NY21-2143 displaying attenuation in rodent-derived cell culture and a fitness advantage in experimentally infected A. americanum. These data suggest genetic and phenotypic diversification of emergent BRBV strains circulating in NYS that could contribute to increased spread of BRBV in the northeastern US.

SARS-CoV-2 Co-Infections and Recombinations Identified by Long-Read Single-Molecule Real-Time Sequencing.

Co-infection with at least 2 strains of virus is the prerequisite for recombination, one of the means of genetic diversification. Little is known about the prevalence of these events in SARS-CoV-2, partly because it is difficult to detect them. We used long-read PacBio single-molecule real-time (SMRT) sequencing technology to sequence whole genomes and targeted regions for haplotyping. We identified 17 co-infections with SARS-CoV-2 strains belonging to different clades in 6829 samples sequenced between January and October, 2022 (prevalence 0.25%). There were 3 Delta/Omicron co-infections and 14 Omicron/Omicron co-infections (4 cases of 21K/21L, 1 case of 21L/22A, 2 cases of 21L/22B, 4 cases of 22A/22B, 2 cases of 22B/22C and 1 case of 22B/22E). Four of these patients (24%) also harbored recombinant minor haplotypes, including one with a recombinant virus that was selected in the viral quasispecies over the course of his chronic infection. While co-infections remain rare among SARS-CoV-2-infected individuals, long-read SMRT sequencing is a useful tool for detecting them as well as recombinant events, providing the basis for assessing their clinical impact, and a precise indicator of epidemic evolution. IMPORTANCE SARS-CoV-2 variants have been responsible for the successive waves of infection over the 3 years of pandemic. While co-infection followed by recombination is one driver of virus evolution, there have been few reports of co-infections, mainly between Delta and Omicron variants or between the first 2 Omicron variants 21K_BA.1 and 21L_BA.2. The 17 co-infections we detected during 2022 included cases with the recent clades of Omicron 22A, 22B, 22C, and 22E; 24% harbored recombinant variants. This study shows that long-read SMRT sequencing is well suited to SARS-CoV-2 genomic surveillance.

MAF functions as a pioneer transcription factor that initiates and sustains myelomagenesis

AbstractDeregulated expression of lineage-affiliated transcription factors (TFs) is a major mechanism of oncogenesis. However, how the deregulation of nonlineage affiliated TF affects chromatin to initiate oncogenic transcriptional programs is not well-known. To address this, we studied the chromatin effects imposed by oncogenic MAF as the cancer-initiating driver in the plasma cell cancer multiple myeloma. We found that the ectopically expressed MAF endows myeloma plasma cells with migratory and proliferative transcriptional potential. This potential is regulated by the activation of enhancers and superenhancers, previously inactive in healthy B cells and plasma cells, and the cooperation of MAF with the plasma cell-defining TF IRF4. Forced ectopic MAF expression confirms the de novo ability of oncogenic MAF to convert transcriptionally inert chromatin to active chromatin with the features of superenhancers, leading to the activation of the MAF-specific oncogenic transcriptome and the acquisition of cancer-related cellular phenotypes such as CCR1-dependent cell migration. These findings establish oncogenic MAF as a pioneer transcription factor that can initiate as well as sustain oncogenic transcriptomes and cancer phenotypes. However, despite its pioneer function, myeloma cells remain MAF-dependent, thus validating oncogenic MAF as a therapeutic target that would be able to circumvent the challenges of subsequent genetic diversification driving disease relapse and drug resistance.

Genetic Diversification of Starch Branching Enzymes during Maize Domestication and Improvement.

Elucidating the genetic basis of starch pasting and gelatinization properties is crucial for enhancing the quality of maize and its utility as feed and industrial raw material. In maize, ZmSBE genes encode important starch branching enzymes in the starch biosynthesis pathway. In this study, we re-sequenced the genomic sequences of ZmSBEI, ZmSBEIIa, ZmSBEIIb, and ZmSBEIII in three lines called 335 inbred lines, 68 landrace lines, and 32 teosinte lines. Analyses of nucleotide polymorphisms and haplotype diversity revealed differences in the selection patterns of ZmSBEI, ZmSBEIIa, ZmSBEIIb, and ZmSBEIII during maize domestication and improvement. A marker-trait association analysis of inbred lines detected 22 significant loci, including 18 SNPs and 4 indels significantly associated with three maize starch physicochemical properties. The allele frequencies of two variants (SNP17249C and SNP5055G) were examined in three lines. The frequency of SNP17249C in ZmSBEIIb was highest in teosinte lines, followed by landrace lines, and inbred lines, whereas there were no significant differences in the frequency of SNP5055G in ZmSBEIII among the three lines. These results suggest that ZmSBE genes play an important role in the phenotypic variations in the starch physicochemical properties in maize. The genetic variants detected in this study may be used to develop functional markers for improving maize starch quality.

Is integrated taxonomy useful to study diversity and ecology? An example from crustacean zooplankton at the Long‐Term ecological research site MareChiara (LTER‐MC)

AbstractMolecular tools increasingly refine and improve the identification of zooplankton organisms based on phenotypic features, providing a more robust and comprehensive species description. Integration of data helps revealing the hidden diversity of zooplankton and facilitating the detection of rare and non‐indigenous species. This approach, merging morphological characters and a diagnostic marker for specific identification, such as the mitochondrial cytochrome c oxidase I (COI), is here used to characterize key taxa from the zooplankton assemblage of the Gulf of Naples at the Long‐Term Ecological Research site MareChiara (LTER‐MC) (Central Tyrrhenian Sea, Western Mediterranean Sea). Zooplankton biodiversity assessment using integrated taxonomy was focused on selected crustacean groups: cyclopod copepods (Agetus typicus, Oithona plumifera, Oncaea mediterranea, Oncaea scottodicarloi); newly records of cladocerans (Evadne nordmanni), euphausiids (Euphausia krohnii, Nematoscelis megalops, Nyctiphanes spp.) and sergestids (Lucifer typus), with the aim to boost the knowledge of real zooplankton biodiversity. The results of our investigation provide new high‐quality molecular references of the analysed taxa and contribute to unveiling the genetic diversification of zooplankton species and their relevant ecological significance for Mediterranean coastal waters.

Progress in Bioengineering of Myotropic Adeno-Associated Viral Gene Therapy Vectors

The ability to specifically, safely, and efficiently transfer therapeutic payloads to the striated musculature via a minimally invasive delivery route remains one of the most important but also most ambitious aims in human gene therapy. Over the past two decades, a flurry of groups have harnessed recombinant adeno-associated viruses (AAVs) for this purpose, carrying cargoes that were packaged either in one of the various wild-type capsids or in a synthetic protein shell derived by molecular bioengineering. In this study, we provide an overview over the most commonly used techniques for the enrichment of muscle-specific (myotropic) AAV capsids, typically starting off with the genetic diversification of one or more extant wild-type sequences, followed by the stratification of the ensuing capsid libraries in different muscle types in small or large animals. These techniques include the shuffling of multiple parental capsid genes, peptide display in exposed capsid loops, mutagenesis of individual capsid residues, creation of chimeras between two viral parents, or combinations thereof. Moreover, we highlight alternative experimental or bioinformatic strategies such as ancestral reconstruction or rational design, all of which have already been employed successfully to derive synthetic AAV capsids or vectors with unprecedented in vivo efficiency and/or specificity in the musculature. Most recently, these efforts have culminated in the isolation of unique clades of myotropic vectors called AAVMYO or MyoAAV that have in common the display of the amino acid motif RGD (arginine-glycine-aspartate) on the capsid surface and that exhibit the highest transduction rate in striated muscles of mice or nonhuman primates reported to date. Finally, we note essential looming improvements that will facilitate and accelerate clinical translation of these latest generations of myotropic AAVs, including the identification and utilization of capsid selection or validation schemes that promise optimal translation in humans, and continued efforts to enhance patient safety by minimizing hepatic off-targeting.

Chiropterans Are a Hotspot for Horizontal Transfer of DNA Transposons in Mammalia.

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.

Phylogenomics of the Olea europaea complex using 15 whole genomes supports recurrent genetic admixture together with differentiation into seven subspecies

BackgroundThe last taxonomic account of Olea recognises six subspecies within Olea europaea L., including the Mediterranean olive tree (subsp. europaea) and five other subspecies (laperrinei, guanchica, maroccana, cerasiformis, and cuspidata) distributed across the Old World, including Macaronesian islands. The evolutionary history of this monophyletic group (O. europaea complex) has revealed a reticulated scenario involving hybridization and polyploidization events, leading to the presence of a polyploid series associated with the subspecies. However, how the polyploids originated, and how the different subspecies contributed to the domestication of the cultivated olive are questions still debated. Tracing the recent evolution and genetic diversification of the species is key for the management and preservation of its genetic resources. To study the recent history of the O. europaea complex, we compared newly sequenced and available genomes for 27 individuals representing the six subspecies.ResultsOur results show discordance between current subspecies distributions and phylogenomic patterns, which support intricate biogeographic patterns. The subspecies guanchica, restricted to the Canary Islands, is closely related to subsp. europaea, and shows a high genetic diversity. The subsp. laperrinei, restricted now to high mountains of the Sahara desert, and the Canarian subsp. guanchica contributed to the formation of the allotetraploid subsp. cerasiformis (Madeira islands) and the allohexaploid subsp. maroccana (western Sahara region). Our phylogenomic data support the recognition of one more taxon (subsp. ferruginea) for the Asian populations, which is clearly segregated from the African subsp. cuspidata.ConclusionsIn sum, the O. europaea complex underwent several processes of hybridization, polyploidy, and geographical isolation resulting in seven independent lineages with certain morphological traits recognised into subspecies.

High Nucleotide Diversity Accompanies Differential DNA Methylation in Naturally Diverging Populations.

Epigenetic mechanisms such as DNA methylation (DNAme) are thought to comprise an invaluable adaptive toolkit in the early stages of local adaptation, especially when genetic diversity is constrained. However, the link between genetic diversity and DNAme has been scarcely examined in natural populations, despite its potential to shed light on the evolutionary forces acting on methylation state. Here, we analyzed reduced-representation bisulfite sequencing and whole-genome pool-seq data from marine and freshwater stickleback populations to examine the relationship between DNAme variation (between- and within-population) and nucleotide diversity in the context of freshwater adaptation. We find that sites that are differentially methylated between populations have higher underlying standing genetic variation, with diversity higher among sites that gained methylation in freshwater than those that lost it. Strikingly, although nucleotide diversity is generally lower in the freshwater population as expected from a population bottleneck, this is not the case for sites that lost methylation, which instead have elevated nucleotide diversity in freshwater compared with marine. Subsequently, we show that nucleotide diversity is higher among sites with ancestrally variable methylation and also positively correlates with the sensitivity to environmentally induced methylation change. The results suggest that as selection on the control of methylation state becomes relaxed, so too does selection against mutations at the sites themselves. Increased epigenetic variance in a population is therefore likely to precede genetic diversification.

Enhancing the quality of staple food crops through CRISPR/Cas-mediated site-directed mutagenesis

The enhancement of CRISPR-Cas gene editing with robust nuclease activity promotes genetic modification of desirable agronomic traits, such as resistance to pathogens, drought tolerance, nutritional value, and yield-related traits in crops. The genetic diversity of food crops has reduced tremendously over the past twelve millennia due to plant domestication. This reduction presents significant challenges for the future especially considering the risks posed by global climate change to food production. While crops with improved phenotypes have been generated through crossbreeding, mutation breeding, and transgenic breeding over the years, improving phenotypic traits through precise genetic diversification has been challenging. The challenges are broadly associated with the randomness of genetic recombination and conventional mutagenesis. This review highlights how emerging gene-editing technologies reduce the burden and time necessary for developing desired traits in plants. Our focus is to provide readers with an overview of the advances in CRISPR-Cas-based genome editing for crop improvement. The use of CRISPR-Cas systems in generating genetic diversity to enhance the quality and nutritional value of staple food crops is discussed. We also outlined recent applications of CRISPR-Cas in developing pest-resistant crops and removing unwanted traits, such as allergenicity from crops. Genome editing tools continue to evolve and present unprecedented opportunities to enhance crop germplasm via precise mutations at the desired loci of the plant genome.

Moderation effects of serotype on dengue severity across pregnancy status in Mexico

BackgroundPregnancy increases a woman’s risk of severe dengue. To the best of our knowledge, the moderation effect of the dengue serotype among pregnant women has not been studied in Mexico. This study explores how pregnancy interacted with the dengue serotype from 2012 to 2020 in Mexico.MethodInformation from 2469 notifying health units in Mexican municipalities was used for this cross-sectional analysis. Multiple logistic regression with interaction effects was chosen as the final model and sensitivity analysis was done to assess potential exposure misclassification of pregnancy status.ResultsPregnant women were found to have higher odds of severe dengue [1.50 (95% CI 1.41, 1.59)]. The odds of dengue severity varied for pregnant women with DENV-1 [1.45, (95% CI 1.21, 1.74)], DENV-2 [1.33, (95% CI 1.18, 1.53)] and DENV-4 [3.78, (95% CI 1.14, 12.59)]. While the odds of severe dengue were generally higher for pregnant women compared with non-pregnant women with DENV-1 and DENV-2, the odds of disease severity were much higher for those infected with the DENV-4 serotype.ConclusionThe effect of pregnancy on severe dengue is moderated by the dengue serotype. Future studies on genetic diversification may potentially elucidate this serotype-specific effect among pregnant women in Mexico.

Molecular Characterization of AZD7442 (Tixagevimab-Cilgavimab) Neutralization of SARS-CoV-2 Omicron Subvariants.

Therapeutic anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) monoclonal antibodies (MAbs) provide immunosuppressed and vulnerable populations with prophylactic and treatment interventions against coronavirus disease 2019 (COVID-19). AZD7442 (tixagevimab-cilgavimab) is a combination of extended-half-life neutralizing MAbs that bind to distinct epitopes on the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. The Omicron variant of concern carries mutations at >35 positions in the spike protein and has undergone further genetic diversification since its emergence in November 2021. Here, we characterize the in vitro neutralization activity of AZD7442 toward major viral subvariants circulating worldwide during the first 9 months of the Omicron wave. BA.2 and its derived subvariants showed the highest susceptibility to AZD7442, while BA.1 and BA.1.1 showed a lower susceptibility. BA.4/BA.5 had a susceptibility level intermediate between BA.1 and BA.2. Mutagenesis of parental Omicron subvariant spike proteins was performed to establish a molecular model to describe the underlying determinants of neutralization by AZD7442 and its component MAbs. The concurrent mutation of residues at positions 446 and 493, located in the tixagevimab and cilgavimab binding sites, was sufficient to enhance in vitro susceptibility of BA.1 to AZD7442 and its component MAbs to levels similar to the Wuhan-Hu-1+D614G virus. AZD7442 maintained neutralization activity against all Omicron subvariants tested up to and including BA.5. The evolving nature of the SARS-CoV-2 pandemic warrants continuing real-time molecular surveillance and assessment of in vitro activity of MAbs used in prophylaxis against and the treatment of COVID-19. IMPORTANCE MAbs are key therapeutic options for COVID-19 prophylaxis and treatment in immunosuppressed and vulnerable populations. Due to the emergence of SARS-CoV-2 variants, including Omicron, it is vital to ensure that neutralization is maintained for MAb-based interventions. We studied the in vitro neutralization of AZD7442 (tixagevimab-cilgavimab), a cocktail of two long-acting MAbs targeting the SARS-CoV-2 spike protein, toward Omicron subvariants circulating from November 2021 to July 2022. AZD7442 neutralized major Omicron subvariants up to and including BA.5. The mechanism of action responsible for the lower in vitro susceptibility of BA.1 to AZD7442 was investigated using in vitro mutagenesis and molecular modeling. A combination of mutations at two spike protein positions, namely, 446 and 493, was sufficient to enhance BA.1 susceptibility to AZD7442 to levels similar to the Wuhan-Hu-1+D614G ancestral virus. The evolving nature of the SARS-CoV-2 pandemic warrants continuing real-time global molecular surveillance and mechanistic studies of therapeutic MAbs for COVID-19.

Luigi Luca Cavalli Sforza and the history of human languages: A linguist’s point of view

L.L. Cavalli Sforza had the great merit of transposing the methods and concepts of a modern natural science, human genetics, into historical research; after this model, it is now possible to transfer the research style and results of another growing discipline, cognitive biolinguistics, to the field of history. It is along this new line that it becomes finally possible to formally pursue Cavalli Sforza’s enterprise of assessing the degree of congruence between genetic and linguistic diversification of human populations.

Interpersonal variability of the human gut virome confounds disease signal detection in IBD

Viruses are increasingly recognised as important components of the human microbiome, fulfilling numerous ecological roles including bacterial predation, immune stimulation, genetic diversification, horizontal gene transfer, microbial interactions, and augmentation of metabolic functions. However, our current view of the human gut virome is tainted by previous sequencing requirements that necessitated the amplification of starting nucleic acids. In this study, we performed an original longitudinal analysis of 40 healthy control, 19 Crohn’s disease, and 20 ulcerative colitis viromes over three time points without an amplification bias, which revealed and highlighted the interpersonal individuality of the human gut virome. In contrast to a 16 S rRNA gene analysis of matched samples, we show that α- and β-diversity metrics of unamplified viromes are not as efficient at discerning controls from patients with inflammatory bowel disease. Additionally, we explored the intrinsic properties of unamplified gut viromes and show there is considerable interpersonal variability in viral taxa, infrequent longitudinal persistence of intrapersonal viruses, and vast fluctuations in the abundance of temporal viruses. Together, these properties of unamplified faecal viromes confound the ability to discern disease associations but significantly advance toward an unbiased and accurate representation of the human gut virome.

Population genetics and phylogeography of alfalfa mosaic virus in China and a comparison with other regional epidemics based on the cp gene.

Alfalfa mosaic virus (AMV) is the most pervasive epidemic virus affecting alfalfa production. However, detailed investigations on the molecular population genetics and evolutionary dynamics of AMV are scarce. This study aimed to report on a large-scale long-term survey of genetic variability in AMV populations from China and perform a comparative analysis of AMV population genetics in the three most thoroughly studied countries to date: China, Iran, and Spain. The study was based on the analysis of the coat protein gene (cp) using two analytical approaches: an analysis of molecular variance (AMOVA) and a Bayesian Markov Chain Monte Carlo approach that investigates the association between geographical origin and phylogeny. Both analytical approaches found significant genetic differentiation within localities, but not among localities nor among provinces. This observation might result from inappropriate agronomical practices involving extensive exchange of plant materials followed by rapid viral diversification within localities. In the Chinese population, both methods found that genetic diversification in AMV was strongly associated with different bioclimatic zones. Rates of molecular evolution were similar in the three countries. The estimated epidemic exponential population size and growth rate suggest that the epidemics grew faster and with higher incidence in Iran, followed by Spain and China. Estimates of the time to the most recent common ancestors suggest that AMV was first seen in Spain by the beginning of the twentieth century and later on in eastern and central Eurasia. After ruling out the existence of recombination breakpoints within the cp gene, a codon-based selection analysis per population was performed and identified many codons under significant negative selection and a few under significant positive selection; the latter varied among countries, suggesting regional differences in selective pressures.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

Contrasting Genetic Diversity of Listeria Pathogenicity Islands 3 and 4 Harbored by Nonpathogenic Listeria spp.

Listeria monocytogenes causes the severe foodborne disease listeriosis. Several clonal groups of L. monocytogenes possess the pathogenicity islands Listeria pathogenicity island 3 (LIPI-3) and LIPI-4. Here, we investigated the prevalence and genetic diversity of LIPI-3 and LIPI-4 among 63 strains of seven nonpathogenic Listeria spp. from the natural environment, i.e., wildlife (black bears [Ursus americanus]) and surface water. Analysis of the whole-genome sequence data suggested that both islands were horizontally acquired but differed considerably in their incidence and genetic diversity. LIPI-3 was identified among half of the L. innocua strains in the same genomic location as in L. monocytogenes (guaA hot spot) in a truncated form, with only three strains harboring full-length LIPI-3, and a highly divergent partial LIPI-3 was observed in three Listeria seeligeri strains, outside the guaA hot spot. Premature stop codons (PMSCs) and frameshifts were frequently noted in the LIPI-3 gene encoding listeriolysin S. On the other hand, full-length LIPI-4 without any PMSCs was found in all Listeria innocua strains, in the same genomic location as L. monocytogenes and with ~85% similarity to the L. monocytogenes counterpart. Our study provides intriguing examples of genetic changes that pathogenicity islands may undergo in nonpathogenic bacterial species, potentially in response to environmental pressures that promote either maintenance or degeneration of the islands. Investigations of the roles that LIPI-3 and LIPI-4 play in nonpathogenic Listeria spp. are warranted to further understand the differential evolution of genetic elements in pathogenic versus nonpathogenic hosts of the same genus. IMPORTANCE Listeria monocytogenes is a serious foodborne pathogen that can harbor the pathogenicity islands Listeria pathogenicity island 3 (LIPI-3) and LIPI-4. Intriguingly, these have also been reported in nonpathogenic L. innocua from food and farm environments, though limited information is available for strains from the natural environment. Here, we analyzed whole-genome sequence data of nonpathogenic Listeria spp. from wildlife and surface water to further elucidate the genetic diversity and evolution of LIPI-3 and LIPI-4 in Listeria. While the full-length islands were found only in L. innocua, LIPI-3 was uncommon and exhibited frequent truncation and genetic diversification, while LIPI-4 was remarkable in being ubiquitous, albeit diversified from L. monocytogenes. These contrasting features demonstrate that pathogenicity islands in nonpathogenic hosts can evolve along different trajectories, leading to either degeneration or maintenance, and highlight the need to examine their physiological roles in nonpathogenic hosts.

Expanded MutaT7 toolkit efficiently and simultaneously accesses all possible transition mutations in bacteria.

Targeted mutagenesis mediated by nucleotide base deaminase-T7 RNA polymerase fusions has recently emerged as a novel and broadly useful strategy to power genetic diversification in the context of in vivo directed evolution campaigns. Here, we expand the utility of this approach by introducing a highly active adenosine deaminase-T7 RNA polymerase fusion protein (eMutaT7A→G), resulting in higher mutation frequencies to enable more rapid directed evolution. We also assess the benefits and potential downsides of using this more active mutator. We go on to show in Escherichia coli that adenosine deaminase-bearing mutators (MutaT7A→G or eMutaT7A→G) can be employed in tandem with a cytidine deaminase-bearing mutator (MutaT7C→T) to introduce all possible transition mutations simultaneously. We illustrate the efficacy of this in vivo mutagenesis approach by exploring mutational routes to antibacterial drug resistance. This work sets the stage for general application of optimized MutaT7 tools able to induceall types oftransition mutations during in vivo directed evolution campaigns across diverse organisms.

Phylogenomic inference and demographic model selection suggest peripatric separation of the cryptic steppe ant species Plagiolepis pyrenaica stat. rev.

The ant Plagiolepis taurica Santschi, 1920 (Hymenoptera, Formicidae) is a typical species of the Eurasian steppes, a large grassland dominated biome that stretches continuously from Central Asia to Eastern Europe and is represented by disjunct outposts also in Central and Western Europe. The extent of this biome has been influenced by the Pleistocene climate, and steppes expanded recurrently during cold stages and contracted in warm stages. Consequently, stenotopic steppe species such as P. taurica repeatedly went through periods of demographic expansion and severe isolation. Here, we explore the impact of these dynamics on the genetic diversification within P. taurica. Delimitation of P. taurica from other Plagiolepis species has been unclear since its initial description, which raised questions on both its classification and its spatiotemporal diversification early on. We re-evaluate species limits and explore underlying mechanisms driving speciation by using an integrative approach based on genomic and morphometric data. We found large intraspecific divergence within P. taurica and resolved geographically coherent western and eastern genetic groups, which likewise differed morphologically. A morphometric survey of type material showed that Plagiolepis from the western group were more similar to P. barbara pyrenaica Emery, 1921 than to P. taurica; we thus lift the former from synonymy and establish it as separate species, P. pyrenaica stat. rev. Explicit evolutionary model testing based on genomic data supported a peripatric speciation for the species pair, probably as a consequence of steppe contraction and isolation during the mid-Pleistocene. We speculate that this scenario could be exemplary for many stenotopic steppe species, given the emphasized dynamics of Eurasian steppes.

Recombination Drives Evolution of Carbapenem-Resistant Klebsiella pneumoniae Sequence Type 11 KL47 to KL64 in China.

Carbapenem-resistant Enterobacteriaceae, especially carbapenemase-producing Klebsiella pneumoniae, is an urgent problem in health care facilities worldwide. K. pneumoniae isolates classified as sequence type 11 (ST11) are largely responsible for the spread of carbapenem-resistant K. pneumoniae (CRKP) in China. Our previous phylogenetic reconstruction suggested that CRKP ST11 capsular locus 64 (KL64) was derived from an ST11-KL47 ancestor through recombination. However, the molecular origin of KL64 remains largely unknown, and our understanding of the recombination is incomplete. Here, we screened a global sample of 22,600 K. pneumoniae genomes and searched for KL64-harboring STs, which were found to be ST1764, ST3685, ST1764-1LV, ST30, ST505, ST147, and ST11, wherein ST1764, ST3685, ST1764-1LV, and ST30 belonged to a clonal complex, CC1764. We compared the genetic structures of representative strains from ST11-KL47, ST11-KL64, CC1764-KL64, ST505-KL64, and ST147-KL64 and further performed phylogenetic analysis and single-nucleotide polymorphism analysis among 248 isolates from all these STs. The results suggested a recombination event has occurred in a homologous ~154-kb region covering KL and the lipopolysaccharide biosynthesis locus (OL) between a recipient ST11-KL47-OL101 and a donor CC1764 (except ST30), giving rise to ST11-KL64-O2v1 strains (recombination I). Furthermore, we also found an infrequent ST11-KL64-O2v1 subclone which was not produced by recombination I but was hybridized from ST11-KL47-OL101 and ST147-KL64-O2v1 strains through recombination of a homologous ~485-kb region covering KL and OL (recombination II). Our findings provide important insights into the role of recombination in the evolution of clinical strains and the diversity of capsule and lipopolysaccharide of widely distributed KPC-associated ST11 K. pneumoniae in China. IMPORTANCE Chromosomal recombination events are considered to contribute to the genetic diversification and ultimate success of many bacterial pathogens. A previous study unravelled the molecular evolution history of ST258 strains, which have been largely responsible for the spread of KPC in the United States. Here, we used comparative genomic analyses to discover two recombination events in ST11 CRKP strains, which is a predominant KPC-associated CRKP clone in China. Two new ST11-CRKP subclones with altered capsule and lipopolysaccharide, which are two primary determinants of antigenicity and antigenic diversity among K. pneumoniae, were produced through these two recombination events, respectively. Horizontal transfer of the KL and OL appears to be a crucial element driving the molecular evolution of K. pneumoniae strains. These findings not only extend our understanding of the molecular evolutionary history of ST11 but also are an important step toward the development of preventive, diagnostic, and therapeutic strategies for CRKP.