Benchmarking next-generation transcriptome sequencing for functional and evolutionary genomics

The sequencing of the key genes and end products in the TLR4 signaling pathway from the kidney of Rana dybowskii exposed to Aeromonas hydrophila

Infectious diseases caused by Aeromonas hydrophila (AH) have reduced the populations of Rana dybowskii). However, little is known about the immune response of R. dybowskii against AH infections. The toll-like receptor (TLR) signaling pathway has been identified as a critical component in innate immunity, responsible for identifying pathogen-associated molecular patterns in pathogens. Our study used the next-generation sequencing technique and single-molecule long-read sequencing to determine the structures of transcript isoforms and functions of genes in the kidneys of R. dybowskii, as well as identify and validate the related genes in the TLR4 signaling pathway. In total, 628,774 reads of inserts were identified, including 300,053 full-length non-chimeric reads and 233,592 non-full-length reads. Among the transcriptome sequences, 124 genes were identified as homologs of known genes in the TLR4 pathway especially inflammatory cytokines and receptors. Our findings shed light on the structures and functions of R. dybowskii genes exposed to AH and confirm the presence of both MyD88-dependent and independent pathways in R. dybowskii. Our work reveals how various functional proteins in amphibians at the initial stage of immune response are activated and complete their corresponding functions in a short time.

TransPi - a comprehensive TRanscriptome ANalysiS PIpeline for de novo transcriptome assembly

The use of RNA-Seq data and the generation of de novo transcriptome assemblies have been pivotal for studies in ecology and evolution. This is distinctly true for non-model organisms, where no genome information is available. In such organisms, studies of differential gene expression, DNA enrichment baits design, and phylogenetics can all be accomplished with de novo transcriptome assemblies. Multiple tools are available for transcriptome assembly, however, no single tool can provide the best assembly for all datasets. Therefore, a multi assembler approach, followed by a reduction step, is often sought to generate an improved representation of the assembly. To reduce errors in these complex analyses while at the same time attaining reproducibility and scalability, automated workflows have been essential in the analysis of RNA-Seq data. However, most of these tools are designed for species where genome data is used as reference for the assembly process, limiting their use in non-model organisms. We present TransPi, a comprehensive pipeline for de novo transcriptome assembly, with minimum user input but without losing the ability of a thorough analysis. A combination of different model organisms, k-mer sets, read lengths, and read quantities were used for assessing the tool. Furthermore, a total of 49 non-model organisms, spanning different phyla, were also analysed. Compared to approaches using single assemblers only, TransPi produces higher BUSCO completeness percentages, and a concurrent significant reduction in duplication rates. TransPi is easy to configure and can be deployed seamlessly using Conda, Docker and Singularity.

Plasticity and modulation of olfactory circuits in insects

Olfactory circuits change structurally and physiologically during development and adult life. This allows insects to respond to olfactory cues in an appropriate and adaptive way according to their physiological and behavioral state, and to adapt to their specific abiotic and biotic natural environment. We highlight here findings on olfactory plasticity and modulation in various model and non-model insects with an emphasis on moths and social Hymenoptera. Different categories of plasticity occur in the olfactory systems of insects. One type relates to the reproductive or feeding state, as well as to adult age. Another type of plasticity is context-dependent and includes influences of the immediate sensory and abiotic environment, but also environmental conditions during postembryonic development, periods of adult behavioral maturation, and short- and long-term sensory experience. Finally, plasticity in olfactory circuits is linked to associative learning and memory formation. The vast majority of the available literature summarized here deals with plasticity in primary and secondary olfactory brain centers, but also peripheral modulation is treated. The described molecular, physiological, and structural neuronal changes occur under the influence of neuromodulators such as biogenic amines, neuropeptides, and hormones, but the mechanisms through which they act are only beginning to be analyzed.

De novo assembly of Auricularia polytricha transcriptome and discovery of genes involved in the degradation of lignocellulose

Auricularia polytricha belonging to Basidiomycota has the ability to degrade lignocellulose. However, there has been no resource in public databases examining the transcriptome of A. polytricha. In this study, high-throughput sequencing platform BGISEQ-500 was used to generate large amount of transcript sequences from A. polytricha for gene discovery and molecular marker development. A total of 28,102 unigenes were discovered from the assembly of clean reads. In addition, functional categorization of the gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) metabolic pathways revealed several important biological processes. GO annotation analysis presented 47 categories, with the major subcategories being catalytic activity, binding, cellular process, metabolic process, and cell. Among the five functional categories and 21 subcategories of processes discovered from KEGG, global and overview maps, carbohydrate metabolism, transport, and catabolism are the main subcategories. Furthermore, among the unigenes related to lignocellulosic degradation discovered by KEGG pathway enrichment analysis, 2, 5, and 16 unigenes in de novo assembly of A. polytricha transcriptome were found to relate to cellulose, hemicellulose, and lignin degradation, respectively. The study provided valuable information on the degradation of lignocellulose to facilitate research on the degradation mechanism, molecular marker, functional research, gene mapping, and other multigenomic studies of species containing lignocellulose.

Species diversity and phylogeography of Cornus kousa (Asian dogwood) captured by genomic and genic microsatellites

Cornus kousa (Asian dogwood), an East Asia native tree, is the most economically important species of the dogwood genus, owing to its desirable horticultural traits and ability to hybridize with North America-native dogwoods. To assess the species genetic diversity and to better inform the ongoing and future breeding efforts, we assembled an herbarium and arboretum collection of 131 noncultivated C. kousa specimens. Genotyping and capillary electrophoresis analyses of our C. kousa collection with the newly developed genic and published nuclear genomic microsatellites permitted assessment of genetic diversity and evolutionary history of the species. Regardless of the microsatellite type used, the study yielded generally similar insights into the C. kousa diversity with subtle differences deriving from and underlining the marker used. The accrued evidence pointed to the species distinct genetic pools related to the plant country of origin. This can be helpful in the development of the commercial cultivars for this important ornamental crop with increased pyramided utility traits. Analyses of the C. kousa evolutionary history using the accrued genotyping datasets pointed to an unsampled ancestor population, possibly now extinct, as per the phylogeography of the region. To our knowledge, there are few studies utilizing the same gDNA collection to compare performance of genomic and genic microsatellites. This is the first detailed report on C. kousa species diversity and evolutionary history inference.

De novo assembly and characterization of transcriptome in the medicinal plant Euphorbia jolkini

BACKGROUND

Euphorbia jolkini, a medicinal herb that grows on the warm beaches in Japan and South Korea, is known to be used for traditional medicines to treat a variety of ailments, including bruises, stiffness, indigestion, toothache, and diabetes.

OBJECTIVE

It is to analyze the whole transcriptome and identify the genes related to the phenylpropanoid biosynthesis in the medicinally important herb E jolkini.

METHODS

Paired-end Illumina HiSeq™ 2500 sequencing technology was employed for cDNA library construction and Illumina sequencing. Public databases like TAIR (The Arabidopsis Information Resource), Swissprot and KEGG (Kyoto Encyclopedia of Genes and Genomes) were used for annotations of unigenes obtained.

RESULTS

The transcriptome of E. jolkini generated 139,215 assembled transcripts with an average length of 868 bp and an N50 value of 1460 bp that were further clustered using CD-HIT into 93,801 unigenes with an average length of 847 bp (N50-1410 bp). Sixty-three percent of the coding sequences (CDS) were annotated from the longest open reading frame (ORF). A remarkable percentage of unigenes were annotated against various databases. The differentially expressed gene analysis revealed that the expression of genes related to the terpenoid backbone biosynthesis pathway was higher in the flowers, whereas that of genes related to the phenylpropanoid biosynthesis pathway was both up- and downregulated in flowers and leaves. A search of against the transcription factor domain found 1023 transcription factors (TFs) that were from 54 TF families.

CONCLUSION

Assembled sequences of the E. jolkini transcriptome are made available for the first time in this study E. jolkini and lay a foundation for the investigation of secondary metabolite biosynthesis.

Metagenomic analysis of relative abundance and diversity of bacterial microbiota in Bemisia tabaci infesting cotton crop in Pakistan

B. tabaci species complex are among the world's most devastating agricultural pests causing economic losses by direct feeding and more importantly by transmitting plant viruses like cotton leaf curl disease (CLCuD) associated viruses to cultivated cotton in Pakistan. Taxonomic diversity of B. tabaci associated bacterial communities using NGS techniques so far is reported from insects grown on artificial diet under lab conditions. In this study 16S rDNA metagenome sequencing analysis was used to characterize bacterial compositions in wild adult B. tabaci infesting cultivated cotton in eight major cotton growing districts of southern Punjab, Pakistan. We have identified 50 known and 7 unknown genera of bacteria belonging to 10 phyla, 20 classes, 30 orders and 40 families. Beta diversity analysis of our data sets reveal that whiteflies infesting cotton in geographically distinct locations had similar bacterial diversity. These results for the first time provide insights into the microbiome diversity of wild type whiteflies infesting a cultivated crop.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm

BACKGROUND

The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity.

RESULTS

We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs.

CONCLUSIONS

Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

Multi-tissue transcriptomes of caecilian amphibians highlight incomplete knowledge of vertebrate gene families

RNA sequencing (RNA-seq) has become one of the most powerful tools to unravel the genomic basis of biological adaptation and diversity. Although challenging, RNA-seq is particularly promising for research on non-model, secretive species that cannot be observed in nature easily and therefore remain comparatively understudied. Among such animals, the caecilians (order Gymnophiona) likely constitute the least known group of vertebrates, despite being an old and remarkably distinct lineage of amphibians. Here, we characterize multi-tissue transcriptomes for five species of caecilians that represent a broad level of diversity across the order. We identified vertebrate homologous elements of caecilian functional genes of varying tissue specificity that reveal a great number of unclassified gene families, especially for the skin. We annotated several protein domains for those unknown candidate gene families to investigate their function. We also conducted supertree analyses of a phylogenomic dataset of 1,955 candidate orthologous genes among five caecilian species and other major lineages of vertebrates, with the inferred tree being in agreement with current views of vertebrate evolution and systematics. Our study provides insights into the evolution of vertebrate protein-coding genes, and a basis for future research on the molecular elements underlying the particular biology and adaptations of caecilian amphibians.

Steroid Hormone Function Controls Non-competitive Plasmodium Development in Anopheles

Transmission of malaria parasites occurs when a female Anopheles mosquito feeds on an infected host to acquire nutrients for egg development. How parasites are affected by oogenetic processes, principally orchestrated by the steroid hormone 20-hydroxyecdysone (20E), remains largely unknown. Here we show that Plasmodium falciparum development is intimately but not competitively linked to processes shaping Anopheles gambiae reproduction. We unveil a 20E-mediated positive correlation between egg and oocyst numbers; impairing oogenesis by multiple 20E manipulations decreases parasite intensities. These manipulations, however, accelerate Plasmodium growth rates, allowing sporozoites to become infectious sooner. Parasites exploit mosquito lipids for faster growth, but they do so without further affecting egg development. These results suggest that P. falciparum has adopted a non-competitive evolutionary strategy of resource exploitation to optimize transmission while minimizing fitness costs to its mosquito vector. Our findings have profound implications for currently proposed control strategies aimed at suppressing mosquito populations.

Usability of reference-free transcriptome assemblies for detection of differential expression: a case study on Aethionema arabicum dimorphic seeds

Background

RNA-sequencing analysis is increasingly utilized to study gene expression in non-model organisms without sequenced genomes. Aethionema arabicum (Brassicaceae) exhibits seed dimorphism as a bet-hedging strategy – producing both a less dormant mucilaginous (M⁺) seed morph and a more dormant non-mucilaginous (NM) seed morph. Here, we compared de novo and reference-genome based transcriptome assemblies to investigate Ae. arabicum seed dimorphism and to evaluate the reference-free versus -dependent approach for identifying differentially expressed genes (DEGs).

Results

A de novo transcriptome assembly was generated using sequences from M⁺ and NM Ae. arabicum dry seed morphs. The transcripts of the de novo assembly contained 63.1% complete Benchmarking Universal Single-Copy Orthologs (BUSCO) compared to 90.9% for the transcripts of the reference genome. DEG detection used the strict consensus of three methods (DESeq2, edgeR and NOISeq). Only 37% of 1533 differentially expressed de novo assembled transcripts paired with 1876 genome-derived DEGs. Gene Ontology (GO) terms distinguished the seed morphs: the terms translation and nucleosome assembly were overrepresented in DEGs higher in abundance in M⁺ dry seeds, whereas terms related to mRNA processing and transcription were overrepresented in DEGs higher in abundance in NM dry seeds. DEGs amongst these GO terms included ribosomal proteins and histones (higher in M⁺), RNA polymerase II subunits and related transcription and elongation factors (higher in NM). Expression of the inferred DEGs and other genes associated with seed maturation (e.g. those encoding late embryogenesis abundant proteins and transcription factors regulating seed development and maturation such as ABI3, FUS3, LEC1 and WRI1 homologs) were put in context with Arabidopsis thaliana seed maturation and indicated that M⁺ seeds may desiccate and mature faster than NM. The 1901 transcriptomic DEG set GO-terms had almost 90% overlap with the 2191 genome-derived DEG GO-terms.

Conclusions

Whilst there was only modest overlap of DEGs identified in reference-free versus -dependent approaches, the resulting GO analysis was concordant in both approaches. The identified differences in dry seed transcriptomes suggest mechanisms underpinning previously identified contrasts between morphology and germination behaviour of M⁺ and NM seeds.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5452-4) contains supplementary material, which is available to authorized users.

Usability of reference-free transcriptome assemblies for detection of differential expression: a case study on Aethionema arabicum dimorphic seeds

BACKGROUND

RNA-sequencing analysis is increasingly utilized to study gene expression in non-model organisms without sequenced genomes. Aethionema arabicum (Brassicaceae) exhibits seed dimorphism as a bet-hedging strategy - producing both a less dormant mucilaginous (M⁺) seed morph and a more dormant non-mucilaginous (NM) seed morph. Here, we compared de novo and reference-genome based transcriptome assemblies to investigate Ae. arabicum seed dimorphism and to evaluate the reference-free versus -dependent approach for identifying differentially expressed genes (DEGs).

RESULTS

A de novo transcriptome assembly was generated using sequences from M⁺ and NM Ae. arabicum dry seed morphs. The transcripts of the de novo assembly contained 63.1% complete Benchmarking Universal Single-Copy Orthologs (BUSCO) compared to 90.9% for the transcripts of the reference genome. DEG detection used the strict consensus of three methods (DESeq2, edgeR and NOISeq). Only 37% of 1533 differentially expressed de novo assembled transcripts paired with 1876 genome-derived DEGs. Gene Ontology (GO) terms distinguished the seed morphs: the terms translation and nucleosome assembly were overrepresented in DEGs higher in abundance in M⁺ dry seeds, whereas terms related to mRNA processing and transcription were overrepresented in DEGs higher in abundance in NM dry seeds. DEGs amongst these GO terms included ribosomal proteins and histones (higher in M⁺), RNA polymerase II subunits and related transcription and elongation factors (higher in NM). Expression of the inferred DEGs and other genes associated with seed maturation (e.g. those encoding late embryogenesis abundant proteins and transcription factors regulating seed development and maturation such as ABI3, FUS3, LEC1 and WRI1 homologs) were put in context with Arabidopsis thaliana seed maturation and indicated that M⁺ seeds may desiccate and mature faster than NM. The 1901 transcriptomic DEG set GO-terms had almost 90% overlap with the 2191 genome-derived DEG GO-terms.

CONCLUSIONS

Whilst there was only modest overlap of DEGs identified in reference-free versus -dependent approaches, the resulting GO analysis was concordant in both approaches. The identified differences in dry seed transcriptomes suggest mechanisms underpinning previously identified contrasts between morphology and germination behaviour of M⁺ and NM seeds.

Transcriptomic Analysis of Aedes aegypti in Response to Mosquitocidal Bacillus thuringiensis LLP29 Toxin

Globally, Aedes aegypti is one of the most dangerous mosquitoes that plays a crucial role as a vector for human diseases, such as yellow fever, dengue, and chikungunya. To identify (1) transcriptomic basis of midgut (2) key genes that are involved in the toxicity process by a comparative transcriptomic analysis between the control and Bacillus thuringiensis (Bt) toxin (LLP29 proteins)-treated groups. Next-generation sequencing technology was used to sequence the midgut transcriptome of A. aegypti. A total of 17130 unigenes, including 574 new unigenes, were identified containing 16358 (95.49%) unigenes that were functionally annotated. According to differentially expressed gene (DEG) analysis, 557 DEGs were annotated, including 226 upregulated and 231 downregulated unigenes in the Bt toxin-treated group. A total of 442 DEGs were functionally annotated; among these, 33 were specific to multidrug resistance, 6 were immune-system-related (Lectin, Defensin, Lysozyme), 28 were related to putative proteases, 7 were lipase-related, 8 were related to phosphatases, and 30 were related to other transporters. In addition, the relative expression of 28 DEGs was further confirmed through quantitative real time polymerase chain reaction. The results provide a transcriptomic basis for the identification and functional authentication of DEGs in A. aegypti.

A comprehensive analysis of the Korean fir (Abies koreana) genes expressed under heat stress using transcriptome analysis

Korean fir (Abies koreana), a rare species endemic to South Korea, is sensitive to climate change. Here, we used next-generation massively parallel sequencing technology and de novo transcriptome assembly to gain a comprehensive overview of the Korean fir transcriptome under heat stress. Sequencing control and heat-treated samples of Korean fir, we obtained more than 194,872,650 clean reads from each sample. After de novo assembly and quantitative assessment, 42,056 unigenes were generated with an average length of 908 bp. In total, 6,401 differentially expressed genes were detected, of which 2,958 were up-regulated and 3,443 down-regulated, between the heat-treated and control samples. A gene ontology analysis of these unigenes revealed heat-stress-related terms, such as "response to stimulus". Further, in depth analysis revealed 204 transcription factors and 189 Hsps as differentially expressed. Finally, 12 regulated candidate genes associated with heat stress were examined using quantitative real-time PCR (qRT-PCR). In this study, we present the first comprehensive characterisation of Korean fir subjected to heat stress using transcriptome analysis. It provides an important resource for future studies of Korean fir with the objective of identifying heat stress tolerant lines.

Omics approaches to study gene regulatory networks for development in echinoderms

Gene regulatory networks (GRNs) describe the interactions for a developmental process at a given time and space. Historically, perturbation experiments represent one of the key methods for analyzing and reconstructing a GRN, and the GRN governing early development in the sea urchin embryo stands as one of the more deeply dissected so far. As technology progresses, so do the methods used to address different biological questions. Next-generation sequencing (NGS) has become a standard experimental technique for genome and transcriptome sequencing and studies of protein-DNA interactions and DNA accessibility. While several efforts have been made toward the integration of different omics approaches for the study of the regulatory genome in many animals, in a few cases, these are applied with the purpose of reconstructing and experimentally testing developmental GRNs. Here, we review emerging approaches integrating multiple NGS technologies for the prediction and validation of gene interactions within echinoderm GRNs. These approaches can be applied to both 'model' and 'non-model' organisms. Although a number of issues still need to be addressed, advances in NGS applications, such as assay for transposase-accessible chromatin sequencing, combined with the availability of embryos belonging to different species, all separated by various evolutionary distances and accessible to experimental regulatory biology, place echinoderms in an unprecedented position for the reconstruction and evolutionary comparison of developmental GRNs. We conclude that sequencing technologies and integrated omics approaches allow the examination of GRNs on a genome-wide scale only if biological perturbation and cis-regulatory analyses are experimentally accessible, as in the case of echinoderm embryos.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

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

Fungal Genomes and Insights into the Evolution of the Kingdom

The kingdom Fungi comprises species that inhabit nearly all ecosystems. Fungi exist as both free-living and symbiotic unicellular and multicellular organisms with diverse morphologies. The genomes of fungi encode genes that enable them to thrive in diverse environments, invade plant and animal cells, and participate in nutrient cycling in terrestrial and aquatic ecosystems. The continuously expanding databases of fungal genome sequences have been generated by individual and large-scale efforts such as Génolevures, Broad Institute's Fungal Genome Initiative, and the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). These efforts have produced a catalog of fungal genes and genomic organization. The genomic datasets can be utilized to better understand how fungi have adapted to their lifestyles and ecological niches. Large datasets of fungal genomic and transcriptomic data have enabled the use of novel methodologies and improved the study of fungal evolution from a molecular sequence perspective. Combined with microscopes, petri dishes, and woodland forays, genome sequencing supports bioinformatics and comparative genomics approaches as important tools in the study of the biology and evolution of fungi.

Generation and characterization of the blood transcriptome of Macaca thibetana and comparative analysis with M. mulatta

Transcriptome profiles provide a large transcript sequence data set for genomic study, particularly in organisms that have no accurate genome data published. The Tibetan macaque (Macaca thibetana) is commonly considered to be an endemic species to China and an important animal in biomedical research in the present day. In the present study, we report the de novo assembly and characterization of the blood transcriptome of the Tibetan macaque from three individuals, and we also sequenced the blood transcriptome of the rhesus macaque (Macaca mulatta) for comparison. Using RNA-seq technology, 138 million sequencing reads of the M. thibetana transcriptome were generated. The assembly included 327 871 transcripts with an N50 of 1571 bp. According to the sequence similarity search, 80 317 (24.5%) transcripts were annotated in the nr protein database. All transcripts from M. thibetana and M. mulatta were functionally classified and compared using GO and KEGG analyses. The two transcriptomes were different in the GO term of nutrient reservoir activity, and in the KEGG subcategories of signaling molecules and interaction, infectious diseases, cell growth and death, and immune system. The transcriptomes in this study would provide a valuable resource for future functional and comparative genomic studies, and even for biological studies of this non-human primate.

Comparative transcriptome analyses reveal the genetic basis underlying the immune function of three amphibians' skin

Skin as the first barrier against external invasions plays an essential role for the survival of amphibians on land. Understanding the genetic basis of skin function is significant in revealing the mechanisms underlying immunity of amphibians. In this study, we de novo sequenced and comparatively analyzed skin transcriptomes from three different amphibian species, Andrias davidianus, Bufo gargarizans, and Rana nigromaculata Hallowell. Functional classification of unigenes in each amphibian showed high accordance, with the most represented GO terms and KEGG pathways related to basic biological processes, such as binding and metabolism and immune system. As for the unigenes, GO and KEGG distributions of conserved orthologs in each species were similar, with the predominantly enriched pathways including RNA polymerase, nucleotide metabolism, and defense. The positively selected orthologs in each amphibian were also similar, which were primarily involved in stimulus response, cell metabolic, membrane, and catalytic activity. Furthermore, a total of 50 antimicrobial peptides from 26 different categories were identified in the three amphibians, and one of these showed high efficiency in inhibiting the growth of different bacteria. Our understanding of innate immune function of amphibian skin has increased basis on the immune-related unigenes, pathways, and antimicrobial peptides in amphibians.

De novo assembly of Phlomis purpurea after challenging with Phytophthora cinnamomi

BACKGROUND

Phlomis plants are a source of biological active substances with potential applications in the control of phytopathogens. Phlomis purpurea (Lamiaceae) is autochthonous of southern Iberian Peninsula and Morocco and was found to be resistant to Phytophthora cinnamomi. Phlomis purpurea has revealed antagonistic effect in the rhizosphere of Quercus suber and Q. ilex against P. cinnamomi. Phlomis purpurea roots produce bioactive compounds exhibiting antitumor and anti-Phytophthora activities with potential to protect susceptible plants. Although these important capacities of P. purpurea have been demonstrated, there is no transcriptomic or genomic information available in public databases that could bring insights on the genes underlying this anti-oomycete activity.

RESULTS

Using Illumina technology we obtained a de novo assembly of P. purpurea transcriptome and differential transcript abundance to identify putative defence related genes in challenged versus non-challenged plants. A total of 1,272,600,000 reads from 18 cDNA libraries were merged and assembled into 215,739 transcript contigs. BLASTX alignment to Nr NCBI database identified 124,386 unique annotated transcripts (57.7%) with significant hits. Functional annotation identified 83,550 out of 124,386 unique transcripts, which were mapped to 141 pathways. 39% of unigenes were assigned GO terms. Their functions cover biological processes, cellular component and molecular functions. Genes associated with response to stimuli, cellular and primary metabolic processes, catalytic and transporter functions were among those identified. Differential transcript abundance analysis using DESeq revealed significant differences among libraries depending on post-challenge times. Comparative cyto-histological studies of P. purpurea roots challenged with P. cinnamomi zoospores and controls revealed specific morphological features (exodermal strips and epi-cuticular layer), that may provide a constitutive efficient barrier against pathogen penetration. Genes involved in cutin biosynthesis and in exodermal Casparian strips formation were up-regulated.

CONCLUSIONS

The de novo assembly of transcriptome using short reads for a non-model plant, P. purpurea, revealed many unique transcripts useful for further gene expression, biological function, genomics and functional genomics studies. The data presented suggest a combination of a constitutive resistance and an increased transcriptional response from P. purpurea when challenged with the pathogen. This knowledge opens new perspectives for the understanding of defence responses underlying pathogenic oomycete/plant interaction upon challenge with P. cinnamomi.

The Euphausia superba transcriptome database, SuperbaSE: An online, open resource for researchers

Antarctic krill (Euphausia superba) is a crucial component of the Southern Ocean ecosystem, acting as the major link between primary production and higher trophic levels with an annual predator demand of up to 470 million tonnes. It also acts as an ecosystem engineer, affecting carbon sequestration and recycling iron and nitrogen, and has increasing importance as a commercial product in the aquaculture and health industries. Here we describe the creation of a de novo assembled head transcriptome for E. superba. As an example of its potential as a molecular resource, we relate its exploitation in identifying and characterizing numerous genes related to the circadian clock in E. superba, including the major components of the central feedback loop. We have made the transcriptome openly accessible for a wider audience of ecologists, molecular biologists, evolutionary geneticists, and others in a user-friendly format at SuperbaSE, hosted at http://www.krill.le.ac.uk.

The Peptidoglycan Recognition Proteins PGRPLA and PGRPLB Regulate Anopheles Immunity to Bacteria and Affect Infection by Plasmodium

Peptidoglycan recognition proteins (PGRPs) form a family of immune regulators that is conserved from insects to mammals. In the malaria vector mosquito Anophelescoluzzii, the peptidoglycan receptor PGRPLC activates the immune-deficiency (Imd) pathway limiting both the microbiota load and Plasmodium infection. Here, we carried out an RNA interference screen to examine the role of all 7 Anopheles PGRPs in infections with Plasmodium berghei and P. falciparum. We show that, in addition to PGRPLC, PGRPLA and PGRPS2/PGRPS3 also participate in antiparasitic defenses, and that PGRPLB promotes mosquito permissiveness to P. falciparum. We also demonstrate that following a mosquito blood feeding, which promotes growth of the gut microbiota, PGRPLA and PGRPLB positively and negatively regulate the activation of the Imd pathway, respectively. Our data demonstrate that PGRPs are important regulators of the mosquito epithelial immunity and vector competence.

Investigating right ovary degeneration in chick embryos by transcriptome sequencing

In asymmetric chick gonads, the left and right female gonads undergo distinct programs during development, generating a functional ovary on the left side only. Despite some progress being made in recent years, the mechanisms of molecular regulation remain incompletely understood, and little genomic information is available regarding the degeneration of the right ovary in the chick embryo testis. In this study, we performed transcriptome sequencing to investigate differentially expressed genes in the left and right ovaries and gene functions at two critical time points; embryonic days 6 (E6) and 10 (E10). Using high-throughput RNA-sequencing technologies, 539 and 1046 genes were identified as being significantly differentially expressed between 6R-VS-6L and 10R-VS-10L. Gene ontology analysis of the differentially expressed genes revealed enrichment in functional pathways. Among these, candidate genes associated with degeneration of the right ovary in the chick embryo were identified. Identification of a pathway involved in ovarian degeneration provides an important resource for the further study of its molecular mechanisms and functions.

Identification of Development-Related Genes in the Ovaries of Adult Harmonia axyridis (Pallas) Lady Beetles Using a Time- Series Analysis by RNA-seq

Adults of the lady beetle species Harmonia axyridis (Pallas) are bred artificially en masse for classic biological control, which requires egg-laying by the H. axyridis ovary. Development-related genes may impact the growth of the H. axyridis adult ovary but have not been reported. Here, we used integrative time-series RNA-seq analysis of the ovary in H. axyridis adults to detect development-related genes. A total of 28,558 unigenes were functionally annotated using seven types of databases to obtain an annotated unigene database for ovaries in H. axyridis adults. We also analysed differentially expressed genes (DEGs) between samples. Based on a combination of the results of this bioinformatics analysis with literature reports and gene expression level changes in four different stages, we focused on the development of oocyte reproductive stem cell and yolk formation process and identified 26 genes with high similarity to development-related genes. 20 DEGs were randomly chosen for quantitative real-time PCR (qRT-PCR) to validate the accuracy of the RNA-seq results. This study establishes a robust pipeline for the discovery of key genes using high-throughput sequencing and the identification of a class of development-related genes for characterization.

Advances in genetics and genomics: use and limitations in achieving malaria elimination goals

Success of the global research agenda towards eradication of malaria will depend on the development of new tools, including drugs, vaccines, insecticides and diagnostics. Genetic and genomic information now available for the malaria parasites, their mosquito vectors and human host, can be harnessed to both develop these tools and monitor their effectiveness. Here we review and provide specific examples of current technological advances and how these genetic and genomic tools have increased our knowledge of host, parasite and vector biology in relation to malaria elimination and in turn enhanced the potential to reach that goal. We then discuss limitations of these tools and future prospects for the successful achievement of global malaria elimination goals.

Pyrosequencing the Midgut Transcriptome of the Banana Weevil Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) Reveals Multiple Protease-Like Transcripts

The banana weevil Cosmopolites sordidus is an important and serious insect pest in most banana and plantain-growing areas of the world. In spite of the economic importance of this insect pest very little genomic and transcriptomic information exists for this species. In the present study, we characterized the midgut transcriptome of C. sordidus using massive 454-pyrosequencing. We generated over 590,000 sequencing reads that assembled into 30,840 contigs with more than 400 bp, representing a significant expansion of existing sequences available for this insect pest. Among them, 16,427 contigs contained one or more GO terms. In addition, 15,263 contigs were assigned an EC number. In-depth transcriptome analysis identified genes potentially involved in insecticide resistance, peritrophic membrane biosynthesis, immunity-related function and defense against pathogens, and Bacillus thuringiensis toxins binding proteins as well as multiple enzymes involved with protein digestion. This transcriptome will provide a valuable resource for understanding larval physiology and for identifying novel target sites and management approaches for this important insect pest.

Transcriptome Analysis and Discovery of Genes Relevant to Development in Bradysia odoriphaga at Three Developmental Stages

Bradysia odoriphaga (Diptera: Sciaridae) is the most important pest of Chinese chive (Allium tuberosum) in Asia; however, the molecular genetics are poorly understood. To explore the molecular biological mechanism of development, Illumina sequencing and de novo assembly were performed in the third-instar, fourth-instar, and pupal B. odoriphaga. The study resulted in 16.2 Gb of clean data and 47,578 unigenes (≥125 bp) contained in 7,632,430 contigs, 46.21% of which were annotated from non-redundant protein (NR), Gene Ontology (GO), Clusters of Orthologous Groups (COG), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. It was found that 19.67% of unigenes matched the homologous species mainly, including Aedes aegypti, Culex quinquefasciatus, Ceratitis capitata, and Anopheles gambiae. According to differentially expressed gene (DEG) analysis, 143, 490, and 309 DEGs were annotated as involved in the developmental process in the GO database respectively, in the comparisons of third-instar and fourth-instar larvae, third-instar larvae and pupae, and fourth-instar larvae and pupae. Twenty-five genes were closely related to these processes, including developmental process, reproduction process, and reproductive organs development and programmed cell death (PCD). The information of unigenes assembled in B. odoriphaga through transcriptome and DEG analyses could provide a detailed genetic basis and regulated information for elaborating the developmental mechanism from the larval, pre-pupal to pupal stages of B. odoriphaga.

Evolution of a Heavy Metal Homeostasis/Resistance Island Reflects Increasing Copper Stress in Enterobacteria

Copper homeostasis in bacteria is challenged by periodic elevation of copper levels in the environment, arising from both natural sources and human inputs. Several mechanisms have evolved to efflux copper from bacterial cells, including thecus(copper sensing copper efflux system), andpco(plasmid-borne copper resistance system) systems. The genes belonging to these two systems can be physically clustered in a Copper Homeostasis and Silver Resistance Island (CHASRI) on both plasmids and chromosomes in Enterobacteria. Increasing use of copper in agricultural and industrial applications raises questions about the role of human activity in the evolution of novel copper resistance mechanisms. Here we present evidence that CHASRI emerged and diversified in response to copper deposition across aerobic and anaerobic environments. An analysis of diversification rates and a molecular clock model suggest that CHASRI experienced repeated episodes of elevated diversification that could correspond to peaks in human copper production. Phylogenetic analyses suggest that CHASRI originated in a relative ofEnterobacter cloacaeas the ultimate product of sequential assembly of several pre-existing two-gene modules. Once assembled, CHASRI dispersed via horizontal gene transfer within Enterobacteriaceae and also to certain members of Shewanellaceae, where the originalpcomodule was replaced by a divergentpcohomolog. Analyses of copper stress mitigation suggest that CHASRI confers increased resistance aerobically, anaerobically, and during shifts between aerobic and anaerobic environments, which could explain its persistence in facultative anaerobes and emergent enteric pathogens.

Deep sequencing revealed molecular signature of horizontal gene transfer of plant like transcripts in the mosquito Anopheles culicifacies: an evolutionary puzzle

In prokaryotes, horizontal gene transfer (HGT) has been regarded as an important evolutionary drive to acquire and retain beneficial genes for their survival in diverse ecologies. However, in eukaryotes, the functional role of HGTs remains questionable, although current genomic tools are providing increased evidence of acquisition of novel traits within non-mating metazoan species. Here, we provide another transcriptomic evidence for the acquisition of massive plant genes in the mosquito, Anopheles culicifacies. Our multiple experimental validations including genomic PCR, RT-PCR, real-time PCR, immuno-blotting and immuno-florescence microscopy, confirmed that plant like transcripts (PLTs) are of mosquito origin and may encode functional proteins. A comprehensive molecular analysis of the PLTs and ongoing metagenomic analysis of salivary microbiome provide initial clues that mosquitoes may have survival benefits through the acquisition of nuclear as well as chloroplast encoded plant genes. Our findings of PLTs further support the similar questionable observation of HGTs in other higher organisms, which is still a controversial and debatable issue in the community of evolutionists. We believe future understanding of the underlying mechanism of the feeding associated molecular responses may shed new insights in the functional role of PLTs in the mosquito.

De novo transcriptome sequencing and analysis of Coccinella septempunctata L. in non-diapause, diapause and diapause-terminated states to identify diapause-associated genes

BACKGROUND

The most common ladybird beetle, Coccinella septempunctata L., is an excellent predator of crop pests such as aphids and white flies, and it shows a wide range of adaptability, a large appetite and a high reproductive ability. Diapause research plays an important role in the artificial propagation and shelf-life extension of insect products. Although this lady beetle's regulatory, physiological and biochemical characteristics in the diapause period are well understood, the molecular mechanism of diapause remains unknown. Therefore, we collected female adults in three different states, i.e., non-diapause, diapause and diapause termination, for transcriptome sequencing.

RESULTS

After transcriptome sequencing using the Illumina HiSeq 2500 platform with pretreatment, a total of 417.6 million clean reads from nine samples were filtered using the program FASTX (version 0.0). Additionally, 106,262 contigs were assembled into 82,820 unigenes with an average length of 921 bp and an N50 of 1,241 bp. All of the unigenes were annotated through BLASTX alignment against the Nr or UniProt database, and 37,872 unigenes were matched. We performed further analysis of these unigenes using the Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Through pairwise comparisons of the non-diapause (ND), diapause (D), and diapause-terminated (DT) groups, 3,501 and 1,427 differentially expressed genes (DEGs) were identified between D and ND and between DT and D, respectively. Moreover, 443 of the DEGs were specifically expressed during the diapause period (i.e., DEGs that were expressed at the highest or lowest levels during diapause compared with the other stages). GO function and KEGG pathway enrichment were performed on all DEGs and showed that RNA-directed DNA polymerase activity and fatty acid metabolism were significantly affected. Furthermore, eight specific expressed genes were selected for validation using qRT-PCR. Among these eight genes, seven genes were up-regulated, and one gene was down-regulated; the change trends of the eight genes were the same between the qRT-PCR and RNA-seq analysis results.

CONCLUSIONS

In this study, a new method for collecting and identifying diapause insects was described. We generated a vast quantity of transcriptome data from C. septempunctata L., providing a resource for gene function research. The diapause-associated genes that we identified establish a foundation for future studies on the molecular mechanisms of diapause.

De novo transcriptome analysis of Inonotus baumii by RNA-seq

Inonotus baumii, a basidiomycete white rot fungus, has been widely used as traditional herbal medicine in China, Korea, Japan and other Asian countries for many years. Its extract is of great medicinal importance and plays a valuable role in the immune response and disease resistance. However, limited genetic resources for I. baumii have hindered exploration of this species. In order to gain a molecular understanding of this fungus, Illumina high-throughput technology was used to sequence and analyze the transcriptome of I. baumii, and 280,691 contigs, 43,890 scaffolds and 30,051 unigenes were obtained. Additionally, based on similarity search with known proteins, unigenes were annotated with gene descriptions, gene ontology (GO), clusters of orthologous group (COG), and database of protein families (Pfam) terms. According to the annotation of unigenes, a total of 12 candidate genes involved in the triterpenoid biosynthesis pathway and 21 putative FOLymes (fungal oxidative lignin enzymes) and 176 CAZymes (carbohydrate-active enzymes) were obtained using homology-based BlastX. Moreover, for better understanding of the transcripts function, the BlastX algorithm was used to search for homologous sequences against the Yeast genome. This is the first study on transcriptome analyses of I. baumii, which provided a dataset for functional gene mining and laid a basis for further functional genomics studies of I. baumii.

Assembly and Analysis of Differential Transcriptome Responses of Hevea brasiliensis on Interaction with Microcyclus ulei

Natural rubber (Hevea brasiliensis) is a tropical tree used commercially for the production of latex, from which 40,000 products are generated. The fungus Microcyclus ulei infects this tree, causing South American leaf blight (SALB) disease. This disease causes developmental delays and significant crop losses, thereby decreasing the production of latex. Currently several groups are working on obtaining clones of rubber tree with durable resistance to SALB through the use of extensive molecular biology techniques. In this study, we used a secondary clone that was resistant to M. ulei isolate GCL012. This clone, FX 3864 was obtained by crossing between clones PB 86 and B 38 (H. brasiliensis x H. brasiliensis). RNA-Seq high-throughput sequencing technology was used to analyze the differential expression of the FX 3864 clone transcriptome at 0 and 48 h post infection (hpi) with the M. ulei isolate GCL012. A total of 158,134,220 reads were assembled using the de novo assembly strategy to generate 90,775 contigs with an N50 of 1672. Using a reference-based assembly, 76,278 contigs were generated with an N50 of 1324. We identified 86 differentially expressed genes associated with the defense response of FX 3864 to GCL012. Seven putative genes members of the AP2/ERF ethylene (ET)-dependent superfamily were found to be down-regulated. An increase in salicylic acid (SA) was associated with the up-regulation of three genes involved in cell wall synthesis and remodeling, as well as in the down-regulation of the putative gene CPR5. The defense response of FX 3864 against the GCL012 isolate was associated with the antagonistic SA, ET and jasmonic acid (JA) pathways. These responses are characteristic of plant resistance to biotrophic pathogens.

Unraveling dual feeding associated molecular complexity of salivary glands in the mosquito Anopheles culicifacies

Mosquito salivary glands are well known to facilitate meal acquisition, however the fundamental question on how adult female salivary gland manages molecular responses during sugar versus blood meal uptake remains unanswered. To investigate these responses, we analyzed a total of 58.5 million raw reads generated from two independent RNAseq libraries of the salivary glands collected from 3–4 day-old sugar and blood fed Anopheles culicifacies mosquitoes. Comprehensive functional annotation analysis of 10,931 contigs unraveled that salivary glands may encode diverse nature of proteins in response to distinct physiological feeding status. Digital gene expression analysis and PCR validation indicated that first blood meal significantly alters the molecular architecture of the salivary glands. Comparative microscopic analysis also revealed that first blood meal uptake not only causes an alteration of at least 12–22% of morphological features of the salivary glands but also results in cellular changes e.g. apoptosis, confirming together that adult female salivary glands are specialized organs to manage meal specific responses. Unraveling the underlying mechanism of mosquito salivary gene expression, controlling dual feeding associated responses may provide a new opportunity to control vector borne diseases.

RNA-Seq based phylogeny recapitulates previous phylogeny of the genus Flaveria (Asteraceae) with some modifications

BACKGROUND

The genus Flaveria has been extensively used as a model to study the evolution of C4 photosynthesis as it contains C3 and C4 species as well as a number of species that exhibit intermediate types of photosynthesis. The current phylogenetic tree of the genus Flaveria contains 21 of the 23 known Flaveria species and has been previously constructed using a combination of morphological data and three non-coding DNA sequences (nuclear encoded ETS, ITS and chloroplast encoded trnL-F).

RESULTS

Here we developed a new strategy to update the phylogenetic tree of 16 Flaveria species based on RNA-Seq data. The updated phylogeny is largely congruent with the previously published tree but with some modifications. We propose that the data collection method provided in this study can be used as a generic method for phylogenetic tree reconstruction if the target species has no genomic information. We also showed that a "F. pringlei" genotype recently used in a number of labs may be a hybrid between F. pringlei (C3) and F. angustifolia (C3-C4).

CONCLUSIONS

We propose that the new strategy of obtaining phylogenetic sequences outlined in this study can be used to construct robust trees in a larger number of taxa. The updated Flaveria phylogenetic tree also supports a hypothesis of stepwise and parallel evolution of C4 photosynthesis in the Flavaria clade.

Transcriptome Analysis Comparison of Lipid Biosynthesis in the Leaves and Developing Seeds of Brassica napus

Brassica napus seed is a lipid storage organ containing approximately 40% oil, while its leaves contain many kinds of lipids for many biological roles, but the overall amounts are less than in seeds. Thus, lipid biosynthesis in the developing seeds and the leaves is strictly regulated which results the final difference of lipids. However, there are few reports about the molecular mechanism controlling the difference in lipid biosynthesis between developing seeds and leaves. In this study, we tried to uncover this mechanism by analyzing the transcriptome data for lipid biosynthesis. The transcriptome data were de novo assembled and a total of 47,216 unigenes were obtained, which had an N50 length and median of 1271 and 755 bp, respectively. Among these unigenes, 36,368 (about 77.02%) were annotated and there were 109 up-regulated unigenes and 72 down-regulated unigenes in the developing seeds lipid synthetic pathway after comparing with leaves. In the oleic acid pathway, 23 unigenes were up-regulated and four unigenes were down-regulated. During triacylglycerol (TAG) synthesis, the key unigenes were all up-regulated, such as phosphatidate phosphatase and diacylglycerol O-acyltransferase. During palmitic acid, palmitoleic acid, stearic acid, linoleic acid and linolenic acid synthesis in leaves, the unigenes were nearly all up-regulated, which indicated that the biosynthesis of these particular fatty acids were more important in leaves. In the developing seeds, almost all the unigenes in the ABI3VP1, RKD, CPP, E2F-DP, GRF, JUMONJI, MYB-related, PHD and REM transcript factor families were up-regulated, which helped us to discern the regulation mechanism underlying lipid biosynthesis. The differential up/down-regulation of the genes and TFs involved in lipid biosynthesis in developing seeds and leaves provided direct evidence that allowed us to map the network that regulates lipid biosynthesis, and the identification of new TFs that are up-regulated in developing seeds will help us to further elucidate the lipids biosynthesis pathway in developing seeds and leaves.

The hub protein loquacious connects the microRNA and short interfering RNA pathways in mosquitoes

Aedes aegypti mosquitoes vector several arboviruses of global health significance, including dengue viruses and chikungunya virus. RNA interference (RNAi) plays an important role in antiviral immunity, gene regulation and protection from transposable elements. Double-stranded RNA binding proteins (dsRBPs) are important for efficient RNAi; in Drosophila functional specialization of the miRNA, endo-siRNA and exo-siRNA pathway is aided by the dsRBPs Loquacious (Loqs-PB, Loqs-PD) and R2D2, respectively. However, this functional specialization has not been investigated in other dipterans. We were unable to detect Loqs-PD in Ae. aegypti; analysis of other dipteran genomes demonstrated that this isoform is not conserved outside of Drosophila. Overexpression experiments and small RNA sequencing following depletion of each dsRBP revealed that R2D2 and Loqs-PA cooperate non-redundantly in siRNA production, and that these proteins exhibit an inhibitory effect on miRNA levels. Conversely, Loqs-PB alone interacted with mosquito dicer-1 and was essential for full miRNA production. Mosquito Loqs interacted with both argonaute 1 and 2 in a manner independent of its interactions with dicer. We conclude that the functional specialization of Loqs-PD in Drosophila is a recently derived trait, and that in other dipterans, including the medically important mosquitoes, Loqs-PA participates in both the miRNA and endo-siRNA based pathways.

Pathways associated with lignin biosynthesis in lignomaniac jute fibres

We generated the bast transcriptomes of a deficient lignified phloem fibre mutant and its wild-type jute (Corchorus capsularis) using Illumina paired-end sequencing. A total of 34,163 wild-type and 29,463 mutant unigenes, with average lengths of 1442 and 1136 bp, respectively, were assembled de novo, ~77-79 % of which were functionally annotated. These annotated unigenes were assigned to COG (~37-40 %) and GO (~22-28 %) classifications and mapped to 189 KEGG pathways (~19-21 %). We discovered 38 and 43 isoforms of 16 and 10 genes of the upstream shikimate-aromatic amino acid and downstream monolignol biosynthetic pathways, respectively, rendered their sequence similarities, confirmed the identities of 22 of these candidate gene families by phylogenetic analyses and reconstructed the pathway leading to lignin biosynthesis in jute fibres. We also identified major genes and bast-related transcription factors involved in secondary cell wall (SCW) formation. The quantitative RT-PCRs revealed that phenylalanine ammonia-lyase 1 (CcPAL1) was co-down-regulated with several genes of the upstream shikimate pathway in mutant bast tissues at an early growth stage, although its expression relapsed to the normal level at the later growth stage. However, cinnamyl alcohol dehydrogenase 7 (CcCAD7) was strongly down-regulated in mutant bast tissues irrespective of growth stages. CcCAD7 disruption at an early growth stage was accompanied by co-up-regulation of SCW-specific genes cellulose synthase A7 (CcCesA7) and fasciclin-like arabinogalactan 6 (CcFLA6), which was predicted to be involved in coordinating the S-layers' deposition in the xylan-type jute fibres. Our results identified CAD as a promising target for developing low-lignin jute fibres using genomics-assisted molecular approaches.

Inferring bona fide transfrags in RNA-Seq derived-transcriptome assemblies of non-model organisms

BACKGROUND

De novo transcriptome assembly of short transcribed fragments (transfrags) produced from sequencing-by-synthesis technologies often results in redundant datasets with differing levels of unassembled, partially assembled or mis-assembled transcripts. Post-assembly processing intended to reduce redundancy typically involves reassembly or clustering of assembled sequences. However, these approaches are mostly based on common word heuristics and often create clusters of biologically unrelated sequences, resulting in loss of unique transfrags annotations and propagation of mis-assemblies.

RESULTS

Here, we propose a structured framework that consists of a few steps in pipeline architecture for Inferring Functionally Relevant Assembly-derived Transcripts (IFRAT). IFRAT combines 1) removal of identical subsequences, 2) error tolerant CDS prediction, 3) identification of coding potential, and 4) complements BLAST with a multiple domain architecture annotation that reduces non-specific domain annotation. We demonstrate that independent of the assembler, IFRAT selects bona fide transfrags (with CDS and coding potential) from the transcriptome assembly of a model organism without relying on post-assembly clustering or reassembly. The robustness of IFRAT is inferred on RNA-Seq data of Neurospora crassa assembled using de Bruijn graph-based assemblers, in single (Trinity and Oases-25) and multiple (Oases-Merge and additive or pooled) k-mer modes. Single k-mer assemblies contained fewer transfrags compared to the multiple k-mer assemblies. However, Trinity identified a comparable number of predicted coding sequence and gene loci to Oases pooled assembly. IFRAT selects bona fide transfrags representing over 94% of cumulative BLAST-derived functional annotations of the unfiltered assemblies. Between 4-6% are lost when orphan transfrags are excluded and this represents only a tiny fraction of annotation derived from functional transference by sequence similarity. The median length of bona fide transfrags ranged from 1.5kb (Trinity) to 2kb (Oases), which is consistent with the average coding sequence length in fungi. The fraction of transfrags that could be associated with gene ontology terms ranged from 33-50%, which is also high for domain based annotation. We showed that unselected transfrags were mostly truncated and represent sequences from intronic, untranslated (5' and 3') regions and non-coding gene loci.

CONCLUSIONS

IFRAT simplifies post-assembly processing providing a reference transcriptome enriched with functionally relevant assembly-derived transcripts for non-model organism.

Development and validation of an rDNA operon based primer walking strategy applicable to de novo bacterial genome finishing

Advances in sequencing technology have drastically increased the depth and feasibility of bacterial genome sequencing. However, little information is available that details the specific techniques and procedures employed during genome sequencing despite the large numbers of published genomes. Shotgun approaches employed by second-generation sequencing platforms has necessitated the development of robust bioinformatics tools for in silico assembly, and complete assembly is limited by the presence of repetitive DNA sequences and multi-copy operons. Typically, re-sequencing with multiple platforms and laborious, targeted Sanger sequencing are employed to finish a draft bacterial genome. Here we describe a novel strategy based on the identification and targeted sequencing of repetitive rDNA operons to expedite bacterial genome assembly and finishing. Our strategy was validated by finishing the genome of Paenibacillus polymyxa strain CR1, a bacterium with potential in sustainable agriculture and bio-based processes. An analysis of the 38 contigs contained in the P. polymyxa strain CR1 draft genome revealed 12 repetitive rDNA operons with varied intragenic and flanking regions of variable length, unanimously located at contig boundaries and within contig gaps. These highly similar but not identical rDNA operons were experimentally verified and sequenced simultaneously with multiple, specially designed primer sets. This approach also identified and corrected significant sequence rearrangement generated during the initial in silico assembly of sequencing reads. Our approach reduces the required effort associated with blind primer walking for contig assembly, increasing both the speed and feasibility of genome finishing. Our study further reinforces the notion that repetitive DNA elements are major limiting factors for genome finishing. Moreover, we provided a step-by-step workflow for genome finishing, which may guide future bacterial genome finishing projects.

De novo assembly and analysis of crow lungs transcriptome

The jungle crow (Corvus macrorhynchos) belongs to the order Passeriformes of bird species and is important for avian ecological and evolutionary genetics studies. However, there is limited information on the transcriptome data of this species. In the present study, we report the characterization of the lung transcriptome of the jungle crow using GS FLX Titanium XLR70. Altogether, 1,510,303 high-quality sequence reads with 581,198,230 bases was de novo assembled into 22,169 isotigs (isotig represents an individual transcript) and 784,009 singletons. Using these isotigs and 581,681 length-filtered (greater than 300 bp) singletons, 20,010 unique protein-coding genes were identified by BLASTx comparison against a nonredundant (nr) protein sequence database. Comparative analysis revealed that 46,604 (70.29%) and 51,642 (72.48%) of the assembled transcripts have significant similarity to zebra finch and chicken RefSeq proteins, respectively. As determined by GO annotation and KEGG pathway mapping, functional annotation of the unigenes recovered diverse biological functions and processes. Transcripts putatively involved in the immune response were identified. Furthermore, 20,599 single nucleotide polymorphisms (SNPs) and 7525 simple sequence repeats (SSRs) were retrieved from the assembled transcript database. This resource should lay an important base for future ecological, evolutionary, and conservation genetic studies on this species and in other related species.

Continental-scale footprint of balancing and positive selection in a small rodent (Microtus arvalis)

Genetic adaptation to different environmental conditions is expected to lead to large differences between populations at selected loci, thus providing a signature of positive selection. Whereas balancing selection can maintain polymorphisms over long evolutionary periods and even geographic scale, thus leads to low levels of divergence between populations at selected loci. However, little is known about the relative importance of these two selective forces in shaping genomic diversity, partly due to difficulties in recognizing balancing selection in species showing low levels of differentiation. Here we address this problem by studying genomic diversity in the European common vole (Microtus arvalis) presenting high levels of differentiation between populations (average F_ST = 0.31). We studied 3,839 Amplified Fragment Length Polymorphism (AFLP) markers genotyped in 444 individuals from 21 populations distributed across the European continent and hence over different environmental conditions. Our statistical approach to detect markers under selection is based on a Bayesian method specifically developed for AFLP markers, which treats AFLPs as a nearly codominant marker system, and therefore has increased power to detect selection. The high number of screened populations allowed us to detect the signature of balancing selection across a large geographic area. We detected 33 markers potentially under balancing selection, hence strong evidence of stabilizing selection in 21 populations across Europe. However, our analyses identified four-times more markers (138) being under positive selection, and geographical patterns suggest that some of these markers are probably associated with alpine regions, which seem to have environmental conditions that favour adaptation. We conclude that despite favourable conditions in this study for the detection of balancing selection, this evolutionary force seems to play a relatively minor role in shaping the genomic diversity of the common vole, which is more influenced by positive selection and neutral processes like drift and demographic history.

Partition enrichment of nucleotide sequences (PINS)--a generally applicable, sequence based method for enrichment of complex DNA samples

The dwindling cost of DNA sequencing is driving transformative changes in various biological disciplines including medicine, thus resulting in an increased need for routine sequencing. Preparation of samples suitable for sequencing is the starting point of any practical application, but enrichment of the target sequence over background DNA is often laborious and of limited sensitivity thereby limiting the usefulness of sequencing. The present paper describes a new method, Probability directed Isolation of Nucleic acid Sequences (PINS), for enrichment of DNA, enabling the sequencing of a large DNA region surrounding a small known sequence. A 275,000 fold enrichment of a target DNA sample containing integrated human papilloma virus is demonstrated. Specifically, a sample containing 0.0028 copies of target sequence per ng of total DNA was enriched to 786 copies per ng. The starting concentration of 0.0028 target copies per ng corresponds to one copy of target in a background of 100,000 complete human genomes. The enriched sample was subsequently amplified using rapid genome walking and the resulting DNA sequence revealed not only the sequence of a the truncated virus, but also 1026 base pairs 5' and 50 base pairs 3' to the integration site in chromosome 8. The demonstrated enrichment method is extremely sensitive and selective and requires only minimal knowledge of the sequence to be enriched and will therefore enable sequencing where the target concentration relative to background is too low to allow the use of other sample preparation methods or where significant parts of the target sequence is unknown.

Assembly and annotation of a non-model gastropod (Nerita melanotragus) transcriptome: a comparison of de novo assemblers

Background

The sequencing, de novo assembly and annotation of transcriptome datasets generated with next generation sequencing (NGS) has enabled biologists to answer genomic questions in non-model species with unprecedented ease. Reliable and accurate de novo assembly and annotation of transcriptomes, however, is a critically important step for transcriptome assemblies generated from short read sequences. Typical benchmarks for assembly and annotation reliability have been performed with model species. To address the reliability and accuracy of de novo transcriptome assembly in non-model species, we generated an RNAseq dataset for an intertidal gastropod mollusc species, Nerita melanotragus, and compared the assembly produced by four different de novo transcriptome assemblers; Velvet, Oases, Geneious and Trinity, for a number of quality metrics and redundancy.

Results

Transcriptome sequencing on the Ion Torrent PGM™ produced 1,883,624 raw reads with a mean length of 133 base pairs (bp). Both the Trinity and Oases de novo assemblers produced the best assemblies based on all quality metrics including fewer contigs, increased N50 and average contig length and contigs of greater length. Overall the BLAST and annotation success of our assemblies was not high with only 15-19% of contigs assigned a putative function.

Conclusions

We believe that any improvement in annotation success of gastropod species will require more gastropod genome sequences, but in particular an increase in mollusc protein sequences in public databases. Overall, this paper demonstrates that reliable and accurate de novo transcriptome assemblies can be generated from short read sequencers with the right assembly algorithms.

In depth annotation of the Anopheles gambiae mosquito midgut transcriptome

Background

Genome sequencing of Anopheles gambiae was completed more than ten years ago and has accelerated research on malaria transmission. However, annotation needs to be refined and verified experimentally, as most predicted transcripts have been identified by comparative analysis with genomes from other species. The mosquito midgut—the first organ to interact with Plasmodium parasites—mounts effective antiplasmodial responses that limit parasite survival and disease transmission. High-throughput Illumina sequencing of the midgut transcriptome was used to identify new genes and transcripts, contributing to the refinement of An. gambiae genome annotation.

Results

We sequenced ~223 million reads from An. gambiae midgut cDNA libraries generated from susceptible (G3) and refractory (L35) mosquito strains. Mosquitoes were infected with either Plasmodium berghei or Plasmodium falciparum, and midguts were collected after the first or second Plasmodium infection. In total, 22,889 unique midgut transcript models were generated from both An. gambiae strain sequences combined, and 76% are potentially novel. Of these novel transcripts, 49.5% aligned with annotated genes and appear to be isoforms or pre-mRNAs of reference transcripts, while 50.5% mapped to regions between annotated genes and represent novel intergenic transcripts (NITs). Predicted models were validated for midgut expression using qRT-PCR and microarray analysis, and novel isoforms were confirmed by sequencing predicted intron-exon boundaries. Coding potential analysis revealed that 43% of total midgut transcripts appear to be long non-coding RNA (lncRNA), and functional annotation of NITs showed that 68% had no homology to current databases from other species. Reads were also analyzed using de novo assembly and predicted transcripts compared with genome mapping-based models. Finally, variant analysis of G3 and L35 midgut transcripts detected 160,742 variants with respect to the An. gambiae PEST genome, and 74% were new variants. Intergenic transcripts had a higher frequency of variation compared with non-intergenic transcripts.

Conclusion

This in-depth Illumina sequencing and assembly of the An. gambiae midgut transcriptome doubled the number of known transcripts and tripled the number of variants known in this mosquito species. It also revealed existence of a large number of lncRNA and opens new possibilities for investigating the biological function of many newly discovered transcripts.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-636) contains supplementary material, which is available to authorized users.

Transcriptome analysis of a petal anthocyanin polymorphism in the arctic mustard, Parrya nudicaulis

Angiosperms are renown for their diversity of flower colors. Often considered adaptations to pollinators, the most common underlying pigments, anthocyanins, are also involved in plants’ stress response. Although the anthocyanin biosynthetic pathway is well characterized across many angiosperms and is composed of a few candidate genes, the consequences of blocking this pathway and producing white flowers has not been investigated at the transcriptome scale. We take a transcriptome-wide approach to compare expression differences between purple and white petal buds in the arctic mustard, Parrya nudicaulis, to determine which genes’ expression are consistently correlated with flower color. Using mRNA-Seq and de novo transcriptome assembly, we assembled an average of 722 bp per gene (49.81% coding sequence based on the A. thaliana homolog) for 12,795 genes from the petal buds of a pair of purple and white samples. Our results correlate strongly with qRT-PCR analysis of nine candidate genes in the anthocyanin biosynthetic pathway where chalcone synthase has the greatest difference in expression between color morphs (P/W = ∼7×). Among the most consistently differentially expressed genes between purple and white samples, we found 3× more genes with higher expression in white petals than in purple petals. These include four unknown genes, two drought-response genes (CDSP32, ERD5), a cold-response gene (GR-RBP2), and a pathogen defense gene (DND1). Gene ontology analysis of the top 2% of genes with greater expression in white relative to purple petals revealed enrichment in genes associated with stress responses including cold, drought and pathogen defense. Unlike the uniform downregulation of chalcone synthase that may be directly involved in the loss of petal anthocyanins, the variable expression of several genes with greater expression in white petals suggest that the physiological and ecological consequences of having white petals may be microenvironment-dependent.

De novo transcriptome sequencing and sequence analysis of the malaria vector Anopheles sinensis (Diptera: Culicidae)

Background

Anopheles sinensis is the major malaria vector in China and Southeast Asia. Vector control is one of the most effective measures to prevent malaria transmission. However, there is little transcriptome information available for the malaria vector. To better understand the biological basis of malaria transmission and to develop novel and effective means of vector control, there is a need to build a transcriptome dataset for functional genomics analysis by large-scale RNA sequencing (RNA-seq).

Methods

To provide a more comprehensive and complete transcriptome of An. sinensis, eggs, larvae, pupae, male adults and female adults RNA were pooled together for cDNA preparation, sequenced using the Illumina paired-end sequencing technology and assembled into unigenes. These unigenes were then analyzed in their genome mapping, functional annotation, homology, codon usage bias and simple sequence repeats (SSRs).

Results

Approximately 51.6 million clean reads were obtained, trimmed, and assembled into 38,504 unigenes with an average length of 571 bp, an N50 of 711 bp, and an average GC content 51.26%. Among them, 98.4% of unigenes could be mapped onto the reference genome, and 69% of unigenes could be annotated with known biological functions. Homology analysis identified certain numbers of An. sinensis unigenes that showed homology or being putative 1:1 orthologues with genomes of other Dipteran species. Codon usage bias was analyzed and 1,904 SSRs were detected, which will provide effective molecular markers for the population genetics of this species.

Conclusions

Our data and analysis provide the most comprehensive transcriptomic resource and characteristics currently available for An. sinensis, and will facilitate genetic, genomic studies, and further vector control of An. sinensis.

Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis

BACKGROUND

Anopheles sinensis is a major malaria vector in China and other Southeast Asian countries, and it is becoming increasingly resistant to the insecticides used for agriculture, net impregnation, and indoor residual spray. Very limited genomic information on this species is available, which has hindered the development of new tools for resistance surveillance and vector control. We used the 454 GS FLX system and generated expressed sequence tag (EST) databases of various life stages of An. sinensis, and we determined the transcriptional differences between deltamethrin resistant and susceptible mosquitoes.

RESULTS

The 454 GS FLX transcriptome sequencing yielded a total of 624,559 reads (average length of 290 bp) with the pooled An. sinensis mosquitoes across various development stages. The de novo assembly generated 33,411 contigs with average length of 493 bp. A total of 8,057 ESTs were generated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation. A total of 2,131 ESTs were differentially expressed between deltamethrin resistant and susceptible mosquitoes collected from the same field site in Jiangsu, China. Among these differentially expressed ESTs, a total of 294 pathways were mapped to the KEGG database, with the predominant ESTs belonging to metabolic pathways. Furthermore, a total of 2,408 microsatellites and 15,496 single nucleotide polymorphisms (SNPs) were identified.

CONCLUSIONS

The annotated EST and transcriptome databases provide a valuable genomic resource for further genetic studies of this important malaria vector species. The differentially expressed ESTs associated with insecticide resistance identified in this study lay an important foundation for further functional analysis. The identified microsatellite and SNP markers will provide useful tools for future population genetic and comparative genomic analyses of malaria vectors.

RNA sequencing analysis of the gametophyte transcriptome from the liverwort, Marchantia polymorpha

The liverwort Marchantia polymorpha is a member of the most basal lineage of land plants (embryophytes) and likely retains many ancestral morphological, physiological and molecular characteristics. Despite its phylogenetic importance and the availability of previous EST studies, M. polymorpha's lack of economic importance limits accessible genomic resources for this species. We employed Illumina RNA-Seq technology to sequence the gametophyte transcriptome of M. polymorpha. cDNA libraries from 6 different male and female developmental tissues were sequenced to delineate a global view of the M. polymorpha transcriptome. Approximately 80 million short reads were obtained and assembled into a non-redundant set of 46,533 transcripts (> =  200 bp) from 46,070 loci. The average length and the N50 length of the transcripts were 757 bp and 471 bp, respectively. Sequence comparison of assembled transcripts with non-redundant proteins from embryophytes resulted in the annotation of 43% of the transcripts. The transcripts were also compared with M. polymorpha expressed sequence tags (ESTs), and approximately 69.5% of the transcripts appeared to be novel. Twenty-one percent of the transcripts were assigned GO terms to improve annotation. In addition, 6,112 simple sequence repeats (SSRs) were identified as potential molecular markers, which may be useful in studies of genetic diversity. A comparative genomics approach revealed that a substantial proportion of the genes (35.5%) expressed in M. polymorpha were conserved across phylogenetically related species, such as Selaginella and Physcomitrella, and identified 580 genes that are potentially unique to liverworts. Our study presents an extensive amount of novel sequence information for M. polymorpha. This information will serve as a valuable genomics resource for further molecular, developmental and comparative evolutionary studies, as well as for the isolation and characterization of functional genes that are involved in sex differentiation and sexual reproduction in this liverwort.