De novo transcriptome assembly for a non-model species, the blood-sucking bug Triatoma brasiliensis, a vector of Chagas disease

Leveraging transcriptome sequence read archives for virus detection in wild and colony populations of triatomines (Hemiptera: Reduviidae: Triatominae)

Triatomines are infamous as vectors of the parasite Trypanosoma cruzi, the causative agent of Chagas disease. However, climate-driven range expansion and urbanization adaptation of triatomine populations, coupled with their highly diverse feeding strategies (vertebrate haematophagy, kleptohaematophagy, and coprophagy), and has elevated interest in triatomines as potential arboviral vectors. Information on the triatomine virome is scant, with prior records including only eight insect-specific viruses: Triatoma virus (TrV) and Rhodnius prolixus viruses 1-7. Here, we leverage publicly available transcriptome datasets to assess viral diversity in 122 wild and colony kissing bugs representing eight species from six countries. In total, six viruses were detected (including Rhodnius prolixus viruses 4-6), and TrV was detected in almost half of all screened triatomines. This is the first report of TrV in Triatoma brasiliensis and in members of the genus Mepraia (M. gajardoi, M. spinolai, and M. parapatrica), and this effort has vastly expanded the publicly available genomic resources of TrV, adding 39 genome sequences to the single genome sequence currently available in the GenBank database. Furthermore, two additional viruses-Meccus longipennis virus 1 and Drosophila melanogaster Nora virus-are herein reported for the first time from kissing bugs. Meccus longipennis virus 1 was detected in Triatoma infestans from Argentina, Brazil, Chile, and Peru, and Drosophila melanogaster Nora virus was found in T. infestans from Argentina. Our results illustrate the advantage and utility of low-cost transcriptome data mining for the discovery of known and novel arboviruses in triatomines and other potential insect vectors.

Population genetics and genomics of Triatoma brasiliensis (Hemiptera, Reduviidae) in an area of high pressure of domiciliary infestation in Northeastern Brazil

Understanding the population dynamics of vectors is crucial for effective control of vector-borne diseases. In the Northeastern Brazilian semi-arid region, Triatoma brasiliensis persists as the most significant Chagas disease vector, frequently displaying recurrent domiciliary infestations. This situation raises relevant public health concerns in the municipality of Currais Novos in the state of Rio Grande do Norte. This area has experienced a high prevalence of peridomiciliary re-infestations by T. brasiliensis, coupled with elevated rates of Trypanosoma cruzi infection. Therefore, we assessed the distribution of genetic variation via mitochondrial Cytochrome b gene (MT-CYB) sequencing (n = 109) and single nucleotide polymorphisms (SNPs, n = 86) to assess the gene flow among distinct populations distributed in varied geographic spots and environments, mainly sylvatic and peridomiciliary. Insects were collected from rural communities at Currais Novos, enclosed within a 16 km radius. Sampling included 13 populations: one intradomiciliary, eight peridomiciliary, and four sylvatic. Furthermore, an external population located 220 km from Currais Novos was also included in the study. The method employed to obtain SNP information relied on ddRAD-seq genotyping-by-sequencing (GBS), enabling a genome-wide analysis to infer genetic variation. Through AMOVA analysis of MT-CYB gene variation, we identified four distinct population groups with statistical significance (FCT= 0.42; p<0.05). We identified a total of 3,013 SNPs through GBS, with 11 loci showing putative signs of being under selection. The variation based on 3,002 neutral loci evidenced low genetic structuration based on low FST values (p>0.05), indicating local panmixia. However, resampling algorithms pointed out that three samples from the external population were assigned (>98 %) in a cluster contrasting from the ones putatively under local panmixia - validating the newly applied genome-wide marker for studies on the population genetics at finer-scale resolution for T. brasiliensis. The presence of population structuring in some of the sampled points, as suggested by the mitochondrial marker, leads us to assume that infestations were probably initiated by small populations of females - demographic event poses a risk for rapid re-infestations. The local panmictic pattern revealed by the GBS marker poses a challenge for vector control measures, as re-infestation foci may be distributed over a wide geographical and ecological range. In such instances, vectors exhibit reduced susceptibility to conventional insecticide spraying operations since sylvatic populations are beyond the reach of these interventions. The pattern of infestation exhibited by T. brasiliensis necessitates integrating innovative strategies into the existing control framework, holding the potential to create a more resilient and adaptive vector control program. In our dataset, the results demonstrated that the genetic signals from both markers were complementary. Therefore, it is essential to consider the nature and inheritance pattern of each marker when inferring the pattern of re-infestations.

Evidence that hematophagous triatomine bugs may eat plants in the wild

Blood feeding is a secondary adaptation in hematophagous bugs. Many proteins are secreted in the saliva that are devoted to coping with the host's defense and to process the blood meal. Digestive enzymes that are no longer required for a blood meal would be expected to be eventually lost. Yet, in many strictly hematophagous arthropods, α-amylase genes, which encode the enzymes that digest starch from plants, are still present and transcribed, including in the kissing bug Rhodnius prolixus (Hemiptera, Reduviidae) and its related species, which transmit the Chagas disease. We hypothesized that retaining α-amylase could be advantageous if the bugs occasionally consume plant tissues. We first checked that the α-amylase protein of Rhodnius robustus retains normal amylolytic activity. Then we surveyed hundreds of gut DNA extracts from the sylvatic R. robustus to detect traces of plants. We found plant DNA in 8% of the samples, mainly identified as Attalea palm trees, where R. robustus are usually found. We suggest that although of secondary importance in the blood-sucking bugs, α-amylase may be needed during occasional plant feeding and thus has been retained.

Phylogenomics for Chagas Disease Vectors of the Rhodnius Genus (Hemiptera, Triatominae): What We Learn From Mito-Nuclear Conflicts and Recommendations

We provide in this study a very large DNA dataset on Rhodnius species including 36 samples representing 16 valid species of the three Rhodnius groups, pictipes, prolixus and pallescens. Samples were sequenced at low-depth with whole-genome shotgun sequencing (Illumina technology). Using phylogenomics including 15 mitochondrial genes (13.3 kb), partial nuclear rDNA (5.2 kb) and 51 nuclear protein-coding genes (36.3 kb), we resolve sticking points in the Rhodnius phylogeny. At the species level, we confirmed the species-specific status of R. montenegrensis and R. marabaensis and we agree with the synonymy of R. taquarussuensis with R. neglectus. We also invite to revisit the species-specific status of R. milesi that is more likely R. nasutus. We proposed to define a robustus species complex that comprises the four close relative species: R. marabaensis, R. montenegrensis, R. prolixus and R. robustus. As Psammolestes tertius was included in the Rhodnius clade, we strongly recommend reclassifying this species as R. tertius. At the Rhodnius group level, molecular data consistently supports the clustering of the pictipes and pallescens groups, more related to each other than they are to the prolixus group. Moreover, comparing mitochondrial and nuclear tree topologies, our results demonstrated that various introgression events occurred in all the three Rhodnius groups, in laboratory strains but also in wild specimens. We demonstrated that introgressions occurred frequently in the prolixus group, involving the related species of the robustus complex but also the pairwise R. nasutus and R. neglectus. A genome wide analysis highlighted an introgression event in the pictipes group between R. stali and R. brethesi and suggested a complex gene flow between the three species of the pallescens group, R. colombiensis, R. pallescens and R. ecuadoriensis. The molecular data supports also a sylvatic distribution of R. prolixus in Brazil (Pará state) and the monophyly of R. robustus. As we detected extensive introgression events and selective pressure on mitochondrial genes, we strongly recommend performing separate mitochondrial and nuclear phylogenies and to take advantages of mito-nuclear conflicts in order to have a comprehensive evolutionary vision of this genus.

Chemosensory Gene Expression for Two Closely Relative Species Rhodnius robustus and R. prolixus (Hemiptera, Reduviidade, Triatominae) Vectors of Chagas Disease

Two closely related species, Rhodnius prolixus and Rhodnius robustus, are the vectors of Trypanosoma cruzi, which is the causative agent of Chagas disease, but clearly exhibit clear-cut differences in their ecological behavior. R. prolixus is considered as a domiciliated species, whereas R. robustus only sporadically visits human houses in Amazonia. We performed a chemosensory gene expression study via RNA-sequencing (RNA-seq) for the two species and also included a laboratory introgressed R. robustus strain. We built an assembled transcriptome for each sample and for both sexes and compiled all in a reference transcriptome for a differential gene expression study. Because the genes specifically expressed in one condition and not expressed in another may also reflect differences in the adaptation of organisms, a comparative study of the presence/absence of transcripts was also performed for the chemosensory transcripts, namely chemosensory proteins (CSPs), odorant-binding proteins (OBPs), odorant receptors (ORs), gustatory receptors (GRs), and ionotropic receptors (IRs), as well as takeout (TO) transcripts because TO proteins have been proposed to be associated with chemosensory perception in both olfactory and taste systems. In this study, 12 novel TO transcripts from the R. prolixus genome were annotated. Among the 199 transcripts, out of interest, annotated in this study, 93% were conserved between R. prolixus and the sylvatic R. robustus. Moreover, 10 transcripts out of interest were specifically expressed in one sex and absent in another. Three chemosensory transcripts were found to be expressed only in the reared R. prolixus (CSP19, OBP9, and OR89) and only one in sylvatic R. robustus (OR22). A large set of transcripts were found to be differentially expressed (DE) between males and females (1,630), with a majority of them (83%) overexpressed in males. Between environmental conditions, 8,596 transcripts were DE, with most (67%) overexpressed in the sylvatic R. robustus samples, including 17 chemosensory transcripts (4 CSPs, 1 OBP, 5 ORs, 1 GR, 4 IR, and 2 TO), but 4 genes (OBP19, OR13, OR40, and OR79) were overexpressed in the reared samples.

Dynamics of food sources, ecotypic distribution and Trypanosoma cruzi infection in Triatoma brasiliensis from the northeast of Brazil

Innovative approaches used to combat Chagas disease transmission tend to combine a set of comprehensive efforts to understand the ecology of local vectors. In this work we identified molecularly the blood meal of 181 Triatoma brasiliensis, distributed in 18 populations (8 sylvatic and 10 peridomestic), which were collected across a range of 240 km (East-West) and 95 km (North-South) in the semi-arid region of northeastern, Brazil. We used the vertebrate mitochondrial gene (cytochrome B) sequencing applied to DNA isolated from bug midgut to identify the insect blood meal sources via the BLAST procedure. The peridomestic populations were classified according to two main hypotheses of site-occupancy for T. brasiliensis: the first says that the infestation is mainly driven by structures that resemble its natural habitat (stony-like ecotopes) and the second assumes that it is associated with key-hosts (rodents and goats). Rodents of the Caviidae family (Galea spixii and Kerodon rupestris) were identified as the key-host of T. brasiliensis, but also the potential Trypanosoma cruzi reservoir-able to connect the sylvatic and domestic T. cruzi cycle. Cats also deserve to be studied better, as potential T. cruzi reservoirs. By modeling the food sources + site-occupancy + T. cruzi natural infection, we identified man-made ecotopes suitable for forming dense triatomine infestations with high rates of T. cruzi natural infection, which may be taken into account for vector control measures.

Evolution, Systematics, and Biogeography of the Triatominae, Vectors of Chagas Disease

In this chapter, we review and update current knowledge about the evolution, systematics, and biogeography of the Triatominae (Hemiptera: Reduviidae)-true bugs that feed primarily on vertebrate blood. In the Americas, triatomines are the vectors of Trypanosoma cruzi, the etiological agent of Chagas disease. Despite declining incidence and prevalence, Chagas disease is still a major public health concern in Latin America. Triatomines occur also in the Old World, where vector-borne T. cruzi transmission has not been recorded. Triatomines evolved from predatory reduviid bugs, most likely in the New World, and diversified extensively across the Americas (including the Caribbean) and in parts of Asia and Oceania. Here, we first discuss our current understanding of how, how many times, and when the blood-feeding habit might have evolved among the Reduviidae. Then we present a summary of recent advances in the systematics of this diverse group of insects, with an emphasis on the contribution of molecular tools to the clarification of taxonomic controversies. Finally, and in the light of both up-to-date phylogenetic hypotheses and a thorough review of distribution records, we propose a global synthesis of the biogeography of the Triatominae. Over 130 triatomine species contribute to maintaining T. cruzi transmission among mammals (sometimes including humans) in almost every terrestrial ecoregion of the Americas. This means that Chagas disease will never be eradicated and underscores the fact that effective disease prevention will perforce require stronger, long-term vector control-surveillance systems.

Coexpression of PalbHLH1 and PalMYB90 Genes From Populus alba Enhances Pathogen Resistance in Poplar by Increasing the Flavonoid Content

Secondary metabolites of the flavonoid pathway participate in plant defense, and bHLH and MYB transcription factors regulate the synthesis of these metabolites. Here, we define the regulatory mechanisms in response to pathogens. Two transcription factors from Populus alba var. pyramidalis, PalbHLH1 and PalMYB90, were overexpressed together in poplar, and transcriptome analysis revealed differences in response to pathogen infection. The transgenic plants showed elevated levels of several key flavonoid pathway components: total phenols, proanthocyanidins (PAs), and anthocyanins and intermediates quercetin and kaempferol. Furthermore, PalbHLH1 and PalMYB90 overexpression in poplar enhanced antioxidase activities and H₂O₂ release and also increased resistance to Botrytis cinerea and Dothiorella gregaria infection. Gene expression profile analysis showed most genes involved in the flavonoid biosynthesis pathway or antioxidant response to be upregulated in MYB90/bHLH1-OE poplar, but significant differential expression occurred in response to pathogen infection. Specifically, expression of PalF3H (flavanone 3-hydroxylase), PalDFR (dihydroflavonol 4-seductase), PalANS (anthocyanin synthase), and PalANR (anthocyanin reductase), which function in initial, middle, and final steps of anthocyanin and PA biosynthesis, respectively, was significantly upregulated in D. gregaria-infected MYB90/bHLH1-OE poplar. Our results highlight that PalbHLH1 and PalMYB90 function as transcriptional activators of flavonoid pathway secondary-metabolite synthesis genes, with differential mechanisms in response to bacterial or fungal infection.

Adaptations in energy metabolism and gene family expansions revealed by comparative transcriptomics of three Chagas disease triatomine vectors

Background

Chagas disease is a parasitic infection caused by Trypanosoma cruzi. It is an important public health problem affecting around seven to eight million people in the Americas. A large number of hematophagous triatomine insect species, occupying diverse natural and human-modified ecological niches transmit this disease. Triatomines are long-living hemipterans that have evolved to explode different habitats to associate with their vertebrate hosts. Understanding the molecular basis of the extreme physiological conditions including starvation tolerance and longevity could provide insights for developing novel control strategies. We describe the normalized cDNA, full body transcriptome analysis of three main vectors in North, Central and South America, Triatoma pallidipennis, T. dimidiata and T. infestans.

Results

Two-thirds of the de novo assembled transcriptomes map to the Rhodnius prolixus genome and proteome. A Triatoma expansion of the calycin family and two types of protease inhibitors, pacifastins and cystatins were identified. A high number of transcriptionally active class I transposable elements was documented in T. infestans, compared with T. dimidiata and T. pallidipennis. Sequence identity in Triatoma-R. prolixus 1:1 orthologs revealed high sequence divergence in four enzymes participating in gluconeogenesis, glycogen synthesis and the pentose phosphate pathway, indicating high evolutionary rates of these genes. Also, molecular evidence suggesting positive selection was found for several genes of the oxidative phosphorylation I, III and V complexes.

Conclusions

Protease inhibitors and calycin-coding gene expansions provide insights into rapidly evolving processes of protease regulation and haematophagy. Higher evolutionary rates in enzymes that exert metabolic flux control towards anabolism and evidence for positive selection in oxidative phosphorylation complexes might represent genetic adaptations, possibly related to prolonged starvation, oxidative stress tolerance, longevity, and hematophagy and flight reduction. Overall, this work generated novel hypothesis related to biological adaptations to extreme physiological conditions and diverse ecological niches that sustain Chagas disease transmission.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4696-8) contains supplementary material, which is available to authorized users.

Tissue-Specific Transcriptome Analysis Reveals Multiple Responses to Salt Stress in Populus euphratica Seedlings

Salt stress is one of the most crucial factors impacting plant growth, development and reproduction. However, information regarding differences in tissue-specific gene expression patterns, which may improve a plant’s tolerance to salt stress, is limited. Here, we investigated the gene expression patterns in tissues of Populus euphratica Oliv. seedlings using RNA sequencing (RNA-Seq) technology. A total of 109.3 million, 125bp paired-end clean reads were generated, and 6428, 4797, 2335 and 3358 differentially expressed genes (DEGs) were identified in leaf, phloem, xylem and root tissues, respectively. While the tissue-specific DEGs under salt stress had diverse functions, “membrane transporter activity” was the most significant leaf function, whereas “oxidation–reduction process” was the most significant function in root tissue. Further analysis of the tissue-specific DEGs showed that the expression patterns or functions of gene families, such as SOS, NHX, GolS, GPX, APX, RBOHF and CBL, were diverse, suggesting that calcium signaling, reactive oxygen species (ROS) and salt overly sensitive (SOS) pathways are all involved in ionic homeostasis in tissues from P. euphratica seedlings. The DEGs, for example the up-regulated antioxidant genes, contribute to ROS-scavenging induced by salt stress but result in decreased Na⁺ concentrations in root vasculature cells and in xylem sap, while the down-regulated rbohF leads to the reverse results. These results suggest that the divergence of DEGs expression patterns contribute to maintenance of ionic and ROS homeostasis in tissues and improve plant salinity tolerance. We comprehensively analyzed the response of P. euphratica seedlings to salt stress and provide helpful genetic resources for studying plant-abiotic stress interactions.

Neuropeptidomics in Triatoma infestans. Comparative transcriptomic analysis among triatomines

Chagas' disease, affecting up to 6-7 million people worldwide, is transmitted to humans through the feces of triatomine kissing bugs. From these, Rhodnius prolixus, Triatoma dimidiata, Triatoma infestans and Triatoma pallidipennis are important vectors distributed throughout the Latin American subcontinent. Resistance to pyrethroids has been developed by some triatomine populations, especially T. infestans, obstructing their control. Given their role in the regulation of physiological processes, neuroendocrine-derived factors have been proposed as a source of molecular targets for new-generation insecticides. However, the involvement of neuropeptides in insecticide metabolism and resistance in insects has been poorly studied. In the present work, the sequences of 20 neuropeptide precursor genes in T. infestans, 16 in T. dimidiata, and 13 in T. pallidipennis detected in transcriptomic databases are reported, and a comparative analysis in triatomines is presented. A total of 59 neuropeptides were validated by liquid chromatography-tandem mass spectrometry in brain and nervous ganglia from T. infestans, revealing the existence of differential post-translational modifications, extended and truncated forms. The results suggest a high sequence conservation in some neuropeptide systems in triatomines, whereas remarkable differences occur in several others within the core domains. Comparisons of the basal expression levels for several neuropeptide precursor genes between pyrethroid sensitive and resistant population of T. infestans are also presented here, in order to introduce a proof of concept to test the involvement of neuropeptides in insecticide resistance. From the precursors tested, NVP and ITG peptides are significantly higher expressed in the resistant population. To our knowledge, this is the first report to associate differential neuropeptide expression with insecticide resistance. The information provided here contributes to creating conditions to widely extend functional and genetic studies involving neuropeptides in triatomines.

Under-Expression of Chemosensory Genes in Domiciliary Bugs of the Chagas Disease Vector Triatoma brasiliensis

BACKGROUND

In Latin America, the bloodsucking bugs Triatominae are vectors of Trypanosoma cruzi, the parasite that causes Chagas disease. Chemical elimination programs have been launched to control Chagas disease vectors. However, the disease persists because native vectors from sylvatic habitats are able to (re)colonize houses-a process called domiciliation. Triatoma brasiliensis is one example. Because the chemosensory system allows insects to interact with their environment and plays a key role in insect adaption, we conducted a descriptive and comparative study of the chemosensory transcriptome of T. brasiliensis samples from different ecotopes.

METHODOLOGY/PRINCIPAL FINDING

In a reference transcriptome built using de novo assembly, we found transcripts encoding 27 odorant-binding proteins (OBPs), 17 chemosensory proteins (CSPs), 3 odorant receptors (ORs), 5 transient receptor potential channel (TRPs), 1 sensory neuron membrane protein (SNMPs), 25 takeout proteins, 72 cytochrome P450s, 5 gluthatione S-transferases, and 49 cuticular proteins. Using protein phylogenies, we showed that most of the OBPs and CSPs for T. brasiliensis had well supported orthologs in the kissing bug Rhodnius prolixus. We also showed a higher number of these genes within the bloodsucking bugs and more generally within all Hemipterans compared to the other species in the super-order Paraneoptera. Using both DESeq2 and EdgeR software, we performed differential expression analyses between samples of T. brasiliensis, taking into account their environment (sylvatic, peridomiciliary and domiciliary) and sex. We also searched clusters of co-expressed contigs using HTSCluster. Among differentially expressed (DE) contigs, most were under-expressed in the chemosensory organs of the domiciliary bugs compared to the other samples and in females compared to males. We clearly identified DE genes that play a role in the chemosensory system.

CONCLUSION/SIGNIFICANCE

Chemosensory genes could be good candidates for genes that contribute to adaptation or plastic rearrangement to an anthropogenic system. The domiciliary environment probably includes less diversity of xenobiotics and probably has more stable abiotic parameters than do sylvatic and peridomiciliary environments. This could explain why both detoxification and cuticle protein genes are less expressed in domiciliary bugs. Understanding the molecular basis for how vectors adapt to human dwellings may reveal new tools to control disease vectors; for example, by disrupting chemical communication.

Global transcriptome changes in perennial ryegrass during early infection by pink snow mould

Lack of resistance to pink snow mould (Microdochium nivale) is a major constraint for adaptation of perennial ryegrass (Lolium perenne L.) to continental regions with long-lasting snow cover at higher latitudes. Almost all investigations of genetic variation in resistance have been performed using cold acclimated plants. However, there may be variation in resistance mechanisms that are functioning independently of cold acclimation. In this study our aim was to identify candidate genes involved in such resistance mechanisms. We first characterized variation in resistance to M. nivale among non-acclimated genotypes from the Norwegian cultivar 'Fagerlin' based on relative regrowth and fungal quantification by real-time qPCR. One resistant and one susceptible genotype were selected for transcriptome analysis using paired-end sequencing by Illumina Hiseq 2000. Transcriptome profiles, GO enrichment and KEGG pathway analysis indicate that defense response related genes are differentially expressed between the resistant and the susceptible genotype. A significant up-regulation of defense related genes, as well as genes involved in cell wall cellulose metabolic processes and aryl-alcohol dehydrogenase (NADP+) activity, was observed in the resistant genotype. The candidate genes identified in this study might be potential molecular marker resources for breeding perennial ryegrass cultivars with improved resistance to pink snow mould.

Disease vectors in the era of next generation sequencing

Almost 20 % of all infectious human diseases are vector borne and, together, are responsible for over one million deaths per annum. Over the past decade, the decreasing costs of massively parallel sequencing technologies have facilitated the agnostic interrogation of insect vector genomes, giving medical entomologists access to an ever-expanding volume of high-quality genomic and transcriptomic data. In this review, we highlight how genomics resources have provided new insights into the physiology, behavior, and evolution of human disease vectors within the context of the global health landscape.

Comparing de novo and reference-based transcriptome assembly strategies by applying them to the blood-sucking bug Rhodnius prolixus

High Throughput Sequencing capabilities have made the process of assembling a transcriptome easier, whether or not there is a reference genome. But the quality of a transcriptome assembly must be good enough to capture the most comprehensive catalog of transcripts and their variations, and to carry out further experiments on transcriptomics. There is currently no consensus on which of the many sequencing technologies and assembly tools are the most effective. Many non-model organisms lack a reference genome to guide the transcriptome assembly. One question, therefore, is whether or not a reference-based genome assembly gives better results than de novo assembly. The blood-sucking insect Rhodnius prolixus-a vector for Chagas disease-has a reference genome. It is therefore a good model on which to compare reference-based and de novo transcriptome assemblies. In this study, we compared de novo and reference-based genome assembly strategies using three datasets (454, Illumina, 454 combined with Illumina) and various assembly software. We developed criteria to compare the resulting assemblies: the size distribution and number of transcripts, the proportion of potentially chimeric transcripts, how complete the assembly was (completeness evaluated both through CEGMA software and R. prolixus proteome fraction retrieved). Moreover, we looked for the presence of two chemosensory gene families (Odorant-Binding Proteins and Chemosensory Proteins) to validate the assembly quality. The reference-based assemblies after genome annotation were clearly better than those generated using de novo strategies alone. Reference-based strategies revealed new transcripts, including new isoforms unpredicted by automatic genome annotation. However, a combination of both de novo and reference-based strategies gave the best result, and allowed us to assemble fragmented transcripts.

Successful Recovery of Nuclear Protein-Coding Genes from Small Insects in Museums Using Illumina Sequencing

In this paper we explore high-throughput Illumina sequencing of nuclear protein-coding, ribosomal, and mitochondrial genes in small, dried insects stored in natural history collections. We sequenced one tenebrionid beetle and 12 carabid beetles ranging in size from 3.7 to 9.7 mm in length that have been stored in various museums for 4 to 84 years. Although we chose a number of old, small specimens for which we expected low sequence recovery, we successfully recovered at least some low-copy nuclear protein-coding genes from all specimens. For example, in one 56-year-old beetle, 4.4 mm in length, our de novo assembly recovered about 63% of approximately 41,900 nucleotides in a target suite of 67 nuclear protein-coding gene fragments, and 70% using a reference-based assembly. Even in the least successfully sequenced carabid specimen, reference-based assembly yielded fragments that were at least 50% of the target length for 34 of 67 nuclear protein-coding gene fragments. Exploration of alternative references for reference-based assembly revealed few signs of bias created by the reference. For all specimens we recovered almost complete copies of ribosomal and mitochondrial genes. We verified the general accuracy of the sequences through comparisons with sequences obtained from PCR and Sanger sequencing, including of conspecific, fresh specimens, and through phylogenetic analysis that tested the placement of sequences in predicted regions. A few possible inaccuracies in the sequences were detected, but these rarely affected the phylogenetic placement of the samples. Although our sample sizes are low, an exploratory regression study suggests that the dominant factor in predicting success at recovering nuclear protein-coding genes is a high number of Illumina reads, with success at PCR of COI and killing by immersion in ethanol being secondary factors; in analyses of only high-read samples, the primary significant explanatory variable was body length, with small beetles being more successfully sequenced.

Deep sequencing of the Mexican avocado transcriptome, an ancient angiosperm with a high content of fatty acids

BACKGROUND

Avocado (Persea americana) is an economically important tropical fruit considered to be a good source of fatty acids. Despite its importance, the molecular and cellular characterization of biochemical and developmental processes in avocado is limited due to the lack of transcriptome and genomic information.

RESULTS

The transcriptomes of seeds, roots, stems, leaves, aerial buds and flowers were determined using different sequencing platforms. Additionally, the transcriptomes of three different stages of fruit ripening (pre-climacteric, climacteric and post-climacteric) were also analyzed. The analysis of the RNAseqatlas presented here reveals strong differences in gene expression patterns between different organs, especially between root and flower, but also reveals similarities among the gene expression patterns in other organs, such as stem, leaves and aerial buds (vegetative organs) or seed and fruit (storage organs). Important regulators, functional categories, and differentially expressed genes involved in avocado fruit ripening were identified. Additionally, to demonstrate the utility of the avocado gene expression atlas, we investigated the expression patterns of genes implicated in fatty acid metabolism and fruit ripening.

CONCLUSIONS

A description of transcriptomic changes occurring during fruit ripening was obtained in Mexican avocado, contributing to a dynamic view of the expression patterns of genes involved in fatty acid biosynthesis and the fruit ripening process.

Next-generation sequencing as a powerful motor for advances in the biological and environmental sciences

Next-generation sequencing (NGS) provides unprecedented insight into (meta)genomes, (meta)transcriptomes (cDNA) and (meta)barcodes of individuals, populations and communities of Archaea, Bacteria and Eukarya, as well as viruses. This special issue combines reviews and original papers reporting technical and scientific advances in genomics and transcriptomics of non-model species, as well as quantification and functional analyses of biodiversity using NGS technologies of the second and third generations. In addition, certain papers also exemplify the transition from Sanger to NGS barcodes in molecular taxonomy.