From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology.

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Gene expression mapping of the neuroectoderm across phyla - conservation and divergence of early brain anlagen between insects and vertebrates.

Gene expression has been employed for homologizing body regions across bilateria. The molecular comparison of vertebrate and fly brains has led to a number of disputed homology hypotheses. Data from the fly Drosophila melanogaster have recently been complemented by extensive data from the red flour beetle Tribolium castaneum with its more insect-typical development. In this review, we revisit the molecular mapping of the neuroectoderm of insects and vertebrates to reconsider homology hypotheses. We claim that the protocerebrum is non-segmental and homologous to the vertebrate fore- and midbrain. The boundary between antennal and ocular regions correspond to the vertebrate mid-hindbrain boundary while the deutocerebrum represents the anterior-most ganglion with serial homology to the trunk. The insect head placode is shares common embryonic origin with the vertebrate adenohypophyseal placode. Intriguingly, vertebrate eyes develop from a different region compared to the insect compound eyes calling organ homology into question. Finally, we suggest a molecular re-definition of the classic concepts of archi- and prosocerebrum.

Tissue dissociation for single-cell and single-nuclei RNA sequencing for low amounts of input material.

BACKGROUND: Recent technological advances opened the opportunity to simultaneously study gene expression for thousands of individual cells on a genome-wide scale. The experimental accessibility of such single-cell RNA sequencing (scRNAseq) approaches allowed gaining insights into the cell type composition of heterogeneous tissue samples of animal model systems and emerging models alike. A major prerequisite for a successful application of the method is the dissociation of complex tissues into individual cells, which often requires large amounts of input material and harsh mechanical, chemical and temperature conditions. However, the availability of tissue material may be limited for small animals, specific organs, certain developmental stages or if samples need to be acquired from collected specimens. Therefore, we evaluated different dissociation protocols to obtain single cells from small tissue samples of Drosophila melanogaster eye-antennal imaginal discs. RESULTS: We show that a combination of mechanical and chemical dissociation resulted in sufficient high-quality cells. As an alternative, we tested protocols for the isolation of single nuclei, which turned out to be highly efficient for fresh and frozen tissue samples. Eventually, we performed scRNAseq and single-nuclei RNA sequencing (snRNAseq) to show that the best protocols for both methods successfully identified relevant cell types. At the same time, snRNAseq resulted in less artificial gene expression that is caused by rather harsh dissociation conditions needed to obtain single cells for scRNAseq. A direct comparison of scRNAseq and snRNAseq data revealed that both datasets share biologically relevant genes among the most variable genes, and we showed differences in the relative contribution of the two approaches to identified cell types. CONCLUSION: We present two dissociation protocols that allow isolating single cells and single nuclei, respectively, from low input material. Both protocols resulted in extraction of high-quality RNA for subsequent scRNAseq or snRNAseq applications. If tissue availability is limited, we recommend the snRNAseq procedure of fresh or frozen tissue samples as it is perfectly suited to obtain thorough insights into cellular diversity of complex tissue.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila.

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Multiple loci linked to inversions are associated with eye size variation in species of the Drosophila virilis phylad.

The size and shape of organs is tightly controlled to achieve optimal function. Natural morphological variations often represent functional adaptations to an ever-changing environment. For instance, variation in head morphology is pervasive in insects and the underlying molecular basis is starting to be revealed in the Drosophila genus for species of the melanogaster group. However, it remains unclear whether similar diversifications are governed by similar or different molecular mechanisms over longer timescales. To address this issue, we used species of the virilis phylad because they have been diverging from D. melanogaster for at least 40 million years. Our comprehensive morphological survey revealed remarkable differences in eye size and head shape among these species with D. novamexicana having the smallest eyes and southern D. americana populations having the largest eyes. We show that the genetic architecture underlying eye size variation is complex with multiple associated genetic variants located on most chromosomes. Our genome wide association study (GWAS) strongly suggests that some of the putative causative variants are associated with the presence of inversions. Indeed, northern populations of D. americana share derived inversions with D. novamexicana and they show smaller eyes compared to southern ones. Intriguingly, we observed a significant enrichment of genes involved in eye development on the 4th chromosome after intersecting chromosomal regions associated with phenotypic differences with those showing high differentiation among D. americana populations. We propose that variants associated with chromosomal inversions contribute to both intra- and interspecific variation in eye size among species of the virilis phylad.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Sex differences in spiders: from phenotype to genomics.

Sexual reproduction is pervasive in animals and has led to the evolution of sexual dimorphism. In most animals, males and females show marked differences in primary and secondary sexual traits. The formation of sex-specific organs and eventually sex-specific behaviors is defined during the development of an organism. Sex determination processes have been extensively studied in a few well-established model organisms. While some key molecular regulators are conserved across animals, the initiation of sex determination is highly diverse. To reveal the mechanisms underlying the development of sexual dimorphism and to identify the evolutionary forces driving the evolution of different sexes, sex determination mechanisms must thus be studied in detail in many different animal species beyond the typical model systems. In this perspective article, we argue that spiders represent an excellent group of animals in which to study sex determination mechanisms. We show that spiders are sexually dimorphic in various morphological, behavioral, and life history traits. The availability of an increasing number of genomic and transcriptomic resources and functional tools provides a great starting point to scrutinize the extensive sexual dimorphism present in spiders on a mechanistic level. We provide an overview of the current knowledge of sex determination in spiders and propose approaches to reveal the molecular and genetic underpinnings of sexual dimorphism in these exciting animals.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies.

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.

Origin and Consequences of Chromosomal Inversions in the virilis Group of Drosophila.

In Drosophila, large variations in rearrangement rate have been reported among different lineages and among Muller's elements. Nevertheless, the mechanisms that are involved in the generation of inversions, their increase in frequency, as well as their impact on the genome are not completely understood. This is in part due to the lack of comparative studies on species distantly related to Drosophila melanogaster. Therefore, we sequenced and assembled the genomes of two species of the virilis phylad (Drosophila novamexicana [15010-1031.00] and Drosophila americana [SF12]), which are diverging from D. melanogaster for more than 40 Myr. Based on these data, we identified the precise location of six novel inversion breakpoints. A molecular characterization provided clear evidence that DAIBAM (a miniature inverted-repeat transposable element) was involved in the generation of eight out of the nine inversions identified. In contrast to what has been previously reported for D. melanogaster and close relatives, ectopic recombination is thus the prevalent mechanism of generating inversions in species of the virilis phylad. Using pool-sequencing data for three populations of D. americana, we also show that common polymorphic inversions create a high degree of genetic differentiation between populations for chromosomes X, 4, and 5 over large physical distances. We did not find statistically significant differences in expression levels between D. americana (SF12) and D. novamexicana (15010-1031.00) strains for the three genes surveyed (CG9588, Fig 4, and fab1) flanking three inversion breakpoints.

Dynamic genome wide expression profiling of Drosophila head development reveals a novel role of Hunchback in retinal glia cell development and blood-brain barrier integrity.

Drosophila melanogaster head development represents a valuable process to study the developmental control of various organs, such as the antennae, the dorsal ocelli and the compound eyes from a common precursor, the eye-antennal imaginal disc. While the gene regulatory network underlying compound eye development has been extensively studied, the key transcription factors regulating the formation of other head structures from the same imaginal disc are largely unknown. We obtained the developmental transcriptome of the eye-antennal discs covering late patterning processes at the late 2nd larval instar stage to the onset and progression of differentiation at the end of larval development. We revealed the expression profiles of all genes expressed during eye-antennal disc development and we determined temporally co-expressed genes by hierarchical clustering. Since co-expressed genes may be regulated by common transcriptional regulators, we combined our transcriptome dataset with publicly available ChIP-seq data to identify central transcription factors that co-regulate genes during head development. Besides the identification of already known and well-described transcription factors, we show that the transcription factor Hunchback (Hb) regulates a significant number of genes that are expressed during late differentiation stages. We confirm that hb is expressed in two polyploid subperineurial glia cells (carpet cells) and a thorough functional analysis shows that loss of Hb function results in a loss of carpet cells in the eye-antennal disc. Additionally, we provide for the first time functional data indicating that carpet cells are an integral part of the blood-brain barrier. Eventually, we combined our expression data with a de novo Hb motif search to reveal stage specific putative target genes of which we find a significant number indeed expressed in carpet cells.

A novel role for Ets4 in axis specification and cell migration in the spider Parasteatoda tepidariorum.

Organizers play important roles during the embryonic development of many animals. The most famous example is the Spemann organizer that sets up embryonic axes in amphibian embryos. In spiders, a group of BMP secreting mesenchymal cells (the cumulus) functions as an organizer of the dorsoventral axis. Similar to experiments performed with the Spemann organizer, transplantation of the cumulus is able to induce a secondary axis in spiders. Despite the importance of this structure, it is unknown which factors are needed to activate cumulus specific gene expression. To address this question, we performed a transcriptomic analysis of early embryonic development in the spider Parasteatoda tepidariorum. Through this work, we found that the transcription factor Pt-Ets4 is needed for cumulus integrity, dorsoventral patterning and for the activation of Pt-hunchback and Pt-twist expression. Furthermore, ectopic expression of Pt-Ets4 is sufficient to induce cell delamination and migration by inducing a mesoderm-like cell fate.

Specific expression and function of the Six3 optix in Drosophila serially homologous organs.

Organ size and pattern results from the integration of two positional information systems. One global information system, encoded by the Hox genes, links organ type with position along the main body axis. Within specific organs, local information is conveyed by signaling molecules that regulate organ growth and pattern. The mesothoracic (T2) wing and the metathoracic (T3) haltere of Drosophila represent a paradigmatic example of this coordination. The Hox gene Ultrabithorax (Ubx), expressed in the developing T3, selects haltere identity by, among other processes, modulating the production and signaling efficiency of Dpp, a BMP2-like molecule that acts as a major regulator of size and pattern. However, the mechanisms of the Hox-signal integration in this well-studied system are incomplete. Here, we have investigated this issue by studying the expression and function of the Six3 transcription factor optix during Drosophila wing and haltere development. We find that in both organs, Dpp defines the expression domain of optix through repression, and that the specific position of this domain in wing and haltere seems to reflect the differential signaling profile among these organs. We show that optix expression in wing and haltere primordia is conserved beyond Drosophila in other higher diptera. In Drosophila, optix is necessary for the growth of wing and haltere. In the wing, optix is required for the growth of the most anterior/proximal region (the 'marginal cell') and for the correct formation of sensory structures along the proximal anterior wing margin; the halteres of optix mutants are also significantly reduced. In addition, in the haltere, optix is necessary for the suppression of sensory bristles.

A robust (re-)annotation approach to generate unbiased mapping references for RNA-seq-based analyses of differential expression across closely related species.

BACKGROUND: RNA-seq based on short reads generated by next generation sequencing technologies has become the main approach to study differential gene expression. Until now, the main applications of this technique have been to study the variation of gene expression in a whole organism, tissue or cell type under different conditions or at different developmental stages. However, RNA-seq also has a great potential to be used in evolutionary studies to investigate gene expression divergence in closely related species. RESULTS: We show that the published genomes and annotations of the three closely related Drosophila species D. melanogaster, D. simulans and D. mauritiana have limitations for inter-specific gene expression studies. This is due to missing gene models in at least one of the genome annotations, unclear orthology assignments and significant gene length differences in the different species. A comprehensive evaluation of four statistical frameworks (DESeq2, DESeq2 with length correction, RPKM-limma and RPKM-voom-limma) shows that none of these methods sufficiently accounts for inter-specific gene length differences, which inevitably results in false positive candidate genes. We propose that published reference genomes should be re-annotated before using them as references for RNA-seq experiments to include as many genes as possible and to account for a potential length bias. We present a straight-forward reciprocal re-annotation pipeline that allows to reliably compare the expression for nearly all genes annotated in D. melanogaster. CONCLUSIONS: We conclude that our reciprocal re-annotation of previously published genomes facilitates the analysis of significantly more genes in an inter-specific differential gene expression study. We propose that the established pipeline can easily be applied to re-annotate other genomes of closely related animals and plants to improve comparative expression analyses.

Size relationships of different body parts in the three dipteran species Drosophila melanogaster, Ceratitis capitata and Musca domestica.

Body size is an integral feature of an organism that influences many aspects of life such as fecundity, life span and mating success. Size of individual organs and the entire body size represent quantitative traits with a large reaction norm, which are influenced by various environmental factors. In the model system Drosophila melanogaster, pupal size and adult traits, such as tibia and thorax length or wing size, accurately estimate the overall body size. However, it is unclear whether these traits can be used in other flies. Therefore, we studied changes in size of pupae and adult organs in response to different rearing temperatures and densities for D. melanogaster, Ceratitis capitata and Musca domestica. We confirm a clear sexual size dimorphism (SSD) for Drosophila and show that the SSD is less uniform in the other species. Moreover, the size response to changing growth conditions is sex dependent. Comparison of static and evolutionary allometries of the studied traits revealed that response to the same environmental variable is genotype specific but has similarities between species of the same order. We conclude that the value of adult traits as estimators of the absolute body size may differ among species and the use of a single trait may result in wrong assumptions. Therefore, we suggest using a body size coefficient computed from several individual measurements. Our data is of special importance for monitoring activities of natural populations of the three dipteran flies, since they are harmful species causing economical damage (Drosophila, Ceratitis) or transferring diseases (Musca).

Molecular characterization and embryonic origin of the eyes in the common house spider Parasteatoda tepidariorum.

BACKGROUND: Two visual systems are present in most arthropod groups: median and lateral eyes. Most of our current knowledge about the developmental and molecular mechanisms involved in eye formation in arthropods comes from research in the model system Drosophila melanogaster. Here, a core set of retinal determination genes, namely, sine-oculis (so), eyes absent (eya), dachshund (dac), and the two pax6 orthologues eyeless (ey) and twin of eyeless (toy) govern early retinal development. By contrast, not much is known about the development of the up-to-eight eyes present in spiders. Therefore, we analyzed the embryonic expression of core retinal determination genes in the common house spider Parasteatoda tepidariorum. RESULTS: We show that the anlagen of the median and lateral eyes in P. tepidariorum originate from different regions of the non-neurogenic ectoderm in the embryonic head. The median eyes are specified as two individual anlagen in an anterior median position in the developing head and subsequently move to their final position following extensive morphogenetic movements of the non-neurogenic ectoderm. The lateral eyes develop from a more lateral position. Intriguingly, they are specified as a unique field of cells that splits into the three individual lateral eyes during late embryonic development. Using gene expression analyses, we identified a unique combination of determination gene expression in the anlagen of the lateral and median eyes, respectively. CONCLUSIONS: This study of retinal determination genes in the common house spider P. tepidariorum represents the first comprehensive analysis of the well-known retinal determination genes in arthropods outside insects. The development of the individual lateral eyes via the subdivision of one single eye primordium might be the vestige of a larger composite eye anlage, and thus supports the notion that the composite eye is the plesiomorphic state of the lateral eyes in arthropods. The molecular distinction of the two visual systems is similar to the one described for compound eyes and ocelli in Drosophila, suggesting that a unique core determination network for median and lateral eyes, respectively, might have been in place already in the last common ancestor of spiders and insects.

Sexual dimorphism and natural variation within and among species in the Drosophila retinal mosaic.

BACKGROUND: Insect compound eyes are composed of ommatidia, which contain photoreceptor cells that are sensitive to different wavelengths of light defined by the specific rhodopsin proteins that they express. The fruit fly Drosophila melanogaster has several different ommatidium types that can be localised to specific retinal regions, such as the dorsal rim area (DRA), or distributed stochastically in a mosaic across the retina, like the 'pale' and 'yellow' types. Variation in these ommatidia patterns very likely has important implications for the vision of insects and could underlie behavioural and environmental adaptations. However, despite the detailed understanding of ommatidia specification in D. melanogaster, the extent to which the frequency and distribution of the different ommatidium types vary between sexes, strains and species of Drosophila is not known. RESULTS: We investigated the frequency and distribution of ommatidium types based on rhodopsin protein expression, and the expression levels of rhodopsin transcripts in the eyes of both sexes of different strains of D. melanogaster, D. simulans and D. mauritiana. We found that while the number of DRA ommatidia was invariant, Rh3 expressing ommatidia were more frequent in the larger eyes of females compared to the males of all species analysed. The frequency and distribution of ommatidium types also differed between strains and species. The D. simulans strain ZOM4 has the highest frequency of Rh3 expressing ommatidia, which is associated with a non-stochastic patch of pale and odd-coupled ommatidia in the dorsal-posterior of their eyes. CONCLUSIONS: Our results show that there is striking variation in the frequency and distribution of ommatidium types between sexes, strains and species of Drosophila. This suggests that evolutionary changes in the underlying regulatory mechanisms can alter the distribution of ommatidium types to promote or restrict their expression in specific regions of the eye within and between species, and that this could cause differences in vision among these flies.

A comprehensive reference transcriptome resource for the common house spider Parasteatoda tepidariorum.

Parasteatoda tepidariorum is an increasingly popular model for the study of spider development and the evolution of development more broadly. However, fully understanding the regulation and evolution of P. tepidariorum development in comparison to other animals requires a genomic perspective. Although research on P. tepidariorum has provided major new insights, gene analysis to date has been limited to candidate gene approaches. Furthermore, the few available EST collections are based on embryonic transcripts, which have not been systematically annotated and are unlikely to contain transcripts specific to post-embryonic stages of development. We therefore generated cDNA from pooled embryos representing all described embryonic stages, as well as post-embryonic stages including nymphs, larvae and adults, and using Illumina HiSeq technology obtained a total of 625,076,514 100-bp paired end reads. We combined these data with 24,360 ESTs available in GenBank, and 1,040,006 reads newly generated from 454 pyrosequencing of a mixed-stage embryo cDNA library. The combined sequence data were assembled using a custom de novo assembly strategy designed to optimize assembly product length, number of predicted transcripts, and proportion of raw reads incorporated into the assembly. The de novo assembly generated 446,427 contigs with an N50 of 1,875 bp. These sequences obtained 62,799 unique BLAST hits against the NCBI non-redundant protein data base, including putative orthologs to 8,917 Drosophila melanogaster genes based on best reciprocal BLAST hit identity compared with the D. melanogaster proteome. Finally, we explored the utility of the transcriptome for RNA-Seq studies, and showed that this resource can be used as a mapping scaffold to detect differential gene expression in different cDNA libraries. This resource will therefore provide a platform for future genomic, gene expression and functional approaches using P. tepidariorum.