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.

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 the 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. We found that some genes were under selection 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.

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.

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.

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.

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.

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.

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.

We previously found the Arabidopsis thaliana PHD-finger protein EDM2 to serve as a chromatin-associated factor controlling expression of the NLR-type immune receptor gene RPP7. EDM2 binds to the transposon-silencing signal H3K9me2 and affects levels of this epigenetic mark at various loci. By genome-wide profiling of transcript- and H3K9me2-levels we now found EDM2 to have a broader role in controlling NLR gene expression. In order to mitigate fitness costs caused by its promoting effects on RPP7 expression and that of several other NLR genes, EDM2 seems to suppress expression of many additional members of this gene family. This observation is in line with multiple reports demonstrating the need for balanced expression of NLRs, which can substantially reduce overall plant fitness, but need to be present at certain minimal levels to confer sufficient immune protection. Our previous results demonstrated that the influence of EDM2 on RPP7 expression was co-opted to this immune receptor gene by the insertion of an EDM2-controlled transposon. Here, we are providing additional examples for transposon-associated effects on NLR gene expression, suggesting that these mobile elements play an important role for NLR genes by equipping members of this rapidly evolving gene family with regulatory mechanisms needed for balanced expression.

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.

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.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1032-y) contains supplementary material, which is available to authorized users.

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.

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.

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.