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.

BRPF3-HBO1 regulates replication origin activation and histone H3K14 acetylation.

During DNA replication, thousands of replication origins are activated across the genome. Chromatin architecture contributes to origin specification and usage, yet it remains unclear which chromatin features impact on DNA replication. Here, we perform a RNAi screen for chromatin regulators implicated in replication control by measuring RPA accumulation upon replication stress. We identify six factors required for normal rates of DNA replication and characterize a function of the bromodomain and PHD finger-containing protein 3 (BRPF3) in replication initiation. BRPF3 forms a complex with HBO1 that specifically acetylates histone H3K14, and genomewide analysis shows high enrichment of BRPF3, HBO1 and H3K14ac at ORC1-binding sites and replication origins found in the vicinity of TSSs. Consistent with this, BRPF3 is necessary for H3K14ac at selected origins and efficient origin activation. CDC45 recruitment, but not MCM2-7 loading, is impaired in BRPF3-depleted cells, identifying a BRPF3-dependent function of HBO1 in origin activation that is complementary to its role in licencing. We thus propose that BRPF3-HBO1 acetylation of histone H3K14 around TSS facilitates efficient activation of nearby replication origins.

Structure of the intergenic spacers in chicken ribosomal DNA.

BACKGROUND: Ribosomal DNA (rDNA) repeats are situated in the nucleolus organizer regions (NOR) of chromosomes and transcribed into rRNA for ribosome biogenesis. Thus, they are an essential component of eukaryotic genomes. rDNA repeat units consist of rRNA gene clusters that are transcribed into single pre-rRNA molecules, each separated by intergenic spacers (IGS) that contain regulatory elements for rRNA gene cluster transcription. Because of their high repeat content, rDNA sequences are usually absent from genome assemblies. In this work, we used the long-read sequencing technology to describe the chicken IGS and fill the knowledge gap on rDNA sequences of one of the key domesticated animals. METHODS: We used the long-read PacBio RSII technique to sequence the BAC clone WAG137G04 (Wageningen BAC library) known to contain chicken NOR elements and the HGAP workflow software suit to assemble the PacBio RSII reads. Whole-genome sequence contigs homologous to the chicken rDNA repetitive unit were identified based on the Gallus_gallus-5.0 assembly with BLAST. We used the Geneious 9.0.5 and Mega software, maximum likelihood method and Chickspress project for sequence evolution analysis, phylogenetic tree construction and analysis of the raw transcriptome data. RESULTS: Three complete IGS sequences in the White Leghorn chicken genome and one IGS sequence in the red junglefowl contig AADN04001305.1 (Gallus_gallus-5.0) were detected. They had various lengths and contained three groups of tandem repeats (some of them being very GC rich) that form highly organized arrays. Initiation and termination sites of rDNA transcription were located within small and large unique regions (SUR and LUR), respectively. No functionally significant sites were detected within the tandem repeat sequences. CONCLUSIONS: Due to the highly organized GC-rich repeats, the structure of the chicken IGS differs from that of IGS in human, apes, Xenopus or fish rDNA. However, the chicken IGS shares some molecular organization features with that of the turtles, which are other representatives of the Sauropsida clade that includes birds and reptiles. Our current results on the structure of chicken IGS together with the previously reported ribosomal gene cluster sequence provide sufficient data to consider that the complete chicken rDNA sequence is assembled with confidence in terms of molecular DNA organization.

The RAD51 paralogs ensure cellular protection against mitotic defects and aneuploidy.

The interplay between homologous DNA recombination and mitotic progression is poorly understood. The five RAD51 paralogs (RAD51B, RAD51C, RAD51D, XRCC2 and XRCC3) are key enzymes for DNA double-strand break repair. In our search for specific functions of the various RAD51 paralogs, we found that inhibition of XRCC3 elicits checkpoint defects, while inhibition of RAD51B or RAD51C induces G2/M cell cycle arrest in HeLa cells. Using live-cell microscopy we show that in XRCC3-knockdown cells the spindle assembly checkpoint persists and there is a higher frequency of chromosome misalignments, anaphase bridges, and aneuploidy. We observed centrosome defects in the absence of XRCC3. While RAD51B and RAD51C act early in homologous recombination, XRCC3 functions jointly with GEN1 later in the pathway at the stage of Holliday junction resolution. Our data demonstrate that Holliday junction resolution has critical functions for preventing aberrant mitosis and aneuploidy in mitotic cells.

Which sub-compartments of fat mass and fat-free mass are related to blood viscosity factors?

The size of body compartments is a determinant of several factors of blood viscosity. Red cell aggregation is proportional to fat mass while hematocrit is proportional to both fat-free mass and abdominal adiposity, but which parts of these body components are involved in this relationship is not known. Segmental bioelectrical impedance analysis (sBIA) provides a possibility to delineate the relationships more precisely between various subdivisions of the body and blood viscosity factors, going farther than preceding studies using non segmental BIA. In this study we investigated in 38 subjects undergoing a standardized breakfast test with mathematical modelling of glucose homeostasis and a segmental bioelectrical impedance analysis (sBIA) the relationships between the various compartments of the body and viscosity factors. Blood and plasma viscosity were measured with the Anton Paar rheometer and analyzed with Quemada's model. The parameters better correlated to hematocrit are fat free mass (r = 0.562) and its two components muscle mass (r = 0.516) and non-muscular fat-free mass (r = 0.452), and also trunk fat mass (r = 0.383) and waist-to hip ratio (r = 0.394). Red cell aggregation measurements were correlated with both truncal and appendicular fat mass (r ranging between 0.603 and 0.728). Weaker correlations of M and M1 are found with waist circumference and hip circumference. This study shows that the correlation between lean mass and hematocrit involves both muscle and non-muscle moieties of lean mass, and that both central and appendicular fat are determinants of red cell aggregation.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range.

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

MRE11-RAD50-NBS1 is a critical regulator of FANCD2 stability and function during DNA double-strand break repair.

Monoubiquitination of the Fanconi anaemia protein FANCD2 is a key event leading to repair of interstrand cross-links. It was reported earlier that FANCD2 co-localizes with NBS1. However, the functional connection between FANCD2 and MRE11 is poorly understood. In this study, we show that inhibition of MRE11, NBS1 or RAD50 leads to a destabilization of FANCD2. FANCD2 accumulated from mid-S to G2 phase within sites containing single-stranded DNA (ssDNA) intermediates, or at sites of DNA damage, such as those created by restriction endonucleases and laser irradiation. Purified FANCD2, a ring-like particle by electron microscopy, preferentially bound ssDNA over various DNA substrates. Inhibition of MRE11 nuclease activity by Mirin decreased the number of FANCD2 foci formed in vivo. We propose that FANCD2 binds to ssDNA arising from MRE11-processed DNA double-strand breaks. Our data establish MRN as a crucial regulator of FANCD2 stability and function in the DNA damage response.

Phytophthora capsici genome assembly for two isolates using long-read Oxford Nanopore Technology sequencing.

The oomycete Phytophthora capsici is a common pathogen of the Solanaceae and Cucurbitaceae families. An improved assembly for the reference isolate LT1534 was constructed using Oxford Nanopore Technologies and Illumina data. Additionally, an unpolished assembly was produced for the European isolate Pc285 collected on chili pepper using Oxford Nanopore reads.

Transcriptomic profiling reveals host-specific evolutionary pathways promoting enhanced fitness in the plant pathogen Ralstonia pseudosolanacearum.

The impact of host diversity on the genotypic and phenotypic evolution of broad-spectrum pathogens is an open issue. Here, we used populations of the plant pathogen Ralstonia pseudosolanacearum that were experimentally evolved on five types of host plants, either belonging to different botanical families or differing in their susceptibility or resistance to the pathogen. We investigated whether changes in transcriptomic profiles, associated with or independent of genetic changes, could occur during the process of host adaptation, and whether transcriptomic reprogramming was dependent on host type. Genomic and transcriptomic variations were established for 31 evolved clones that showed better fitness in their experimental host than the ancestral clone. Few genomic polymorphisms were detected in these clones, but significant transcriptomic variations were observed, with a large number of differentially expressed genes (DEGs). In a very clear way, a group of genes belonging to the network of regulation of the bacterial virulence such as efpR, efpH or hrpB, among others, were deregulated in several independent evolutionary lineages and appeared to play a key role in the transcriptomic rewiring observed in evolved clones. A double hierarchical clustering based on the 400 top DEGs for each clone revealed 2 major patterns of gene deregulation that depend on host genotype, but not on host susceptibility or resistance to the pathogen. This work therefore highlights the existence of two major evolutionary paths that result in a significant reorganization of gene expression during adaptive evolution and underscore clusters of co-regulated genes associated with bacterial adaptation on different host lines.

Comparative Genome Analysis of Enterococcus cecorum Reveals Intercontinental Spread of a Lineage of Clinical Poultry Isolates.

Enterococcus cecorum is an emerging pathogen responsible for osteomyelitis, spondylitis, and femoral head necrosis causing animal suffering and mortality and requiring antimicrobial use in poultry. Paradoxically, E. cecorum is a common inhabitant of the intestinal microbiota of adult chickens. Despite evidence suggesting the existence of clones with pathogenic potential, the genetic and phenotypic relatedness of disease-associated isolates remains little investigated. Here, we sequenced and analyzed the genomes and characterized the phenotypes of more than 100 isolates, the majority of which were collected over the last 10 years from 16 French broiler farms. Comparative genomics, genome-wide association studies, and the measured susceptibility to serum, biofilm-forming capacity, and adhesion to chicken type II collagen were used to identify features associated with clinical isolates. We found that none of the tested phenotypes could discriminate the origin of the isolates or the phylogenetic group. Instead, we found that most clinical isolates are grouped phylogenetically, and our analyses selected six genes that discriminate 94% of isolates associated with disease from those that are not. Analysis of the resistome and the mobilome revealed that multidrug-resistant clones of E. cecorum cluster into a few clades and that integrative conjugative elements and genomic islands are the main carriers of antimicrobial resistance. This comprehensive genomic analysis shows that disease-associated clones of E. cecorum belong mainly to one phylogenetic clade. IMPORTANCE Enterococcus cecorum is an important pathogen of poultry worldwide. It causes a number of locomotor disorders and septicemia, particularly in fast-growing broilers. Animal suffering, antimicrobial use, and associated economic losses require a better understanding of disease-associated E. cecorum isolates. To address this need, we performed whole-genome sequencing and analysis of a large collection of isolates responsible for outbreaks in France. By providing the first data set on the genetic diversity and resistome of E. cecorum strains circulating in France, we pinpoint an epidemic lineage that is probably also circulating elsewhere that should be targeted preferentially by preventive strategies in order to reduce the burden of E. cecorum-related diseases.

Genome structures resolve the early diversification of teleost fishes.

Accurate species phylogenies are a prerequisite for all evolutionary research. Teleosts are the largest and most diversified group of extant vertebrates, but relationships among their three oldest extant lineages remain unresolved. On the basis of seven high-quality new genome assemblies in Elopomorpha (tarpons, eels), we revisited the topology of the deepest branches of the teleost phylogeny using independent gene sequence and chromosomal rearrangement phylogenomic approaches. These analyses converged to a single scenario that unambiguously places the Elopomorpha and Osteoglossomorpha (arapaima, elephantnose fish) in a monophyletic sister group to all other teleosts, i.e., the Clupeocephala lineage (zebrafish, medaka). This finding resolves more than 50 years of controversy on the evolutionary relationships of these lineages and highlights the power of combining different levels of genome-wide information to solve complex phylogenies.

Suspicions of two bridgehead invasions of Xylella fastidiosa subsp. multiplex in France.

Of American origin, a wide diversity of Xylella fastidiosa strains belonging to different subspecies have been reported in Europe since 2013 and its discovery in Italian olive groves. Strains from the subspecies multiplex (ST6 and ST7) were first identified in France in 2015 in urban and natural areas. To trace back the most probable scenario of introduction in France, the molecular evolution rate of this subspecies was estimated at 3.2165 × 10-7 substitutions per site per year, based on heterochronous genome sequences collected worldwide. This rate allowed the dating of the divergence between French and American strains in 1987 for ST6 and in 1971 for ST7. The development of a new VNTR-13 scheme allowed tracing the spread of the bacterium in France, hypothesizing an American origin. Our results suggest that both sequence types were initially introduced and spread in Provence-Alpes-Côte d'Azur (PACA); then they were introduced in Corsica in two waves from the PACA bridgehead populations.