One's trash is someone else's treasure: sequence read archives from Lepidoptera genomes provide material for genome reconstruction of their endosymbionts.

BACKGROUND: Maternally inherited bacterial symbionts are extremely widespread in insects. They owe their success to their ability to promote their own transmission through various manipulations of their hosts' life-histories. Many symbionts however very often go undetected. Consequently, we have only a restricted idea of the true symbiont diversity in insects, which may hinder our understanding of even bigger questions in the field such as the evolution or establishment of symbiosis. RESULTS: In this study, we screened publicly available Lepidoptera genomic material for two of the most common insect endosymbionts, namely Wolbachia and Spiroplasma, in 1904 entries, encompassing 106 distinct species. We compared the performance of two screening software, Kraken2 and MetaPhlAn2, to identify the bacterial infections and using a baiting approach we reconstruct endosymbiont genome assemblies. Of the 106 species screened, 20 (19%) and nine (8.5%) were found to be infected with either Wolbachia or Spiroplasma, respectively. Construction of partial symbiotic genomes and phylogenetic analyses suggested the Wolbachia strains from the supergroup B were the most prevalent type of symbionts, while Spiroplasma infections were scarce in the Lepidoptera species screened here. CONCLUSIONS: Our results indicate that many of the host-symbiont associations remain largely unexplored, with the majority of associations we identify never being recorded before. This highlights the usefulness of public databases to explore the hidden diversity of symbiotic entities, allowing the development of hypotheses regarding host-symbiont associations. The ever-expanding genomic databases provide a diverse databank from which one can characterize and explore the true diversity of symbiotic entities.

Pediococcus pentosaceus IMI 507025 genome sequencing data.

The genome sequence data for the pickled cucumbers isolate, Pediococcus pentosaceus IMI 507025, is reported. The raw reads and analysed genome reads were deposited at NCBI under Bioproject with the accession number PRJNA814992. The number of contigs before and after trimming were 17 and 12 contigs, respectively. The total size of the genome was 1,795,439 bp containing 1,811 total genes, of which 1,751 were coding sequences. IMI 507025 identity was determined via average nucleotide identity (ANI), obtaining an identity value of 99.5994% between IMI 507025 and the type strain P. pentosaceus ATCC 33316, identifying the strain as P. pentosaceus. Screening for the antimicrobial resistance (AMR) and virulence genes in the genome of IMI 507025 showed no hits, confirming the safety of the tested strain. Presence of plasmids was not found.

Genome sequence data of Lactiplantibacillus plantarum IMI 507028.

Here, we report the genome sequencing data for the fermented milk isolate, Lactiplantibacillus plantarum (L. plantarum) IMI 507028. The Bioproject, SRA, and GenBank data were deposited at NCBI under accession numbers PRJNA801616, SRR18323693, and JAKMAX000000000, respectively. The size of the genome was 3,231,321 bp, with a GC% of 44.52. Before sequence trimming, the genome contained 40 contigs, in which 35 contigs were annotated, revealing 2937 coding sequences out of 3052 total genes. The strain was identified as L. plantarum with an average nucleotide identity (ANI) value of 99.9922% between IMI 507028 and L. plantarum JDM1. Genes related to antimicrobial resistance or pathogenic factors were not found during screening.

Draft Genome Sequence of Lactiplantibacillus plantarum IMI 507026.

Here, we announce the draft genome sequence of Lactiplantibacillus plantarum isolated from corn silage in Nicholasville, KY. L. plantarum IMI 507026 is deposited in the Centre for Agriculture and Bioscience International (CABI) Culture Collection with the accession number IMI 507026.

Draft Genome Sequence of Lacticaseibacillus rhamnosus IMI 507023.

We report here the draft genome sequence of Lacticaseibacillus rhamnosus strain IMI 507023, a lactic acid bacterium, isolated from corn silage in Nicholasville, Kentucky, USA. The strain is deposited in the Centre for Agriculture and Bioscience International (CABI) Culture Collection with the accession number IMI 507023.

Whole genome sequence data of Lactiplantibacillus plantarum IMI 507027.

Here we report the draft genome sequence of the Lactiplantibacillus plantarum IMI 507027 strain. The genome consists of 37 contigs with a total size of 3,235,614 bp and a GC% of 44.51. After sequence trimming, 31 contigs were annotated, revealing 3,126 genes, of which 3,030 were coding sequences. The Average Nucleotide Identity (ANI) gave a value of 99.9926% between IMI 507027 and L. plantarum JDM1, identifying the strain as L. plantarum. No genes of concern for safety-related traits such as antimicrobial resistance or virulence factors were found. The annotated genome and raw sequence reads were deposited at NCBI under Bioproject with the accession number PRJNA791753.

Draft Genome Sequence of Pediococcus pentosaceus IMI 507024.

We report here the draft genome sequence of Pediococcus pentosaceus strain IMI 507024, a lactic acid bacterium isolated from fermented sausage in Kentucky (Nicholasville, KY, USA). The strain is deposited in the Centre for Agriculture and Bioscience International (CABI) Culture Collection with the accession number IMI 507024.

Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map.

BACKGROUND: The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses. RESULTS: The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models. CONCLUSIONS: The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

Expression of the GAF Sensor, Carbohydrate-Active Enzymes, Elicitins, and RXLRs Differs Markedly Between Two Phytophthora cactorum Isolates.

The phytopathogen Phytophthora cactorum infects economically important herbaceous and woody plant species. P. cactorum isolates differ in host specificity; for example, strawberry crown rot is often caused by a specialized pathotype. Here we compared the transcriptomes of two P. cactorum isolates that differ in their virulence to garden strawberry (Pc407: high virulence; Pc440: low virulence). De novo transcriptome assembly and clustering of contigs resulted in 19,372 gene clusters. Two days after inoculation of Fragaria vesca roots, 3,995 genes were differently expressed between the P. cactorum isolates. One of the genes that were highly expressed only in Pc407 encodes a GAF sensor protein potentially involved in membrane trafficking processes. Two days after inoculation, elicitins were highly expressed in Pc407 and lipid catabolism appeared to be more active than in Pc440. Of the carbohydrate-active enzymes, those that degrade pectin were often more highly expressed in Pc440, whereas members of glycosyl hydrolase family 1, potentially involved in the metabolism of glycosylated secondary metabolites, were more highly expressed in Pc407 at the time point studied. Differences were also observed among the RXLR effectors: Pc407 appears to rely on a smaller set of key RXLR effectors, whereas Pc440 expresses a greater number of RXLRs. This study is the first step toward improving understanding of the molecular basis of differences in the virulence of P. cactorum isolates. Identification of the key effectors is important, as it enables effector-assisted breeding strategies toward crown rot-resistant strawberry cultivars.

Transcriptional effects of cadmium on iron homeostasis differ in calamine accessions of Noccaea caerulescens.

Calamine accessions of the zinc/cadmium/nickel hyperaccumulator, Noccaea caerulescens, exhibit striking variation in foliar cadmium accumulation in nature. The Ganges accession (GA) from Southern France displays foliar cadmium hyperaccumulation (>1000 µg g-1 DW), whereas the accession La Calamine (LC) from Belgium, with similar local soil metal composition, does not (<100 µg g-1 DW). All calamine accessions are cadmium hypertolerant. To find out the differences between LC and GA in their basic adaptation mechanisms, we bypassed the cadmium excluding phenotype of LC by exposing the plants to 50 µm cadmium in hydroponics, achieving equal cadmium accumulation in the shoots. The iron content increased in the roots of both accessions. GA exhibited significant decreases in manganese and zinc contents in the roots and shoots, approaching those in LC. Altogether 702 genes responded differently to cadmium exposure between the accessions, 157 and 545 in the roots and shoots, respectively. Cadmium-exposed LC showed a stress response and had decreased levels of a wide range of photosynthesis-related transcripts. GA showed less changes, mainly exhibiting an iron deficiency-like response. This included increased expression of genes encoding five iron deficiency-regulated bHLH transcription factors, ferric reduction oxidase FRO2, iron transporters IRT1 and OPT3, and nicotianamine synthase NAS1, and decreased expression of genes encoding ferritins and NEET (a NEET family iron-sulfur protein), which is possibly involved in iron transfer, distribution and/or management. The function of the IRT1 gene in the accessions was compared. We conclude that the major difference between the two accessions is in the way they cope with iron under cadmium exposure.

Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch.

Silver birch (Betula pendula) is a pioneer boreal tree that can be induced to flower within 1 year. Its rapid life cycle, small (440-Mb) genome, and advanced germplasm resources make birch an attractive model for forest biotechnology. We assembled and chromosomally anchored the nuclear genome of an inbred B. pendula individual. Gene duplicates from the paleohexaploid event were enriched for transcriptional regulation, whereas tandem duplicates were overrepresented by environmental responses. Population resequencing of 80 individuals showed effective population size crashes at major points of climatic upheaval. Selective sweeps were enriched among polyploid duplicates encoding key developmental and physiological triggering functions, suggesting that local adaptation has tuned the timing of and cross-talk between fundamental plant processes. Variation around the tightly-linked light response genes PHYC and FRS10 correlated with latitude and longitude and temperature, and with precipitation for PHYC. Similar associations characterized the growth-promoting cytokinin response regulator ARR1, and the wood development genes KAK and MED5A.

De novo transcriptome assemblies of four accessions of the metal hyperaccumulator plant Noccaea caerulescens.

Noccaea caerulescens of the Brassicaceae family has become the key model plant among the metal hyperaccumulator plants. Populations/accessions of N. caerulescens from geographic locations with different soil metal concentrations differ in their ability to hyperaccumulate and hypertolerate metals. Comparison of transcriptomes in several accessions provides candidates for detailed exploration of the mechanisms of metal accumulation and tolerance and local adaptation. This can have implications in the development of plants for phytoremediation and improved mineral nutrition. Transcriptomes from root and shoot tissues of four N. caerulescens accessions with contrasting Zn, Cd and Ni hyperaccumulation and tolerance traits were sequenced with Illumina Hiseq2000. Transcriptomes were assembled using the Trinity de novo assembler and were annotated and the protein sequences predicted. The comparison against the BUSCO plant early release dataset indicated high-quality assemblies. The predicted protein sequences have been clustered into ortholog groups with closely related species. The data serve as important reference sequences in whole transcriptome studies, in analyses of genetic differences between the accessions and other species, and for primer design.

Reprogramming of Strawberry (Fragaria vesca) Root Transcriptome in Response to Phytophthora cactorum.

Crown rot (Phytophthora cactorum) causes significant economic losses in strawberry production. The best control strategy would be to use resistant cultivars, but polygenically inherited resistance makes the breeding of the garden strawberry (Fragaria × ananassa) challenging. The diploid wild strawberry Fragaria vesca Hawaii 4 genotype was shown previously to have resistance against crown rot. To explore the resistance mechanisms, we inoculated the roots of Hawaii 4 with P. cactorum in a novel in vitro hydroponic system to minimize interference caused by other microbes. Major reprogramming of the root transcriptome occurred, involving 30% of the genes. The surveillance system of the plant shifted from the development mode to the defense mode. Furthermore, the immune responses as well as many genes involved in the biosynthesis of the defense hormones jasmonic acid, ethylene and salicylic acid were up-regulated. Several major allergen-like genes encoding PR-10 proteins were highly expressed in the inoculated plants, suggesting that they also have a crucial role in the defense responses against P. cactorum. Additionally, flavonoids and terpenoids may be of vital importance, as several genes involved in their biosynthesis were up-regulated. The cell wall biosynthesis and developmental processes were down-regulated, possibly as a result of the down-regulation of the key genes involved in the biosynthesis of growth-promoting hormones brassinosteroids and auxin. Of particular interest was the expression of potential resistance genes in the recently identified P. cactorum resistance locus RPc-1. These new findings help to target the breeding efforts aiming at more resistant strawberry cultivars.

De Novo assembly and annotation of the freshwater crayfish Astacus astacus transcriptome.

We generated RNA-seq data to assemble the transcriptome of the noble crayfish (Astacus astacus) from four combined tissues (abdominal muscle, hepatopancreas, ovaries, green glands). A total of 194 million read pairs with a length of 100 bp were generated. The transcriptome was assembled de novo using Trinity software, producing 158,649 non-redundant transcripts. Lowly expressed transcripts were filtered out leaving 45,415 transcripts of which 14,559 were found to contain open reading frames with predicted gene function. The Transrate software revealed that 91% of the total reads were realigned to the assembly. Furthermore, BUSCO analysis indicated that our assembly is 64% complete. A total of 13,770 transcripts were assigned at least one GO term. This first de novo transcriptome assembly is an important foundation for future genomic research on the noble crayfish and adds to the general knowledge and further characterization of transcriptomes of non-model organisms.

Comparative transcriptome analysis of the metal hyperaccumulator Noccaea caerulescens.

The metal hyperaccumulator Noccaea caerulescens is an established model to study the adaptation of plants to metalliferous soils. Various comparators have been used in these studies. The choice of suitable comparators is important and depends on the hypothesis to be tested and methods to be used. In high-throughput analyses such as microarray, N. caerulescens has been compared to non-tolerant, non-accumulator plants like Arabidopsis thaliana or Thlaspi arvense rather than to the related hypertolerant or hyperaccumulator plants. An underutilized source is N. caerulescens populations with considerable variation in their capacity to accumulate and tolerate metals. Whole transcriptome sequencing (RNA-Seq) is revealing interesting variation in their gene expression profiles. Combining physiological characteristics of N. caerulescens accessions with their RNA-Seq has a great potential to provide detailed insight into the underlying molecular mechanisms, including entirely new gene products. In this review we will critically consider comparative transcriptome analyses carried out to explore metal hyperaccumulation and hypertolerance of N. caerulescens, and demonstrate the potential of RNA-Seq analysis as a tool in evolutionary genomics.

Gene expression differences between Noccaea caerulescens ecotypes help to identify candidate genes for metal phytoremediation.

Populations of Noccaea caerulescens show tremendous differences in their capacity to hyperaccumulate and hypertolerate metals. To explore the differences that could contribute to these traits, we undertook SOLiD high-throughput sequencing of the root transcriptomes of three phenotypically well-characterized N. caerulescens accessions, i.e., Ganges, La Calamine, and Monte Prinzera. Genes with possible contribution to zinc, cadmium, and nickel hyperaccumulation and hypertolerance were predicted. The most significant differences between the accessions were related to metal ion (di-, trivalent inorganic cation) transmembrane transporter activity, iron and calcium ion binding, (inorganic) anion transmembrane transporter activity, and antioxidant activity. Analysis of correlation between the expression profile of each gene and the metal-related characteristics of the accessions disclosed both previously characterized (HMA4, HMA3) and new candidate genes (e.g., for nickel IRT1, ZIP10, and PDF2.3) as possible contributors to the hyperaccumulation/tolerance phenotype. A number of unknown Noccaea-specific transcripts also showed correlation with Zn(2+), Cd(2+), or Ni(2+) hyperaccumulation/tolerance. This study shows that N. caerulescens populations have evolved great diversity in the expression of metal-related genes, facilitating adaptation to various metalliferous soils. The information will be helpful in the development of improved plants for metal phytoremediation.