Cooperation and cheating orchestrate Vibrio assemblages and polymicrobial synergy in oysters infected with OsHV-1 virus.

Polymicrobial infections threaten the health of humans and animals but remain understudied in natural systems. We recently described the Pacific Oyster Mortality Syndrome (POMS), a polymicrobial disease affecting oyster production worldwide. In the French Atlantic coast, the disease involves coinfection with ostreid herpesvirus 1 (OsHV-1) and virulent Vibrio. However, it is unknown whether consistent Vibrio populations are associated with POMS in different regions, how Vibrio contribute to POMS, and how they interact with OsHV-1 during pathogenesis. By connecting field-based approaches in a Mediterranean ecosystem, laboratory infection assays and functional genomics, we uncovered a web of interdependencies that shape the structure and function of the POMS pathobiota. We show that Vibrio harveyi and Vibrio rotiferianus are predominant in OsHV-1-diseased oysters and that OsHV-1 drives the partition of the Vibrio community observed in the field. However only V. harveyi synergizes with OsHV-1 by promoting mutual growth and accelerating oyster death. V. harveyi shows high-virulence potential and dampens oyster cellular defenses through a type 3 secretion system, making oysters a more favorable niche for microbe colonization. In addition, V. harveyi produces a key siderophore called vibrioferrin. This important resource promotes the growth of V. rotiferianus, which cooccurs with V. harveyi in diseased oysters, and behaves as a cheater by benefiting from V. harveyi metabolite sharing. Our data show that cooperative behaviors contribute to synergy between bacterial and viral coinfecting partners. Additional cheating behaviors further shape the polymicrobial consortium. Controlling cooperative behaviors or countering their effects opens avenues for mitigating polymicrobial diseases.

Transcriptomic responses of Antarctic clam Laternula elliptica to nanoparticles, at single and combined exposures reveal ecologically relevant biomarkers.

In recent years micro- and nanoplastics and metal-oxide nanomaterials have been found in several environmental compartments. The Antarctic soft clam Laternula elliptica is an endemic Antarctic species having a wide distribution in the Southern Ocean. Being a filter-feeder, it could act as suitable bioindicator of pollution from nanoparticles also considering its sensitivity to various sources of stress. The present study aims to assess the impact of polystyrene nanoparticles (PS-NP) and the nanometal titanium-dioxide (n-TiO2) on genome-wide transcript expression of L. elliptica either alone and in combination and at two toxicological relevant concentrations (5 and 50 µg/L) during 96 h exposure. Transcript-target qRT-PCR was performed with the aim to identify suitable biomarkers of exposure and effects. As expected, at the highest concentration tested, the clustering was clearer between control and exposed clams. A total of 221 genes resulted differentially expressed in exposed clams and control ones, and 21 of them had functional annotation such as ribosomal proteins, antioxidant, ion transport (osmoregulation), acid-base balance, immunity, lipid metabolism, cell adhesion, cytoskeleton, apoptosis, chromatin condensation and cell signaling. At functional level, relevant transcripts were shared among some treatments and could be considered as general stress due to nanoparticle exposure. After applying transcript-target approach duplicating the number of clam samples, four ecologically relevant transcripts were revealed as biomarkers for PS-NP, n-TiO2 and their combination at 50 µg/L, that could be used for monitoring clams' health status in different Antarctic localities.

Diverging patterns of introgression from Schistosoma bovis across S. haematobium African lineages.

Hybridization is a fascinating evolutionary phenomenon that raises the question of how species maintain their integrity. Inter-species hybridization occurs between certain Schistosoma species that can cause important public health and veterinary issues. In particular hybrids between Schistosoma haematobium and S. bovis associated with humans and animals respectively are frequently identified in Africa. Recent genomic evidence indicates that some S. haematobium populations show signatures of genomic introgression from S. bovis. Here, we conducted a genomic comparative study and investigated the genomic relationships between S. haematobium, S. bovis and their hybrids using 19 isolates originating from a wide geographical range over Africa, including samples initially classified as S. haematobium (n = 11), S. bovis (n = 6) and S. haematobium x S. bovis hybrids (n = 2). Based on a whole genomic sequencing approach, we developed 56,181 SNPs that allowed a clear differentiation of S. bovis isolates from a genomic cluster including all S. haematobium isolates and a natural S. haematobium-bovis hybrid. All the isolates from the S. haematobium cluster except the isolate from Madagascar harbored signatures of genomic introgression from S. bovis. Isolates from Corsica, Mali and Egypt harbored the S. bovis-like Invadolysin gene, an introgressed tract that has been previously detected in some introgressed S. haematobium populations from Niger. Together our results highlight the fact that introgression from S. bovis is widespread across S. haematobium and that the observed introgression is unidirectional.

Species-specific mechanisms of cytotoxicity toward immune cells determine the successful outcome of Vibrio infections.

Vibrio species cause infectious diseases in humans and animals, but they can also live as commensals within their host tissues. How Vibrio subverts the host defenses to mount a successful infection remains poorly understood, and this knowledge is critical for predicting and managing disease. Here, we have investigated the cellular and molecular mechanisms underpinning infection and colonization of 2 virulent Vibrio species in an ecologically relevant host model, oyster, to study interactions with marine Vibrio species. All Vibrio strains were recognized by the immune system, but only nonvirulent strains were controlled. We showed that virulent strains were cytotoxic to hemocytes, oyster immune cells. By analyzing host and bacterial transcriptional responses to infection, together with Vibrio gene knock-outs, we discovered that Vibrio crassostreae and Vibrio tasmaniensis use distinct mechanisms to cause hemocyte lysis. Whereas V. crassostreae cytotoxicity is dependent on a direct contact with hemocytes and requires an ancestral gene encoding a protein of unknown function, r5.7, V. tasmaniensis cytotoxicity is dependent on phagocytosis and requires intracellular secretion of T6SS effectors. We conclude that proliferation of commensal vibrios is controlled by the host immune system, preventing systemic infections in oysters, whereas the successful infection of virulent strains relies on Vibrio species-specific molecular determinants that converge to compromise host immune cell function, allowing evasion of the host immune system.

Sympatric versus allopatric evolutionary contexts shape differential immune response in Biomphalaria / Schistosoma interaction.

Selective pressures between hosts and their parasites can result in reciprocal evolution or adaptation of specific life history traits. Local adaptation of resident hosts and parasites should lead to increase parasite infectivity/virulence (higher compatibility) when infecting hosts from the same location (in sympatry) than from a foreign location (in allopatry). Analysis of geographic variations in compatibility phenotypes is the most common proxy used to infer local adaptation. However, in some cases, allopatric host-parasite systems demonstrate similar or greater compatibility than in sympatry. In such cases, the potential for local adaptation remains unclear. Here, we study the interaction between Schistosoma and its vector snail Biomphalaria in which such discrepancy in local versus foreign compatibility phenotype has been reported. Herein, we aim at bridging this gap of knowledge by comparing life history traits (immune cellular response, host mortality, and parasite growth) and molecular responses in highly compatible sympatric and allopatric Schistosoma/Biomphalaria interactions originating from different geographic localities (Brazil, Venezuela and Burundi). We found that despite displaying similar prevalence phenotypes, sympatric schistosomes triggered a rapid immune suppression (dual-RNAseq analyses) in the snails within 24h post infection, whereas infection by allopatric schistosomes (regardless of the species) was associated with immune cell proliferation and triggered a non-specific generalized immune response after 96h. We observed that, sympatric schistosomes grow more rapidly. Finally, we identify miRNAs differentially expressed by Schistosoma mansoni that target host immune genes and could be responsible for hijacking the host immune response during the sympatric interaction. We show that despite having similar prevalence phenotypes, sympatric and allopatric snail-Schistosoma interactions displayed strong differences in their immunobiological molecular dialogue. Understanding the mechanisms allowing parasites to adapt rapidly and efficiently to new hosts is critical to control disease emergence and risks of Schistosomiasis outbreaks.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar.

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Histone methylation changes are required for life cycle progression in the human parasite Schistosoma mansoni.

Epigenetic mechanisms and chromatin structure play an important role in development. Their impact is therefore expected to be strong in parasites with complex life cycles and multiple, strikingly different, developmental stages, i.e. developmental plasticity. Some studies have already described how the chromatin structure, through histone modifications, varies from a developmental stage to another in a few unicellular parasites. While H3K4me3 profiles remain relatively constant, H3K27 trimethylation and bivalent methylation show strong variation. Inhibitors (A366 and GSK343) of H3K27 histone methyltransferase activity in S. mansoni efficiently blocked miracidium to sporocyst transition indicating that H3K27 trimethylation is required for life cycle progression. As S. mansoni is a multicellular parasite that significantly affects both the health and economy of endemic areas, a better understanding of fluke developmental processes within the definitive host will likely highlight novel disease control strategies. Towards this goal, we also studied H4K20me1 in female cercariae and adults. In particular, we found that bivalent trimethylation of H3K4 and H3K27 at the transcription start site of genes is a landmark of the cercarial stage. In cercariae, H3K27me3 presence and strong enrichment in H4K20me1 over long regions (10-100 kb) is associated with development related genes. Here, we provide a broad overview of the chromatin structure of a metazoan parasite throughout its most important lifecycle stages. The five developmental stages studied here present distinct chromatin structures, indicating that histone methylation plays an important role during development. Hence, components of the histone methylation (and demethylation) machinery may provide suitable Schistosomiasis control targets.

Developmental Methylome of the Medicinal Plant Catharanthus roseus Unravels the Tissue-Specific Control of the Monoterpene Indole Alkaloid Pathway by DNA Methylation.

Catharanthus roseus produces a wide spectrum of monoterpene indole alkaloids (MIAs). MIA biosynthesis requires a tightly coordinated pathway involving more than 30 enzymatic steps that are spatio-temporally and environmentally regulated so that some MIAs specifically accumulate in restricted plant parts. The first regulatory layer involves a complex network of transcription factors from the basic Helix Loop Helix (bHLH) or AP2 families. In the present manuscript, we investigated whether an additional epigenetic layer could control the organ-, developmental- and environmental-specificity of MIA accumulation. We used Whole-Genome Bisulfite Sequencing (WGBS) together with RNA-seq to identify differentially methylated and expressed genes among nine samples reflecting different plant organs and experimental conditions. Tissue specific gene expression was associated with specific methylation signatures depending on cytosine contexts and gene parts. Some genes encoding key enzymatic steps from the MIA pathway were found to be simultaneously differentially expressed and methylated in agreement with the corresponding MIA accumulation. In addition, we found that transcription factors were strikingly concerned by DNA methylation variations. Altogether, our integrative analysis supports an epigenetic regulation of specialized metabolisms in plants and more likely targeting transcription factors which in turn may control the expression of enzyme-encoding genes.

Notos - a galaxy tool to analyze CpN observed expected ratios for inferring DNA methylation types.

BACKGROUND: DNA methylation patterns store epigenetic information in the vast majority of eukaryotic species. The relatively high costs and technical challenges associated with the detection of DNA methylation however have created a bias in the number of methylation studies towards model organisms. Consequently, it remains challenging to infer kingdom-wide general rules about the functions and evolutionary conservation of DNA methylation. Methylated cytosine is often found in specific CpN dinucleotides, and the frequency distributions of, for instance, CpG observed/expected (CpG o/e) ratios have been used to infer DNA methylation types based on higher mutability of methylated CpG. RESULTS: Predominantly model-based approaches essentially founded on mixtures of Gaussian distributions are currently used to investigate questions related to the number and position of modes of CpG o/e ratios. These approaches require the selection of an appropriate criterion for determining the best model and will fail if empirical distributions are complex or even merely moderately skewed. We use a kernel density estimation (KDE) based technique for robust and precise characterization of complex CpN o/e distributions without a priori assumptions about the underlying distributions. CONCLUSIONS: We show that KDE delivers robust descriptions of CpN o/e distributions. For straightforward processing, we have developed a Galaxy tool, called Notos and available at the ToolShed, that calculates these ratios of input FASTA files and fits a density to their empirical distribution. Based on the estimated density the number and shape of modes of the distribution is determined, providing a rational for the prediction of the number and the types of different methylation classes. Notos is written in R and Perl.

Frequency and mitotic heritability of epimutations in Schistosoma mansoni.

Schistosoma mansoni is a parasitic platyhelminth responsible for intestinal bilharzia. It has a complex life cycle, infecting a freshwater snail of the Biomphalaria genus, and then a mammalian host. Schistosoma mansoni adapts rapidly to new (allopatric) strains of its intermediate host. To study the importance of epimutations in this process, we infected sympatric and allopatric mollusc strains with parasite clones. ChIP-Seq was carried out on four histone modifications (H3K4me3, H3K27me3, H3K27ac and H4K20me1) in parallel with genomewide DNA resequencing (i) on parasite larvae shed by the infected snails and (ii) on adult worms that had developed from the larvae. No change in single nucleotide polymorphisms and no mobilization of transposable elements were observed, but 58-105 copy number variations (CNVs) within the parasite clones in different molluscs were detected. We also observed that the allopatric environment induces three types of chromatin structure changes: (i) host-induced changes on larvae epigenomes in 51 regions of the genome that are independent of the parasites' genetic background, (ii) spontaneous changes (not related to experimental condition or genotype of the parasite) at 64 locations and (iii) 64 chromatin structure differences dependent on the parasite genotype. Up to 45% of the spontaneous, but none of the host-induced chromatin structure changes were transmitted to adults. In our model, the environment induces epigenetic changes at specific loci but only spontaneous epimutations are mitotically heritable and have therefore the potential to contribute to transgenerational inheritance. We also show that CNVs are the only source of genetic variation and occur at the same order of magnitude as epimutations.

Genome structure and metabolic features in the red seaweed Chondrus crispus shed light on evolution of the Archaeplastida.

Red seaweeds are key components of coastal ecosystems and are economically important as food and as a source of gelling agents, but their genes and genomes have received little attention. Here we report the sequencing of the 105-Mbp genome of the florideophyte Chondrus crispus (Irish moss) and the annotation of the 9,606 genes. The genome features an unusual structure characterized by gene-dense regions surrounded by repeat-rich regions dominated by transposable elements. Despite its fairly large size, this genome shows features typical of compact genomes, e.g., on average only 0.3 introns per gene, short introns, low median distance between genes, small gene families, and no indication of large-scale genome duplication. The genome also gives insights into the metabolism of marine red algae and adaptations to the marine environment, including genes related to halogen metabolism, oxylipins, and multicellularity (microRNA processing and transcription factors). Particularly interesting are features related to carbohydrate metabolism, which include a minimalistic gene set for starch biosynthesis, the presence of cellulose synthases acquired before the primary endosymbiosis showing the polyphyly of cellulose synthesis in Archaeplastida, and cellulases absent in terrestrial plants as well as the occurrence of a mannosylglycerate synthase potentially originating from a marine bacterium. To explain the observations on genome structure and gene content, we propose an evolutionary scenario involving an ancestral red alga that was driven by early ecological forces to lose genes, introns, and intergenetic DNA; this loss was followed by an expansion of genome size as a consequence of activity of transposable elements.

Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis.

When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.

Combining a transcriptomic approach and a targeted metabolomics approach for deciphering the molecular bases of compatibility phenotype in the snail Biomphalaria glabrata toward Schistosoma mansoni.

Biomphalaria glabrata is a freshwater snail and the obligatory intermediate host of Schistosoma mansoni parasite, the etiologic agent of intestinal Schistosomiasis, in South America and Caribbean. Interestingly in such host-parasite interactions, compatibility varies between populations, strains or individuals. This observed compatibility polymorphism is based on a complex molecular-matching-phenotype, the molecular bases of which have been investigated in numerous studies, notably by comparing between different strains or geographical isolates or clonal selected snail lines. Herein we propose to decipher the constitutive molecular support of this interaction in selected non-clonal resistant and susceptible snail strain originating from the same natural population from Brazil and thus having the same genetic background. Thanks to a global RNAseq transcriptomic approach on whole snail, we identified a total of 328 differentially expressed genes between resistant and susceptible phenotypes among which 129 were up-regulated and 199 down-regulated. Metabolomic studies were used to corroborate the RNAseq results. The activation of immune genes and specific metabolic pathways in resistant snails might provide them with the capacity to better respond to parasite infection.

Double-stranded RNA binding proteins DRB2 and DRB4 have an antagonistic impact on polymerase IV-dependent siRNA levels in Arabidopsis.

Biogenesis of the vast majority of plant siRNAs depends on the activity of the plant-specific RNA polymerase IV (PolIV) enzyme. As part of the RNA-dependent DNA methylation (RdDM) process, PolIV-dependent siRNAs (p4-siRNAs) are loaded onto an ARGONAUTE4-containing complex and guide de novo DNA methyltransferases to target loci. Here we show that the double-stranded RNA binding proteins DRB2 and DRB4 are required for proper accumulation of p4-siRNAs. In flowers, loss of DRB2 results in increased accumulation of p4-siRNAs but not ta-siRNAs, inverted repeat (IR)-derived siRNAs, or miRNA. Loss of DRB2 does not impair uniparental expression of p4-dependent siRNAs in developing endosperm, indicating that p4-siRNA increased accumulation is not the result of the activation of the polIV pathway in the male gametophyte. In contrast to drb2, drb4 mutants exhibit reduced p4-siRNA levels, but the extent of this reduction is variable, according to the nature and size of the p4-siRNAs. Loss of DRB4 also leads to a spectacular increase of p4-independent IR-derived 24-nt siRNAs, suggesting a reallocation of factors from p4-dependent to p4-independent siRNA pathways in drb4. Opposite effects of drb2 and drb4 mutations on the accumulation of p4-siRNAs were also observed in vegetative tissues. Moreover, transgenic plants overexpressing DRB2 mimicked drb4 mutants at the morphological and molecular levels, confirming the antagonistic roles of DRB2 and DRB4.

Targeted insertion and reporter transgene activity at a gene safe harbor of the human blood fluke, Schistosoma mansoni.

The identification and characterization of genomic safe harbor sites (GSHs) can facilitate consistent transgene activity with minimal disruption to the host cell genome. We combined computational genome annotation and chromatin structure analysis to predict the location of four GSHs in the human blood fluke, Schistosoma mansoni, a major infectious pathogen of the tropics. A transgene was introduced via CRISPR-Cas-assisted homology-directed repair into one of the GSHs in the egg of the parasite. Gene editing efficiencies of 24% and transgene-encoded fluorescence of 75% of gene-edited schistosome eggs were observed. The approach advances functional genomics for schistosomes by providing a tractable path for generating transgenics using homology-directed, repair-catalyzed transgene insertion. We also suggest that this work will serve as a roadmap for the development of similar approaches in helminths more broadly.

Transpositional landscape of the rice genome revealed by paired-end mapping of high-throughput re-sequencing data.

Transposable elements (TEs) are mobile entities that densely populate most eukaryotic genomes and contribute to both their structural and functional dynamics. However, most TE-related sequences in both plant and animal genomes correspond to inactive, degenerated elements, due to the combined effect of silencing pathways and elimination through deletions. One of the major difficulties in fully characterizing the molecular basis of genetic diversity of a given species lies in establishing its genome-wide transpositional activity. Here, we provide an extensive survey of the transpositional landscape of a plant genome using a deep sequencing strategy. This was achieved through paired-end mapping of a fourfold coverage of the genome of rice mutant line derived from an in vitro callus culture using Illumina technology. Our study shows that at least 13 TE families are active in this genotype, causing 34 new insertions. This next-generation sequencing-based strategy provides new opportunities to quantify the impact of TEs on the genome dynamics of the species.

Exceptional diversity, non-random distribution, and rapid evolution of retroelements in the B73 maize genome.

Recent comprehensive sequence analysis of the maize genome now permits detailed discovery and description of all transposable elements (TEs) in this complex nuclear environment. Reiteratively optimized structural and homology criteria were used in the computer-assisted search for retroelements, TEs that transpose by reverse transcription of an RNA intermediate, with the final results verified by manual inspection. Retroelements were found to occupy the majority (>75%) of the nuclear genome in maize inbred B73. Unprecedented genetic diversity was discovered in the long terminal repeat (LTR) retrotransposon class of retroelements, with >400 families (>350 newly discovered) contributing >31,000 intact elements. The two other classes of retroelements, SINEs (four families) and LINEs (at least 30 families), were observed to contribute 1,991 and approximately 35,000 copies, respectively, or a combined approximately 1% of the B73 nuclear genome. With regard to fully intact elements, median copy numbers for all retroelement families in maize was 2 because >250 LTR retrotransposon families contained only one or two intact members that could be detected in the B73 draft sequence. The majority, perhaps all, of the investigated retroelement families exhibited non-random dispersal across the maize genome, with LINEs, SINEs, and many low-copy-number LTR retrotransposons exhibiting a bias for accumulation in gene-rich regions. In contrast, most (but not all) medium- and high-copy-number LTR retrotransposons were found to preferentially accumulate in gene-poor regions like pericentromeric heterochromatin, while a few high-copy-number families exhibited the opposite bias. Regions of the genome with the highest LTR retrotransposon density contained the lowest LTR retrotransposon diversity. These results indicate that the maize genome provides a great number of different niches for the survival and procreation of a great variety of retroelements that have evolved to differentially occupy and exploit this genomic diversity.

Detailed analysis of a contiguous 22-Mb region of the maize genome.

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on approximately 1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.

From Methylome to Integrative Analysis of Tissue Specificity.

DNA methylation is the most studied epigenetic mark in both plants and animals. The gold standard for assaying genome-wide DNA methylation at single-base resolution is whole-genome bisulfite sequencing (WGBS). Here, we describe an improved procedure for WGBS and original bioinformatic workflows applied to unravel tissue-specific variations of the methylome in relation to gene expression and accumulation of secondary metabolites in the medicinal plant Catharanthus roseus.

Single cell RNA sequencing reveals hemocyte heterogeneity in Biomphalaria glabrata: Plasticity over diversity.

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni, the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites.