Presence of phage-plasmids in multiple serovars of Salmonella enterica.

Evidence is accumulating in the literature that the horizontal spread of antimicrobial resistance (AMR) genes mediated by bacteriophages and bacteriophage-like plasmid (phage-plasmid) elements is much more common than previously envisioned. For instance, we recently identified and characterized a circular P1-like phage-plasmid harbouring a bla CTX-M-15 gene conferring extended-spectrum beta-lactamase (ESBL) resistance in Salmonella enterica serovar Typhi. As the prevalence and epidemiological relevance of such mechanisms has never been systematically assessed in Enterobacterales, in this study we carried out a follow-up retrospective analysis of UK Salmonella isolates previously sequenced as part of routine surveillance protocols between 2016 and 2021. Using a high-throughput bioinformatics pipeline we screened 47 784 isolates for the presence of the P1 lytic replication gene repL, identifying 226 positive isolates from 25 serovars and demonstrating that phage-plasmid elements are more frequent than previously thought. The affinity for phage-plasmids appears highly serovar-dependent, with several serovars being more likely hosts than others; most of the positive isolates (170/226) belonged to S. Typhimurium ST34 and ST19. The phage-plasmids ranged between 85.8 and 98.2 kb in size, with an average length of 92.1 kb; detailed analysis indicated a high amount of diversity in gene content and genomic architecture. In total, 132 phage-plasmids had the p0111 plasmid replication type, and 94 the IncY type; phylogenetic analysis indicated that both horizontal and vertical gene transmission mechanisms are likely to be involved in phage-plasmid propagation. Finally, phage-plasmids were present in isolates that were resistant and non-resistant to antimicrobials. In addition to providing a first comprehensive view of the presence of phage-plasmids in Salmonella, our work highlights the need for a better surveillance and understanding of phage-plasmids as AMR carriers, especially through their characterization with long-read sequencing.

Sequencing and Analysis of Lumpy Skin Disease Virus Whole Genomes Reveals a New Viral Subgroup in West and Central Africa.

Lumpy skin disease virus (LSDV) is a member of the capripoxvirus (CPPV) genus of the Poxviridae family. LSDV is a rapidly emerging, high-consequence pathogen of cattle, recently spreading from Africa and the Middle East into Europe and Asia. We have sequenced the whole genome of historical LSDV isolates from the Pirbright Institute virus archive, and field isolates from recent disease outbreaks in Sri Lanka, Mongolia, Nigeria and Ethiopia. These genome sequences were compared to published genomes and classified into different subgroups. Two subgroups contained vaccine or vaccine-like samples ("Neethling-like" clade 1.1 and "Kenya-like" subgroup, clade 1.2.2). One subgroup was associated with outbreaks of LSD in the Middle East/Europe (clade 1.2.1) and a previously unreported subgroup originated from cases of LSD in west and central Africa (clade 1.2.3). Isolates were also identified that contained a mix of genes from both wildtype and vaccine samples (vaccine-like recombinants, grouped in clade 2). Whole genome sequencing and analysis of LSDV strains isolated from different regions of Africa, Europe and Asia have provided new knowledge of the drivers of LSDV emergence, and will inform future disease control strategies.

KmerAperture: Retaining k-mer synteny for alignment-free extraction of core and accessory differences between bacterial genomes.

By decomposing genome sequences into k-mers, it is possible to estimate genome differences without alignment. Techniques such as k-mer minimisers, for example MinHash, have been developed and are often accurate approximations of distances based on full k-mer sets. These and other alignment-free methods avoid the large temporal and computational expense of alignment. However, these k-mer set comparisons are not entirely accurate within-species and can be completely inaccurate within-lineage. This is due, in part, to their inability to distinguish core polymorphism from accessory differences. Here we present a new approach, KmerAperture, which uses information on the k-mer relative genomic positions to determine the type of polymorphism causing differences in k-mer presence and absence between pairs of genomes. Single SNPs are expected to result in k unique contiguous k-mers per genome. On the other hand, contiguous series > k may be caused by accessory differences of length S-k+1; when the start and end of the sequence are contiguous with homologous sequence. Alternatively, they may be caused by multiple SNPs within k bp from each other and KmerAperture can determine whether that is the case. To demonstrate use cases KmerAperture was benchmarked using datasets including a very low diversity simulated population with accessory content independent from the number of SNPs, a simulated population where SNPs are spatially dense, a moderately diverse real cluster of genomes (Escherichia coli ST1193) with a large accessory genome and a low diversity real genome cluster (Salmonella Typhimurium ST34). We show that KmerAperture can accurately distinguish both core and accessory sequence diversity without alignment, outperforming other k-mer based tools.

Inference of Infectious Disease Transmission through a Relaxed Bottleneck Using Multiple Genomes Per Host.

In recent times, pathogen genome sequencing has become increasingly used to investigate infectious disease outbreaks. When genomic data is sampled densely enough amongst infected individuals, it can help resolve who infected whom. However, transmission analysis cannot rely solely on a phylogeny of the genomes but must account for the within-host evolution of the pathogen, which blurs the relationship between phylogenetic and transmission trees. When only a single genome is sampled for each host, the uncertainty about who infected whom can be quite high. Consequently, transmission analysis based on multiple genomes of the same pathogen per host has a clear potential for delivering more precise results, even though it is more laborious to achieve. Here, we present a new methodology that can use any number of genomes sampled from a set of individuals to reconstruct their transmission network. Furthermore, we remove the need for the assumption of a complete transmission bottleneck. We use simulated data to show that our method becomes more accurate as more genomes per host are provided, and that it can infer key infectious disease parameters such as the size of the transmission bottleneck, within-host growth rate, basic reproduction number, and sampling fraction. We demonstrate the usefulness of our method in applications to real datasets from an outbreak of Pseudomonas aeruginosa amongst cystic fibrosis patients and a nosocomial outbreak of Klebsiella pneumoniae.

Allelic Entropy and Times until Loss or Fixation of Neutral Polymorphisms.

In the diffusion approximation of the neutral Wright-Fisher model, the expected time until fixation or loss of a neutral allele is proportional to the initial entropy of the distribution of the allele in the population. No explanation is known for this coincidence. In this paper, we show that the rate of entropy dissipation is proportional to the number of segregating alleles. Since the final fixed state has zero entropy, the expected lifetime of segregating alleles is proportional to the initial entropy in the system. We show that classical formulae on the average time to loss of segregating alleles and the expected time to fixation of the last segregating allele stem from these properties of the diffusion process. We also extend our results to the case of population size changing in time. The dissipation of heterozygosity and entropy shows that superlinear population growth leads to infinite expected fixation times, i.e., neutral alleles in fast-growing populations could segregate forever without ever becoming fixed or disappearing by genetic drift.

A novel factor modulating X chromosome dosage compensation in Anopheles.

Dosage compensation (DC), a process countering chromosomal imbalance in individuals with heteromorphic sex chromosomes, has been molecularly characterized only in mammals, Caenorhabditis elegans, and fruit flies.1 In Drosophila melanogaster males, it is achieved by an approximately 2-fold hypertranscription of the monosomic X chromosome mediated by the MSL complex.2,3 The complex is not assembled on female X chromosomes because production of its key protein MSL-2 is prevented due to intron retention and inhibition of translation by Sex-lethal, a female-specific protein operating at the top of the sex determination pathway.4 It remains unclear how DC is mechanistically regulated in other insects. In the malaria mosquito Anopheles gambiae, an approximately 2-fold hypertranscription of the male X also occurs5 by a yet-unknown molecular mechanism distinct from that in D. melanogaster.6 Here we show that a male-specifically spliced gene we call 007, which arose by a tandem duplication in the Anopheles ancestral lineage, is involved in the control of DC in males. Homozygous 007 knockouts lead to a global downregulation of the male X, phenotypically manifested by a slower development compared to wild-type mosquitoes or mutant females-however, without loss of viability or fertility. In females, a 007 intron retention promoted by the sex determination protein Femaleless, known to prevent hypertranscription from both X chromosomes,7 introduces a premature termination codon apparently rendering the female transcripts non-productive. In addition to providing a unique perspective on DC evolution, the 007, with its conserved properties, may represent an important addition to a genetic toolbox for malaria vector control.

Distinct monkeypox virus lineages co-circulating in humans before 2022.

The 2022 global mpox outbreak raises questions about how this zoonotic disease established effective human-to-human transmission and its potential for further adaptation. The 2022 outbreak virus is related to an ongoing outbreak in Nigeria originally reported in 2017, but the evolutionary path linking the two remains unclear due to a lack of genomic data between 2018, when virus exportations from Nigeria were first recorded, and 2022, when the global mpox outbreak began. Here, 18 viral genomes obtained from patients across southern Nigeria in 2019-2020 reveal multiple lineages of monkeypox virus (MPXV) co-circulated in humans for several years before 2022, with progressive accumulation of mutations consistent with APOBEC3 activity over time. We identify Nigerian A.2 lineage isolates, confirming the lineage that has been multiply exported to North America independently of the 2022 outbreak originated in Nigeria, and that it has persisted by human-to-human transmission in Nigeria for more than 2 years before its latest exportation. Finally, we identify a lineage-defining APOBEC3-style mutation in all A.2 isolates that disrupts gene A46R, encoding a viral innate immune modulator. Collectively, our data demonstrate MPXV capacity for sustained diversification within humans, including mutations that may be consistent with established mechanisms of poxvirus adaptation.

Distinguishing imported cases from locally acquired cases within a geographically limited genomic sample of an infectious disease.

MOTIVATION: The ability to distinguish imported cases from locally acquired cases has important consequences for the selection of public health control strategies. Genomic data can be useful for this, for example, using a phylogeographic analysis in which genomic data from multiple locations are compared to determine likely migration events between locations. However, these methods typically require good samples of genomes from all locations, which is rarely available. RESULTS: Here, we propose an alternative approach that only uses genomic data from a location of interest. By comparing each new case with previous cases from the same location, we are able to detect imported cases, as they have a different genealogical distribution than that of locally acquired cases. We show that, when variations in the size of the local population are accounted for, our method has good sensitivity and excellent specificity for the detection of imports. We applied our method to data simulated under the structured coalescent model and demonstrate relatively good performance even when the local population has the same size as the external population. Finally, we applied our method to several recent genomic datasets from both bacterial and viral pathogens, and show that it can, in a matter of seconds or minutes, deliver important insights on the number of imports to a geographically limited sample of a pathogen population. AVAILABILITY AND IMPLEMENTATION: The R package DetectImports is freely available from https://github.com/xavierdidelot/DetectImports. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of a P1-bacteriophage-like plasmid (phage-plasmid) harbouring bla CTX-M-15 in Salmonella enterica serovar Typhi.

Antimicrobial-resistance (AMR) genes can be transferred between microbial cells via horizontal gene transfer (HGT), which involves mobile and integrative elements such as plasmids, bacteriophages, transposons, integrons and pathogenicity islands. Bacteriophages are found in abundance in the microbial world, but their role in virulence and AMR has not fully been elucidated in the Enterobacterales. With short-read sequencing paving the way to systematic high-throughput AMR gene detection, long-read sequencing technologies now enable us to establish how such genes are structurally connected into meaningful genomic units, raising questions about how they might cooperate to achieve their biological function. Here, we describe a novel ~98 kbp circular P1-bacteriophage-like plasmid termed ph681355 isolated from a clinical Salmonella enterica serovar Typhi isolate. It carries bla CTX-M-15, an IncY plasmid replicon (repY gene) and the ISEcP1 mobile element and is, to our knowledge, the first reported P1-bacteriophage-like plasmid (phage-plasmid) in S. enterica Typhi. We compared ph681355 to two previously described phage-plasmids, pSJ46 from S. enterica serovar Indiana and pMCR-1-P3 from Escherichia coli, and found high nucleotide similarity across the backbone. However, we saw low ph681355 backbone similarity to plasmid p60006 associated with the extensively drug-resistant S. enterica Typhi outbreak isolate in Pakistan, providing evidence of an alternative route for bla CTX-M-15 transmission. Our discovery highlights the importance of utilizing long-read sequencing in interrogating bacterial genomic architecture to fully understand AMR mechanisms and their clinical relevance. It also raises questions regarding how widespread bacteriophage-mediated HGT might be, suggesting that the resulting genomic plasticity might be higher than previously thought.

Transcriptomic Profiling of Dromedary Camels Immunised with a MERS Vaccine Candidate.

Middle East Respiratory Syndrome coronavirus (MERS-CoV) infects dromedary camels and zoonotically infects humans, causing a respiratory disease with severe pneumonia and death. With no approved antiviral or vaccine interventions for MERS, vaccines are being developed for camels to prevent virus transmission into humans. We have previously developed a chimpanzee adenoviral vector-based vaccine for MERS-CoV (ChAdOx1 MERS) and reported its strong humoral immunogenicity in dromedary camels. Here, we looked back at total RNA isolated from whole blood of three immunised dromedaries pre and post-vaccination during the first day; and performed RNA sequencing and bioinformatic analysis in order to shed light on the molecular immune responses following a ChAdOx1 MERS vaccination. Our finding shows that a number of transcripts were differentially regulated as an effect of the vaccination, including genes that are involved in innate and adaptive immunity, such as type I and II interferon responses. The camel Bcl-3 and Bcl-6 transcripts were significantly upregulated, indicating a strong activation of Tfh cell, B cell, and NF-κB pathways. In conclusion, this study gives an overall view of the first changes in the immune transcriptome of dromedaries after vaccination; it supports the potency of ChAdOx1 MERS as a potential camel vaccine to block transmission and prevent new human cases and outbreaks.

Mutagenesis Mapping of RNA Structures within the Foot-and-Mouth Disease Virus Genome Reveals Functional Elements Localized in the Polymerase (3Dpol)-Encoding Region.

RNA structures can form functional elements that play crucial roles in the replication of positive-sense RNA viruses. While RNA structures in the untranslated regions (UTRs) of several picornaviruses have been functionally characterized, the roles of putative RNA structures predicted for protein coding sequences (or open reading frames [ORFs]) remain largely undefined. Here, we have undertaken a bioinformatic analysis of the foot-and-mouth disease virus (FMDV) genome to predict 53 conserved RNA structures within the ORF. Forty-six of these structures were located in the regions encoding the nonstructural proteins (nsps). To investigate whether structures located in the regions encoding the nsps are required for FMDV replication, we used a mutagenesis method, CDLR mapping, where sequential coding segments were shuffled to minimize RNA secondary structures while preserving protein coding, native dinucleotide frequencies, and codon usage. To examine the impact of these changes on replicative fitness, mutated sequences were inserted into an FMDV subgenomic replicon. We found that three of the RNA structures, all at the 3' termini of the FMDV ORF, were critical for replicon replication. In contrast, disruption of the other 43 conserved RNA structures that lie within the regions encoding the nsps had no effect on replicon replication, suggesting that these structures are not required for initiating translation or replication of viral RNA. Conserved RNA structures that are not essential for virus replication could provide ideal targets for the rational attenuation of a wide range of FMDV strains. IMPORTANCE Some RNA structures formed by the genomes of RNA viruses are critical for viral replication. Our study shows that of 46 conserved RNA structures located within the regions of the foot-and-mouth disease virus (FMDV) genome that encode the nonstructural proteins, only three are essential for replication of an FMDV subgenomic replicon. Replicon replication is dependent on RNA translation and synthesis; thus, our results suggest that the three RNA structures are critical for either initiation of viral RNA translation and/or viral RNA synthesis. Although further studies are required to identify whether the remaining 43 RNA structures have other roles in virus replication, they may provide targets for the rational large-scale attenuation of a wide range of FMDV strains. FMDV causes a highly contagious disease, posing a constant threat to global livestock industries. Such weakened FMDV strains could be investigated as live-attenuated vaccines or could enhance biosecurity of conventional inactivated vaccine production.

Engineered Promoter-Switched Viruses Reveal the Role of Poxvirus Maturation Protein A26 as a Negative Regulator of Viral Spread.

Vaccinia virus produces two types of virions known as single-membraned intracellular mature virus (MV) and double-membraned extracellular enveloped virus (EV). EV production peaks earlier when initial MVs are further wrapped and secreted to spread infection within the host. However, late during infection, MVs accumulate intracellularly and become important for host-to-host transmission. The process that regulates this switch remains elusive and is thought to be influenced by host factors. Here, we examined the hypothesis that EV and MV production are regulated by the virus through expression of F13 and the MV-specific protein A26. By switching the promoters and altering the expression kinetics of F13 and A26, we demonstrate that A26 expression downregulates EV production and plaque size, thus limiting viral spread. This process correlates with A26 association with the MV surface protein A27 and exclusion of F13, thus reducing EV titers. Thus, MV maturation is controlled by the abundance of the viral A26 protein, independently of other factors, and is rate limiting for EV production. The A26 gene is conserved within vertebrate poxviruses but is strikingly lost in poxviruses known to be transmitted exclusively by biting arthropods. A26-mediated virus maturation thus has the appearance to be an ancient evolutionary adaptation to enhance transmission of poxviruses that has subsequently been lost from vector-adapted species, for which it may serve as a genetic signature. The existence of virus-regulated mechanisms to produce virions adapted to fulfill different functions represents a novel level of complexity in mammalian viruses with major impacts on evolution, adaptation, and transmission. IMPORTANCE Chordopoxviruses are mammalian viruses that uniquely produce a first type of virion adapted to spread within the host and a second type that enhances transmission between hosts, which can take place by multiple ways, including direct contact, respiratory droplets, oral/fecal routes, or via vectors. Both virion types are important to balance intrahost dissemination and interhost transmission, so virus maturation pathways must be tightly controlled. Here, we provide evidence that the abundance and kinetics of expression of the viral protein A26 regulates this process by preventing formation of the first form and shifting maturation toward the second form. A26 is expressed late after the initial wave of progeny virions is produced, so sufficient viral dissemination is ensured, and A26 provides virions with enhanced environmental stability. Conservation of A26 in all vertebrate poxviruses, but not in those transmitted exclusively via biting arthropods, reveals the importance of A26-controlled virus maturation for transmission routes involving environmental exposure.

Whole genome de novo sequencing and comparative genomic analyses suggests that Chlamydia psittaci strain 84/2334 should be reclassified as Chlamydia abortus species.

BACKGROUND: Chlamydia abortus and Chlamydia psittaci are important pathogens of livestock and avian species, respectively. While C. abortus is recognized as descended from C. psittaci species, there is emerging evidence of strains that are intermediary between the two species, suggesting they are recent evolutionary ancestors of C. abortus. Such strains include C. psittaci strain 84/2334 that was isolated from a parrot. Our aim was to classify this strain by sequencing its genome and explore its evolutionary relationship to both C. abortus and C. psittaci. RESULTS: In this study, methods based on multi-locus sequence typing (MLST) of seven housekeeping genes and on typing of five species discriminant proteins showed that strain 84/2334 clustered with C. abortus species. Furthermore, whole genome de novo sequencing of the strain revealed greater similarity to C. abortus in terms of GC content, while 16S rRNA and whole genome phylogenetic analysis, as well as network and recombination analysis showed that the strain clusters more closely with C. abortus strains. The analysis also suggested a closer evolutionary relationship between this strain and the major C. abortus clade, than to two other intermediary avian C. abortus strains or C. psittaci strains. Molecular analyses of genes (polymorphic membrane protein and transmembrane head protein genes) and loci (plasticity zone), found in key virulence-associated regions that exhibit greatest diversity within and between chlamydial species, reveal greater diversity than present in sequenced C. abortus genomes as well as similar features to both C. abortus and C. psittaci species. The strain also possesses an extrachromosomal plasmid, as found in most C. psittaci species but absent from all sequenced classical C. abortus strains. CONCLUSION: Overall, the results show that C. psittaci strain 84/2334 clusters very closely with C. abortus strains, and are consistent with the strain being a recent C. abortus ancestral species. This suggests that the strain should be reclassified as C. abortus. Furthermore, the identification of a C. abortus strain bearing an extra-chromosomal plasmid has implications for plasmid-based transformation studies to investigate gene function as well as providing a potential route for the development of a next generation vaccine to protect livestock from C. abortus infection.

Correction: Pervasive within-host recombination and epistasis as major determinants of the molecular evolution of the foot-and-mouth disease virus capsid.

[This corrects the article DOI: 10.1371/journal.ppat.1008235.].

A Systematic Evaluation of High-Throughput Sequencing Approaches to Identify Low-Frequency Single Nucleotide Variants in Viral Populations.

High-throughput sequencing such as those provided by Illumina are an efficient way to understand sequence variation within viral populations. However, challenges exist in distinguishing process-introduced error from biological variance, which significantly impacts our ability to identify sub-consensus single-nucleotide variants (SNVs). Here we have taken a systematic approach to evaluate laboratory and bioinformatic pipelines to accurately identify low-frequency SNVs in viral populations. Artificial DNA and RNA "populations" were created by introducing known SNVs at predetermined frequencies into template nucleic acid before being sequenced on an Illumina MiSeq platform. These were used to assess the effects of abundance and starting input material type, technical replicates, read length and quality, short-read aligner, and percentage frequency thresholds on the ability to accurately call variants. Analyses revealed that the abundance and type of input nucleic acid had the greatest impact on the accuracy of SNV calling as measured by a micro-averaged Matthews correlation coefficient score, with DNA and high RNA inputs (107 copies) allowing for variants to be called at a 0.2% frequency. Reduced input RNA (105 copies) required more technical replicates to maintain accuracy, while low RNA inputs (103 copies) suffered from consensus-level errors. Base errors identified at specific motifs identified in all technical replicates were also identified which can be excluded to further increase SNV calling accuracy. These findings indicate that samples with low RNA inputs should be excluded for SNV calling and reinforce the importance of optimising the technical and bioinformatics steps in pipelines that are used to accurately identify sequence variants.

Full genome sequencing of archived wild type and vaccine rinderpest virus isolates prior to their destruction.

When rinderpest virus (RPV) was declared eradicated in 2011, the only remaining samples of this once much-feared livestock virus were those held in various laboratories. In order to allow the destruction of our institute's stocks of RPV while maintaining the ability to recover the various viruses if ever required, we have determined the full genome sequence of all our distinct samples of RPV, including 51 wild type viruses and examples of three different types of vaccine strain. Examination of the sequences of these virus isolates has shown that the African isolates form a single disparate clade, rather than two separate clades, which is more in accord with the known history of the virus in Africa. We have also identified two groups of goat-passaged viruses which have acquired an extra 6 bases in the long untranslated region between the M and F protein coding sequences, and shown that, for more than half the genomes sequenced, translation of the F protein requires translational frameshift or non-standard translation initiation. Curiously, the clade containing the lapinised vaccine viruses that were developed originally in Korea appears to be more similar to the known African viruses than to any other Asian viruses.

Pervasive Differential Splicing in Marek's Disease Virus can Discriminate CVI-988 Vaccine Strain from RB-1B Very Virulent Strain in Chicken Embryonic Fibroblasts.

Marek's disease is a major scourge challenging poultry health worldwide. It is caused by the highly contagious Marek's disease virus (MDV), an alphaherpesvirus. Here, we showed that, similar to other members of its Herpesviridae family, MDV also presents a complex landscape of splicing events, most of which are uncharacterised and/or not annotated. Quite strikingly, and although the biological relevance of this fact is unknown, we found that a number of viral splicing isoforms are strain-specific, despite the close sequence similarity of the strains considered: very virulent RB-1B and vaccine CVI-988. We validated our findings by devising an assay that discriminated infections caused by the two strains in chicken embryonic fibroblasts on the basis of the presence of some RNA species. To our knowledge, this study is the first to accomplish such a result, emphasizing how relevant a comprehensive picture of the viral transcriptome is to fully understand viral pathogenesis.

Pervasive within-host recombination and epistasis as major determinants of the molecular evolution of the foot-and-mouth disease virus capsid.

Although recombination is known to occur in foot-and-mouth disease virus (FMDV), it is considered only a minor determinant of virus sequence diversity. Analysis at phylogenetic scales shows inter-serotypic recombination events are rare, whereby recombination occurs almost exclusively in non-structural proteins. In this study we have estimated recombination rates within a natural host in an experimental setting. African buffaloes were inoculated with a SAT-1 FMDV strain containing two major viral sub-populations differing in their capsid sequence. This population structure enabled the detection of extensive within-host recombination in the genomic region coding for structural proteins and allowed recombination rates between the two sub-populations to be estimated. Quite surprisingly, the effective recombination rate in VP1 during the acute infection phase turns out to be about 0.1 per base per year, i.e. comparable to the mutation/substitution rate. Using a high-resolution map of effective within-host recombination in the capsid-coding region, we identified a linkage disequilibrium pattern in VP1 that is consistent with a mosaic structure with two main genetic blocks. Positive epistatic interactions between co-evolved variants appear to be present both within and between blocks. These interactions are due to intra-host selection both at the RNA and protein level. Overall our findings show that during FMDV co-infections by closely related strains, capsid-coding genes recombine within the host at a much higher rate than expected, despite the presence of strong constraints dictated by the capsid structure. Although these intra-host results are not immediately translatable to a phylogenetic setting, recombination and epistasis must play a major and so far underappreciated role in the molecular evolution of the virus at all scales.

Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence.

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

Persistent Infection of African Buffalo (Syncerus caffer) with Foot-and-Mouth Disease Virus: Limited Viral Evolution and No Evidence of Antibody Neutralization Escape.

African buffaloes (Syncerus caffer) are the principal "carrier" hosts of foot-and-mouth disease virus (FMDV). Currently, the epithelia and lymphoid germinal centers of the oropharynx have been identified as sites for FMDV persistence. We carried out studies in FMDV SAT1 persistently infected buffaloes to characterize the diversity of viruses in oropharyngeal epithelia, germinal centers, probang samples (oropharyngeal scrapings), and tonsil swabs to determine if sufficient virus variation is generated during persistence for immune escape. Most sequencing reads of the VP1 coding region of the SAT1 virus inoculum clustered around 2 subpopulations differing by 22 single-nucleotide variants of intermediate frequency. Similarly, most sequences from oropharynx tissue clustered into two subpopulations, albeit with different proportions, depending on the day postinfection (dpi). There was a significant difference between the populations of viruses in the inoculum and in lymphoid tissue taken at 35 dpi. Thereafter, until 400 dpi, no significant variation was detected in the viral populations in samples from individual animals, germinal centers, and epithelial tissues. Deep sequencing of virus from probang or tonsil swab samples harvested prior to postmortem showed less within-sample variability of VP1 than that of tissue sample sequences analyzed at the same time. Importantly, there was no significant difference in the ability of sera collected between 14 and 400 dpi to neutralize the inoculum or viruses isolated at later time points in the study from the same animal. Therefore, based on this study, there is no evidence of escape from antibody neutralization contributing to FMDV persistent infection in African buffalo.IMPORTANCE Foot-and-mouth disease virus (FMDV) is a highly contagious virus of cloven-hoofed animals and is recognized as the most important constraint to international trade in animals and animal products. African buffaloes (Syncerus caffer) are efficient carriers of FMDV, and it has been proposed that new virus variants are produced in buffalo during the prolonged carriage after acute infection, which may spread to cause disease in livestock populations. Here, we show that despite an accumulation of low-frequency sequence variants over time, there is no evidence of significant antigenic variation leading to immune escape. Therefore, carrier buffalo are unlikely to be a major source of new virus variants.