Normalization of full-length enriched cDNA

A multi-omics integrative network map of maize

Networks are powerful tools to uncover functional roles of genes in phenotypic variation at a system-wide scale. Here, we constructed a maize network map that contains the genomic, transcriptomic, translatomic and proteomic networks across maize development. This map comprises over 2.8 million edges in more than 1,400 functional subnetworks, demonstrating an extensive network divergence of duplicated genes. We applied this map to identify factors regulating flowering time and identified 2,651 genes enriched in eight subnetworks. We validated the functions of 20 genes, including 18 with previously unknown connections to flowering time in maize. Furthermore, we uncovered a flowering pathway involving histone modification. The multi-omics integrative network map illustrates the principles of how molecular networks connect different types of genes and potential pathways to map a genome-wide functional landscape in maize, which should be applicable in a wide range of species.

Identification of Olfactory Genes in Monochamus saltuarius and Effects of Bursaphelenchus xylophilus Infestation on Their Expression

The pine wood nematode (PWN) Bursaphelenchus xylophilus has caused disastrous losses of pine forests in many countries, and the success of PWN depends strongly on interactions with its insect vectors. Monochamus saltuarius is a newly recorded vector in Northeast China. Feeding (i.e., immature) and egg-laying (i.e., mature) Monochamus spp. target different host plants, and olfactory cues play important roles regarding host choice. Whether infestation with PWN affects olfactory mechanisms in M. saltuarius related to feeding and oviposition is of interest as this may affect the spread of nematodes to new healthy hosts. However, little is known about molecular mechanisms of the olfactory system of M. saltuarius. We identified chemosensory-related genes in adult M. saltuarius and examined the influence of B. xylophilus on the respective expression patterns. Fifty-three odorant-binding proteins (OBPs), 15 chemosensory proteins, 15 olfactory receptors (ORs), 10 gustatory receptors, 22 ionotropic receptors (IRs), and two sensory neuron membrane proteins were identified, and sex bias among non-infested beetles was mainly found with respect to expression of OBPs. Interestingly, OBPs and ORs were markedly down-regulated in male M. saltuarius infested with B. xylophilus, which may reduce olfactory sensitivity of male M. saltuarius and affect the spreading of B. xylophilus to new hosts. Our results will help understand the interactions between B. xylophilus and M. saltuarius, which may lead to the identification of new control targets in the olfactory system of M. saltuarius.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

Distinct Transcriptomic Reprogramming in the Wheat Stripe Rust Fungus During the Initial Infection of Wheat and Barberry

Puccinia striiformis f. sp. tritici is the causal agent of wheat stripe rust that causes severe yield losses all over the world. As a macrocyclic heteroecious rust fungus, it is able to infect two unrelated host plants, wheat and barberry. Its urediniospores infect wheat and cause disease epidemic, while its basidiospores parasitize barberry to fulfill the sexual reproduction. This complex life cycle poses interesting questions on the different mechanisms of pathogenesis underlying the infection of the two different hosts. In the present study, transcriptomes of P. striiformis f. sp. tritici during the initial infection of wheat and barberry leaves were qualitatively and quantitatively compared. As a result, 142 wheat-specifically expressed genes (WEGs) were identified, which was far less than the 2,677 barberry-specifically expressed genes (BEGs). A larger proportion of evolutionarily conserved genes were observed in BEGs than that in WEGs, implying a longer history of the interaction between P. striiformis f. sp. tritici and barberry. Additionally, P. striiformis f. sp. tritici differentially expressed genes (DEGs) between wheat at 1 and 2 days postinoculation (dpi) and barberry at 3 and 4 dpi were identified by quantitative analysis. Gene Ontology analysis of these DEGs and expression patterns of P. striiformis f. sp. tritici pathogenic genes, including those encoding candidate secreted effectors, cell wall-degrading enzymes, and nutrient transporters, demonstrated that urediniospores and basidiospores exploited distinct strategies to overcome host defense systems. These results represent the first analysis of the P. striiformis f. sp. tritici transcriptome in barberry and contribute to a better understanding of the evolutionary processes and strategies of different types of rust spores during the infection process on different hosts.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Characteristics of Chemosensory and Opsin Genes in Newly Emerged and Sexually Mature Agrilus planipennis, an Important Quarantine Forest Beetle

The emerald ash borer (EAB), Agrilus planipennis, is a highly destructive quarantine pest. The olfactory and visual systems of A. planipennis play different but critical roles at newly emerged and sexually mature stages; however, the molecular basis underlying these differences remain unclear. Consequently, based on deep transcriptome sequencing, we evaluated the expression levels of chemosensory-related proteins and opsins at the two developmental stages of A. planipennis. We found 15 new chemosensory-related genes in our transcriptome assembly compared with the previous genome assembly, including 6 that code for odorant-binding proteins (OBPs) and 9 for chemosensory proteins (CSPs). The expression of several chemosensory-related genes (OBP7, OBP10, CSP1, and CSP12) differed markedly between newly emerged and sexually mature A. planipennis. We also found that the expression of UV opsin 2 and LW opsin 1 was higher in sexually mature male A. planipennis, which may be associated with their strong visual mate detection ability. This study forms the basis for further investigation of the chemosensory and visual system of A. planipennis, and these differentially expressed genes between newly emerged and sexually mature stages may serve as targets for the management of this destructive forest pest after sexual maturity.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Dynamic Changes in Chemosensory Gene Expression during the Dendrolimus punctatus Mating Process

The insect chemosensory system is pivotal for interactions with their environments, and moths have especially sensitive olfaction. Exploration of the connection between the plasticity of olfactory-guided and molecular level pathways in insects is important for understanding the olfactory recognition mechanisms of insects. The pine caterpillar moth, Dendrolimus punctatus Walker, is a dominant conifer defoliator in China, and mating is the priority for adults of this species, during which sex pheromone recognition and oviposition site location are the main activities; these activities are all closely related to chemosensory genes. Thus, we aimed to identify chemosensory related genes and monitor the spectrum of their dynamic expression during the entire mating process in D. punctatus. In this study, we generated transcriptome data from male and female adult D. punctatus specimens at four mating stages: eclosion, calling, copulation, and post-coitum. These data were analyzed using bioinformatics tools to identify the major olfactory-related gene families and determine their expression patterns during mating. Levels of odorant binding proteins (OBPs), chemosensory proteins (CSPs), and odorant receptors (ORs) were closely correlated with mating behavior. Comparison with ORs from other Dendrolimus and Lepidoptera species led to the discovery of a group of ORs specific to Dendrolimus. Furthermore, we identified several genes encoding OBPs and ORs that were upregulated after mating in females; these genes may mediate the location of host plants for oviposition via plant-emitted volatiles. This work will facilitate functional research into D. punctatus chemosensory genes, provide information about the relationship between chemosensory genes and important physiological activities, and promote research into the mechanisms underlying insect olfactory recognition.

Identification and Expression Profiling of Chemosensory Genes in Dendrolimus punctatus Walker

Dendrolimus punctatus Walker is a serious pest affecting conifers in southern China. As extensive pesticide spraying is currently required to control D. punctatus, new control strategies are urgently needed. Chemosensory genes represent potential molecular targets for development of alternative pest control strategies, and the expression characteristics of these genes provide an indication of their function. To date, little information is available regarding chemosensory genes in D. punctatus or their expression profiles at different development stages and in various tissues. Here, we assembled and analyzed the transcriptomes of D. punctatus collected at different developmental stages and in a range of organs, using next-generation sequencing. A total of 171 putative chemosensory genes were identified, encoding 53 odorant binding proteins, 26 chemosensory proteins, 60 odorant receptors (OR), 12 gustatory receptors (GR), 18 ionotropic receptors (IR), and 2 sensory neuron membrane proteins (SNMPs). Expression analysis indicated that the antennae possess the largest number of highly expressed olfactory genes and that olfactory gene expression patterns in the eggs, larvae, and head were similar to one another, with each having moderate numbers of highly expressed olfactory genes. Fat body, ovary, midgut, and testis tissues also had similar olfactory gene expression patterns, including few highly expressed olfactory genes. Of particular note, we identified only two pheromone binding proteins and no pheromone receptors in D. punctatus, similar to our previous findings in Dendrolimus houi and Dendrolimus kikuchii, suggesting that insects of the Dendrolimus genus have different pheromone recognition characteristics to other Lepidopteran insects. Overall, this extensive expression profile analysis provides a clear map of D. punctatus chemosensory genes, and will facilitate functional studies and the development of new pest control methods in the future.

Automated multiplex genome-scale engineering in yeast

Genome-scale engineering is indispensable in understanding and engineering microorganisms, but the current tools are mainly limited to bacterial systems. Here we report an automated platform for multiplex genome-scale engineering in Saccharomyces cerevisiae, an important eukaryotic model and widely used microbial cell factory. Standardized genetic parts encoding overexpression and knockdown mutations of >90% yeast genes are created in a single step from a full-length cDNA library. With the aid of CRISPR-Cas, these genetic parts are iteratively integrated into the repetitive genomic sequences in a modular manner using robotic automation. This system allows functional mapping and multiplex optimization on a genome scale for diverse phenotypes including cellulase expression, isobutanol production, glycerol utilization and acetic acid tolerance, and may greatly accelerate future genome-scale engineering endeavours in yeast.

Targeted and Untargeted Approaches Unravel Novel Candidate Genes and Diagnostic SNPs for Quantitative Resistance of the Potato (Solanum tuberosum L.) to Phytophthora infestans Causing the Late Blight Disease

The oomycete Phytophthora infestans causes late blight of potato, which can completely destroy the crop. Therefore, for the past 160 years, late blight has been the most important potato disease worldwide. The identification of cultivars with high and durable field resistance to P. infestans is an objective of most potato breeding programs. This type of resistance is polygenic and therefore quantitative. Its evaluation requires multi-year and location trials. Furthermore, quantitative resistance to late blight correlates with late plant maturity, a negative agricultural trait. Knowledge of the molecular genetic basis of quantitative resistance to late blight not compromised by late maturity is very limited. It is however essential for developing diagnostic DNA markers that facilitate the efficient combination of superior resistance alleles in improved cultivars. We used association genetics in a population of 184 tetraploid potato cultivars in order to identify single nucleotide polymorphisms (SNPs) that are associated with maturity corrected resistance (MCR) to late blight. The population was genotyped for almost 9000 SNPs from three different sources. The first source was candidate genes specifically selected for their function in the jasmonate pathway. The second source was novel candidate genes selected based on comparative transcript profiling (RNA-Seq) of groups of genotypes with contrasting levels of quantitative resistance to P. infestans. The third source was the first generation 8.3k SolCAP SNP genotyping array available in potato for genome wide association studies (GWAS). Twenty seven SNPs from all three sources showed robust association with MCR. Some of those were located in genes that are strong candidates for directly controlling quantitative resistance, based on functional annotation. Most important were: a lipoxygenase (jasmonate pathway), a 3-hydroxy-3-methylglutaryl coenzyme A reductase (mevalonate pathway), a P450 protein (terpene biosynthesis), a transcription factor and a homolog of a major gene for resistance to P. infestans from the wild potato species Solanum venturii. The candidate gene approach and GWAS complemented each other as they identified different genes. The results of this study provide new insight in the molecular genetic basis of quantitative resistance in potato and a toolbox of diagnostic SNP markers for breeding applications.

Sensory and immune genes identification and analysis in a widely used parasitoid wasp Trichogramma dendrolimi (Hymenoptera: Trichogrammatidae)

Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae) is one of the preponderant egg parasitoids of Dendrolimus spp., which are important defoliators of coniferous forests. This parasitoid wasp has been widely released to control pine caterpillar and other lepidopteran pests, but its control efficiency needs to be improved. Sensory systems are crucial for T. dendrolimi to locate hosts, and immunity is probably involved after egg deposition in the host cavity; however, few reports have focused on the molecular mechanism of olfactory detection and survival of T. dendrolimi. It is necessary to identify these genes before further functional research is conducted. In this study, we assembled and analyzed the transcriptome of T. dendrolimi using next-generation sequencing technology. The sequencing and assembly resulted in 38 565 contigs with N50 of 3422 bp. Sequence comparison indicate that T. dendrolimi sequences are very similar to those of another parasitoid Nasonia vitripennis. Then the olfactory, vision, and immune-related gene families were identified, and phylogenetic analyses were performed with these genes from T. dendrolimi and other model insect species. Furthermore, phylogenetic tree with odorant binding proteins of T. dendrolimi and their host Dendrolimus was constructed to determine whether convergent evolution exists. These genes can be valid targets for further gene function research. The present study may help us to understand host location and survival mechanisms of T. dendrolimi and to use them more efficiently for pest control in the future.

Large-Scale Transcriptome Analysis in Faba Bean (Vicia faba L.) under Ascochyta fabae Infection

Faba bean is an important food crop worldwide. However, progress in faba bean genomics lags far behind that of model systems due to limited availability of genetic and genomic information. Using the Illumina platform the faba bean transcriptome from leaves of two lines (29H and Vf136) subjected to Ascochyta fabae infection have been characterized. De novo transcriptome assembly provided a total of 39,185 different transcripts that were functionally annotated, and among these, 13,266 were assigned to gene ontology against Arabidopsis. Quality of the assembly was validated by RT-qPCR amplification of selected transcripts differentially expressed. Comparison of faba bean transcripts with those of better-characterized plant genomes such as Arabidopsis thaliana, Medicago truncatula and Cicer arietinum revealed a sequence similarity of 68.3%, 72.8% and 81.27%, respectively. Moreover, 39,060 single nucleotide polymorphism (SNP) and 3,669 InDels were identified for genotyping applications. Mapping of the sequence reads generated onto the assembled transcripts showed that 393 and 457 transcripts were overexpressed in the resistant (29H) and susceptible genotype (Vf136), respectively. Transcripts involved in plant-pathogen interactions such as leucine rich proteins (LRR) or plant growth regulators involved in plant adaptation to abiotic and biotic stresses were found to be differently expressed in the resistant line. The results reported here represent the most comprehensive transcript database developed so far in faba bean, providing valuable information that could be used to gain insight into the pathways involved in the resistance mechanism against A. fabae and to identify potential resistance genes to be further used in marker assisted selection.

Normalization of environmental metagenomic DNA enhances the discovery of under-represented microbial community members

Normalization is a procedure classically employed to detect rare sequences in cellular expression profiles (i.e. cDNA libraries). Here, we present a normalization protocol involving the direct treatment of extracted environmental metagenomic DNA with S1 nuclease, referred to as normalization of metagenomic DNA: NmDNA. We demonstrate that NmDNA, prior to post hoc PCR-based experiments (16S rRNA gene T-RFLP fingerprinting and clone library), increased the diversity of sequences retrieved from environmental microbial communities by detection of rarer sequences. This approach could be used to enhance the resolution of detection of ecologically relevant rare members in environmental microbial assemblages and therefore is promising in enabling a better understanding of ecosystem functioning.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study is the first testing 'normalization' on environmental metagenomic DNA (mDNA). The aim of this procedure was to improve the identification of rare phylotypes in environmental communities. Using hypoliths as model systems, we present evidence that this post-mDNA extraction molecular procedure substantially enhances the detection of less common phylotypes and could even lead to the discovery of novel microbial genotypes within a given environment.

Antennal transcriptome analysis and comparison of olfactory genes in two sympatric defoliators, Dendrolimus houi and Dendrolimus kikuchii (Lepidoptera: Lasiocampidae)

The Yunnan pine and Simao pine caterpillar moths, Dendrolimus houi Lajonquière and Dendrolimus kikuchii Matsumura (Lepidoptera: Lasiocampidae), are two closely related and sympatric pests of coniferous forests in southwestern China, and olfactory communication systems of these two insects have received considerable attention because of their economic importance. However, there is little information on the molecular aspect of odor detection about these insects. Furthermore, although lepidopteran species have been widely used in studies of insect olfaction, few work made comparison between sister moths on the olfactory recognition mechanisms. In this study, next-generation sequencing of the antennal transcriptome of these two moths were performed to identify the major olfactory genes. After comparing the antennal transcriptome of these two moths, we found that they exhibit highly similar transcripts-associated GO terms. Chemosensory gene families were further analyzed in both species. We identified 23 putative odorant binding proteins (OBP), 17 chemosensory proteins (CSP), two sensory neuron membrane proteins (SNMP), 33 odorant receptors (OR), and 10 ionotropic receptors (IR) in D. houi; and 27 putative OBPs, 17 CSPs, two SNMPs, 33 ORs, and nine IRs in D. kikuchii. All these transcripts were full-length or almost full-length. The predicted protein sequences were compared with orthologs in other species of Lepidoptera and model insects, including Bombyx mori, Manduca sexta, Heliothis virescens, Danaus plexippus, Sesamia inferens, Cydia pomonella, and Drosophila melanogaster. The sequence homologies of the orthologous genes in D. houi and D. kikuchii are very high. Furthermore, the olfactory genes were classed according to their expression level, and the highly expressed genes are our target for further function investigation. Interestingly, many highly expressed genes are ortholog gene of D. houi and D. kikuchii. We also found that the Classic OBPs were further separated into three groups according to their motifs, which will help future functional researches. Surprisingly, no pheromone receptor was identified in the two Dendrolimus species, which may indicate a special pheromone identification mechanism in Dendrolimus. Our work allows for further functional studies of pheromones and host volatile recognition genes, and give novel candidate targets for pest management.

Gene discovery through transcriptome sequencing for the invasive mussel Limnoperna fortunei

The success of the Asian bivalve Limnoperna fortunei as an invader in South America is related to its high acclimation capability. It can inhabit waters with a wide range of temperatures and salinity and handle long-term periods of air exposure. We describe the transcriptome of L. fortunei aiming to give a first insight into the phenotypic plasticity that allows non-native taxa to become established and widespread. We sequenced 95,219 reads from five main tissues of the mussel L. fortunei using Roche’s 454 and assembled them to form a set of 84,063 unigenes (contigs and singletons) representing partial or complete gene sequences. We annotated 24,816 unigenes using a BLAST sequence similarity search against a NCBI nr database. Unigenes were divided into 20 eggNOG functional categories and 292 KEGG metabolic pathways. From the total unigenes, 1,351 represented putative full-length genes of which 73.2% were functionally annotated. We described the first partial and complete gene sequences in order to start understanding bivalve invasiveness. An expansion of the hsp70 gene family, seen also in other bivalves, is present in L. fortunei and could be involved in its adaptation to extreme environments, e.g. during intertidal periods. The presence of toll-like receptors gives a first insight into an immune system that could be more complex than previously assumed and may be involved in the prevention of disease and extinction when population densities are high. Finally, the apparent lack of special adaptations to extremely low O₂ levels is a target worth pursuing for the development of a molecular control approach.

Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae)

Background

The European spruce bark beetle, Ips typographus, and the North American mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae), are severe pests of coniferous forests. Both bark beetle species utilize aggregation pheromones to coordinate mass-attacks on host trees, while odorants from host and non-host trees modulate the pheromone response. Thus, the bark beetle olfactory sense is of utmost importance for fitness. However, information on the genes underlying olfactory detection has been lacking in bark beetles and is limited in Coleoptera. We assembled antennal transcriptomes from next-generation sequencing of I. typographus and D. ponderosae to identify members of the major chemosensory multi-gene families.

Results

Gene ontology (GO) annotation indicated that the relative abundance of transcripts associated with specific GO terms was highly similar in the two species. Transcripts with terms related to olfactory function were found in both species. Focusing on the chemosensory gene families, we identified 15 putative odorant binding proteins (OBP), 6 chemosensory proteins (CSP), 3 sensory neuron membrane proteins (SNMP), 43 odorant receptors (OR), 6 gustatory receptors (GR), and 7 ionotropic receptors (IR) in I. typographus; and 31 putative OBPs, 11 CSPs, 3 SNMPs, 49 ORs, 2 GRs, and 15 IRs in D. ponderosae. Predicted protein sequences were compared with counterparts in the flour beetle, Tribolium castaneum, the cerambycid beetle, Megacyllene caryae, and the fruit fly, Drosophila melanogaster. The most notable result was found among the ORs, for which large bark beetle-specific expansions were found. However, some clades contained receptors from all four beetle species, indicating a degree of conservation among some coleopteran OR lineages. Putative GRs for carbon dioxide and orthologues for the conserved antennal IRs were included in the identified receptor sets.

Conclusions

The protein families important for chemoreception have now been identified in three coleopteran species (four species for the ORs). Thus, this study allows for improved evolutionary analyses of coleopteran olfaction. Identification of these proteins in two of the most destructive forest pests, sharing many semiochemicals, is especially important as they might represent novel targets for population control.

Consequences of normalizing transcriptomic and genomic libraries of plant genomes using a duplex-specific nuclease and tetramethylammonium chloride

Several applications of high throughput genome and transcriptome sequencing would benefit from a reduction of the high-copy-number sequences in the libraries being sequenced and analyzed, particularly when applied to species with large genomes. We adapted and analyzed the consequences of a method that utilizes a thermostable duplex-specific nuclease for reducing the high-copy components in transcriptomic and genomic libraries prior to sequencing. This reduces the time, cost, and computational effort of obtaining informative transcriptomic and genomic sequence data for both fully sequenced and non-sequenced genomes. It also reduces contamination from organellar DNA in preparations of nuclear DNA. Hybridization in the presence of 3 M tetramethylammonium chloride (TMAC), which equalizes the rates of hybridization of GC and AT nucleotide pairs, reduced the bias against sequences with high GC content. Consequences of this method on the reduction of high-copy and enrichment of low-copy sequences are reported for Arabidopsis and lettuce.

Construction and evaluation of normalized cDNA libraries enriched with full-length sequences for rapid discovery of new genes from Sisal (Agave sisalana Perr.) different developmental stages

To provide a resource of sisal-specific expressed sequence data and facilitate this powerful approach in new gene research, the preparation of normalized cDNA libraries enriched with full-length sequences is necessary. Four libraries were produced with RNA pooled from Agave sisalana multiple tissues to increase efficiency of normalization and maximize the number of independent genes by SMART™ method and the duplex-specific nuclease (DSN). This procedure kept the proportion of full-length cDNAs in the subtracted/normalized libraries and dramatically enhanced the discovery of new genes. Sequencing of 3875 cDNA clones of libraries revealed 3320 unigenes with an average insert length about 1.2 kb, indicating that the non-redundancy of libraries was about 85.7%. These unigene functions were predicted by comparing their sequences to functional domain databases and extensively annotated with Gene Ontology (GO) terms. Comparative analysis of sisal unigenes and other plant genomes revealed that four putative MADS-box genes and knotted-like homeobox (knox) gene were obtained from a total of 1162 full-length transcripts. Furthermore, real-time PCR showed that the characteristics of their transcripts mainly depended on the tight expression regulation of a number of genes during the leaf and flower development. Analysis of individual library sequence data indicated that the pooled-tissue approach was highly effective in discovering new genes and preparing libraries for efficient deep sequencing.

Development of a porcine (Sus scofa) embryo-specific microarray: array annotation and validation

BACKGROUND

The domestic pig is an important livestock species and there is strong interest in the factors that affect the development of viable embryos and offspring in this species. A limited understanding of the molecular mechanisms involved in early embryonic development has inhibited our ability to fully elucidate these factors. Next generation deep sequencing and microarray technologies are powerful tools for delineation of molecular pathways involved in the developing embryo.

RESULTS

Here we present the development of a porcine-embryo-specific microarray platform created from a large expressed sequence tag (EST) analysis generated by Roche/454 next-generation sequencing of cDNAs constructed from critical stages of in vivo or in vitro porcine preimplantation embryos. Two cDNA libraries constructed from in vitro and in vivo produced preimplantation porcine embryos were normalized and sequenced using 454 Titanium pyrosequencing technology. Over one million high-quality EST sequences were obtained and used to develop the EMbryogene Porcine Version 1 (EMPV1) microarray composed of 43,795 probes. Based on an initial probe sequence annotation, the EMPV1 features 17,409 protein-coding, 473 pseudogenes, 46 retrotransposed, 2,359 non-coding RNA, 4,121 splice variants in 2,862 genes and a total of 12,324 Novel Transcript Regions (NTR). After re-annotation, the total unique genes increased from 11,961 to 16,281 and 1.9% of them belonged to a large olfactory receptor (OR) gene family. Quality control on the EMPV1 was performed and revealed an even distribution of ten clusters of spiked-in control spots and array to array (dye-swap) correlation was 0.97.

CONCLUSIONS

Using next-generation deep sequencing we have produced a large EST dataset to allow for the selection of probe sequences for the development of the EMPV1 microarray platform. The quality of this embryo-specific array was confirmed with a high-level of reproducibility using current Agilent microarray technology. With more than an estimated 20,000 unique genes represented on the EMPV1, this platform will provide the foundation for future research into the in vivo and in vitro factors that affect the viability of porcine embryos, as well as the effects of these factors on the live offspring that result from these embryos.

Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration

Background

The naid annelids contain a number of species that vary in their ability to regenerate lost body parts, making them excellent candidates for evolution of regeneration studies. However, scant sequence data exists to facilitate such studies. We constructed a cDNA library from the naid Pristina leidyi, a species that is highly regenerative and also reproduces asexually by fission, using material from a range of regeneration and fission stages for our library. We then sequenced the transcriptome of P. leidyi using 454 technology.

Results

454 sequencing produced 1,550,174 reads with an average read length of 376 nucleotides. Assembly of 454 sequence reads resulted in 64,522 isogroups and 46,679 singletons for a total of 111,201 unigenes in this transcriptome. We estimate that over 95% of the transcripts in our library are present in our transcriptome. 17.7% of isogroups had significant BLAST hits to the UniProt database and these include putative homologs of a number of genes relevant to regeneration research. Although many sequences are incomplete, the mean sequence length of transcripts (isotigs) is 707 nucleotides. Thus, many sequences are large enough to be immediately useful for downstream applications such as gene expression analyses. Using in situ hybridization, we show that two Wnt/β-catenin pathway genes (homologs of frizzled and β-catenin) present in our transcriptome are expressed in the regeneration blastema of P. leidyi, demonstrating the usefulness of this resource for regeneration research.

Conclusions

454 sequencing is a rapid and efficient approach for identifying large numbers of genes in an organism that lacks a sequenced genome. This transcriptome dataset will be a valuable resource for molecular analyses of regeneration in P. leidyi and will serve as a starting point for comparisons to non-regenerating naids. It also contributes significantly to the still limited genomic resources available for annelids and lophotrochozoans more generally.

[Preparation of prokaryotic cDNA for high-throughput transcriptome analysis]

High contents of non-coding RNA in total bacteria RNA complicates considerably transcriptome analysis using standard approaches like high-throughput sequencing, gene expression profiles, subtractive hybridization. We suggest a procedure of preparation of bacterial cDNA for transcriptomics that includes rRNA and tRNA depletion with preservation of relative abundance of coding sequences. The method is based on the second order hybridization kinetics and unique properties of Kanchatka crab duplex-specific nuclease. The method efficacy was demonstrated on a model experiments.

Increased sensitivity of next generation sequencing-based expression profiling after globin reduction in human blood RNA

Background

Transcriptome analysis is of great interest in clinical research, where significant differences between individuals can be translated into biomarkers of disease. Although next generation sequencing provides robust, comparable and highly informative expression profiling data, with several million of tags per blood sample, reticulocyte globin transcripts can constitute up to 76% of total mRNA compromising the detection of low abundant transcripts. We have removed globin transcripts from 6 human whole blood RNA samples with a human globin reduction kit and compared them with the same non-reduced samples using deep Serial Analysis of Gene Expression.

Results

Globin tags comprised 52-76% of total tags in our samples. Out of 21,633 genes only 87 genes were detected at significantly lower levels in the globin reduced samples. In contrast, 11,338 genes were detected at significantly higher levels in the globin reduced samples. Removing globin transcripts allowed us to also identify 2112 genes that could not be detected in the non-globin reduced samples, with roles in cell surface receptor signal transduction, G-protein coupled receptor protein signalling pathways and neurological processes.

Conclusions

The reduction of globin transcripts in whole blood samples constitutes a reproducible and reliable method that can enrich data obtained from next generation sequencing-based expression profiling.

Single cell genome sequencing

Whole genome amplification and next-generation sequencing of single cells have become a powerful approach for studying uncultivated microorganisms that represent 90-99% of all environmental microbes. Single cell sequencing enables not only the identification of microbes but also linking of functions to species, a feat not achievable by metagenomic techniques. Moreover, it allows the analysis of low abundance species that may be missed in community-based analyses. It has also proved very useful in complementing metagenomics in the assembly and binning of single genomes. With the advent of drastically cheaper and higher throughput sequencing technologies, it is expected that single cell sequencing will become a standard tool in studying the genome and transcriptome of microbial communities.

[Application of the duplex-specific nuclease for fast analysis of single nucleotide polymorphisms and detection of target DNA in complex PCR products]

We have developed a simple method for fast analysis of single nucleotide polymorphisms and identification of target clones from cloned complex PCR products. The method utilizes Kamchatka crab duplex-specific nuclease and universal fluorescent probe and is alternative to laborious screening procedures using radioactive probes, restriction analysis followed by gel electrophoresis or expensive sequencing. The method efficacy was demonstrated in several model experiments.

Duplex-specific nuclease efficiently removes rRNA for prokaryotic RNA-seq

Next-generation sequencing has great potential for application in bacterial transcriptomics. However, unlike eukaryotes, bacteria have no clear mechanism to select mRNAs over rRNAs; therefore, rRNA removal is a critical step in sequencing-based transcriptomics. Duplex-specific nuclease (DSN) is an enzyme that, at high temperatures, degrades duplex DNA in preference to single-stranded DNA. DSN treatment has been successfully used to normalize the relative transcript abundance in mRNA-enriched cDNA libraries from eukaryotic organisms. In this study, we demonstrate the utility of this method to remove rRNA from prokaryotic total RNA. We evaluated the efficacy of DSN to remove rRNA by comparing it with the conventional subtractive hybridization (Hyb) method. Illumina deep sequencing was performed to obtain transcriptomes from Escherichia coli grown under four growth conditions. The results clearly showed that our DSN treatment was more efficient at removing rRNA than the Hyb method was, while preserving the original relative abundance of mRNA species in bacterial cells. Therefore, we propose that, for bacterial mRNA-seq experiments, DSN treatment should be preferred to Hyb-based methods.

Synthetic spike-in standards for RNA-seq experiments

High-throughput sequencing of cDNA (RNA-seq) is a widely deployed transcriptome profiling and annotation technique, but questions about the performance of different protocols and platforms remain. We used a newly developed pool of 96 synthetic RNAs with various lengths, and GC content covering a 2(20) concentration range as spike-in controls to measure sensitivity, accuracy, and biases in RNA-seq experiments as well as to derive standard curves for quantifying the abundance of transcripts. We observed linearity between read density and RNA input over the entire detection range and excellent agreement between replicates, but we observed significantly larger imprecision than expected under pure Poisson sampling errors. We use the control RNAs to directly measure reproducible protocol-dependent biases due to GC content and transcript length as well as stereotypic heterogeneity in coverage across transcripts correlated with position relative to RNA termini and priming sequence bias. These effects lead to biased quantification for short transcripts and individual exons, which is a serious problem for measurements of isoform abundances, but that can partially be corrected using appropriate models of bias. By using the control RNAs, we derive limits for the discovery and detection of rare transcripts in RNA-seq experiments. By using data collected as part of the model organism and human Encyclopedia of DNA Elements projects (ENCODE and modENCODE), we demonstrate that external RNA controls are a useful resource for evaluating sensitivity and accuracy of RNA-seq experiments for transcriptome discovery and quantification. These quality metrics facilitate comparable analysis across different samples, protocols, and platforms.

The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus

BACKGROUND

Most evolutionary developmental biology ("evo-devo") studies of emerging model organisms focus on small numbers of candidate genes cloned individually using degenerate PCR. However, newly available sequencing technologies such as 454 pyrosequencing have recently begun to allow for massive gene discovery in animals without sequenced genomes. Within insects, although large volumes of sequence data are available for holometabolous insects, developmental studies of basally branching hemimetabolous insects typically suffer from low rates of gene discovery.

RESULTS

We used 454 pyrosequencing to sequence over 500 million bases of cDNA from the ovaries and embryos of the milkweed bug Oncopeltus fasciatus, which lacks a sequenced genome. This indirectly developing insect occupies an important phylogenetic position, branching basal to Diptera (including fruit flies) and Hymenoptera (including honeybees), and is an experimentally tractable model for short-germ development. 2,087,410 reads from both normalized and non-normalized cDNA assembled into 21,097 sequences (isotigs) and 112,531 singletons. The assembled sequences fell into 16,617 unique gene models, and included predictions of splicing isoforms, which we examined experimentally. Discovery of new genes plateaued after assembly of ~1.5 million reads, suggesting that we have sequenced nearly all transcripts present in the cDNA sampled. Many transcripts have been assembled at close to full length, and there is a net gain of sequence data for over half of the pre-existing O. fasciatus accessions for developmental genes in GenBank. We identified 10,775 unique genes, including members of all major conserved metazoan signaling pathways and genes involved in several major categories of early developmental processes. We also specifically address the effects of cDNA normalization on gene discovery in de novo transcriptome analyses.

CONCLUSIONS

Our sequencing, assembly and annotation framework provide a simple and effective way to achieve high-throughput gene discovery for organisms lacking a sequenced genome. These data will have applications to the study of the evolution of arthropod genes and genetic pathways, and to the wider evolution, development and genomics communities working with emerging model organisms.[The sequence data from this study have been submitted to GenBank under study accession number SRP002610 (http://www.ncbi.nlm.nih.gov/sra?term=SRP002610). Custom scripts generated are available at http://www.extavourlab.com/protocols/index.html. Seven Additional files are available.].

Normalization of full-length-enriched cDNA

A well-recognized obstacle to efficient high-throughput analysis of cDNA libraries is the differential abundance of various transcripts in any particular cell type. Decreasing the prevalence of clones representing abundant transcripts before sequencing, using cDNA normalization, may significantly increase the efficacy of random sequencing and is essential for rare gene discovery. Duplex-specific nuclease (DSN) normalization allows the generation of normalized full-length-enriched cDNA libraries to permit a high gene discovery rate. The method is based on the unique properties of DSN from the Kamchatka crab and involves denaturation-reassociation of cDNA, degradation of the ds-fraction formed by abundant transcripts by DSN, and PCR amplification of the remaining ss-DNA fraction. The method has been evaluated in various plant and animal models.

Normalization of genomic DNA using duplex-specific nuclease

An application of duplex-specific nuclease (DSN) normalization technology to whole-genome shotgun sequencing of genomes with a large proportion of repetitive DNA is described. The method uses a thermostable DSN from the Kamchatka crab that specifically hydrolyzes dsDNA. In model experiments on human genomic DNA, we demonstrated that DSN normalization of double-stranded DNA formed during C0t analysis is effective against abundant repetitive sequences with high sequence identity, while retaining highly divergent repeats and coding regions at base-line levels. Thus, DSN normalization applied to C0t analysis can be used to eliminate evolutionarily young repetitive elements from genomic DNA before sequencing, and should prove invaluable in studies of large eukaryotic genomes, such as those of higher plants.

A novel expression profile of the Loxosceles intermedia spider venomous gland revealed by transcriptome analysis

Spiders of the Loxosceles genus are cosmopolitan, and their venom components possess remarkable biological properties associated with their ability to act upon different molecules and receptors. Accidents with Loxosceles intermedia specimens are recognized as a public health problem in the south of Brazil. To describe the transcriptional profile of the L. intermedia venom gland, we generated a wide cDNA library, and its transcripts were functionally and structurally analyzed. After initial analyses, 1843 expressed sequence tags (ESTs) produced readable sequences that were grouped into 538 clusters, 281 of which were singletons. 985 reads (53% of total ESTs) matched to known proteins. Similarity searches showed that toxin-encoding transcripts account for 43% of the total library and comprise a great number of ESTs. The most frequent toxins were from the LiTx family, which are known for their insecticidal activity. Both phospholipase D and astacin-like metalloproteases toxins account for approximately 9% of total transcripts. Toxins components such as serine proteases, hyaluronidases and venom allergens were also found but with minor representation. Almost 10% of the ESTs encode for proteins involved in cellular processes. These data provide an important overview of the L. intermedia venom gland expression scenario and revealed significant differences from profiles of other spiders from the Loxosceles genus. Furthermore, our results also confirm that this venom constitutes an amazing source of novel compounds with potential agrochemical, industrial and pharmacological applications.

Heart transcriptome of the bank vole (Myodes glareolus): towards understanding the evolutionary variation in metabolic rate

BACKGROUND

Understanding the genetic basis of adaptive changes has been a major goal of evolutionary biology. In complex organisms without sequenced genomes, de novo transcriptome assembly using a longer read sequencing technology followed by expression profiling using short reads is likely to provide comprehensive identification of adaptive variation at the expression level and sequence polymorphisms in coding regions. We performed sequencing and de novo assembly of the bank vole heart transcriptome in lines selected for high metabolism and unselected controls.

RESULTS

A single 454 Titanium run produced over million reads, which were assembled into 63,581 contigs. Searches against the SwissProt protein database and the ENSEMBL collection of mouse transcripts detected similarity to 11,181 and 14,051 genes, respectively. As judged by the representation of genes from the heart-related Gene Ontology categories and UniGenes detected in the mouse heart, our detection of the genes expressed in the heart was nearly complete (> 95% and almost 90% respectively). On average, 38.7% of the transcript length was covered by our sequences, with notably higher (45.0%) coverage of coding regions than of untranslated regions (24.5% of 5' and 32.7% of 3'UTRs). Lower sequence conservation between mouse and bank vole in untranslated regions was found to be partially responsible for poorer UTR representation. Our data might suggest a widespread transcription from noncoding genomic regions, a finding not reported in previous studies regarding transcriptomes in non-model organisms. We also identified over 19 thousand putative single nucleotide polymorphisms (SNPs). A much higher fraction of the SNPs than expected by chance exhibited variant frequency differences between selection regimes.

CONCLUSION

Longer reads and higher sequence yield per run provided by the 454 Titanium technology in comparison to earlier generations of pyrosequencing proved beneficial for the quality of assembly. An almost full representation of genes known to be expressed in the mouse heart was identified. Usage of the extensive genomic resources available for the house mouse, a moderately (20-40 mln years) divergent relative of the voles, enabled a comprehensive assessment of the transcript completeness. Transcript sequences generated in the present study allowed the identification of candidate SNPs associated with divergence of selection lines and constitute a valuable permanent resource forming a foundation for RNAseq experiments aiming at detection of adaptive changes both at the level of gene expression and sequence variants, that would facilitate studies of the genetic basis of evolutionary divergence.

Shedding light on an extremophile lifestyle through transcriptomics

The tropical intertidal ecosystem is defined by trees - mangroves - which are adapted to an extreme and extremely variable environment. The genetic basis underlying these adaptations is, however, virtually unknown. Based on advances in pyrosequencing, we present here the first transcriptome analysis for plants for which no prior genomic information was available. We selected the mangroves Rhizophora mangle (Rhizophoraceae) and Heritiera littoralis (Malvaceae) as ecologically important extremophiles employing markedly different physiological and life-history strategies for survival and dominance in this extreme environment. For maximal representation of conditional transcripts, mRNA was obtained from a variety of developmental stages, tissues types, and habitats. For each species, a normalized cDNA library of pooled mRNAs was analysed using GSFLX pyrosequencing. A total of 537,635 sequences were assembled de novo and annotated as > 13,000 distinct gene models for each species. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology annotations highlighted remarkable similarities in the mangrove transcriptome profiles, which differed substantially from the model plants Arabidopsis and Populus. Similarities in the two species suggest a unique mangrove lifestyle overarching the effects of transcriptome size, habitat, tissue type, developmental stage, and biogeographic and phylogenetic differences between them.

EST and EST-SSR marker resources for Iris

BACKGROUND

Limited DNA sequence and DNA marker resources have been developed for Iris (Iridaceae), a monocot genus of 200-300 species in the Asparagales, several of which are horticulturally important. We mined an I. brevicaulis-I. fulva EST database for simple sequence repeats (SSRs) and developed ortholog-specific EST-SSR markers for genetic mapping and other genotyping applications in Iris. Here, we describe the abundance and other characteristics of SSRs identified in the transcript assembly (EST database) and the cross-species utility and polymorphisms of I. brevicaulis-I. fulva EST-SSR markers among wild collected ecotypes and horticulturally important cultivars.

RESULTS

Collectively, 6,530 ESTs were produced from normalized leaf and root cDNA libraries of I. brevicaulis (IB72) and I. fulva (IF174), and assembled into 4,917 unigenes (1,066 contigs and 3,851 singletons). We identified 1,447 SSRs in 1,162 unigenes and developed 526 EST-SSR markers, each tracing a different unigene. Three-fourths of the EST-SSR markers (399/526) amplified alleles from IB72 and IF174 and 84% (335/399) were polymorphic between IB25 and IF174, the parents of I. brevicaulis x I. fulva mapping populations. Forty EST-SSR markers were screened for polymorphisms among 39 ecotypes or cultivars of seven species - 100% amplified alleles from wild collected ecotypes of Louisiana Iris (I.brevicaulis, I.fulva, I. nelsonii, and I. hexagona), whereas 42-52% amplified alleles from cultivars of three horticulturally important species (I. pseudacorus, I. germanica, and I. sibirica). Ecotypes and cultivars were genetically diverse - the number of alleles/locus ranged from two to 18 and mean heterozygosity was 0.76.

CONCLUSION

Nearly 400 ortholog-specific EST-SSR markers were developed for comparative genetic mapping and other genotyping applications in Iris, were highly polymorphic among ecotypes and cultivars, and have broad utility for genotyping applications within the genus.