Changes in DNA methylation contribute to rapid adaptation in bacterial plant pathogen evolution.

Adaptation is usually explained by beneficial genetic mutations that are transmitted from parents to offspring and become fixed in the adapted population. However, genetic mutation analysis alone is not sufficient to fully explain the adaptive processes, and several studies report the existence of nongenetic (or epigenetic) inheritance that can enable adaptation to new environments. In the present work, we tested the hypothesis of the role of DNA methylation, a form of epigenetic modification, in adaptation of the plant pathogen Ralstonia pseudosolanacearum to the host during experimental evolution. Using SMRT-seq technology, we analyzed the methylomes of 31 experimentally evolved clones obtained after serial passages on 5 different plant species during 300 generations. Comparison with the methylome of the ancestral clone revealed a list of 50 differential methylated sites (DMSs) at the GTWWAC motif. Gene expression analysis of the 39 genes targeted by these DMSs revealed limited correlation between differential methylation and differential expression of the corresponding genes. Only 1 gene showed a correlation, the RSp0338 gene encoding the EpsR regulator protein. The MSRE-qPCR technology, used as an alternative approach for DNA methylation analysis, also found the 2 DMSs upstream RSp0338. Using site-directed mutagenesis, we demonstrated the contribution of these 2 DMSs in host adaptation. As these DMSs appeared very early in the experimental evolution, we hypothesize that such fast epigenetic changes can allow rapid adaptation to the plant stem environment. In addition, we found that the change in DNA methylation upstream RSp0338 remains stable at least for 100 generations outside the host and thus can contribute to long-term adaptation to the host plant. To our knowledge, this is the first study showing a direct link between bacterial epigenetic variation and adaptation to a new environment.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae.

BACKGROUND: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. RESULTS: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. CONCLUSIONS: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes.

BACKGROUND: The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses. RESULTS: Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants. CONCLUSIONS: Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

The TIP60 Complex Regulates Bivalent Chromatin Recognition by 53BP1 through Direct H4K20me Binding and H2AK15 Acetylation.

The NuA4/TIP60 acetyltransferase complex is a key regulator of genome expression and stability. Here we identified MBTD1 as a stable subunit of the complex, and we reveal that, via a histone reader domain for H4K20me1/2, MBTD1 allows TIP60 to associate with specific gene promoters and to promote the repair of DNA double-strand breaks by homologous recombination. It was previously suggested that TIP60-dependent acetylation of H4 regulates binding of the non-homologous end joining factor 53BP1, which engages chromatin through simultaneous binding of H4K20me2 and H2AK15ub. We find that the TIP60 complex regulates association of 53BP1 partly by competing for H4K20me2 and by regulating H2AK15ub. Ubiquitylation of H2AK15 by RNF168 inhibits chromatin acetylation by TIP60, while this residue can be acetylated by TIP60 in vivo, blocking its ubiquitylation. Altogether, these results uncover an intricate mechanism orchestrated by the TIP60 complex to regulate 53BP1-dependent repair through competitive bivalent binding and modification of chromatin.

MRG Proteins Are Shared by Multiple Protein Complexes With Distinct Functions.

MRG15/MORF4L1 is a highly conserved protein in eukaryotes that contains a chromodomain (CHD) recognizing methylation of lysine 36 on histone H3 (H3K36me3) in chromatin. Intriguingly, it has been reported in the literature to interact with several different factors involved in chromatin modifications, gene regulation, alternative mRNA splicing, and DNA repair by homologous recombination. To get a complete and reliable picture of associations in physiological conditions, we used genome editing and tandem affinity purification to analyze the stable native interactome of human MRG15, its paralog MRGX/MORF4L2 that lacks the CHD, and MRGBP (MRG-binding protein) in isogenic K562 cells. We found stable interchangeable association of MRG15 and MRGX with the NuA4/TIP60 histone acetyltransferase/chromatin remodeler, Sin3B histone deacetylase/demethylase, ASH1L histone methyltransferase, and PALB2-BRCA2 DNA repair protein complexes. These associations were further confirmed and analyzed by CRISPR tagging of endogenous proteins and comparison of expressed isoforms. Importantly, based on structural information, point mutations could be introduced that specifically disrupt MRG15 association with some complexes but not others. Most interestingly, we also identified a new abundant native complex formed by MRG15/X-MRGBP-BRD8-EP400NL (EP400 N-terminal like) that is functionally similar to the yeast TINTIN (Trimer Independent of NuA4 for Transcription Interactions with Nucleosomes) complex. Our results show that EP400NL, being homologous to the N-terminal region of NuA4/TIP60 subunit EP400, creates TINTIN by competing for BRD8 association. Functional genomics indicate that human TINTIN plays a role in transcription of specific genes. This is most likely linked to the H4ac-binding bromodomain of BRD8 along the H3K36me3-binding CHD of MRG15 on the coding region of transcribed genes. Taken together, our data provide a complete detailed picture of human MRG proteins-associated protein complexes, which are essential to understand and correlate their diverse biological functions in chromatin-based nuclear processes.

Convergent Rewiring of the Virulence Regulatory Network Promotes Adaptation of Ralstonia solanacearum on Resistant Tomato.

The evolutionary and adaptive potential of a pathogen is a key determinant for successful host colonization and proliferation but remains poorly known for most of the pathogens. Here, we used experimental evolution combined with phenotyping, genomics, and transcriptomics to estimate the adaptive potential of the bacterial plant pathogen Ralstonia solanacearum to overcome the quantitative resistance of the tomato cultivar Hawaii 7996. After serial passaging over 300 generations, we observed pathogen adaptation to within-plant environment of the resistant cultivar but no plant resistance breakdown. Genomic sequence analysis of the adapted clones revealed few genetic alterations, but we provide evidence that all but one were gain of function mutations. Transcriptomic analyses revealed that even if different adaptive events occurred in independently evolved clones, there is convergence toward a global rewiring of the virulence regulatory network as evidenced by largely overlapping gene expression profiles. A subset of four transcription regulators, including HrpB, the activator of the type 3 secretion system regulon and EfpR, a global regulator of virulence and metabolic functions, emerged as key nodes of this regulatory network that are frequently targeted to redirect the pathogen's physiology and improve its fitness in adverse conditions. Significant transcriptomic variations were also detected in evolved clones showing no genomic polymorphism, suggesting that epigenetic modifications regulate expression of some of the virulence network components and play a major role in adaptation as well.

DNAModAnnot: a R toolbox for DNA modification filtering and annotation.

MOTIVATION: Long-read sequencing technologies can be employed to detect and map DNA modifications at the nucleotide resolution on a genome-wide scale. However, published software packages neglect the integration of genomic annotation and comprehensive filtering when analyzing patterns of modified bases detected using Pacific Biosciences (PacBio) or Oxford Nanopore Technologies (ONT) data. Here, we present DNA Modification Annotation (DNAModAnnot), a R package designed for the global analysis of DNA modification patterns using adapted filtering and visualization tools. RESULTS: We tested our package using PacBio sequencing data to analyze patterns of the 6-methyladenine (6mA) in the ciliate Paramecium tetraurelia, in which high 6mA amounts were previously reported. We found P. tetraurelia 6mA genome-wide distribution to be similar to other ciliates. We also performed 5-methylcytosine (5mC) analysis in human lymphoblastoid cells using ONT data and confirmed previously known patterns of 5mC. DNAModAnnot provides a toolbox for the genome-wide analysis of different DNA modifications using PacBio and ONT long-read sequencing data. AVAILABILITY AND IMPLEMENTATION: DNAModAnnot is distributed as a R package available via GitHub (https://github.com/AlexisHardy/DNAModAnnot). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of the master sex determining gene in Northern pike (Esox lucius) reveals restricted sex chromosome differentiation.

Teleost fishes, thanks to their rapid evolution of sex determination mechanisms, provide remarkable opportunities to study the formation of sex chromosomes and the mechanisms driving the birth of new master sex determining (MSD) genes. However, the evolutionary interplay between the sex chromosomes and the MSD genes they harbor is rather unexplored. We characterized a male-specific duplicate of the anti-Müllerian hormone (amh) as the MSD gene in Northern Pike (Esox lucius), using genomic and expression evidence as well as by loss-of-function and gain-of-function experiments. Using RAD-Sequencing from a family panel, we identified Linkage Group (LG) 24 as the sex chromosome and positioned the sex locus in its sub-telomeric region. Furthermore, we demonstrated that this MSD originated from an ancient duplication of the autosomal amh gene, which was subsequently translocated to LG24. Using sex-specific pooled genome sequencing and a new male genome sequence assembled using Nanopore long reads, we also characterized the differentiation of the X and Y chromosomes, revealing a small male-specific insertion containing the MSD gene and a limited region with reduced recombination. Our study reveals an unexpectedly low level of differentiation between a pair of sex chromosomes harboring an old MSD gene in a wild teleost fish population, and highlights both the pivotal role of genes from the amh pathway in sex determination, as well as the importance of gene duplication as a mechanism driving the turnover of sex chromosomes in this clade.

Characterisation of hydrocarbon degradation, biosurfactant production, and biofilm formation in Serratia sp. Tan611: a new strain isolated from industrially contaminated environment in Algeria.

A novel bacterial strain was isolated from industrially contaminated waste water. In the presence of crude oil, this strain was shown to reduce the rate of total petroleum hydrocarbons (TPH) up to 97.10% in 24 h. This bacterium was subsequently identified by 16S rRNA gene sequence analysis and affiliated to the Serratia genus by the RDP classifier. Its genome was sequenced and annotated, and genes coding for catechol 1,2 dioxygenase and naphthalene 1,2-dioxygenase system involved in aromatic hydrocarbon catabolism, and LadA-type monooxygenases involved in alkane degradation, were identified. Gas Chromatography-Mass Spectrometry (GC-MS) analysis of crude oil after biological treatment showed that Serratia sp. Tan611 strain was able to degrade n-alkanes (from C13 to C25). This bacterium was also shown to produce a biosurfactant, the emulsification index (E24) reaching 43.47% and 65.22%, against vegetable and crude oil, respectively. Finally, the formation of a biofilm was increased in the presence of crude oil. These observations make Serratia sp. Tan611 a good candidate for hydrocarbon bioremediation.

A Chromosome-Scale Genome Assembly Resource for Myriosclerotinia sulcatula Infecting Sedge Grass (Carex sp.).

The fungus Myriosclerotinia sulcatula is a close relative of the notorious polyphagous plant pathogens Botrytis cinerea and Sclerotinia sclerotiorum but exhibits a host range restricted to plants from the Carex genus (Cyperaceae family). To date, there are no genomic resources available for fungi in the Myriosclerotinia genus. Here, we present a chromosome-scale reference genome assembly for M. sulcatula. The assembly contains 24 contigs with a total length of 43.53 Mbp, with scaffold N50 of 2,649.7 kbp and N90 of 1,133.1 kbp. BRAKER-predicted gene models were manually curated using WebApollo, resulting in 11,275 protein-coding genes that we functionally annotated. We provide a high-quality reference genome assembly and annotation for M. sulcatula as a resource for studying evolution and pathogenicity in fungi from the Sclerotiniaceae family.