Population Genomic Evidence of Adaptive Response during the Invasion History of Plasmodium falciparum in the Americas.

Plasmodium falciparum, the most virulent agent of human malaria, spread from Africa to all continents following the out-of-Africa human migrations. During the transatlantic slave trade between the 16th and 19th centuries, it was introduced twice independently to the Americas where it adapted to new environmental conditions (new human populations and mosquito species). Here, we analyzed the genome-wide polymorphisms of 2,635 isolates across the current P. falciparum distribution range in Africa, Asia, Oceania, and the Americas to investigate its genetic structure, invasion history, and selective pressures associated with its adaptation to the American environment. We confirmed that American populations originated from Africa with at least two independent introductions that led to two genetically distinct clusters, one in the North (Haiti and Colombia) and one in the South (French Guiana and Brazil), and an admixed Peruvian group. Genome scans revealed recent and more ancient signals of positive selection in the American populations. Particularly, we detected positive selection signals in genes involved in interactions with hosts (human and mosquito) cells and in genes involved in resistance to malaria drugs in both clusters. Analyses suggested that for five genes, adaptive introgression between clusters or selection on standing variation was at the origin of this repeated evolution. This study provides new genetic evidence on P. falciparum colonization history and on its local adaptation in the Americas.

The Arabidopsis PHD-finger protein EDM2 has multiple roles in balancing NLR immune receptor gene expression.

Plant NLR-type receptors serve as sensitive triggers of host immunity. Their expression has to be well-balanced, due to their interference with various cellular processes and dose-dependency of their defense-inducing activity. A genetic "arms race" with fast-evolving pathogenic microbes requires plants to constantly innovate their NLR repertoires. We previously showed that insertion of the COPIA-R7 retrotransposon into RPP7 co-opted the epigenetic transposon silencing signal H3K9me2 to a new function promoting expression of this Arabidopsis thaliana NLR gene. Recruitment of the histone binding protein EDM2 to COPIA-R7-associated H3K9me2 is required for optimal expression of RPP7. By profiling of genome-wide effects of EDM2, we now uncovered additional examples illustrating effects of transposons on NLR gene expression, strongly suggesting that these mobile elements can play critical roles in the rapid evolution of plant NLR genes by providing the "raw material" for gene expression mechanisms. We further found EDM2 to have a global role in NLR expression control. Besides serving as a positive regulator of RPP7 and a small number of other NLR genes, EDM2 acts as a suppressor of a multitude of additional NLR genes. We speculate that the dual functionality of EDM2 in NLR expression control arose from the need to compensate for fitness penalties caused by high expression of some NLR genes by suppression of others. Moreover, we are providing new insights into functional relationships of EDM2 with its interaction partner, the RNA binding protein EDM3/AIPP1, and its target gene IBM1, encoding an H3K9-demethylase.

Genome-Wide Reinforcement of DNA Methylation Occurs during Somatic Embryogenesis in Soybean.

Somatic embryogenesis is an important tissue culture technique that sometimes leads to phenotypic variation via genetic and/or epigenetic changes. To understand the genomic and epigenomic impacts of somatic embryogenesis, we characterized soybean (Glycine max) epigenomes sampled from embryos at 10 different stages ranging from 6 weeks to 13 years of continuous culture. We identified genome-wide increases in DNA methylation from cultured samples, especially at CHH sites. The hypermethylation almost exclusively occurred in regions previously possessing non-CG methylation and was accompanied by increases in the expression of genes encoding the RNA-directed DNA methylation (RdDM) machinery. The epigenomic changes were similar between somatic and zygotic embryogenesis. Following the initial global wave of hypermethylation, rare decay events of maintenance methylation were observed, and the extent of the decay increased with time in culture. These losses in DNA methylation were accompanied by downregulation of genes encoding the RdDM machinery and transcriptome reprogramming reminiscent of transcriptomes during late-stage seed development. These results reveal a process for reinforcing already silenced regions to maintain genome integrity during somatic embryogenesis over the short term, which eventually decays at certain loci over longer time scales.

Evolutionary history of Plasmodium vivax and Plasmodium simium in the Americas.

Malaria is a vector-borne disease caused by protozoan parasites of the genus Plasmodium. Plasmodium vivax is the most prevalent human-infecting species in the Americas. However, the origins of this parasite in this continent are still debated. Similarly, it is now accepted that the existence of Plasmodium simium is explained by a P. vivax transfer from humans to monkey in America. However, many uncertainties still exist concerning the origin of the transfer and whether several transfers occurred. In this review, the most recent studies that addressed these questions using genetic and genomic approaches are presented.

Biosynthesis and defensive function of Nδ-acetylornithine, a jasmonate-induced Arabidopsis metabolite.

Since research on plant interactions with herbivores and pathogens is often constrained by the analysis of already known compounds, there is a need to identify new defense-related plant metabolites. The uncommon nonprotein amino acid N(δ)-acetylornithine was discovered in a targeted search for Arabidopsis thaliana metabolites that are strongly induced by the phytohormone methyl jasmonate (MeJA). Stable isotope labeling experiments show that, after MeJA elicitation, Arg, Pro, and Glu are converted to Orn, which is acetylated by NATA1 to produce N(δ)-acetylornithine. MeJA-induced N(δ)-acetylornithine accumulation occurs in all tested Arabidopsis accessions, other Arabidopsis species, Capsella rubella, and Boechera stricta, but not in less closely related Brassicaceae. Both insect feeding and Pseudomonas syringae infection increase NATA1 expression and N(δ)-acetylornithine accumulation. NATA1 transient expression in Nicotiana tabacum and the addition of N(δ)-acetylornithine to an artificial diet both decrease Myzus persicae (green peach aphid) reproduction, suggesting a direct toxic or deterrent effect. However, since broad metabolic changes that are induced by MeJA in wild-type Arabidopsis are attenuated in a nata1 mutant strain, there may also be indirect effects on herbivores and pathogens. In the case of P. syringae, growth on a nata1 mutant is reduced compared with wild-type Arabidopsis, but growth in vitro is unaffected by N(δ)-acetylornithine addition.

Host preference patterns in domestic and wild settings: Insights into Anopheles feeding behavior.

The adaptation of Anopheles malaria vectors to domestic settings is directly linked to their ability to feed on humans. The strength of this species-habitat association is unequal across the species within the genus, with the major vectors being particularly dependent on humans. However, our understanding of how blood-feeding behavior interacts with and adapts to environmental settings, including the presence of humans, remains limited. Using a field-based approach, we first investigated Anopheles community structure and feeding behavior patterns in domestic and sylvatic settings in La Lopé National Park in Gabon, Central Africa. We characterized the preference indices using a dual-host choice sampling approach across mosquito species, habitats, and seasons. We then quantified the plastic biting behavior of mosquito species in each habitat. We collected individuals from 16 Anopheles species that exhibited significant differences in species composition and abundance between sylvatic and domestic settings. The host-seeking behavior also varied among the seven most abundant species. The general attractiveness to each host, human or animal, remained relatively constant for each species, but with significant variations between habitats across species. These variations, to more generalist and to more anthropophilic behavior, were related to seasonal changes and distance from the village, respectively. Finally, we pointed out that the host choice of major malaria vectors changed in the absence of humans, revealing a plastic feeding behavior of these species. This study highlights the effect of humans on Anopheles distribution and feeding evolution. The characterization of feeding behavior in wild and domestic settings provides opportunities to better understand the interplay between genetic determinants of host preference and ecological factors. Our findings suggest that protected areas may offer alternative thriving conditions to major malaria vectors.

Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon.

Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human. Significance: Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression.

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Limited association between Wolbachia and Plasmodium falciparum infections in natural populations of the major malaria mosquito Anopheles moucheti.

Since the discovery of natural malaria vector populations infected by the endosymbiont bacterium Wolbachia, a renewed interest has arisen for using this bacterium as an alternative for malaria control. Among naturally infected mosquitoes, Anopheles moucheti, a major malaria mosquito in Central Africa, exhibits one of the highest prevalences of Wolbachia infection. To better understand whether this maternally inherited bacterium could be used for malaria control, we investigated Wolbachia influence in An. moucheti populations naturally infected by the malaria parasite Plasmodium falciparum. To this end, we collected mosquitoes in a village from Cameroon, Central Africa, where this mosquito is the main malaria vector. We found that the prevalence of Wolbachia bacterium was almost fixed in the studied mosquito population, and was higher than previously recorded. We also quantified Wolbachia in whole mosquitoes and dissected abdomens, confirming that the bacterium is also elsewhere than in the abdomen, but at lower density. Finally, we analyzed the association of Wolbachia presence and density on P. falciparum infection. Wolbachia density was slightly higher in mosquitoes infected with the malaria parasite than in uninfected mosquitoes. However, we observed no correlation between the P. falciparum and Wolbachia densities. In conclusion, our study indicates that naturally occurring Wolbachia infection is not associated to P. falciparum development within An. moucheti mosquitoes.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide.

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Population genomic evidence of Plasmodium vivax Southeast Asian origin.

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

Soybean aphids adapted to host-plant resistance by down regulating putative effectors and up regulating transposable elements.

In agricultural systems, crops equipped with host-plant resistance (HPR) have enhanced protection against pests, and are used as a safe and sustainable tool in pest management. In soybean, HPR can control the soybean aphid (Aphis glycines), but certain aphid populations have overcome this resistance (i.e., virulence). The molecular mechanisms underlying aphid virulence to HPR are unknown, but likely involve effector proteins that are secreted by aphids to modulate plant defenses. Another mechanism to facilitate adaptation is through the activity of transposable elements, which can become activated by stress. In this study, we performed RNA sequencing of virulent and avirulent soybean aphids fed susceptible or resistant (Rag1 + Rag2) soybean. Our goal was to better understand the molecular mechanisms underlying soybean aphid virulence. Our data showed that virulent aphids mostly down regulate putative effector genes relative to avirulent aphids, especially when aphids were fed susceptible soybean. Decreased expression of effectors may help evade HPR plant defenses. Virulent aphids also transcriptionally up regulate a diverse set of transposable elements and nearby genes, which is consistent with stress adaptation. Our work demonstrates two mechanisms of pest adaptation to resistance, and identifies effector gene targets for future functional testing.

A new role for histone demethylases in the maintenance of plant genome integrity.

Histone modifications deposited by the Polycomb repressive complex 2 (PRC2) play a critical role in the control of growth, development, and adaptation to environmental fluctuations of most multicellular eukaryotes. The catalytic activity of PRC2 is counteracted by Jumonji-type (JMJ) histone demethylases, which shapes the genomic distribution of H3K27me3. Here, we show that two JMJ histone demethylases in Arabidopsis, EARLY FLOWERING 6 (ELF6) and RELATIVE OF EARLY FLOWERING 6 (REF6), play distinct roles in H3K27me3 and H3K27me1 homeostasis. We show that failure to reset these chromatin marks during sexual reproduction results in the transgenerational inheritance of histone marks, which cause a loss of DNA methylation at heterochromatic loci and transposon activation. Thus, Jumonji-type histone demethylases play a dual role in plants by helping to maintain transcriptional states through development and safeguard genome integrity during sexual reproduction.

EpiTEome: Simultaneous detection of transposable element insertion sites and their DNA methylation levels.

The genome-wide investigation of DNA methylation levels has been limited to reference transposable element positions. The methylation analysis of non-reference and mobile transposable elements has only recently been performed, but required both genome resequencing and MethylC-seq datasets. We have created epiTEome, a program that detects both new transposable element insertion sites and their methylation states from a single MethylC-seq dataset. EpiTEome outperforms other split-read insertion site detection programs, even while functioning on bisulfite-converted reads. EpiTEome characterizes the previously discarded fraction of DNA methylation at sites of new insertions, enabling future investigation into the epigenetic regulation of non-reference and transposed elements.

Full-length autonomous transposable elements are preferentially targeted by expression-dependent forms of RNA-directed DNA methylation.

BACKGROUND: Chromatin modifications such as DNA methylation are targeted to transposable elements by small RNAs in a process termed RNA-directed DNA methylation (RdDM). In plants, canonical RdDM functions through RNA polymerase IV to reinforce pre-existing transposable element silencing. Recent investigations have identified a "non-canonical" form of RdDM dependent on RNA polymerase II expression to initiate and re-establish silencing of active transposable elements. This expression-dependent RdDM mechanism functions through RNAi degradation of transposable element mRNAs into small RNAs guided by the RNA-dependent RNA polymerase 6 (RDR6) protein and is therefore referred to as RDR6-RdDM. RESULTS: We performed whole-genome MethylC-seq in 20 mutants that distinguish RdDM mechanisms when transposable elements are either transcriptionally silent or active. We identified a new mechanism of expression-dependent RdDM, which functions through DICER-LIKE3 (DCL3) but bypasses the requirement of both RNA polymerase IV and RDR6 (termed DCL3-RdDM). We found that RNA polymerase II expression-dependent forms of RdDM function on over 20 % of transcribed transposable elements, including the majority of full-length elements with all of the domains required for autonomous transposition. Lastly, we find that RDR6-RdDM preferentially targets long transposable elements due to the specificity of primary small RNAs to cleave full-length mRNAs. CONCLUSIONS: Expression-dependent forms of RdDM function to critically target DNA methylation to full-length and transcriptionally active transposable elements, suggesting that these pathways are key to suppressing mobilization. This targeting specificity is initiated on the mRNA cleavage-level, yet manifested as chromatin-level silencing that in plants is epigenetically inherited from generation to generation.

Small-scale gene duplications played a major role in the recent evolution of wheat chromosome 3B.

BACKGROUND: Bread wheat is not only an important crop, but its large (17 Gb), highly repetitive, and hexaploid genome makes it a good model to study the organization and evolution of complex genomes. Recently, we produced a high quality reference sequence of wheat chromosome 3B (774 Mb), which provides an excellent opportunity to study the evolutionary dynamics of a large and polyploid genome, specifically the impact of single gene duplications. RESULTS: We find that 27 % of the 3B predicted genes are non-syntenic with the orthologous chromosomes of Brachypodium distachyon, Oryza sativa, and Sorghum bicolor, whereas, by applying the same criteria, non-syntenic genes represent on average only 10 % of the predicted genes in these three model grasses. These non-syntenic genes on 3B have high sequence similarity to at least one other gene in the wheat genome, indicating that hexaploid wheat has undergone massive small-scale interchromosomal gene duplications compared to other grasses. Insertions of non-syntenic genes occurred at a similar rate along the chromosome, but these genes tend to be retained at a higher frequency in the distal, recombinogenic regions. The ratio of non-synonymous to synonymous substitution rates showed a more relaxed selection pressure for non-syntenic genes compared to syntenic genes, and gene ontology analysis indicated that non-syntenic genes may be enriched in functions involved in disease resistance. CONCLUSION: Our results highlight the major impact of single gene duplications on the wheat gene complement and confirm the accelerated evolution of the Triticeae lineage among grasses.

Organization and evolution of transposable elements along the bread wheat chromosome 3B.

BACKGROUND: The 17 Gb bread wheat genome has massively expanded through the proliferation of transposable elements (TEs) and two recent rounds of polyploidization. The assembly of a 774 Mb reference sequence of wheat chromosome 3B provided us with the opportunity to explore the impact of TEs on the complex wheat genome structure and evolution at a resolution and scale not reached so far. RESULTS: We develop an automated workflow, CLARI-TE, for TE modeling in complex genomes. We delineate precisely 56,488 intact and 196,391 fragmented TEs along the 3B pseudomolecule, accounting for 85% of the sequence, and reconstruct 30,199 nested insertions. TEs have been mostly silent for the last one million years, and the 3B chromosome has been shaped by a succession of bursts that occurred between 1 to 3 million years ago. Accelerated TE elimination in the high-recombination distal regions is a driving force towards chromosome partitioning. CACTAs overrepresented in the high-recombination distal regions are significantly associated with recently duplicated genes. In addition, we identify 140 CACTA-mediated gene capture events with 17 genes potentially created by exon shuffling and show that 19 captured genes are transcribed and under selection pressure, suggesting the important role of CACTAs in the recent wheat adaptation. CONCLUSION: Accurate TE modeling uncovers the dynamics of TEs in a highly complex and polyploid genome. It provides novel insights into chromosome partitioning and highlights the role of CACTA transposons in the high level of gene duplication in wheat.

Structural and functional partitioning of bread wheat chromosome 3B.

We produced a reference sequence of the 1-gigabase chromosome 3B of hexaploid bread wheat. By sequencing 8452 bacterial artificial chromosomes in pools, we assembled a sequence of 774 megabases carrying 5326 protein-coding genes, 1938 pseudogenes, and 85% of transposable elements. The distribution of structural and functional features along the chromosome revealed partitioning correlated with meiotic recombination. Comparative analyses indicated high wheat-specific inter- and intrachromosomal gene duplication activities that are potential sources of variability for adaption. In addition to providing a better understanding of the organization, function, and evolution of a large and polyploid genome, the availability of a high-quality sequence anchored to genetic maps will accelerate the identification of genes underlying important agronomic traits.