Transcription analysis of neonicotinoid resistance in Mediterranean (MED) populations of B. tabaci reveal novel cytochrome P450s, but no nAChR mutations associated with the phenotype.

BACKGROUND: Bemisia tabaci is one of the most damaging agricultural pests world-wide. Although its control is based on insecticides, B. tabaci has developed resistance against almost all classes of insecticides, including neonicotinoids. RESULTS: We employed an RNA-seq approach to generate genome wide expression data and identify genes associated with neonicotinoid resistance in Mediterranean (MED) B. tabaci (Q1 biotype). Twelve libraries from insecticide resistant and susceptible whitefly populations were sequenced on an Illumina Next-generation sequencing platform, and genomic sequence information of approximately 73 Gbp was generated. A reference transcriptome was built by de novo assembly and functionally annotated. A total of 146 P450s, 18 GSTs and 23 CCEs enzymes (unigenes) potentially involved in the detoxification of xenobiotics were identified, along with 78 contigs encoding putative target proteins of six different insecticide classes. Ten unigenes encoding nicotinic Acetylcholine Receptors (nAChR), the target of neoinicotinoids, were identified and phylogenetically classified. No nAChR polymorphism potentially related with the resistant phenotypes, was observed among the studied strains. DE analysis revealed that among the 550 differentially (logFC > 1) over-transcribed unigenes, 52 detoxification enzymes were over expressed including unigenes with orthologues in P450s, GSTs, CCE and UDP-glucuronosyltransferases. Eight P450 unigenes belonging to clades CYP2, CYP3 and CYP4 were highly up-regulated (logFC > 2) including CYP6CM1, a gene already known to confer imidacloprid resistance in B. tabaci. Using quantitative qPCRs, a larger screening of field MED B. tabaci from Crete with known neonicotinoid phenotype was performed to associate expression levels of P450s with resistance levels. Expression levels of five P450s, including CYP6CM1, were found associated with neonicotinoid resistance. However, a significant correlation was found only in CYP303 and CYP6CX3, with imidacloprid and acetamiprid respectively. CONCLUSION: Our work has generated new toxicological data and genomic resources which will significantly enrich the available dataset and substantially facilitate the molecular studies in MED B. tabaci. No evidence of target site neonicotinoid resistance has been found. Eight P450 unigenes, including CYP6CM1, were found significantly over-expressed in resistant B. tabaci. This study suggests at least two novel P450s (CYP303 and CYP6CX3) as candidates for their functional characterization as detoxification mechanisms of neonicotinoid resistance in B. tabaci.

The sex-specific transcriptome of the hermaphrodite sparid sharpsnout seabream (Diplodus puntazzo).

BACKGROUND: Teleosts are characterized by a remarkable breadth of sexual mechanisms including various forms of hermaphroditism. Sparidae is a fish family exhibiting gonochorism or hermaphroditism even in closely related species. The sparid Diplodus puntazzo (sharpsnout seabream), exhibits rudimentary hermaphroditism characterized by intersexual immature gonads but single-sex mature ones. Apart from the intriguing reproductive biology, it is economically important with a continuously growing aquaculture in the Mediterranean Sea, but limited available genetic resources. Our aim was to characterize the expressed transcriptome of gonads and brains through RNA-Sequencing and explore the properties of genes that exhibit sex-biased expression profiles. RESULTS: Through RNA-Sequencing we obtained an assembled transcriptome of 82,331 loci. The expression analysis uncovered remarkable differences between male and female gonads, while male and female brains were almost identical. Focused search for known targets of sex determination and differentiation in vertebrates built the sex-specific expression profile of sharpsnout seabream. Finally, a thorough genetic marker discovery pipeline led to the retrieval of 85,189 SNPs and 29,076 microsatellites enriching the available genetic markers for this species. CONCLUSIONS: We obtained a nearly complete source of transcriptomic sequence as well as marker information for sharpsnout seabream, laying the ground for understanding the complex process of sex differentiation of this economically valuable species. The genes involved include known candidates from other vertebrate species, suggesting a conservation of the toolkit between gonochorists and hermaphrodites.

Microbial diversity in four Mediterranean irciniid sponges.

This paper describes a dataset of microbial communities from four different sponge species: Irciniaoros (Schmidt, 1864), Irciniavariabilis (Schmidt, 1862), Sarcotragusspinosulus Schmidt, 1862 and Sarcotragusfasciculatus (Pallas, 1766). The examined sponges all belong to Demospongiae (Class); Keratosa (Subclass); Dictyoceratida (Order); Irciniidae (Family). Samples were collected by scuba diving at depths between 6-14 m from two sampling sites of rocky formations at the northern coast of Crete (Cretan Sea, eastern Mediterranean) and were subjected to metabarcoding for the V5-V6 region of the 16S rRNA gene.

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.

SnakeCube: containerized and automated pipeline for de novo genome assembly in HPC environments.

OBJECTIVE: The rapid progress in sequencing technology and related bioinformatics tools aims at disentangling diversity and conservation issues through genome analyses. The foremost challenges of the field involve coping with questions emerging from the swift development and application of new algorithms, as well as the establishment of standardized analysis approaches that promote transparency and transferability in research. RESULTS: Here, we present SnakeCube, an automated and containerized whole de novo genome assembly pipeline that runs within isolated, secured environments and scales for use in High Performance Computing (HPC) domains. SnakeCube was optimized for its performance and tested for its effectiveness with various inputs, highlighting its safe and robust universal use in the field.

The Rm1 and Rm2 Resistance Genes to Green Peach Aphid (Myzus persicae) Encode the Same TNL Proteins in Peach (Prunus persica L.).

The green peach aphid (GPA), Myzus persicae, is an important pest of the peach crop. Three major dominant resistance genes have already been detected, Rm1 in the Weeping Flower Peach (WFP) clone, Rm2 in the Rubira clone, and Rm3 in the Fen Shouxing clone. In this study, after NGS resequencing of WFP and Rubira, we found that their genomic sequences in the Rm1 and Rm2 region were similar but very different from that of the susceptible reference peach Lovell. We constructed a BAC library for the GPA-resistant WFP and screened four BAC clones to sequence the target region. The new sequence was 61.7 Kb longer than Lovell and was annotated with four different TIR_NBS_LRR genes. Among them, the TNL1 gene was very overexpressed in WFP leaves 24 h after GPA infestation. This gene was also present and expressed in the Rubira clone and had the same sequence as the candidate Rm3 gene, supporting the hypothesis that the three genes share the same origin. In addition, we identified a second TNL, TNL2, located at 35.4 Kb from TNL1 and slightly overexpressed after GPA infestation. Kasp and size molecular markers were designed for use in marker-assisted selection and were validated in a peach segregating population.

Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits.

Most molecularly characterized plant resistance genes (R genes) belong to the nucleotide-binding-site-leucine-rich-repeat (NLR) receptor family and are prone to duplication and transposition with high sequence diversity. In this family, the Vat gene in melon is one of the few R genes known for conferring resistance to insect, i.e., Aphis gossypii, but it has been misassembled and/or mispredicted in the whole genomes of Cucurbits. We examined 14 genomic regions (about 400 kb) derived from long-read assemblies spanning Vat-related genes in Cucumis melo, Cucumis sativus, Citrullus lanatus, Benincasa hispida, Cucurbita argyrosperma, and Momordica charantia. We built the phylogeny of those genes. Investigating the paleohistory of the Vat gene cluster, we revealed a step by step process beginning from a common ancestry in cucurbits older than 50 my. We highlighted Vat exclusively in the Cucumis genera, which diverged about 20 my ago. We then focused on melon, evaluating a minimum duplication rate of Vat in 80 wild and cultivated melon lines using generalist primers; our results suggested that duplication started before melon domestication. The phylogeny of 44 Vat-CDS obtained from 21 melon lines revealed gain and loss of leucine-rich-repeat domains along diversification. Altogether, we revealed the high putative recognition scale offered in melon based on a combination of SNPs, number of leucine-rich-repeat domains within each homolog and number of homologs within each cluster that might jointly confer resistance to a large pest and pathogen spectrum. Based on our findings, we propose possible avenues for breeding programs.

0s and 1s in marine molecular research: a regional HPC perspective.

High-performance computing (HPC) systems have become indispensable for modern marine research, providing support to an increasing number and diversity of users. Pairing with the impetus offered by high-throughput methods to key areas such as non-model organism studies, their operation continuously evolves to meet the corresponding computational challenges. Here, we present a Tier 2 (regional) HPC facility, operating for over a decade at the Institute of Marine Biology, Biotechnology, and Aquaculture of the Hellenic Centre for Marine Research in Greece. Strategic choices made in design and upgrades aimed to strike a balance between depth (the need for a few high-memory nodes) and breadth (a number of slimmer nodes), as dictated by the idiosyncrasy of the supported research. Qualitative computational requirement analysis of the latter revealed the diversity of marine fields, methods, and approaches adopted to translate data into knowledge. In addition, hardware and software architectures, usage statistics, policy, and user management aspects of the facility are presented. Drawing upon the last decade's experience from the different levels of operation of the Institute of Marine Biology, Biotechnology, and Aquaculture HPC facility, a number of lessons are presented; these have contributed to the facility's future directions in light of emerging distribution technologies (e.g., containers) and Research Infrastructure evolution. In combination with detailed knowledge of the facility usage and its upcoming upgrade, future collaborations in marine research and beyond are envisioned.

NOBLAST and JAMBLAST: New Options for BLAST and a Java Application Manager for BLAST results.

UNLABELLED: NOBLAST (New Options for BLAST) is an open source program that provides a new user-friendly tabular output format for various NCBI BLAST programs (Blastn, Blastp, Blastx, Tblastn, Tblastx, Mega BLAST and Psi BLAST) without any use of a parser and provides E-value correction in case of use of segmented BLAST database. JAMBLAST using the NOBLAST output allows the user to manage, view and filter the BLAST hits using a number of selection criteria. AVAILABILITY: A distribution package of NOBLAST and JAMBLAST including detailed installation procedure is freely available from http://sourceforge.net/projects/JAMBLAST/ and http://sourceforge.net/projects/NOBLAST. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Unravelling paralogous gene expression dynamics during three-spined stickleback embryogenesis.

Development requires the implementation of a plethora of molecular mechanisms, involving a large set of genes to ensure proper cell differentiation, morphogenesis of tissues and organs as well as the growth of the organism. Genome duplication and resulting paralogs are considered to provide the raw genetic materials important for new adaptation opportunities and boosting evolutionary innovation. The present study investigated paralogous genes, involved in three-spined stickleback (Gasterosteus aculeatus) development. Therefore, the transcriptomes of five early stages comprising developmental leaps were explored. Obtained expression profiles reflected the embryo's needs at different stages. Early stages, such as the morula stage comprised transcripts mainly involved in energy requirements while later stages were mostly associated with GO terms relevant to organ development and morphogenesis. The generated transcriptome profiles were further explored for differential expression of known and new paralogous genes. Special attention was given to hox genes, with hoxa13a being of particular interest and to pigmentation genes where itgb1, involved in the melanophore development, displayed a complementary expression pattern throughout studied stages. Knowledge obtained by untangling specific paralogous gene functions during development might not only significantly contribute to the understanding of teleost ontogenesis but might also shed light on paralogous gene evolution.

A de novo transcriptome assembly for the bath sponge Spongia officinalis, adjusting for microsymbionts.

OBJECTIVES: We report a transcriptome acquisition for the bath sponge Spongia officinalis, a non-model marine organism that hosts rich symbiotic microbial communities. To this end, a pipeline was developed to efficiently separate between bacterial expressed genes from those of eukaryotic origin. The transcriptome was produced to support the assessment of gene expression and, thus, the response of the sponge, to elevated temperatures, replicating conditions currently occurring in its native habitat. DATA DESCRIPTION: We describe the assembled transcriptome along with the bioinformatic pipeline used to discriminate between signals of metazoan and prokaryotic origin. The pipeline involves standard read pre-processing steps and incorporates extra analyses to identify and filter prokaryotic reads out of the analysis. The proposed pipeline can be followed to overcome the technical RNASeq problems characteristic for symbiont-rich metazoan organisms with low or non-existent tissue differentiation, such as sponges and cnidarians. At the same time, it can be valuable towards the development of approaches for parallel transcriptomic studies of symbiotic communities and the host.

The Gene Toolkit Implicated in Functional Sex in Sparidae Hermaphrodites: Inferences From Comparative Transcriptomics.

Sex-biased gene expression is the mode through which sex dimorphism arises from a nearly identical genome, especially in organisms without genetic sex determination. Teleost fishes show great variations in the way the sex phenotype forms. Among them, Sparidae, that might be considered as a model family displays a remarkable diversity of reproductive modes. In this study, we sequenced and analyzed the sex-biased transcriptome in gonads and brain (the tissues with the most profound role in sexual development and reproduction) of two sparids with different reproductive modes: the gonochoristic common dentex, Dentex dentex, and the protandrous hermaphrodite gilthead seabream, Sparus aurata. Through comparative analysis with other protogynous and rudimentary protandrous sparid transcriptomes already available, we put forward common male and female-specific genes and pathways that are probably implicated in sex-maintenance in this fish family. Our results contribute to the understanding of the complex processes behind the establishment of the functional sex, especially in hermaphrodite species and set the groundwork for future experiments by providing a gene toolkit that can improve efforts to control phenotypic sex in finfish in the ever-increasingly important field of aquaculture.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird.

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.

Transcriptome Profiling and Genetic Study Reveal Amplified Carboxylesterase Genes Implicated in Temephos Resistance, in the Asian Tiger Mosquito Aedes albopictus.

BACKGROUND: The control of Aedes albopictus, a major vector for viral diseases, such as dengue fever and chikungunya, has been largely reliant on the use of the larvicide temephos for many decades. This insecticide remains a primary control tool for several countries and it is a potential reliable reserve, for emergency epidemics or new invasion cases, in regions such as Europe which have banned its use. Resistance to temephos has been detected in some regions, but the mechanism responsible for the trait has not been investigated. PRINCIPAL FINDINGS: Temephos resistance was identified in an Aedes albopictus population isolated from Greece, and subsequently selected in the laboratory for a few generations. Biochemical assays suggested the association of elevated carboxylesterases (CCE), but not target site resistance (altered AChE), with this phenotype. Illumina transcriptomic analysis revealed the up-regulation of three transcripts encoding CCE genes in the temephos resistant strain. CCEae3a and CCEae6a showed the most striking up-regulation (27- and 12-folds respectively, compared to the reference susceptible strain); these genes have been previously shown to be involved in temephos resistance also in Ae. aegypti. Gene amplification was associated with elevated transcription levels of both CCEae6a and CCEae3a genes. Genetic crosses confirmed the genetic link between CCEae6a and CCEae3a amplification and temephos resistance, by demonstrating a strong association between survival to temephos exposure and gene copy numbers in the F2 generation. Other transcripts, encoding cytochrome P450s, UDP-glycosyltransferases (UGTs), cuticle and lipid biosynthesis proteins, were upregulated in resistant mosquitoes, indicating that the co-evolution of multiple mechanisms might contribute to resistance. SIGNIFICANCE: The identification of specific genes associated with insecticide resistance in Ae. albopictus for the first time is an important pre-requirement for insecticide resistance management. The genomic resources that were produced will be useful to the community, to study relevant aspects of Ae. albopictus biology.

Exploring a Nonmodel Teleost Genome Through RAD Sequencing-Linkage Mapping in Common Pandora, Pagellus erythrinus and Comparative Genomic Analysis.

Common pandora (Pagellus erythrinus) is a benthopelagic marine fish belonging to the teleost family Sparidae, and a newly recruited species in Mediterranean aquaculture. The paucity of genetic information relating to sparids, despite their growing economic value for aquaculture, provides the impetus for exploring the genomics of this fish group. Genomic tool development, such as genetic linkage maps provision, lays the groundwork for linking genotype to phenotype, allowing fine-mapping of loci responsible for beneficial traits. In this study, we applied ddRAD methodology to identify polymorphic markers in a full-sib family of common pandora. Employing the Illumina MiSeq platform, we sampled and sequenced a size-selected genomic fraction of 99 individuals, which led to the identification of 920 polymorphic loci. Downstream mapping analysis resulted in the construction of 24 robust linkage groups, corresponding to the karyotype of the species. The common pandora linkage map showed varying degrees of conserved synteny with four other teleost genomes, namely the European seabass (Dicentrarchus labrax), Nile tilapia (Oreochromis niloticus), stickleback (Gasterosteus aculeatus), and medaka (Oryzias latipes), suggesting a conserved genomic evolution in Sparidae. Our work exploits the possibilities of genotyping by sequencing to gain novel insights into genome structure and evolution. Such information will boost the study of cultured species and will set the foundation for a deeper understanding of the complex evolutionary history of teleosts.

Second generation genetic linkage map for the gilthead sea bream Sparus aurata L.

An updated second linkage map was constructed for the gilthead sea bream, Sparus aurata L., a fish species of great economic importance for the Mediterranean aquaculture industry. In contrast to the first linkage map which mainly consisted of genomic microsatellites (SSRs), the new linkage map is highly enriched with SSRs found in Expressed Sequence Tags (EST-SSRs), which greatly facilitates comparative mapping with other teleosts. The new map consists of 321 genetic markers in 27 linkage groups (LGs): 232 genomic microsatellites, 85 EST-SSRs and 4 SNPs; of those, 13 markers were linked to LGs but were not ordered. Eleven markers (5 SSRs, 5 EST-SSRs and 1 SNP) are not assigned to any LG. The total length of the sex-averaged map is 1769.7cM, 42% longer than the previously published one, and the number of markers in each LG ranges from 2 to 30. The inter-marker distance varies from 0 to 75.6cM, with an average of 5.75cM. The male and female maps have a length of 1349.2 and 2172.1cM, respectively, and the average distance between markers is 4.38 and 7.05cM, respectively. Comparative mapping with the three-spined stickleback (Gasterosteus acuulatus) chromosomes and scaffolds showed conserved synteny with 132 S. aurata markers (42.9% of those mapped) having a hit on the stickleback genome.

In silico mining and characterization of simple sequence repeats from gilthead sea bream (Sparus aurata) expressed sequence tags (EST-SSRs); PCR amplification, polymorphism evaluation and multiplexing and cross-species assays.

We screened for simple sequence repeats (SSRs) found in ESTs derived from an EST-database development project ('Marine Genomics Europe' Network of Excellence). Different motifs of di-, tri-, tetra-, penta- and hexanucleotide SSRs were evaluated for variation in length and position in the expressed sequences, relative abundance and distribution in gilthead sea bream (Sparus aurata). We found 899 ESTs that harbor 997 SSRs (4.94%). On average, one SSR was found per 2.95 kb of EST sequence and the dinucleotide SSRs are the most abundant accounting for 47.6% of the total number. EST-SSRs were used as template for primer design. 664 primer pairs could be successfully identified and a subset of 206 pairs of primers was synthesized, PCR-tested and visualized on ethidium bromide stained agarose gels. The main objective was to further assess the potential of EST-SSRs as informative markers and investigate their cross-species amplification in sixteen teleost fish species: seven sparid species and nine other species from different families. Approximately 78% of the primer pairs gave PCR products of expected size in gilthead sea bream, and as expected, the rate of successful amplification of sea bream EST-SSRs was higher in sparids, lower in other perciforms and even lower in species of the Clupeiform and Gadiform orders. We finally determined the polymorphism and the heterozygosity of 63 markers in a wild gilthead sea bream population; fifty-eight loci were found to be polymorphic with the expected heterozygosity and the number of alleles ranging from 0.089 to 0.946 and from 2 to 27, respectively. These tools and markers are expected to enhance the available genetic linkage map in gilthead sea bream, to assist comparative mapping and genome analyses for this species and further with other model fish species and finally to help advance genetic analysis for cultivated and wild populations and accelerate breeding programs.

Gilthead sea bream (Sparus auratus) and European sea bass (Dicentrarchus labrax) expressed sequence tags: Characterization, tissue-specific expression and gene markers.

The gilthead sea bream, Sparus auratus, and the European sea bass, Dicentrarchus labrax, are two of the most important marine species cultivated in Southern Europe. This study aimed at increasing genomic resources for the two species and produced and annotated two sets of 30,000 expressed sequence tags (EST) each from 14 normalized tissue-specific cDNA libraries from sea bream and sea bass. Clustering and assembly of the ESTs formed 5268 contigs and 12,928 singletons for sea bream and 4573 contigs and 13,143 singletons for sea bass, representing 18,196 and 17,716 putative unigenes, respectively. Assuming a similar number of genes in sea bass, sea bream and in the model fish Gasterosteus aculeatus genomes, it was estimated that approximately two thirds of the sea bream and the sea bass transcriptomes were covered by the unigene collections. BLAST sequence similarity searches (using a cut off of e-value <10(-5)) against fully the curated SwissProt (and TrEMBL) databases produced matches of 28%(37%) and 43%(53%) of the sea bream and sea bass unigene datasets respectively, allowing some putative designation of function. A comparative approach is described using human Ensembl peptide ID homolog's for functional annotation, which increased the number of unigenes with GO terms assigned and resulted in more GO terms assigned per unigene. This allowed the identification of tissue-specific genes using enrichment analysis for GO pathways and protein domains. The comparative annotation approach represents a good strategy for transferring more relevant biological information from highly studied species to genomic resource poorer species. It was possible to confirm by interspecies mRNA-to-genomic alignments 25 and 21 alternative splice events in sea bream and sea bass genes, respectively. Even using normalized cDNA from relatively few pooled individuals it was possible to identify 1145 SNPs and 1748 microsatellites loci for genetic marker development. The EST data are being applied to a range of projects, including the development microarrays, genetic and radiation hybrid maps and QTL genome scans. This highlights the important role of ESTs for generating genetic and genomic resources of aquaculture species.

FDC-specific functions of p55TNFR and IKK2 in the development of FDC networks and of antibody responses.

The FDC-specific molecular signals required in the formation of FDC networks, B cell follicles, and germinal centers (GCs) have remained poorly understood. We used FDC-specific gene targeting to investigate the function of p55TNFR and IKK2 in lymphoid organ structure and function. Here we show that FDC-specific expression of p55TNFR is necessary and sufficient to promote FDC network and B cell follicle formation, restore the expression of CXCL13 and VCAM-1/ICAM-1 in FDCs, and lead to productive GCs. Notably, FDC-specific disruption of IKK2 does not affect formation of FDC networks. Yet, after antigen engagement or immune complex (IC) deposition, FDCs lacking IKK2 fail to upregulate VCAM-1 and ICAM-1, and GCs remain sterile. These findings demonstrate that IKK2-independent function of p55TNFR on FDCs is sufficient to support the development of FDC networks and GCs, while FDC-specific IKK2 is indispensable for the generation of efficient humoral immune responses.