Increased evolutionary rate in the Z-chromosome of sympatric and allopatric species of Morpho butterflies.

Divergent evolution of genomes among closely related species is shaped by both neutral processes and ecological forces, such as local adaptation and reinforcement. These factors can drive accelerated evolution of sex chromosomes relative to autosomes. Comparative genomic analyses between allopatric and sympatric species with overlapping or divergent ecological niches offer insights into reinforcement and ecological specialization on genome evolution. In the butterfly genus Morpho, several species coexist in sympatry, with specialization across forest strata and temporal niches. We analyzed the genomes of eight Morpho species, along with previously published genomes of three others, to compare chromosomal rearrangements and signs of positive selection in the Z chromosome vs autosomes. We found extensive chromosomal rearrangements in Z chromosome, particularly in sympatric species with similar ecological niches, suggesting a role for inversions in restricting gene flow at a postzygotic level. Z-linked genes also exhibited significantly higher dN/dS ratios than autosomal genes across the genus, with pronounced differences in closely related species living in sympatry. Additionally, we examined the evolution of eight circadian clock genes, detecting positive selection in Period, located on the Z chromosome. Our findings suggest that the Z chromosome evolves more rapidly than autosomes, particularly among closely related species, raising questions about its role in pre- and post-zygotic isolation mechanisms.

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.

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.

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.

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.

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.

A chromosome-level, haplotype-phased Vanilla planifolia genome highlights the challenge of partial endoreplication for accurate whole-genome assembly.

Vanilla planifolia, the species cultivated to produce one of the world's most popular flavors, is highly prone to partial genome endoreplication, which leads to highly unbalanced DNA content in cells. We report here the first molecular evidence of partial endoreplication at the chromosome scale by the assembly and annotation of an accurate haplotype-phased genome of V. planifolia. Cytogenetic data demonstrated that the diploid genome size is 4.09 Gb, with 16 chromosome pairs, although aneuploid cells are frequently observed. Using PacBio HiFi and optical mapping, we assembled and phased a diploid genome of 3.4 Gb with a scaffold N50 of 1.2 Mb and 59 128 predicted protein-coding genes. The atypical k-mer frequencies and the uneven sequencing depth observed agreed with our expectation of unbalanced genome representation. Sixty-seven percent of the genes were scattered over only 30% of the genome, putatively linking gene-rich regions and the endoreplication phenomenon. By contrast, low-coverage regions (non-endoreplicated) were rich in repeated elements but also contained 33% of the annotated genes. Furthermore, this assembly showed distinct haplotype-specific sequencing depth variation patterns, suggesting complex molecular regulation of endoreplication along the chromosomes. This high-quality, anchored assembly represents 83% of the estimated V. planifolia genome. It provides a significant step toward the elucidation of this complex genome. To support post-genomics efforts, we developed the Vanilla Genome Hub, a user-friendly integrated web portal that enables centralized access to high-throughput genomic and other omics data and interoperable use of bioinformatics tools.

Complete Genome Sequence of Stenotrophomonas maltophilia 1800, a New Bacterial Strain with Potential for Bioremediation of Oil-Contaminated Environments.

Stenotrophomonas maltophilia strain 1800 was isolated from the effluent of an industrial oil refinery in Algeria. Its genome was sequenced using Illumina MiSeq (2 × 150-bp read pairs) and Oxford Nanopore (long reads) technologies and assembled using Unicycler. It is composed of one chromosome of 4.83 Mb.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage.

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.

A duplicated copy of id2b is an unusual sex-determining candidate gene on the Y chromosome of arapaima (Arapaima gigas).

Arapaima gigas is one of the largest freshwater fish species of high ecological and economic importance. Overfishing and habitat destruction are severe threats to the remaining wild populations. By incorporating a chromosomal Hi-C contact map, we improved the arapaima genome assembly to chromosome-level, revealing an unexpected high degree of chromosome rearrangements during evolution of the bonytongues (Osteoglossiformes). Combining this new assembly with pool-sequencing of male and female genomes, we identified id2bbY, a duplicated copy of the inhibitor of DNA binding 2b (id2b) gene on the Y chromosome as candidate male sex-determining gene. A PCR-test for id2bbY was developed, demonstrating that this gene is a reliable male-specific marker for genotyping. Expression analyses showed that this gene is expressed in juvenile male gonads. Its paralog, id2ba, exhibits a male-biased expression in immature gonads. Transcriptome analyses and protein structure predictions confirm id2bbY as a prime candidate for the master sex-determiner. Acting through the TGFß signaling pathway, id2bbY from arapaima would provide the first evidence for a link of this family of transcriptional regulators to sex determination. Our study broadens our current understanding about the evolution of sex determination genetic networks and provide a tool for improving arapaima aquaculture for commercial and conservation purposes.

Detection and Typing of a Fowl Adenovirus Type 1 Agent of Pancreatitis in Guinea Fowl.

Adenoviral pancreatitis has been amply described for decades in guinea fowl. Although its pathologic picture has been characterized fairly well, its etiology still remains only partially clarified. Based on several outbreaks diagnosed on commercial guinea flocks raised in France since 2017, we performed direct whole-genome sequencing from pancreatic lesional tissue by using the Oxford Nanopore Technologies (ONT) sequencing method. We generated 4781 viral reads and assembled a whole genome of 43,509 bp, clustering within fowl adenovirus type 1 (FAdV-1). A phylogenetic analysis based on a partial sequence of the hexon and short fiber genes on viruses collected in France showed 98.7% and 99.8% nucleotide identity, respectively. Altogether, these results confirm that an FAdV-1 closely related to chicken and other avian strains is the agent of pancreatitis in guinea fowl. This study illustrates the potential of ONT sequencing method to achieve rapid whole-genome sequencing directly from pathologic material.

Detección y tipificación de un adenovirus aviar tipo 1 (FAdV-1), agente de pancreatitis en gallinas de Guinea. La pancreatitis adenoviral se ha descrito ampliamente durante décadas en gallinas de Guinea. Aunque su cuadro patológico se ha caracterizado bastante bien, su etiología todavía permanece sólo parcialmente aclarada. Sobre la base de varios brotes diagnosticados en parvadas comerciales de guineas criadas en Francia desde el año 2017, se realizó una secuenciación directa del genoma completo a partir del tejido de la lesión pancreática mediante el método de secuenciación desarrollado por Oxford Nanopore Technologies. Se generaron 4781 lecturas virales y se ensambló un genoma completo de 43,509 pb, que se agrupó dentro del adenovirus aviar tipo 1 (FAdV-1). Un análisis filogenético basado en una secuencia parcial de los genes hexón y de fibra corta de virus recolectados en Francia mostró identidades de nucleótidos de 98.7% y 99.8%, respectivamente. En conjunto, estos resultados confirman que un adenovirus aviar tipo 1 estrechamente relacionado con el pollo y otras cepas aviares es el agente de la pancreatitis en la gallina de Guinea. Este estudio ilustra el potencial de las tecnologías desarrolladas por Oxford Nanopore Thechnologies para lograr una secuenciación rápida de todo el genoma directamente a partir de material patológico.

A supernumerary "B-sex" chromosome drives male sex determination in the Pachón cavefish, Astyanax mexicanus.

Sex chromosomes are generally derived from a pair of classical type-A chromosomes, and relatively few alternative models have been proposed up to now.1,2 B chromosomes (Bs) are supernumerary and dispensable chromosomes with non-Mendelian inheritance found in many plant and animal species3,4 that have often been considered as selfish genetic elements that behave as genome parasites.5,6 The observation that in some species Bs can be either restricted or predominant in one sex7-14 raised the interesting hypothesis that Bs could play a role in sex determination.15 The characterization of putative B master sex-determining (MSD) genes, however, has not yet been provided to support this hypothesis. Here, in Astyanax mexicanus cavefish originating from Pachón cave, we show that Bs are strongly male predominant. Based on a high-quality genome assembly of a B-carrying male, we characterized the Pachón cavefish B sequence and found that it contains two duplicated loci of the putative MSD gene growth differentiation factor 6b (gdf6b). Supporting its role as an MSD gene, we found that the Pachón cavefish gdf6b gene is expressed specifically in differentiating male gonads, and that its knockout induces male-to-female sex reversal in B-carrying males. This demonstrates that gdf6b is necessary for triggering male sex determination in Pachón cavefish. Altogether these results bring multiple and independent lines of evidence supporting the conclusion that the Pachón cavefish B is a "B-sex" chromosome that contains duplicated copies of the gdf6b gene, which can promote male sex determination in this species.

A 180 Myr-old female-specific genome region in sturgeon reveals the oldest known vertebrate sex determining system with undifferentiated sex chromosomes.

Several hypotheses explain the prevalence of undifferentiated sex chromosomes in poikilothermic vertebrates. Turnovers change the master sex determination gene, the sex chromosome or the sex determination system (e.g. XY to WZ). Jumping master genes stay main triggers but translocate to other chromosomes. Occasional recombination (e.g. in sex-reversed females) prevents sex chromosome degeneration. Recent research has uncovered conserved heteromorphic or even homomorphic sex chromosomes in several clades of non-avian and non-mammalian vertebrates. Sex determination in sturgeons (Acipenseridae) has been a long-standing basic biological question, linked to economical demands by the caviar-producing aquaculture. Here, we report the discovery of a sex-specific sequence from sterlet (Acipenser ruthenus). Using chromosome-scale assemblies and pool-sequencing, we first identified an approximately 16 kb female-specific region. We developed a PCR-genotyping test, yielding female-specific products in six species, spanning the entire phylogeny with the most divergent extant lineages (A. sturio, A. oxyrinchus versus A. ruthenus, Huso huso), stemming from an ancient tetraploidization. Similar results were obtained in two octoploid species (A. gueldenstaedtii, A. baerii). Conservation of a female-specific sequence for a long period, representing 180 Myr of sturgeon evolution, and across at least one polyploidization event, raises many interesting biological questions. We discuss a conserved undifferentiated sex chromosome system with a ZZ/ZW-mode of sex determination and potential alternatives. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Population genomics of apricots unravels domestication history and adaptive events.

Among crop fruit trees, the apricot (Prunus armeniaca) provides an excellent model to study divergence and adaptation processes. Here, we obtain nearly 600 Armeniaca apricot genomes and four high-quality assemblies anchored on genetic maps. Chinese and European apricots form two differentiated gene pools with high genetic diversity, resulting from independent domestication events from distinct wild Central Asian populations, and with subsequent gene flow. A relatively low proportion of the genome is affected by selection. Different genomic regions show footprints of selection in European and Chinese cultivated apricots, despite convergent phenotypic traits, with predicted functions in both groups involved in the perennial life cycle, fruit quality and disease resistance. Selection footprints appear more abundant in European apricots, with a hotspot on chromosome 4, while admixture is more pervasive in Chinese cultivated apricots. Our study provides clues to the biology of selected traits and targets for fruit tree research and breeding.

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.

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.

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.

Complete Genome Sequence of Microbacterium sp. Strain Nx66, Isolated from Waters Contaminated with Petrochemicals in El Saf-Saf Valley, Algeria.

Microbacterium sp. strain Nx66 was isolated from waters contaminated by petrochemical effluents collected in Algeria. Its genome was sequenced using Illumina MiSeq (2 × 150-bp read pairs) and Oxford Nanopore (long reads) technologies and was assembled using Unicycler. It is composed of one chromosome of 3.42 Mb and one plasmid of 34.22 kb.

Genome structure and content of the rice root-knot nematode (Meloidogyne graminicola).

Discovered in the 1960s, Meloidogyne graminicola is a root-knot nematode species considered as a major threat to rice production. Yet, its origin, genomic structure, and intraspecific diversity are poorly understood. So far, such studies have been limited by the unavailability of a sufficiently complete and well-assembled genome. In this study, using a combination of Oxford Nanopore Technologies and Illumina sequencing data, we generated a highly contiguous reference genome (283 scaffolds with an N50 length of 294 kb, totaling 41.5 Mb). The completeness scores of our assembly are among the highest currently published for Meloidogyne genomes. We predicted 10,284 protein-coding genes spanning 75.5% of the genome. Among them, 67 are identified as possibly originating from horizontal gene transfers (mostly from bacteria), which supposedly contribute to nematode infection, nutrient processing, and plant defense manipulation. Besides, we detected 575 canonical transposable elements (TEs) belonging to seven orders and spanning 2.61% of the genome. These TEs might promote genomic plasticity putatively related to the evolution of M. graminicola parasitism. This high-quality genome assembly constitutes a major improvement regarding previously available versions and represents a valuable molecular resource for future phylogenomic studies of Meloidogyne species. In particular, this will foster comparative genomic studies to trace back the evolutionary history of M. graminicola and its closest relatives.