De novo assembly of RNA-Seq data allows us to study transcriptomes without the need for a genome sequence, such as in non-model organisms of ecological and evolutionary importance, cancer samples, or the microbiome. In this protocol, we describe the use of the Trinity platform for de novo transcriptome assembly from RNA-Seq data in non-model organisms. We also present Trinity’s supported companion utilities for downstream applications, including RSEM for transcript abundance estimation, R/Bioconductor packages for identifying differentially expressed transcripts across samples, and approaches to identify protein coding genes. In an included tutorial we provide a workflow for genome-independent transcriptome analysis leveraging the Trinity platform. The software, documentation and demonstrations are freely available from http://trinityrnaseq.sf.net.

Gene silencing through RNA interference (RNAi) has revolutionized the study of gene 98 function, particularly in non-model insects. However, in Lepidoptera (moths and butterflies) 99 RNAi has many times proven to be difficult to achieve. Most of the negative results have been 100 anecdotal and the positive experiments have not been collected in such a way that they are 101 possible to analyze. In this review, we have collected detailed data from more than 150 102 experiments including all to date published and many unpublished experiments. Despite a 103 large variation in the data, trends that are found are that RNAi is particularly successful in the 104 family Saturniidae and in genes involved in immunity. On the contrary, gene expression in 105 epidermal tissues seems to be most difficult to silence. In addition, gene silencing by feeding 106 dsRNA requires high concentrations for success. Possible causes for the variability of success 107 in RNAi experiments in Lepidoptera are discussed. The review also points to a need to further 108 investigate the mechanism of RNAi in lepidopteran insects and its possible connection to the 109 innate immune response. Our general understanding of RNAi in Lepidoptera will be further 110 aided in the future as our public database at http://insectacentral.org/RNAi will continue to 111 gather information on RNAi experiments.

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

BackgroundThe Mediterranean fruit fly (medfly), Ceratitis capitata, is a major destructive insect pest due to its broad host range, which includes hundreds of fruits and vegetables. It exhibits a unique ability to invade and adapt to ecological niches throughout tropical and subtropical regions of the world, though medfly infestations have been prevented and controlled by the sterile insect technique (SIT) as part of integrated pest management programs (IPMs). The genetic analysis and manipulation of medfly has been subject to intensive study in an effort to improve SIT efficacy and other aspects of IPM control.ResultsThe 479 Mb medfly genome is sequenced from adult flies from lines inbred for 20 generations. A high-quality assembly is achieved having a contig N50 of 45.7 kb and scaffold N50 of 4.06 Mb. In-depth curation of more than 1800 messenger RNAs shows specific gene expansions that can be related to invasiveness and host adaptation, including gene families for chemoreception, toxin and insecticide metabolism, cuticle proteins, opsins, and aquaporins. We identify genes relevant to IPM control, including those required to improve SIT.ConclusionsThe medfly genome sequence provides critical insights into the biology of one of the most serious and widespread agricultural pests. This knowledge should significantly advance the means of controlling the size and invasive potential of medfly populations. Its close relationship to Drosophila, and other insect species important to agriculture and human health, will further comparative functional and structural studies of insect genomes that should broaden our understanding of gene family evolution.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1049-2) contains supplementary material, which is available to authorized users.

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.

Insects and their arthropod relatives including mites, spiders, and crustaceans play major roles in the world's terrestrial, aquatic, and marine ecosystems. Arthropods compete with humans for food and transmit devastating diseases. They also comprise the most diverse and successful branch of metazoan evolution, with millions of extant species. Here, we describe an international effort to guide arthropod genomic efforts, from species prioritization to methodology and informatics. The 5000 arthropod genomes initiative (i5K) community met formally in 2012 to discuss a roadmap for sequencing and analyzing 5000 high-priority arthropods and is continuing this effort via pilot projects, the development of standard operating procedures, and training of students and career scientists. With university, governmental, and industry support, the i5K Consortium aspires to deliver sequences and analytical tools for each of the arthropod branches and each of the species having beneficial and negative effects on humankind.

The mimetic wing patterns of Heliconius butterflies are an excellent example of both adaptive radiation and convergent evolution. Alleles at the HmYb and HmSb loci control the presence ⁄ absence of hindwing bar and hindwing margin phenotypes respectively between divergent races of Heliconius melpomene, and also between sister species. Here, we used fine-scale linkage mapping to identify and sequence a BAC tilepath across the HmYb ⁄ Sb loci. We also generated transcriptome sequence data for two wing pattern forms of H. melpomene that differed in HmYb ⁄ Sb alleles using 454 sequencing technology. Custom scripts were used to process the sequence traces and generate transcriptome assemblies. Genomic sequence for the HmYb ⁄ Sb candidate region was annotated both using the MAKER pipeline and manually using transcriptome sequence reads. In total, 28 genes were identified in the HmYb ⁄ Sb candidate region, six of which have alternative splice forms. None of these are orthologues of genes previously identified as being expressed in butterfly wing pattern development, implying previously undescribed molecular mechanisms of pattern determination on Heliconius wings. The use of next-generation sequencing has therefore facilitated DNA annotation of a poorly characterized genome, and generated hypotheses regarding the identity of wing pattern at the HmYb ⁄ Sb loci.