Search

Abstract Background Aphid adaptation to harsh winter conditions is illustrated by an alternation of their reproductive mode. Aphids detect photoperiod shortening by sensing the length of the night and switch from viviparous parthenogenesis in spring and summer, to oviparous sexual reproduction in autumn. The photoperiodic signal is transduced from the head to the reproductive tract to change the fate of the future oocytes from mitotic diploid embryogenesis to haploid formation of gametes. This process takes place in three consecutive generations due to viviparous parthenogenesis. To understand the molecular basis of the switch in the reproductive mode, transcriptomic and proteomic approaches were used to detect significantly regulated transcripts and polypeptides in the heads of the pea aphid Acyrthosiphon pisum. Results The transcriptomic profiles of the heads of the first generation were slightly affected by photoperiod shortening. This suggests that trans-generation signalling between the grand-mothers and the viviparous embryos they contain is not essential. By analogy, many of the genes and some of the proteins regulated in the heads of the second generation are implicated in visual functions, photoreception and cuticle structure. The modification of the cuticle could be accompanied by a down-regulation of the N-β-alanyldopamine pathway and desclerotization. In Drosophila, modification of the insulin pathway could cause a decrease of juvenile hormones in short-day reared aphids. Conclusion This work led to the construction of hypotheses for photoperiodic regulation of the switch of the reproductive mode in aphids.

Males and females share essentially the same genome but differ in their optimal values for many phenotypic traits, which can result in intra-locus conflict between the sexes. Aphids display XX/X0 sex chromosomes and combine unusual X chromosome inheritance with cyclical parthenogenesis. Theoretical and empirical works support the hypothesis that the large excess of male-biased genes observed on the aphid X chromosome compared to autosomes evolved in response to sexual conflicts, by restricting the products of sexually antagonistic alleles to the sex they benefits. However, whether such masculinization of the X affects all tissues (as expected if it evolved in response to sexual conflicts) or is limited to specific tissues remains an open question. Here, we measured gene expression in three different somatic and gonadic tissues of males, sexual females and parthenogenetic females of the pea aphid. We observed a masculinization of the X in each of the studied tissues, with male-biased genes being 2.5 to 3.5 more frequent on the X than expected. We also tested the hypothesis that gene duplication can facilitate the attenuation of conflicts by allowing gene copies to neo- or sub-functionalize and reach sex-specific optima. As predicted, X-linked copies of duplicated genes having their other copies on autosomes were more frequently male-biased (40.5% of the genes) than duplicated autosomal genes (6.6%) or X-linked single-copy genes (32.5%). These results highlight a peculiar pattern of expression of X-linked genes in aphids at the tissue level and provide further support for sex-biased expression as a mechanism to attenuate intra-locus sexual conflicts.

Improving resistance durability involves to be able to predict the adaptation speed of pathogen populations. Identifying the genetic bases of pathogen adaptation to plant resistances is a useful step to better understand and anticipate this phenomenon. Globodera pallida is a major pest of potato crop for which a resistance QTL, GpaV vrn, has been identified in Solanum vernei. However, its durability is threatened as G. pallida populations are able to adapt to the resistance in few generations. The aim of this study was to investigate the genomic regions involved in the resistance breakdown by coupling experimental evolution and high-density genome scan. We performed a whole-genome resequencing of pools of individuals (Pool-Seq) belonging to G. pallida lineages derived from two independent populations having experimentally evolved on susceptible and resistant potato cultivars. About 1.6 million SNPs were used to perform the genome scan using a recent model testing for adaptive differentiation and association to population-specific covariables. We identified 275 outliers and 31 of them, which also showed a significant reduction in diversity in adapted lineages, were investigated for their genic environment. Some candidate genomic regions contained genes putatively encoding effectors and were enriched in SPRYSECs, known in cyst nematodes to be involved in pathogenicity and in (a)virulence. Validated candidate SNPs will provide a useful molecular tool to follow frequencies of virulence alleles in natural G. pallida populations and define efficient strategies of use of potato resistances maximizing their durability. [ABSTRACT FROM AUTHOR]

Unmapped reads are often discarded from the analysis of whole-genome re-sequencing, but new biological information and insights can be uncovered through their analysis. In this paper, we investigate unmapped reads from the re-sequencing data of 33 pea aphid genomes from individuals specialized on different host plants. The unmapped reads for each individual were retrieved following mapping to the Acyrthosiphon pisum reference genome and its mitochondrial and symbiont genomes. These sets of unmapped reads were then cross-compared, revealing that a significant number of these unmapped sequences were conserved across individuals. Interestingly, sequences were most commonly shared between individuals adapted to the same host plant, suggesting that these sequences may contribute to the divergence between host plant specialized biotypes. Analysis of the contigs obtained from assembling the unmapped reads pooled by biotype allowed us to recover some divergent genomic regions previously excluded from analysis and to discover putative novel sequences of A. pisum and its symbionts. In conclusion, this study emphasizes the interest of the unmapped component of re-sequencing data sets and the potential loss of important information. We here propose strategies to aid the capture and interpretation of this information. [ABSTRACT FROM AUTHOR]

The recent assembly of the silkworm Bombyx mon genome with 432 Mb on 28 holocentric chromosomes has become a reference in the genomic analysis of the very diverse Order of Lepidoptera. We sequenced BAC5 from two major pests, the noctuid moths Helicoverpa armigera and Spodoptera frugiperda. corresponding to 15 regions distributed on 11 B. mon chromosomes, each BAC/region being anchored by known orthologous gene(s) to analyze syntenic relationships and genome rearrangements among the three species. Nearly 300 genes and numerous transposable elements were identified, with long interspersed nuclear elements and terminal inverted repeats the most abundant transposable element classes. There was a high degree of synteny conservation between B. mon and the two noctuid species. Conserved syntenic blocks of identified genes were very small, however, approximately 1.3 genes per block between B. mon and the two noctuid species and 2.0 genes per block between S. frugiperda and H. armigera. This corresponds to approximately two chromosome breaks per Mb DNA per My. This is a much higher evolution rate than among species of the Drosophila genus and may be related to the holocentric nature of the lepidopteran genomes. We report a large cluster of eight members of the aminopeptidase N gene family that we estimate to have been present since the Jurassic. In contrast, several clusters of cytochrome P450 genes showed multiple lineage-specific duplication events, in particular in the lepidopteran CYP9A subfamily. Our study highlights the value of the silkworm genome as a reference in lepidopteran comparative genomics. [ABSTRACT FROM AUTHOR]

The expansion of multi-omics datasets raises significant challenges for data integration and querying. To overcome these challenges, we developed a generic RDF-based integration schema that connects various types of differential -omics data, epigenomics, and regulatory information. This schema employs the FALDO ontology to enable querying based on genomic locations. It is designed to be fully or partially populated, providing both flexibility and extensibility while supporting complex queries. We validated the schema by reproducing two recently published studies, one in biomedicine and the other in environmental science, proving its genericity and its ability to integrate data efficiently. This schema serves as an effective tool for managing and querying a wide range of multi-omics datasets., Competing Interests: The authors have no conflict of interest to disclose regarding this study., (© 2024 The Authors.)

We present the first nuclear genome assembly and a complete mitogenome for Hylesia metabus (Arthropoda; Insecta; Lepidoptera; Saturniidae). The assembled nuclear genome sequence is 1,271 Mb long, which is among the 10 largest lepidopteran genome assemblies published to date. It is scaffolded in 31 pseudo chromosomes, has a BUSCO score of 99.5%, and has a highly conserved synteny compared to phylogenetically close species. Repetitive elements make up 67% of the nuclear genome and are mainly located in intergenic regions, among which LINEs were predominant, with CR1-Zenon being the most abundant. Phylogenetic and comparative analyses of H. metabus assembly and 17 additional Saturniidae and Sphingidae assemblies suggested that an accumulation of repetitive elements likely led to the increased size of H. metabus' genome. Gene annotation using Helixer identified 26,122 transcripts. The Z scaffold was identified using both a synteny analysis and variations of coverage for two resequenced male and female H. metabus. The H. metabus nuclear genome and mitogenome assemblies can be found and browsed on the BIPAA website and constitute useful resources for future population and comparative genomics studies., (© The American Genetic Association. 2024.)

Grapevine downy mildew, caused by the oomycete Plasmopara viticola (P. viticola, Berk. & M. A. Curtis; Berl. & De Toni), is a global threat to Eurasian wine grapes Vitis vinifera. Although resistant grapevine varieties are becoming more accessible, P. viticola populations are rapidly evolving to overcome these resistances. We aimed to uncover avirulence genes related to Rpv3.1-mediated grapevine resistance. We sequenced the genomes and characterized the development of 136 P. viticola strains on resistant and sensitive grapevine cultivars. A genome-wide association study was conducted to identify genomic variations associated with resistant-breaking phenotypes. We identified a genomic region associated with the breakdown of Rpv3.1 grapevine resistance (avrRpv3.1 locus). A diploid-aware reassembly of the P. viticola INRA-Pv221 genome revealed structural variations in this locus, including a 30 kbp deletion. Virulent P. viticola strains displayed multiple deletions on both haplotypes at the avrRpv3.1 locus. These deletions involve two paralog genes coding for proteins with 800-900 amino acids and signal peptides. These proteins exhibited a structure featuring LWY-fold structural modules, common among oomycete effectors. When transiently expressed, these proteins induced cell death in grapevines carrying Rpv3.1 resistance, confirming their avirulence nature. This discovery sheds light on the genetic mechanisms enabling P. viticola to adapt to grapevine resistance, laying a foundation for developing strategies to manage this destructive crop pathogen., (© 2024 The Authors. New Phytologist © 2024 New Phytologist Foundation.)

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposoter didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Lorenzi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

We present the first chromosome-level genome assembly and annotation of the pearly heath Coenonympha arcania, generated with a PacBio HiFi sequencing approach and complemented with Hi-C data. We additionally compare synteny, gene, and repeat content between C. arcania and other Lepidopteran genomes. This reference genome will enable future population genomics studies with Coenonympha butterflies, a species-rich genus that encompasses some of the most highly endangered butterfly taxa in Europe., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The endosymbiotic bacteria Spiroplasma (Mollicutes) infect diverse plants and arthropods, and some of which induce male killing, where male hosts are killed during development. Male-killing Spiroplasma strains belong to either the phylogenetically distant Citri-Poulsonii or Ixodetis groups. In Drosophila flies, Spiroplasma poulsonii induces male killing via the Spaid toxin. While Spiroplasma ixodetis infects a wide range of insects and arachnids, little is known about the genetic basis of S. ixodetis-induced male killing. Here, we analyzed the genome of S. ixodetis strains in the pea aphid Acyrthosiphon pisum (Aphididae, Hemiptera). Genome sequencing constructed a complete genome of a male-killing strain, sAp269, consisting of a 1.5 Mb circular chromosome and an 80 Kb plasmid. sAp269 encoded putative virulence factors containing either ankyrin repeat, ovarian tumor-like deubiquitinase, or ribosome inactivating protein domains, but lacked the Spaid toxin. Further comparative genomics of Spiroplasma strains in A. pisum biotypes adapted to different host plants revealed their phylogenetic associations and the diversity of putative virulence factors. Although the mechanisms of S. ixodetis-induced male killing in pea aphids remain elusive, this study underlines the dynamic genome evolution of S. ixodetis and proposes independent acquisition events of male-killing mechanisms in insects., (© The Author(s) 2024. Published by Oxford University Press on behalf of FEMS.)

In order to capture the drought impacts on seed quality acquisition in Brassica napus and its potential interaction with early biotic stress, seeds of the 'Express' genotype of oilseed rape were characterized from late embryogenesis to full maturity from plants submitted to reduced watering (WS) with or without pre-occurring inoculation by the telluric pathogen Plasmodiophora brassicae (Pb + WS or Pb, respectively), and compared to control conditions (C). Drought as a single constraint led to significantly lower accumulation of lipids, higher protein content and reduced longevity of the WS-treated seeds. In contrast, when water shortage was preceded by clubroot infection, these phenotypic differences were completely abolished despite the upregulation of the drought sensor RD20. A weighted gene co-expression network of seed development in oilseed rape was generated using 72 transcriptomes from developing seeds from the four treatments and identified 33 modules. Module 29 was highly enriched in heat shock proteins and chaperones that showed a stronger upregulation in Pb + WS compared to the WS condition, pointing to a possible priming effect by the early P. brassicae infection on seed quality acquisition. Module 13 was enriched with genes encoding 12S and 2S seed storage proteins, with the latter being strongly upregulated under WS conditions. Cis-element promotor enrichment identified PEI1/TZF6, FUS3 and bZIP68 as putative regulators significantly upregulated upon WS compared to Pb + WS. Our results provide a temporal co-expression atlas of seed development in oilseed rape and will serve as a resource to characterize the plant response towards combinations of biotic and abiotic stresses., (© 2024 The Authors. Physiologia Plantarum published by John Wiley & Sons Ltd on behalf of Scandinavian Plant Physiology Society.)

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Background: Local assembly with short and long reads has proven to be very useful in many applications: reconstruction of the sequence of a locus of interest, gap-filling in draft assemblies, as well as alternative allele reconstruction of large Structural Variants. Whereas linked-read technologies have a great potential to assemble specific loci as they provide long-range information while maintaining the power and accuracy of short-read sequencing, there is a lack of local assembly tools for linked-read data., Results: We present MTG-Link, a novel local assembly tool dedicated to linked-reads. The originality of the method lies in its read subsampling step which takes advantage of the barcode information contained in linked-reads mapped in flanking regions. We validated our approach on several datasets from different linked-read technologies. We show that MTG-Link is able to assemble successfully large sequences, up to dozens of Kb. We also demonstrate that the read subsampling step of MTG-Link considerably improves the local assembly of specific loci compared to other existing short-read local assembly tools. Furthermore, MTG-Link was able to fully characterize large insertion variants and deletion breakpoints in a human genome and to reconstruct dark regions in clinically-relevant human genes. It also improved the contiguity of a 1.3 Mb locus of biological interest in several individual genomes of the mimetic butterfly Heliconius numata., Conclusions: MTG-Link is an efficient local assembly tool designed for different linked-read sequencing technologies. MTG-Link source code is available at https://github.com/anne-gcd/MTG-Link and as a Bioconda package., (© 2023. The Author(s).)

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies., (© 2023 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The role of lncRNAs in immune defence has been demonstrated in many multicellular and unicellular organisms. However, investigation of the identification and characterization of long non-coding RNAs (lncRNAs) involved in the insect immune response is still limited. In this study, we used RNA sequencing (RNA-seq) to investigate the expression profiles of lncRNAs and mRNAs in the fall armyworm Spodoptera frugiperda in response to virus infection. To assess the tissue- and virus-specificity of lncRNAs, we analysed and compared their expression profiles in haemocytes and fat body of larvae infected with two entomopathogenic viruses with different lifestyles, i.e. the polydnavirus HdIV ( Hyposoter didymator IchnoVirus) and the densovirus JcDV ( Junonia coenia densovirus). We identified 1883 candidate lncRNAs, of which 529 showed differential expression following viral infection. Expression profiles differed considerably between samples, indicating that many differentially expressed (DE) lncRNAs showed virus- and tissue-specific expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and target prediction analyses indicated that DE-LncRNAs were mainly enriched in metabolic process, DNA replication and repair, immune response, metabolism of insect hormone and cell adhesion. In addition, we identified three DE-lncRNAs potentially acting as microRNA host genes, suggesting that they participate in gene regulation by producing miRNAs in response to virus infection. This study provides a catalogue of lncRNAs expressed in two important immune tissues and potential insight into their roles in the antiviral defence in S. frugiperda . The results may help future in-depth functional studies to better understand the biological function of lncRNAs in interaction between viruses and the fall armyworm.

The fall armyworm (FAW; Spodoptera frugiperda) is one of the major agricultural pest insects. FAW is native to the Americas, and its invasion was first reported in West Africa in 2016. Then it quickly spread through Africa, Asia, and Oceania, becoming one of the main threats to corn production. We analyzed whole genome sequences of 177 FAW individuals from 12 locations on four continents to infer evolutionary processes of invasion. Principal component analysis from the TPI gene and whole genome sequences shows that invasive FAW populations originated from the corn strain. Ancestry coefficient and phylogenetic analyses from the nuclear genome indicate that invasive populations are derived from a single ancestry, distinct from native populations, while the mitochondrial phylogenetic tree supports the hypothesis of multiple introductions. Adaptive evolution specific to invasive populations was observed in detoxification, chemosensory, and digestion genes. We concluded that extant invasive FAW populations originated from the corn strain with potential contributions of adaptive evolution., (© 2022. The Author(s).)

Background: Divergent selection on host-plants is one of the main evolutionary forces driving ecological speciation in phytophagous insects. The ecological speciation might be challenging in the presence of gene flow and assortative mating because the direction of divergence is not necessarily the same between ecological selection (through host-plant adaptation) and assortative mating. The fall armyworm (FAW), a major lepidopteran pest species, is composed of two sympatric strains, corn and rice strains, named after two of their preferred host-plants. These two strains have been hypothesized to undergo incipient speciation, based on (i) several lines of evidence encompassing both pre- and post-zygotic reproductive isolation, and (ii) the presence of a substantial level of genetic differentiation. Even though the status of these two strains has been established a long time ago, it is still yet to be found whether these two strains indeed exhibit a marked level of genetic differentiation from a large number of genomic loci. Here, we analyzed whole genome sequences from 56 FAW individuals either collected from pasture grasses (a part of the favored host range of the rice strain) or corn to assess the role of host-plant adaptation in incipient speciation., Results: Principal component analysis of whole genome data shows that the pattern of divergence in the fall armyworm is predominantly explained by the genetic differentiation associated with host-plants. The level of genetic differentiation between corn and rice strains is particularly marked in the Z chromosome. We identified one autosomal locus and two Z chromosome loci targeted by selective sweeps specific to rice strain and corn strain, respectively. The autosomal locus has both increased D XY and F ST while the Z chromosome loci had decreased D XY and increased F ST ., Conclusion: These results show that the FAW population structure is dominated by the genetic differentiation between corn and rice strains. This differentiation involves divergent selection targeting at least three loci, which include a locus potentially causing reproductive isolation. Taken together, these results suggest the evolutionary scenario that host-plant speciation is a driver of incipient speciation in the fall armyworm., (© 2022. The Author(s).)

The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Nitrogen fertilization has been reported to influence the development of clubroot, a root disease of Brassicaceae species, caused by the obligate protist Plasmodiophora brassicae . Our previous works highlighted that low-nitrogen fertilization induced a strong reduction of clubroot symptoms in some oilseed rape genotypes. To further understand the underlying mechanisms, the response to P. brassicae infection was investigated in two genotypes "Yudal" and HD018 harboring sharply contrasted nitrogen-driven modulation of resistance toward P. brassicae . Targeted hormone and metabolic profiling, as well as RNA-seq analysis, were performed in inoculated and non-inoculated roots at 14 and 27 days post-inoculation, under high and low-nitrogen conditions. Clubroot infection triggered a large increase of SA concentration and an induction of the SA gene markers expression whatever the genotype and nitrogen conditions. Overall, metabolic profiles suggested that N-driven induction of resistance was independent of SA signaling, soluble carbohydrate and amino acid concentrations. Low-nitrogen-driven resistance in "Yudal" was associated with the transcriptional regulation of a small set of genes, among which the induction of NRT2 - and NR -encoding genes. Altogether, our results indicate a possible role of nitrate transporters and auxin signaling in the crosstalk between plant nutrition and partial resistance to pathogens., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Aigu, Daval, Gazengel, Marnet, Lariagon, Laperche, Legeai, Manzanares-Dauleux and Gravot.)

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera., (© 2021 European Society for Evolutionary Biology.)

Motivation: Linked-Reads technologies combine both the high quality and low cost of short-reads sequencing and long-range information, through the use of barcodes tagging reads which originate from a common long DNA molecule. This technology has been employed in a broad range of applications including genome assembly, phasing and scaffolding, as well as structural variant calling. However, to date, no tool or API dedicated to the manipulation of Linked-Reads data exist., Results: We introduce LRez, a C++ API and toolkit that allows easy management of Linked-Reads data. LRez includes various functionalities, for computing numbers of common barcodes between genomic regions, extracting barcodes from BAM files, as well as indexing and querying BAM, FASTQ and gzipped FASTQ files to quickly fetch all reads or alignments containing a given barcode. LRez is compatible with a wide range of Linked-Reads sequencing technologies, and can thus be used in any tool or pipeline requiring barcode processing or indexing, in order to improve their performances., Availability and Implementation: LRez is implemented in C++, supported on Unix-based platforms and available under AGPL-3.0 License at https://github.com/morispi/LRez, and as a bioconda module., Supplementary Information: Supplementary data are available at Bioinformatics Advances online., (© The Author(s) 2021. Published by Oxford University Press.)

Oilseed rape ( Brassica napus L.) is the third largest oil crop worldwide. Like other crops, oilseed rape faces unfavorable environmental conditions resulting from multiple and combined actions of abiotic and biotic constraints that occur throughout the growing season. In particular drought severely reduces seed yield but also impacts seed quality in oilseed rape. In addition, clubroot disease, caused by the pathogen Plasmodiophora brassicae , limits the yield of the oilseed rape crops grown in infected areas. Clubroot induces swellings or galls on the roots that decrease the flow of water and nutrients within the plant. Furthermore, combinations of different stresses lead to complex plant responses that can not be predicted by the simple addition of individual stress responses. Indeed, an abiotic constraint can either reduce or stimulate the plant response to a pathogen or pest. Transcriptome datasets from different conditions are key resources to improve our knowledge of environmental stress-resistance mechanisms in plant organs. Here, we describe a RNA-seq dataset consisting of 72 samples of immature B. napus seeds from plants grown either under drought, infected with P. brassicae , or a combination of both stresses. A total of 67.6 Gb of transcriptome paired-end reads were filtered, mapped onto the B. napus reference genome Darmor- bzh and used for identification of differentially expressed genes and gene ontology enrichment. The raw reads are available under accession PRJNA738318 at NCBI Sequence Read Archive (SRA) repository. The dataset is a resource for the scientific community exploring seed plasticity., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships which have or could be perceived to have influenced the work reported in this article., (© 2021 The Author(s). Published by Elsevier Inc.)