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.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide.

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

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.

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.

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.

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.

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.

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.

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).

Asymmetrically expressed axin required for anterior development in Tribolium.

Canonical Wnt signaling has been implicated in an AP axis polarizing mechanism in most animals, despite limited evidence from arthropods. In the long-germ insect, Drosophila, Wnt signaling is not required for global AP patterning, but in short-germ insects including Tribolium castaneum, loss of Wnt signaling affects development of segments in the growth zone but not those defined in the blastoderm. To determine the effects of ectopic Wnt signaling, we analyzed the expression and function of axin, which encodes a highly conserved negative regulator of the pathway. We found Tc-axin transcripts maternally localized to the anterior pole in freshly laid eggs. Expression spread toward the posterior pole during the early cleavage stages, becoming ubiquitous by the time the germ rudiment formed. Tc-axin RNAi produced progeny phenotypes that ranged from mildly affected embryos with cuticles displaying a graded loss of anterior structures, to defective embryos that condensed at the posterior pole in the absence of serosa. Altered expression domains of several blastodermal markers indicated anterior expansion of posterior fates. Analysis of other canonical Wnt pathway components and the expansion of Tc-caudal expression, a Wnt target, suggest that the effects of Tc-axin depletion are mediated through this pathway and that Wnt signaling must be inhibited for proper anterior development in Tribolium. These studies provide unique evidence that canonical Wnt signaling must be carefully regulated along the AP axis in an arthropod, and support an ancestral role for Wnt activity in defining AP polarity and patterning in metazoan development.

Changes in anterior head patterning underlie the evolution of long germ embryogenesis.

Early embryonic stages differ significantly among related animal taxa while subsequent development converges at the conserved phylotypic stage before again diverging. Although this phenomenon has long been observed, its underlying genetic mechanisms remain enigmatic. The dipteran Drosophila melanogaster develops as a long germ embryo where the head anlagen form a cap at the anterior pole of the blastoderm. Consequently, the anterior and terminal maternal systems give crucial input for head patterning. However, in the short germ beetle Tribolium castaneum, as in most insects, the head anlagen is located at a ventral position distant from the anterior pole of the blastoderm. In line with these divergent embryonic anlagen, several differences in the axis formation between the insects have been discovered. We now ask to what extent patterning and morphogenesis of the anterior median region (AMR) of the head, including clypeolabral and stomodeal anlagen, differ among these insects. Unexpectedly, we find that Tc-huckebein is not a terminal gap gene and, unlike its Drosophila ortholog, is not involved in Tribolium head development. Instead, Tc-six3 acts upstream of Tc-crocodile and Tc-cap'n'collar to pattern posterior and anterior parts of the AMR, respectively. We further find that instead of huckebein, Tc-crocodile is required for stomodeum development by activating Tc-forkhead. Finally, a morphogenetic movement not found in Drosophila shapes the embryonic head of Tribolium. Apparently, with anterior displacement of the head anlagen during long germ evolution of Drosophila, the ancestral regulation by the bilaterian anterior control gene six3 was replaced by the anterior and terminal maternal systems, which were further elaborated by adding bicoid, tailless and huckebein as anterior regionalization genes.

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.

Candidate gene screen in the red flour beetle Tribolium reveals six3 as ancient regulator of anterior median head and central complex development.

Several highly conserved genes play a role in anterior neural plate patterning of vertebrates and in head and brain patterning of insects. However, head involution in Drosophila has impeded a systematic identification of genes required for insect head formation. Therefore, we use the red flour beetle Tribolium castaneum in order to comprehensively test the function of orthologs of vertebrate neural plate patterning genes for a function in insect head development. RNAi analysis reveals that most of these genes are indeed required for insect head capsule patterning, and we also identified several genes that had not been implicated in this process before. Furthermore, we show that Tc-six3/optix acts upstream of Tc-wingless, Tc-orthodenticle1, and Tc-eyeless to control anterior median development. Finally, we demonstrate that Tc-six3/optix is the first gene known to be required for the embryonic formation of the central complex, a midline-spanning brain part connected to the neuroendocrine pars intercerebralis. These functions are very likely conserved among bilaterians since vertebrate six3 is required for neuroendocrine and median brain development with certain mutations leading to holoprosencephaly.

Genetic and developmental analysis of differences in eye and face morphology between Drosophila simulans and Drosophila mauritiana.

Eye and head morphology vary considerably among insects and even between closely related species of Drosophila. Species of the D. melanogaster subgroup, and other Drosophila species, exhibit a negative correlation between eye size and face width (FW); for example, D. mauritiana generally has bigger eyes composed of larger ommatidia and conversely a narrower face than its sibling species. To better understand the evolution of eye and head morphology, we investigated the genetic and developmental basis of differences in eye size and FW between male D. mauritiana and D. simulans. QTL mapping of eye size and FW showed that the major loci responsible for the interspecific variation in these traits are localized to different genomic regions. Introgression of the largest effect QTL underlying the difference in eye size resulted in flies with larger eyes but no significant difference in FW. Moreover,introgression of a QTL region on the third chromosome that contributes to the FW difference between these species affected FW, but not eye size. We also observed that this difference in FW is detectable earlier in the development of the eye­antennal disc than the difference in the size of the retinal field. Our results suggest that different loci that act at different developmental stages underlie changes in eye size and FW. Therefore, while there is a negative correlation between these traits in Drosophila, we show genetically that they also have the potential to evolve independently and this may help to explain the evolution of these traits in other insects.

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.

Insect Tc-six4 marks a unit with similarity to vertebrate placodes.

Cranial placodes are specialized ectodermal regions in the developing vertebrate head that give rise to both neural and non-neural cell types of the neuroendocrine system and the sense organs of the visual, olfactory and acoustic systems. The cranial placodes develop from a panplacodal region which is specifically marked by genes of the eyes absent/eya and two "six homeobox" family members (sine oculis/six1 and six4). It had been believed that cranial placodes are evolutionary novelties of vertebrates. However, data from non-vertebrate chordates suggest that placode-like structures evolved in the chordate ancestor already. Here, we identify a morphological structure in the embryonic head of the beetle Tribolium castaneum with placode-like features. It is marked by the orthologs of the panplacodal markers Tc-six4, Tc-eya and Tc-sine oculis/six1 (Tc-six1) and expresses several genes known to be involved in adenohypophyseal placode development in vertebrates. Moreover, it contributes to both epidermal and neural tissues. We identify Tc-six4 as a specific marker for this structure that we term the insect head placode. Finally, we reveal the regulatory gene network of the panplacodal genes Tc-six4, Tc-eya and Tc-six1 and identify them as head epidermis patterning genes. Our finding of a placode-like structure in an insect suggests that a placode precursor was already present in the last common ancestor of bilaterian animals.

Formation of the insect head involves lateral contribution of the intercalary segment, which depends on Tc-labial function.

The insect head is composed of several segments. During embryonic development, the segments fuse to form a rigid head capsule where obvious segmental boundaries are lacking. Hence, the assignment of regions of the insect head to specific segments is hampered, especially with respect to dorsal (vertex) and lateral (gena) parts. We show that upon Tribolium labial (Tc-lab) knock down, the intercalary segment is deleted but not transformed. Furthermore, we find that the intercalary segment contributes to lateral parts of the head cuticle in Tribolium. Based on several additional mutant and RNAi phenotypes that interfere with gnathal segment development, we show that these segments do not contribute to the dorsal head capsule apart from the dorsal ridge. Opposing the classical view but in line with findings in the vinegar fly Drosophila melanogaster and the milkweed bug Oncopeltus fasciatus, we propose a "bend and zipper" model for insect head capsule formation.