The many definitions of multiplicity of infection

Soil Giant Phage: Genome and Biological Characteristics of Sinorhizobium Jumbo Phage

This paper presents the first in-depth research on the biological and genomic properties of lytic rhizobiophage AP-J-162 isolated from the soils of the mountainous region of Dagestan (North Caucasus), which belongs to the centers of origin of cultivated plants, according to Vavilov N.I. The rhizobiophage host strains are nitrogen-fixing bacteria of the genus Sinorhizobium spp., symbionts of leguminous forage grasses. The phage particles have a myovirus virion structure. The genome of rhizobiophage AP-J-162 is double-stranded DNA of 471.5 kb in length; 711 ORFs are annotated and 41 types of tRNAs are detected. The closest phylogenetic relative of phage AP-J-162 is Agrobacterium phage Atu-ph07, but no rhizobiophages are known. The replicative machinery, capsid, and baseplate proteins of phage AP-J-162 are structurally similar to those of Escherichia phage T4, but there is no similarity between their tail protein subunits. Amino acid sequence analysis shows that 339 of the ORFs encode hypothetical or functionally relevant products, while the remaining 304 ORFs are unique. Additionally, 153 ORFs are similar to those of Atu_ph07, with one-third of the ORFs encoding different enzymes. The biological properties and genomic characteristics of phage AP-J-162 distinguish it as a unique model for exploring phage-microbe interactions with nitrogen-fixing symbiotic microorganisms.

SNP-slice resolves mixed infections: simultaneously unveiling strain haplotypes and linking them to hosts

MOTIVATION

Multi-strain infection is a common yet under-investigated phenomenon of many pathogens. Currently, biologists analyzing SNP information sometimes have to discard mixed infection samples as many downstream analyses require monogenomic inputs. Such a protocol impedes our understanding of the underlying genetic diversity, co-infection patterns, and genomic relatedness of pathogens. A scalable tool to learn and resolve the SNP-haplotypes from polygenomic data is an urgent need in molecular epidemiology.

RESULTS

We develop a slice sampling Markov Chain Monte Carlo algorithm, named SNP-Slice, to learn not only the SNP-haplotypes of all strains in the populations but also which strains infect which hosts. Our method reconstructs SNP-haplotypes and individual heterozygosities accurately without reference panels and outperforms the state-of-the-art methods at estimating the multiplicity of infections and allele frequencies. Thus, SNP-Slice introduces a novel approach to address polygenomic data and opens a new avenue for resolving complex infection patterns in molecular surveillance. We illustrate the performance of SNP-Slice on empirical malaria and HIV datasets and provide recommendations for using our method on empirical datasets.

AVAILABILITY AND IMPLEMENTATION

The implementation of the SNP-Slice algorithm, as well as scripts to analyze SNP-Slice outputs, are available at https://github.com/nianqiaoju/snp-slice.

Estimating multiplicity of infection, allele frequencies, and prevalences accounting for incomplete data

BACKGROUND

Molecular surveillance of infectious diseases allows the monitoring of pathogens beyond the granularity of traditional epidemiological approaches and is well-established for some of the most relevant infectious diseases such as malaria. The presence of genetically distinct pathogenic variants within an infection, referred to as multiplicity of infection (MOI) or complexity of infection (COI) is common in malaria and similar infectious diseases. It is an important metric that scales with transmission intensities, potentially affects the clinical pathogenesis, and a confounding factor when monitoring the frequency and prevalence of pathogenic variants. Several statistical methods exist to estimate MOI and the frequency distribution of pathogen variants. However, a common problem is the quality of the underlying molecular data. If molecular assays fail not randomly, it is likely to underestimate MOI and the prevalence of pathogen variants.

METHODS AND FINDINGS

A statistical model is introduced, which explicitly addresses data quality, by assuming a probability by which a pathogen variant remains undetected in a molecular assay. This is different from the assumption of missing at random, for which a molecular assay either performs perfectly or fails completely. The method is applicable to a single molecular marker and allows to estimate allele-frequency spectra, the distribution of MOI, and the probability of variants to remain undetected (incomplete information). Based on the statistical model, expressions for the prevalence of pathogen variants are derived and differences between frequency and prevalence are discussed. The usual desirable asymptotic properties of the maximum-likelihood estimator (MLE) are established by rewriting the model into an exponential family. The MLE has promising finite sample properties in terms of bias and variance. The covariance matrix of the estimator is close to the Cramér-Rao lower bound (inverse Fisher information). Importantly, the estimator's variance is larger than that of a similar method which disregards incomplete information, but its bias is smaller.

CONCLUSIONS

Although the model introduced here has convenient properties, in terms of the mean squared error it does not outperform a simple standard method that neglects missing information. Thus, the new method is recommendable only for data sets in which the molecular assays produced poor-quality results. This will be particularly true if the model is extended to accommodate information from multiple molecular markers at the same time, and incomplete information at one or more markers leads to a strong depletion of sample size.

Evolutionary genetics of malaria

Many standard-textbook population-genetic results apply to a wide range of species. Sometimes, however, population-genetic models and principles need to be tailored to a particular species. This is particularly true for malaria, which next to tuberculosis and HIV/AIDS ranks among the economically most relevant infectious diseases. Importantly, malaria is not one disease—five human-pathogenic species of Plasmodium exist. P. falciparum is not only the most severe form of human malaria, but it also causes the majority of infections. The second most relevant species, P. vivax, is already considered a neglected disease in several endemic areas. All human-pathogenic species have distinct characteristics that are not only crucial for control and eradication efforts, but also for the population-genetics of the disease. This is particularly true in the context of selection. Namely, fitness is determined by so-called fitness components, which are determined by the parasites live-history, which differs between malaria species. The presence of hypnozoites, i.e., dormant liver-stage parasites, which can cause disease relapses, is a distinct feature of P. vivax and P. ovale sp. In P. malariae inactivated blood-stage parasites can cause a recrudescence years after the infection was clinically cured. To properly describe population-genetic processes, such as the spread of anti-malarial drug resistance, these features must be accounted for appropriately. Here, we introduce and extend a population-genetic framework for the evolutionary dynamics of malaria, which applies to all human-pathogenic malaria species. The model focuses on, but is not limited to, the spread of drug resistance. The framework elucidates how the presence of dormant liver stage or inactivated blood stage parasites that act like seed banks delay evolutionary processes. It is shown that, contrary to standard population-genetic theory, the process of selection and recombination cannot be decoupled in malaria. Furthermore, we discuss the connection between haplotype frequencies, haplotype prevalence, transmission dynamics, and relapses or recrudescence in malaria.