High-Throughput Genotyping of Plasmodium vivax in the Peruvian Amazon via Molecular Inversion Probes.

Plasmodium vivax transmission occurs throughout the tropics and is an emerging threat in areas of Plasmodium falciparum decline, causing relapse infections that complicate treatment and control. Targeted sequencing for P. falciparum has been widely deployed to detect population structure and the geographic spread of antimalarial and diagnostic resistance. However, there are fewer such tools for P. vivax . Leveraging global variation data, we designed four molecular inversion probe (MIP) genotyping panels targeting geographically differentiating SNPs, neutral SNPs, putative antimalarial resistance genes, and vaccine candidate genes. We deployed these MIP panels on 866 infections from the Peruvian Amazon and identified transmission networks with clonality (IBD>0.99), copy number variation in Pvdbp and multiple Pvrbps , fixation of putative antimalarial resistance, and balancing selection in 13 vaccine candidate genes. Our MIP panels are the broadest genotyping panel currently available and are poised for successful deployment in other regions of P. vivax transmission.

Country wide surveillance reveals prevalent artemisinin partial resistance mutations with evidence for multiple origins and expansion of high level sulfadoxine-pyrimethamine resistance mutations in northwest Tanzania.

Background: Emergence of artemisinin partial resistance (ART-R) in Plasmodium falciparum is a growing threat to the efficacy of artemisinin combination therapies (ACT) and the efforts for malaria elimination. The emergence of Plasmodium falciparum Kelch13 (K13) R561H in Rwanda raised concern about the impact in neighboring Tanzania. In addition, regional concern over resistance affecting sulfadoxine-pyrimethamine (SP), which is used for chemoprevention strategies, is high. Methods: To enhance longitudinal monitoring, the Molecular Surveillance of Malaria in Tanzania (MSMT) project was launched in 2020 with the goal of assessing and mapping antimalarial resistance. Community and clinic samples were assessed for resistance polymorphisms using a molecular inversion probe platform. Findings: Genotyping of 6,278 samples collected countrywide in 2021 revealed a focus of K13 561H mutants in northwestern Tanzania (Kagera) with prevalence of 7.7% (50/649). A small number of 561H mutants (about 1%) were found as far as 800 km away in Tabora, Manyara, and Njombe. Genomic analysis suggests some of these parasites are highly related to isolates collected in Rwanda in 2015, supporting regional spread of 561H. However, a novel haplotype was also observed, likely indicating a second origin in the region. Other validated resistance polymorphisms (622I and 675V) were also identified. A focus of high sulfadoxine-pyrimethamine drug resistance was also identified in Kagera with a prevalence of dihydrofolate reductase 164L of 15% (80/526). Interpretation: These findings demonstrate the K13 561H mutation is entrenched in the region and that multiple origins of ART-R, similar as to what was seen in Southeast Asia, have occurred. Mutations associated with high levels of SP resistance are increasing. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region. Funding: This study was funded by the Bill and Melinda Gates Foundation and the National Institutes of Health.

Evolutionary dynamics of microRNA target sites across vertebrate evolution.

MicroRNAs (miRNAs) control the abundance of the majority of the vertebrate transcriptome. The recognition sequences, or target sites, for bilaterian miRNAs are found predominantly in the 3' untranslated regions (3'UTRs) of mRNAs, and are amongst the most highly conserved motifs within 3'UTRs. However, little is known regarding the evolutionary pressures that lead to loss and gain of such target sites. Here, we quantify the selective pressures that act upon miRNA target sites. Notably, selective pressure extends beyond deeply conserved binding sites to those that have undergone recent substitutions. Our approach reveals that even amongst ancient animal miRNAs, which exert the strongest selective pressures on 3'UTR sequences, there are striking differences in patterns of target site evolution between miRNAs. Considering only ancient animal miRNAs, we find three distinct miRNA groups, each exhibiting characteristic rates of target site gain and loss during mammalian evolution. The first group both loses and gains sites rarely. The second group shows selection only against site loss, with site gains occurring at a neutral rate, whereas the third loses and gains sites at neutral or above expected rates. Furthermore, mutations that alter the strength of existing target sites are disfavored. Applying our approach to individual transcripts reveals variation in the distribution of selective pressure across the transcriptome and between miRNAs, ranging from strong selection acting on a small subset of targets of some miRNAs, to weak selection on many targets for other miRNAs. miR-20 and miR-30, and many other miRNAs, exhibit broad, deeply conserved targeting, while several other comparably ancient miRNAs show a lack of selective constraint, and a small number, including mir-146, exhibit evidence of rapidly evolving target sites. Our approach adds valuable perspective on the evolution of miRNAs and their targets, and can also be applied to characterize other 3'UTR regulatory motifs.

A widespread sequence-specific mRNA decay pathway mediated by hnRNPs A1 and A2/B1.

3'-untranslated regions (UTRs) specify post-transcriptional fates of mammalian messenger RNAs (mRNAs), yet knowledge of the underlying sequences and mechanisms is largely incomplete. Here, we identify two related novel 3' UTR motifs in mammals that specify transcript degradation. These motifs are interchangeable and active only within 3' UTRs, where they are often preferentially conserved; furthermore, they are found in hundreds of transcripts, many encoding regulatory proteins. We found that degradation occurs via mRNA deadenylation, mediated by the CCR4-NOT complex. We purified trans factors that recognize the motifs and identified heterogeneous nuclear ribonucleoproteins (hnRNPs) A1 and A2/B1, which are required for transcript degradation, acting in a previously unknown manner. We used RNA sequencing (RNA-seq) to confirm hnRNP A1 and A2/B1 motif-dependent roles genome-wide, profiling cells depleted of these factors singly and in combination. Interestingly, the motifs are most active within the distal portion of 3' UTRs, suggesting that their role in gene regulation can be modulated by alternative processing, resulting in shorter 3' UTRs.

Inferring the evolutionary history of primate microRNA binding sites: overcoming motif counting biases.

The first microRNAs (miRNAs) were identified as essential, conserved regulators of gene expression, targeting the same genes across nearly all bilaterians. However, there are also prominent examples of conserved miRNAs whose functions appear to have shifted dramatically, sometimes over very brief periods of evolutionary time. To determine whether the functions of conserved miRNAs are stable or dynamic over evolutionary time scales, we have here defined the neutral turnover rates of short sequence motifs in predicted primate 3'-UTRs. We find that commonly used approaches to quantify motif turnover rates, which use a presence/absence scoring in extant lineages to infer ancestral states, are inherently biased to infer the accumulation of new motifs, leading to the false inference of continually increasing regulatory complexity over time. Using a maximum likelihood approach to reconstruct individual ancestral nucleotides, we observe that binding sites of conserved miRNAs in fact have roughly equal numbers of gain and loss events relative to ancestral states and turnover extremely slowly relative to nearly identical permutations of the same motif. Contrary to case studies showing examples of functional turnover, our systematic study of miRNA binding sites suggests that in primates, the regulatory roles of conserved miRNAs are strongly conserved. Our revised methodology may be used to quantify the mechanism by which regulatory networks evolve.

Recurrent and recent selective sweeps in the piRNA pathway.

Uncontrolled transposable element (TE) insertions and excisions can cause chromosome breaks and mutations with dramatic deleterious effects. The PIWI interacting RNA (piRNA) pathway functions as an adaptive TE silencing system during germline development. Several essential piRNA pathway proteins appear to be rapidly evolving, suggesting that TEs and the silencing machinery may be engaged in a classical "evolutionary arms race." Using a variety of molecular evolutionary and population genetic approaches, we find that the piRNA pathway genes rhino, krimper, and aubergine show patterns suggestive of extensive recurrent positive selection across Drosophila species. We speculate that selection on these proteins reflects crucial roles in silencing unfamiliar elements during vertical and horizontal transmission of TEs into naïve populations and species, respectively.

Recent progress in polymorphism-based population genetic inference.

The recent availability of whole-genome sequencing data affords tremendous power for statistical inference. With this, there has been great interest in the development of polymorphism-based approaches for the estimation of population genetic parameters. These approaches seek to estimate, for example, recently fixed or sweeping beneficial mutations, the rate of recurrent positive selection, the distribution of selection coefficients, and the demographic history of the population. Yet despite estimating similar parameters using similar data sets, results between methodologies are far from consistent. We here summarize the current state of the field, compare existing approaches, and attempt to reconcile emerging discrepancies. We also discuss the biases in selection estimators introduced by ignoring the demographic history of the population, discuss the biases in demographic estimators introduced by assuming neutrality, and highlight the important challenge to the field of achieving a true joint estimation procedure to circumvent these confounding effects.

MicroRNA-9: functional evolution of a conserved small regulatory RNA.

The functional significance of microRNA-9 (miR-9) during evolution is evidenced by its conservation at the nucleotide level from flies to humans but not its diverse expression patterns. Recent studies in several model systems reveal that miR-9 can regulate neurogenesis through its actions in neural or non-neural cell lineages. In vertebrates, miR-9 exerts diverse cell-autonomous effects on the proliferation, migration, and differentiation of neural progenitor cells by modulating different mRNA targets. In some developmental contexts, miR-9 suppresses apoptosis and is misregulated in several types of cancer cells, influencing proliferation or metastasis formation. Moreover, downregulation of miR-9 in postmitotic neurons is also implicated in some neurodegenerative diseases. Thus, miR-9 is emerging as an important regulator in development and disease through its ability to modulate different targets in a manner dependent on the developmental stage and the cellular context.