The iBeetle large-scale RNAi screen reveals gene functions for insect development and physiology

Abstract

Genetic screens are powerful tools to identify the genes required for a given biological process. However, for technical reasons, comprehensive screens have been restricted to very few model organisms. Therefore, although deep sequencing is revealing the genes of ever more insect species, the functional studies predominantly focus on candidate genes previously identified in Drosophila, which is biasing research towards conserved gene functions. RNAi screens in other organisms promise to reduce this bias. Here we present the results of the iBeetle screen, a large-scale, unbiased RNAi screen in the red flour beetle, Tribolium castaneum, which identifies gene functions in embryonic and postembryonic development, physiology and cell biology. The utility of Tribolium as a screening platform is demonstrated by the identification of genes involved in insect epithelial adhesion. This work transcends the restrictions of the candidate gene approach and opens fields of research not accessible in Drosophila.

ARTICLE

Received 4 Dec 2014 | Accepted 15 Jun 2015 | Published 28 Jul 2015

Christian Schmitt-Engel1,2, Dorothea Schultheis2,*, Jonas Schwirz1,*, Nadi Strhlein2,*, Nicole Troelenberg2,*, Upalparna Majumdar2,*, Van Anh Dao3,*, Daniela Grossmann1, Tobias Richter2, Maike Tech4, Jrgen Dnitz1,4,5, Lizzy Gerischer6, Mirko Theis7,8, Inga Schild1, Jochen Trauner2, Nikolaus D.B. Koniszewski1, Elke Kster1, Sebastian Kittelmann1, Yonggang Hu1, Sabrina Lehmann1, Janna Siemanowski1, Julia Ulrich1, Kristen A. Panlio3, Reinhard Schrder9, Burkhard Morgenstern4, Mario Stanke6, Frank Buchhholz7, Manfred Frasch2,Siegfried Roth3, Ernst A. Wimmer1, Michael Schoppmeier2,**, Martin Klingler2,** & Gregor Bucher1,**

Genetic screens are powerful tools to identify the genes required for a given biological process. However, for technical reasons, comprehensive screens have been restricted to very few model organisms. Therefore, although deep sequencing is revealing the genes of ever more insect species, the functional studies predominantly focus on candidate genes previously identied in Drosophila, which is biasing research towards conserved gene functions. RNAi screens in other organisms promise to reduce this bias. Here we present the results of the iBeetle screen, a large-scale, unbiased RNAi screen in the red our beetle, Tribolium castaneum, which identies gene functions in embryonic and postembryonic development, physiology and cell biology. The utility of Tribolium as a screening platform is demonstrated by the identication of genes involved in insect epithelial adhesion. This work transcends the restrictions of the candidate gene approach and opens elds of research not accessible in Drosophila.

1 Johann-Friedrich-Blumenbach-Institut, GZMB, Georg-August-Universitat Gttingen, Justus-von-Liebig-Weg 11, 37077 Gttingen, Germany. 2 Department Biologie, Entwicklungsbiologie, Friedrich-Alexander-Universitat Erlangen-Nrnberg, Staudtstra 5, 91058 Erlangen, Germany. 3 Institut fr

Entwicklungsbiologie, Universitat zu Kln, Zlpicher Stra 47b, 50674 Cologne, Germany. 4 Abteilung Bioinformatik, Institut fr Mikrobiologie und Genetik, Georg-August-Universitat Gttingen, Goldschmidtstr. 1, 37077 Gttingen, Germany. 5 Abteilung fr Bioinformatik, Universitatsmedizin Gttingen, Goldschmidtstr. 1, 37077 Gttingen, Germany. 6 Institut fr Mathematik und Informatik, Ernst Moritz Arndt Universitat Greifswald, Walther-Rathenau-Stra 47, 17487 Greifswald, Germany. 7 TU Dresden, Medical Faculty Carl Gustav Carus, Medical Systems Biology, Fetscherstr. 74, 01307 Dresden, Germany.

8 Eupheria Biotech GmbH, Dresden, Germany. 9 Abteilung fr Genetik, Universitat Rostock, Albert-Einsteinstr. 3, 18059 Rostock, Germany. * These authors contributed equally to this work. ** These authors jointly supervised this work. Correspondence and requests for materials should be addressed toM.K. (email: mailto:[email protected]

Web End [email protected] ) or to G.B. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2015 Macmillan Publishers Limited. All rights reserved.

DOI: 10.1038/ncomms8822 OPEN

The iBeetle large-scale RNAi screen reveals gene functions for insect development and physiology

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822

The importance of Drosophila melanogaster as a model system is in large part due to its amenability to elegant genetic screens, which allow the comprehensive identica

tion of genes required for a given biological process1,2. Recently, unbiased genetic screens have also been performed in a few other insects, for example, the hymenopteran Nasonia vitripennis and the beetle Tribolium castaneum36. However, technical constraints prohibit saturation screens in these species. Hence, most of what we know about insect gene function remains based on Drosophila work. RNA interference (RNAi) has emerged as an alternative tool to knockdown gene function and has thus far been used for genome-wide screens in the nematode Caenorhabditis elegans, in Drosophila and in cell culture711.

In recent years, reverse genetics based on deep sequencing and RNAi has enabled functional investigations in ever more insect species, broadening the range of biological phenomena that can be analysed. However, the candidate gene approach still prevails: genes are selected based on previous ndings in Drosophila or other model systems1221. Consequently, the eld of insect functional genetics suffers from several limitations. First, the candidate gene approach leads to a bias towards the study of conserved gene functions. Second, it has remained difcult to identify genes required for processes that are not represented in Drosophila. Finally, technical limitations and lineage specic gene losses or duplications prohibit the identication of comprehensive gene sets for a particular process in any single species. An unbiased, large-scale RNAi screen in a non-dipteran insect species should overcome many of these limitations.

The red our beetle, T. castaneum, is well suited for this aim. In many respects its biology is more representative for insects than Drosophila, as its segmentation proceeds from a posterior growth zone, the larval head is not involuted and its extraembryonic tissues are well developed. Further, its modes of oogenesis and metamorphosis resemble those of other non-dipteran insects2228. Further, Tribolium is a representative of the most species-rich animal taxon on earth, the coleopterans (beetles), including many devastating pests29. Finally, Tribolium research builds on an expanding transgenic toolkit6,3032 and a particularly strong systemic RNAi response. Knockdown phenotypes can be induced at all life stages through dsRNA injection into the body cavity, and the RNAi effect spreads throughout the animal and is transferred to the offspring, often phenocopying null mutants3336.

Here we present the results of the iBeetle screen, where we used RNAi targeting about 5,000 genes to identify novel gene functions during oogenesis, embryogenesis and metamorphosis in an unbiased way. We show that this screen has the power of overcoming the current limitations imposed by the candidate gene approach by the identication of unexpected novel players required for long studied processes. For instance, we describe the rst bicaudal phenotype in Tribolium elicited by knockdown of the homeobox gene Tc-homeobrain; a gene so far not related to anteriorposterior axis formation. Further, we show that the unbiased detection of gene function in Tribolium allows opening new elds of research. For instance, many insects have odoriferous glands used for communication and defence but such glands are missing in Drosophila. In the screen, a set of genes were identied, which are involved in producing the defensive chemicals of the Tribolium odoriferous glands. Importantly, many of these genes were not identied in a recent RNA-seq approach conrming the power of a phenotypic screen. Finally, we show that Tribolium is an excellent alternative screening platform, where insect gene functions are efciently identied. One example is the gene Tc-Rbm24, which we found to be required for muscle development in Tribolium. Drosophila does not have a respective ortholog while the vertebrate ortholog has

recently been shown to be involved in muscle formation. As second example, we identied novel genes required for epithelial adhesion in Tribolium, the orthologs of which are required in Drosophila as well but had not been discovered in Drosophila screens.

ResultsDesign of the iBeetle screen. We developed a procedure that allowed efcient screening of several biological processes. Two screens were performed in parallel by injection of different developmental stages. In the pupal injection screen, injected pig19 transgenic female pupae (somatic muscles marked with EGFP) were scored for late metamorphosis phenotypes and, upon maturation to adults, their offspring embryos were analysed for muscle and cuticle phenotypes as readouts for defects in embryogenesis. This treatment knocked down both maternal and zygotic transcripts in developing embryos. In case of reduced egg production, ovaries of the injected females were analysed for oogenesis defects.

In the larval injection screen, penultimate instar larvae (L6) were injected. Female larvae were derived from a cross between D17Xred (adult ight muscles marked with EGFP; X-linked DsRed marker allowed sexing of larvae) and pearl (white eyed) strains. Muscle phenotypes were scored during the pupal stage, and general morphological defects both at pupal and at adult stages. Ovaries were dissected and analysed whenever egg production was found to be strongly reduced. Finally, adult odoriferous glands were scored for alterations in size or colouring and dissected for closer inspection. Importantly, the larval injection screen allowed the identication of gene functions during metamorphosis without affecting essential functions during embryogenesis that would prevent analysis of later stages. We screened 5,300 genes in the (ongoing) pupal and 4,480 genes in the larval injection screen and present analyses of 3,400 genes included in both screens.

The iBeetle screen was designed as a rst pass screen wherein each experiment was performed once, and off-target controls were done only for selected genes. We aimed at minimizing false negative annotations with the trade-off of an increased false positive rate. The genes to be knocked down were selected randomly, except that their annotations were based on RNA sequence data, which may have led to some enrichment of highly expressed genes (see Supplementary Table 1). Using the DEQOR prediction algorithm37, templates were selected for high RNAi efciency and a low number of possible off-target sites.dsRNA fragments with an average length of 479 bp were injected at a concentration of 1 mg ml 1. Phenotypes were annotated in an online database according to the EQM system (entity, quality, modier)38 and using a controlled vocabulary based on the Tribolium morphological ontology39. In addition, the penetrance of phenotypes was recorded, and pictures and free text elds were used for further documentation (see Methods section for details). All abovementioned data sets of the larval and pupal injection screen in addition to sequence and orthology information of the entire Tribolium gene set are available at http://ibeetle-base.uni-goettingen.de

Web End =http://ibeetle-base.uni-goettingen.de 40

Tests for Sensitivity and Reproducibility. To assess the sensitivity of our screen, roughly 5% of screened genes were positive controls from a set of 41 different genes (Supplementary Fig. 1). In addition, 48 previously published genes that happened to be within our gene set served as additional blind positive controls. Ninety-three per cent of the selected controls and 95% of the previously published genes were identied during the screen (Fig. 1a,b; Supplementary Table 2), which is similar to the

2 NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822 ARTICLE

a

b

Recognition of positive controls Recognition of genes with known phenotypes

7%

18%

71%

4%

4%

4% 4%3%

100%

80%

60%

40%

20%

0%

1%

2%

17%

3%

100%

80%

60%

40%

20%

0%

Missed Missed

Technical lethality

Technical lethality

10%

16%

84%

90% 83%

Partially recognized

Recognized

78%

Reproducibly different phenotype

Recognized

Larval (196)

Pupal (201)

Sum (397)

Larval(28)

Pupal(37)

Sum(65)

c d

Negative control injections

Reproducibility

1%

2%

2%

2%

100%

80%

60%

40%

20%

0%

4%

15%

77%

4%

1%

1%

12%

False positive

No phenotype

14%

98% 96% 97%

Technical lethality

6% 68%

Reproduced

Strain-specific

Fragment-specific

Not reproduced

(n =158)

Larval (124)

Pupal (127)

Sum (251)

Figure 1 | Sensitivity and reproducibility. (a) Recognition rates of 41 different positive controls shown separately for the larval and pupal injection screens (left and middle bars) and for both together (right bar). About 80% of the positive controls were fully recognized while another 10% were partially recognized (that is, not all phenotypic aspects were annotated). Only 4% of the positive controls were missed. technical lethality: Expected phenotype not recognized owing to lethality of the animals for example, by injection. (b) Recognition rates for dsRNAs targeting 48 genes with published phenotypes, which had by chance been included in the screen. Of all, 78% were recognized with the published phenotype while 17% were annotated with a reproducibly different phenotype; that is, the differing phenotype was reproduced in independent experiments under iBeetle conditions. Hence, these different phenotypes are biologically meaningful and reect that the timing and the degree of gene knockdown inuences the phenotype. See Supplementary Note 1 for discussion of these cases. (c) Only 2% of all buffer injections led to false positive annotations. (d) A total of 158 dsRNAs were tested in independent injections with non-overlapping fragments. When the phenotype differed from the screening result, we analysed whether it was a false positive (not reproduced), or whether the genetic background was the reason for the difference (strain specic). Finally, we tested whether the outcome depended on the dsRNA fragment used (fragment specic), which indicated off-target effects or splice variant specic knockdown.

ndings in a genome-wide RNAi screen in C. elegans11. Interestingly, we found reproducibly different or additional phenotypes for 17% of the published genes, which likely reects the dependence of the RNAi effect on injection time and dsRNA concentration (see Supplementary Note 1). As expected, negative control injections usually produced no phenotype (Fig. 1c; Supplementary Note 1; Supplementary Fig. 2).

To test for reproducibility, we repeated the screening procedure for 158 genes with high penetrance phenotypes (Fig. 1d; see Supplementary Data 1 for details). Seventy-four per cent of the phenotypes were reproduced (Fig. 1d). Notably, 6% of those phenotypes turned out to depend on the genetic background. This is in line with the emerging view that strain specicity of phenotypes may be more prevalent than appreciated previously41,42. Twelve per cent of the re-tested phenotypes turned out to be false positives, while 14% could be reproduced by the original dsRNA fragment but not with a non-overlapping fragment. The latter nding probably reects off-target effects but in some cases may reect biologically meaningful differences due to the isoform specic nature of RNAi knockdown in Tribolium43. Overall, the proportion of off-target effects using systemic RNAi with long dsRNA fragments in Tribolium was similar to that observed with short transgenic hairpin constructs in Drosophila, where 13% (n 9) of the constructs induced

unexpected lethality44, but it was signicantly lower than the 28%

(n 18) or 24% (n 65) found with long hairpin constructs8,45.

Considering 12% of not reproduced phenotypes and a maximum of 14% off-target effects, the maximum false positive rate in the iBeetle screen is 26%. Hence, the phenotypes detected in our rst pass screen need to be conrmed by a second assay using non-overlapping fragments.

Reproducibility highly depended on the biological process scored. For lethality or wing blister phenotypes the reproducibility was 495%, while for L1 cuticle phenotypes it was about 60% (see Supplementary Fig. 3 for numbers for different phenotype classes).

To test in how far we missed the strong phenotypes due to incomplete knockdown, we injected dsRNAs targeting 98 genes at concentrations of 1 and 3 mg ml 1, and compared the phenotypic quality and strength. We found that phenotype strength was comparable for 86% of these genes.

Essential genes identied in the larval and pupal screens. Of the 3,400 genes tested in both the larval and pupal screens, 56.3% gave any phenotype (Fig. 2a), with 49.6% being lethal for at least one developmental stage (Fig. 2b). In all, 22.9% of the genes displayed a phenotype in both screens, while almost twice as many genes showed a phenotype exclusively in the pupal screen(21.1%) compared with phenotypes restricted to the larval screen

NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822

a

b

Screening summary (3,400)

No phenotype

Any phenotype pupal screen

Any phenotype in both screens

Any phenotype in larval screen

a b

c' d '56.0%

43.7%

56.3%

Injections leading tolethality (1,686) 50.4%

49.6%

717 (21.1%)

510 (15.0%)

175 (5.1 %)

794 (23.5%)

780 (22.9%)

418 (12.3%)

1485 (43.7%)

199 (5.9%)

(0.2%)

%)

%)

8(0.2%)

Lethal after larval injection

Early lethal after larval injection

Parental lethal after pupal injection

Embryonic lethal after pupal injection

(884, 26.0%)

(317, 9.3%)

(207, 6.1%)

(969, 28.5%)

35.2%

c

d

Phenotypes after pupal injection (1,497)

44.0%

Phenotypes after larvalinjection (1,198) 64.8%

1 (0.02%)

24 (0.7%)

198 (5.8%)

206 (6.1%) 866 (25.5%)

20 (0.6%)

94 (2.8%)

90 (2.7%)

11 (0.3%)

539 (15.9%)

197 (5.6%)

%

%)

%)

80 (2.4%)

15 (0.4%)

161 (4.7 %)

5 (0.2%)

78 (2.3%)

Ovary defects (227, 6.7%)

Embryonic lethal (969, 28.5%)

Parental lethal (207, 6.1%)

Embryonic muscle defects

Lethal after larval injection (884, 26.0%)

Defects during the process of metamorphosis

(267, 7.9%)

Defects of adult structures (406, 11.9%)

Ovary defects (100, 2.9%)

No epidermal patterning defects

Morphological defects of L1 larval cuticle

Embryonic lethal before cuticle secretion

(186, 5.5%)

(281, 8.3%)

(502, 14.8%)

Figure 2 | Essential and lethal genes. (a) For more than 56% of the injected genes, phenotypes were observed. The pupal injection screen revealed phenotypes for a larger portion of genes compared with the larval injection screen. (b) Death of the injected animals was scored 22 days post injection (larval injection; blue circle) and 11 days post injection (pupal and larval injection; dark green and hatched blue circles). Note that embryonic lethality is based on maternal and zygotic gene knockdown. Parental lethal: death of the injected animal. (c) Selected phenotypic categories after pupal injection. Embryonic lethal injections are further categorized showing that more than half of the embryonic lethal genes lead to abortion of embryogenesis before cuticle secretion. (d) Phenotypic categories after larval injection. Defects during the process of metamorphosis: metamorphosis not completed or entered precociously. Insets: relations to the entire data set.

(12.3%; Fig. 2a). Thirteen per cent of the genes showing any phenotype (including lethality) were beetle specic genes (Supplementary Table 1), illustrating the importance of screens in additional model organisms (see Supplementary Table 3 for denition of phenotype classes and Supplementary Data 2 for lists of treatments in the respective classes).

Basically, all genes required for survival to adulthood after pupal injection (5.9% of all genes; Fig. 2b) were also required after larval injection, likely due to housekeeping functions. The set of larval lethal genes appeared to be much larger (26% of all genes;

Fig. 2b). However, this difference is in large parts due to lowpenetrance lethal phenotypes, which were more likely to reach the threshold in the larval screen, where lethality was checked 22 days post injection instead of 11 days post injection in the pupal screen.

Embryonic lethality was found for 28.5% of all genes following knockdown of maternal and zygotic gene function (Fig. 2b,c). Altogether, 5.5% did not show obvious cuticle defects, while 8.3% displayed cuticle aberrations (Fig. 2c); 14.8% of all genes did not develop a cuticle, leading to a so called empty egg phenotype.

4 NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822 ARTICLE

This phenotypic class comprises genes with diverse essential functions, including housekeeping, cuticle formation itself or fertilization; however, some early patterning genes are also known to result in death before cuticle formation in Tribolium4648. Further analyses using molecular markers are required to determine which process was affected for a given empty egg phenotype. Interestingly, 5.6% of all genes in the larval injection screen led to alterations of adult morphology without affecting larval survival or fertility, making them interesting candidates for understanding how adult morphologies develop and how they evolved (Fig. 2d).

Comparison of embryonic and postembryonic patterning. For the rst time our data allow the systematic comparison of the gene sets required for embryonic and postembryonic development in an insect with typical metamorphosis. Importantly, larval cells are largely re-used to form the adult epidermis in most insects, instead of being replaced by imaginal cells as is the case in Drosophila24. Nevertheless, the gene sets involved in embryonic

and postembryonic patterning turned out to be largely non-overlapping (Fig. 3a). This is true for processes as different as leg or muscle development (Fig. 3b,c). In the case of oogenesis, the respective numbers are probably an overestimation because most genes leading to reduced egg production in the pupal screen were lethal in the larval screen (Fig. 3d). We assume that many of these apparent oogenesis phenotypes reect incomplete knockdown of genes with basic physiological function because we found that many animals with reduced oogenesis showed a strongly reduced fat body. Subtracting these genes (those outside the dashed line in Fig. 3d), the overlap of genes with defects in both screening parts increases substantially. Together, these data reveal that development during typical insect metamorphosis partially relies on different mechanisms than during embryogenesis.

Essential genes of Tribolium and Drosophila. The classical genetic screens for embryonic phenotypes in Drosophila revealed that about 5,000 genes were lethal when mutated (36% of the Drosophila protein coding genes), and for an additional 1,000

a

b

Morphological defects of external structures (549)

Leg defects (104)

16.2%83.8%

3.1%96.9%

209 (6.1%) 49 (1.4%) 49 (1.4%)

6 (0.2%)

292 (8.6%)

48 (1.4%)

Morphological defects of pupae or adults (larval screen, 340)

Morphological defects of pupae or adults legs (larval screen, 55)

Morphological defects of larval cuticle (pupal screen, 257)

Morphological defects of larval legs (pupal screen, 55)

Muscle defects (150) Ovary defects (304)

4.4%

95.6%

c

d

8.9%91.1%

35 (1.0%)

77 (2.3%)

23 (0.7%)

55 (1.6%)

6 (0.2%)

89 (2.6%)

204 (6.0%)

Defects of adult thoracic musculature (larval screen, 61)

Defects of larval musculature (pupal screen, 95)

Ovary defects (larval screen, 100)

Ovary defects (pupal screen, 227)

Ovary defects in pupal screen and not lethal in larval screen (58)

Figure 3 | Comparison of gene sets involved in embryonic versus postembryonic development. (a) The gene sets required for cuticle morphology (that is, epidermal patterning) during embryogenesis and typical insect metamorphosis are largely non-overlapping. This indicates that patterning principles may differ to quite some extent between these two stages of major morphological change. (b,c) This observation also holds true for the subsets affecting leg morphology (b) and GFP marked somatic muscles (c) indicating that both ectodermal and other patterning processes differ. (d) Gene sets required for ovary function. Many genes required for egg production in the pupal injection screen (green circle) were lethal in the larval injection screen. Hence, reduced egg production for these genes was probably due to starvation (green area outside hatched line). When comparing the non-lethal treatments (blue circle and green circle with hatched blue outline) the number of genes with an ovary phenotype in the pupal and larval injection screen are more similar. Insets: relations to the entire data set.

NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822

genes non-lethal phenotypes were identied. Hence, mutations in 43% of all loci revealed a phenotype of some kind4951. The respective numbers for Tribolium are very similar (37 and 42%; all numbers in this paragraph are corrected assuming 26% false positives; see Methods section for calculations). However, in the iBeetle screen a much larger portion of the lethal genes is embryonic lethal compared with the Drosophila genetic screens (58% versus 20%). This is mainly due to a much higher number of embryonic lethal genes that show cuticle defects in Tribolium (81% versus 15% in Drosophila)49. One likely reason for this increased lethality is that in our screen we knocked down both maternal and zygotic gene functions, while in most Drosophila genetic screens only the maternal or zygotic contribution was affected52. As a consequence, in many cases zygotic functions were not detected because maternal contribution of gene products rescued zygotic mutations throughout Drosophila embryogenesis, leading to death after cuticle formation.

Overcoming the candidate gene approach. Most Drosophila segmentation genes have already been tested in Tribolium, revealing that a different gene set is involved in axis formation46,53,54. Overcoming this exhausted candidate gene approach, the iBeetle screen indeed identied novel players. For instance, knockdown of the homeobox gene Tc-homeobrain elicited a mirror image duplication of the abdomen similar to the

Drosophila bicaudal phenotype and was found to be one of the earliest anteriorly expressed zygotic genes (Fig. 4a,e,f). No cuticle phenotype has been described for Drosophila homeobrain, and so far no Tribolium gene had been found to elicit a bicaudal phenotype. SoxNeuro is required for Drosophila central nervous system development55,56, but was not reported to be involved in cuticle patterning. In the iBeetle screen, Tc-SoxNeuro knockdown resulted in a dorsalized phenotype, suggesting that Tc-SoxNeuro inuences early dorsoventral patterning (Fig. 4c). Another example is Drosophila Dscam, which is a cell adhesion gene with extensive alternative splicing and which is known to act in neurogenesis and immunity57,58. The iBeetle screen identied potential additional essential roles of Tc-Dscam in sensory organ formation (Fig. 4d).

New elds of research. Odoriferous stink glands play crucial roles in insect defence and communication but are not present in Drosophila59,60. In the iBeetle screen, we identied 32 genes with relevant phenotypes, including the absence of the gland contents, altered colour and composition of secretions, or melanosis. Interestingly, only 5 among these 32 genes showed an enrichment of greater than fourfold in odoriferous gland transcriptomes compared with mid-abdominal tissues60, illustrating that a phenotypic screen can only partially be replaced by transcriptomic approaches (see Supplementary

a b

T1 A1 A8

u

P

wt wt

T

c d

Tc-SoxN (iB_05549) RNAi

Tc-DSCAM (iB_05264) RNAi

Tc-homeobrain (iB_04564)

e f

A1

A

A

A A A

A8

Tc-homeobrain (iB_04564) RNAi

wt Tc-Rbm24 (iB_00289) RNAi

g h

Figure 4 | Embryonic phenotypes. (a,b) Wild-type L1 cuticles, with head setae marked by circles (b). T1: rst thoracic segment; A1: rst abdominal segment; A8: eighth abdominal segment; U: urogomphi; P: pygopods. (c) Unexpectedly, Tc-SoxN RNAi led to a strongly dorsalized cuticle phenotype without clear axes (embryo: lled arrowheads; vitelline membrane: open arrowhead). (d) Tc-DSCAM RNAi induced the deletion of head setae. (e) Tc-homeobrain RNAi caused a bicaudal phenotype (mirror image abdomina). This function is not known from Drosophila homeobrain and in Tribolium no bicaudal phenotype has been described before. (f) Early anterior zygotic expression of Tc-homeobrain. (g,h) Tc-Rbm24 RNAi led to detached and shortened body wall muscles (wild type pattern in g). A muscle function is conserved in vertebrates while the ortholog was lost in Drosophila. Scale bars indicate 100 mm.

6 NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822 ARTICLE

a

b

3 2 1

c

I II

IV

3 2 1

**

* *

**

III

III

wt

wt

wt

Tc-retained (iB_04537) RNAi

d

e

f

IV

IV

II

II

Tc-MED24 (iB_00555) RNAi

iB_04887 RNAi wt

wt

Tc-ATP7 (iB_02517) RNAi

i

j

g

h

CG8078 RNAi

Figure 5 | Postembryonic phenotypes. (a) Wild-type ovary stained for F-actin (red) and DNA (blue). Pro-oocytes (asterisks) become encapsulated by somatic follicle cells and separated by stalk cells (arrow). (d) Upon Tc-MED24 RNAi, egg chambers are misarranged, not separated by stalk cells (arrow) and subsequently they fuse (IV). (b,c) After Tc-retained-RNAi the three most distal antennomeres (13) of the adult antenna are fused. (e,f) RNAi against Tc-ATP7 led to strongly reduced odoriferous gland content and partially melanized secretions (white arrowheads; remnant of posterior abdominal cuticle marked by open arrowhead:). (gj) The knockdown of iB_04887 led to wing blisters in Tribolium pupae (arrowhead in h). Transgenic RNAi against the Drosophila ortholog showed the same phenotype, revealing a novel candidate for integrin mediated adhesion (arrowhead in j). Scale bars indicate 100 mm in (af) and 1 mm in (gj).

Fig. 4). For instance, the Tc-copper-transporting-ATPase-I (Tc-ATP7) is neither upregulated nor downregulated (Supplementary Fig. 5) but RNAi mediated knockdown (iB_02517) caused a reduced gland content and melanosis phenotype (Fig. 5f) and a loss of benzoquinones (Supplementary Fig. 6). Another emerging eld is the shaping of the adult body during typical holometabolous insect metamorphosis were larval epidermal cells are reused24. For instance, Tc-retained led to rounded female genitalia and the fusion of distal antennal segments (Fig. 5c). Interestingly, these anatomical features vary in Tenebrionids61, making this gene a good candidate for morphological evolutionary studies. Finally, one difference between Tribolium telotrophic oogenesis and Drosophila meroistic oogenesis is that in Tribolium germ line stem cells stop proliferating at larval stages while the somatic stem cell lineage remains active throughout life23. In Drosophila, both lineages remain active and dependent on each other, making it difcult to study the somatic lineage independently. In the iBeetle screen we identied several genes probably required for the somatic lineage, such as Tc-MED24 whose knockdown led to incomplete separation of egg chambers and a reduced number of follicle cells (Fig. 5d).

Tribolium as a screening platform. Wing blisters indicate the loss of adhesion between the two epithelial sheets of the wing blade. This phenotype has been used in Drosophila to identify components of integrin mediated adhesion62,63. We found 49 genes associated with wing blisters in our screen. Thirty-four were re-tested with non-overlapping fragments. All were conrmed in the screening strain but one led to lethality in another genetic background. Seventeen of these were previously annotated with GO terms connected to cytoskeletal function or adhesion, but 14 had unrelated functional annotations and 5 did not have any

functional annotation at all64. We tested 19 genes in Drosophila by transgenic RNAi knockdown, using two different wing disc driver lines. Out of seven genes without previous annotations with respect to cell adhesion or cytoskeletal function in Drosophila, four uncovered a wing blister or crumbled wing phenotype, indicating an involvement in Drosophila epithelial adhesion (See Supplementary Table 4). One example is the Tribolium ortholog of CG8078, which is predicted to be involved in tRNA thiomodication and was not implicated in cell adhesion before (Fig. 5h,j). Notably, some wing blister genes known from Drosophila were not recovered in the iBeetle screen owing to larval lethality before wing development (see Supplementary Table 5). Therefore, we injected at a later stage (L7) and indeed, two additional genes showed the wing blister phenotype. Hence, a new screen focused on injections at later larval stages is likely to reveal additional wing blister genes. In summary, the Drosophila and Tribolium screening platforms appear to reveal different subsets of genes involved in a common cellular process, suggesting that the use of Tribolium as alternative screening platform may reveal novel players relevant for general insect biology.

We identied many known and novel genes required for muscle development. In Tc-Rbm24 RNAi embryos, for instance, the muscles form small, round syncytia, which eventually seem to decay (Fig. 4h). The gene codes for an RNA binding protein of the RRM superfamily. Of note, it lacks an ortholog in Drosophila but vertebrate orthologs are active in myogenesis65,66. Hence, our data show the conservation of myogenic function of Rbm24 in Bilateria.

DiscussionHere we showed that the iBeetle RNAi screen achieved efcient and sensitive detection of novel gene functions at four developmental stages and in several processes in the red our

NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822

beetle Tribolium. The reproducibility varied depending on the phenotypic class. Highly penetrant phenotypes with simple readout (for example, wing blister, lethality) reproduced with more than 95% while the reproduction rate of embryonic developmental defects was signicantly lower (about 60%). One reason is that iBeetle was designed as rst pass screen without replicates while in cell based RNAi screens replicates and different reagents targeting the same gene function are usually used to minimize the portion of false positive data. Another important reason is that in Tribolium moderate knockdown of housekeeping genes during embryogenesis leads to abortion of embryogenesis at different stages leading to diverse cuticle defects. Hence, any off-target effect affecting a housekeeping gene is likely to result in diverse cuticle defects. Hence, it is advisable to focus on phenotypes annotated with high penetrance (450%) and conrmation of phenotypes using non-overlapping dsRNA fragments is essential.

Our data show that a signicant portion of insect gene functions becomes apparent when using additional model systems for several reasons. First, some aspects of insect biology are more representative in Tribolium compared with Drosophila like for instance metamorphosis based on the re-use of larval epidermal cells. Indeed, we identify novel genes associated with this process, the orthologs of which might have lost their function in Drosophila due to its derived mode of metamorphosis. Second, the different characteristics of the screening procedures may lead to detection of different subsets of the gene set involved in a conserved process. For instance, we identied novel genes leading to a wing blister phenotype but not all genes known from Drosophila were recovered in our screen. This suggests that neither system appears to have the power of detecting the comprehensive gene set of a given process. Finally, 13% of the genes associated with phenotypes in our screen do not have orthologs in Drosophila and therefore need to be studied in other model systems.

In summary, the iBeetle screen helps to overcome several current limitations in insect functional genetics. First, the candidate gene approach prevailing in emerging model organisms cannot reveal novel gene functions. This is superseded by the hypothesis independent identication of genes. Second, new biological processes can now be investigated, the genetic bases of which remained obscure because they are too derived or not present in the y, or not amenable in the y due to technical constraints. Third, we show that Tribolium is a powerful complementary screening platform for basic processes that are being studied in other organisms, such as epithelial adhesion. Importantly, the dsRNA template library generated in this project facilitates future screens focused on additional topics. Taken together, the iBeetle screen aids in bridging the gap between large-scale gene discovery by next generation sequencing and the small-scale approaches used for uncovering the function of novel genes.

Methods

Selection of genes and dsRNA production. The genes to be knocked down were selected at random. However, we required that the predictions were well supported by cDNA sequence data. Further, we excluded genes close to neighbouring same-strand gene predictions to avoid the double injection of genes, which were erroneously annotated as two genes. As the coverage by RNA-seq was low at project start, the gene set may have been enriched in genes expressed above average. This may explain the higher portion of conserved genes in the set of screened genes (82% in the iBeetle screen versus 59% in the ofcial gene set; see Supplementary Table 1).

The target transcript sequences were analysed using the DEQOR algorithm to identify stretches with highest amount of sequences predicted to be efciently recognized by the RNAi machinery and the lowest amount of potential off-target sites. Fragments within these stretches were amplied by PCR with gene-specic primers, which were tagged with parts of the T7 (30 primer) and SP6 (50 primer)

promoter sequences (T7 tag 50-CTCACTATAGGGAGA-30; SP6 tag 50-TGACAC TATAGAAGTG-30). The sequences are available at the iBeetle-Base (http://ibeetle-base.uni-goettingen.de/

Web End =http://ibeetle- http://ibeetle-base.uni-goettingen.de/

Web End =base.uni-goettingen.de/ ) and in Supplementary Data 3. The products from these PCR reactions were used as templates in a second PCR using T7 and SP6-T7 primers to generate templates for bidirectional in vitro transcription using T7-RNA polymerase. (Note that due to the different tags on 50 and 30 end of the fragments, the iBeetle library can be used for the generation of in situ probes as well.) For quality control, all products of the second PCR were checked on a gel and sequenced. All dsRNAs were assessed for purity by gel electrophoresis, the concentration was measured by the Ribogreen Assay (Life Technologies), and adjusted to 1mg ml 1 using injection buffer. Template production was performed by Eupheria Biotech GmbH (Dresden). Overall, 7,200 in silico dened templates were channelled into the pipeline leading to the production of 6,147 templates of which 5,670 dsRNAs were produced and sent to the screening centres. The success rate was about 85% for template generation and 92% for dsRNA production, resulting in an overall success rate of about 79%. The average length of the templates was 479 bp.

The screening procedure. Animals for 24 experiments were reared and processed in parallel using the equipment developed by Berghammer et al.67. For the injections two transgenic lines were used, where EGFP enhancer traps marked larval (pig19) or adult muscles (D17) in the pupal and larval injection screen, respectively (see Supplementary Fig. 7). An X-linked transgene insertion (Xred, expressing DsRed in larval eyes) was used for sexing larvae. dsRNA solution in the volume 1 mg ml 1 was injected into 10 animals per experiment using a FemtoJet express device (Eppendorf). It was injected as much as possible without interfering with survival. For injection and morphological inspection, larvae and adult beetles were anaesthetized with ice or carbon dioxide. Pupae were afxed to microscope slides for injection using a double sided sticky tape or rubber based cement (Fixogum). Inspection of pupae and adult morphology and phenotype documentation was performed with epiuorescence stereomicroscopes (Leica M205 FA). The data were documented electronically during analysis using an online interface. The interface allowed the documentation of technical remarks and offered dropdown lists with controlled vocabularies for documentation. In addition, pictures were uploaded for the documentation of the annotations and remained linked to the respective annotation. This allowed displaying the relevant pictures in the search results. For analysis of embryonic muscles, living embryos were dechorionated, embedded in Halocarbon oil (Voltalef 10S) and analysed using upright uorescence microscopes (Zeiss, Jena). Fourty-one dsRNAs targeting genes with known phenotypes were randomly introduced as positive controls with a frequency of 1 or 2 per 24 injections. The rst and last of the 24 injections performed on 1 day were negative control injections (buffer). The throughput of the screen was 21 and 28 genes per week in the pupal and larval injection screen, respectively. The workow is shown in Supplementary Table 6; the schedules are shown in Supplementary Figs 8 and 9. All phenotypes shown in this paper were reproduced with non-overlapping fragments using the SB strain as a different genetic background (see Supplementary Table 7 for sequences). The data are available at http://ibeetle-base.uni-goettingen.de/

Web End =http://ibeetle-base.uni-goettingen.de/ 40.

GCMS of odoriferous stink gland contents. dsRNA was injected in animals at mid-pupa stage (SB strain). Injected pupae as well as uninjected control pupae were kept on whole-grain our at 32 C. Ten days after hatching both the prothoracic and the abdominal glands of one beetle were dissected and crushed pairwise in 50 ml methanol (Merck Millipore: SupraSolv). The samples were stored at 20 C

and subjected to gas chromatography mass spectrometry (GCMS) measurements within 48 h. A volume of 1 ml was loaded per sample by a split injector into a gas chromatograph (Agilent Technologies 6890N Network GC System) connected to a mass spectrometer (Agilent Technologies 5973 Network Mass Selective Detector). For chromatogram analysis the software MSD ChemStation D.02.00.275 (Agilent Technologies) was used.

Comparison of essential genes. The portions of essential (that is, lethal) genes in Drosophila that had been published previously vary with the assumed total number of genes at the time. For our comparison, we related the published absolute numbers of essential genes to the number of protein coding genes of the current Drosophila genome release (6.02). For Tribolium we reduced the numbers observed in the iBeetle screen by the estimated portion of false positives in our data set (26%; see above). See Supplementary Table 8 for numbers and calculations.

References

1. Nsslein-Volhard, C. & Wieschaus, E. Mutations affecting segment number and polarity in Drosophila. Nature 287, 795801 (1980).

2. St Johnston, D. The art and design of genetic screens: drosophila melanogaster. Nat. Rev. Genet. 3, 176188 (2002).

3. Sulston, I. A. & Anderson, K. V. Embryonic patterning mutants in Tribolium castaneum. Development 122, 805814 (1996).

4. Maderspacher, F., Bucher, G. & Klingler, M. Pair-rule and gap gene mutants in the our beetle Tribolium castaneum. Dev. Genes Evol. 208, 558568 (1998).

8 NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822 ARTICLE

5. Pultz, M. A. et al. A Genetic screen for zygotic embryonic lethal mutations affecting cuticular morphology in the wasp Nasonia vitripennis. Genetics 154, 12131229 (2000).

6. Trauner, J. et al. Large-scale insertional mutagenesis of a coleopteran stored grain pest, the red our beetle Tribolium castaneum, identies embryonic lethal mutations and enhancer traps. BMC Biol. 7, 73 (2009).

7. Boutros, M. et al. Genome-wide RNAi analysis of growth and viability in Drosophila cells. Science 303, 832835 (2004).

8. Dietzel, G. et al. A genome-wide transgenic RNAi library for conditional gene inactivation in Drosophila. Nature 448, 151157 (2007).

9. Fraser, A. G. et al. Functional genomic analysis of C. elegans chromosome I by systematic RNA interference. Nature 408, 325330 (2000).

10. Gonczy, P. et al. Functional genomic analysis of cell division in C. elegans using RNAi of genes on chromosome III. Nature 408, 331336 (2000).

11. Sonnichsen, B. et al. Full-genome RNAi proling of early embryogenesis in Caenorhabditis elegans. Nature 434, 462469 (2005).

12. Bucher, G. & Klingler, M. Divergent segmentation mechanism in the short germ insect Tribolium revealed by giant expression and function. Development 131, 17291740 (2004).

13. Cerny, A. C., Bucher, G., Schroder, R. & Klingler, M. Breakdown of abdominal patterning in the Tribolium Kruppel mutant jaws. Development 132, 53535363 (2005).

14. Cerny, A. C., Grossmann, D., Bucher, G. & Klingler, M. The Tribolium ortholog of knirps and knirps-related is crucial for head segmentation but plays a minor role during abdominal patterning. Dev. Biol. 321, 284294 (2008).

15. Chen, G., Handel, K. & Roth, S. The maternal NF-kappaB/dorsal gradient of Tribolium castaneum: dynamics of early dorsoventral patterning in a short-germ beetle. Development 127, 51455156 (2000).

16. Kittelmann, S., Ulrich, J., Posnien, N. & Bucher, G. Changes in anterior head patterning underlie the evolution of long germ embryogenesis. Dev. Biol. 374, 174184 (2013).

17. Choe, C. P., Miller, S. C. & Brown, S. J. A pair-rule gene circuit denes segments sequentially in the short-germ insect Tribolium castaneum. Proc. Natl Acad. Sci. USA 103, 65606564 (2006).

18. Bolognesi, R., Farzana, L., Fischer, T. D. & Brown, S. J. Multiple Wnt genes are required for segmentation in the short-germ embryo of Tribolium castaneum. Curr. Biol. 18, 16241629 (2008).

19. Brown, S. et al. Implications of the Tribolium Deformed mutant phenotype for the evolution of Hox gene function. Proc. Natl Acad. Sci. USA 97, 45104514 (2000).

20. Posnien, N., Koniszewski, N. D. B., Hein, H. & Bucher, G. Candidate gene screen in the red our beetle Tribolium reveals six3 as ancient regulator of anterior median head and central complex development. PLoS Genet. 7, e1002416 (2011).

21. Schaeper, N. D., Pechmann, M., Damen, W. G., Prpic, N. M. & Wimmer, E. A. Evolutionary plasticity of collier function in head development of diverse arthropods. Dev. Biol. 344, 363376 (2010).

22. Klingler, M. Tribolium. Curr. Biol. 14, R639R640 (2004).23. Trauner, J. & Buning, J. Germ-cell cluster formation in the telotrophic meroistic ovary of Tribolium castaneum (Coleoptera, Polyphaga, Tenebrionidae) and its implication on insect phylogeny. Dev. Genes Evol. 217, 1327 (2007).

24. Snodgrass, R. E. Smithsonian Miscellaneous Collections, vol. 122 (Literary Licensing, Washington, 1954).

25. Panlio, K. A. & Roth, S. Epithelial reorganization events during late extraembryonic development in a hemimetabolous insect. Dev. Biol. 340, 100115 (2010).

26. Baumer, D., Strohlein, N. M. & Schoppmeier, M. Opposing effects of Notch-signaling in maintaining the proliferative state of follicle cells in the telotrophic ovary of the beetle Tribolium. Front Zool 9, 15 (2012).

27. Posnien, N., Schinko, J. B., Kittelmann, S. & Bucher, G. Genetics, development and composition of the insect head - A beetles view. Arthropod Struct. Dev. 39, 399410 (2010).

28. Panlio, K. A., Oberhofer, G. & Roth, S. High plasticity in epithelial morphogenesis during insect dorsal closure. Biol. Open 2, 11081118 (2013).

29. Richards, S. et al. The genome of the model beetle and pest Tribolium castaneum. Nature 452, 949955 (2008).

30. Berghammer, A. J., Klingler, M. & Wimmer, E. A. A universal marker for transgenic insects. Nature 402, 370371 (1999).

31. Schinko, J. B. et al. Functionality of the GAL4/UAS system in Tribolium requires the use of endogenous core promoters. BMC Dev. Biol. 10, 53 (2010).

32. Sarrazin, A. F., Peel, A. D. & Averof, M. A segmentation clock with two-segment periodicity in insects. Science 336, 338341 (2012).

33. Brown, S. J., Mahaffey, J. P., Lorenzen, M. D., Denell, R. E. & Mahaffey, J. W. Using RNAi to investigate orthologous homeotic gene function during development of distantly related insects. Evol. Dev. 1, 1115 (1999).

34. Bucher, G., Scholten, J. & Klingler, M. Parental RNAi in Tribolium (Coleoptera). Curr. Biol. 12, R85R86 (2002).

35. Peel, A. D. et al. Tc-knirps plays different roles in the specicationof antennal and mandibular parasegment boundaries and is regulated bya pair-rule gene in the beetle Tribolium castaneum. BMC Dev. Biol. 13, 25 (2013).

36. Tomoyasu, Y. & Denell, R. E. Larval RNAi in Tribolium (Coleoptera) for analyzing adult development. Dev. Genes Evol. 214, 575578 (2004).

37. Henschel, A., Buchholz, F. & Habermann, B. DEQOR: a web-based tool for the design and quality control of siRNAs. Nucleic Acids Res. 32, W113W120 (2004).

38. Mungall, C. J. et al. Integrating phenotype ontologies across multiple species. Genome Biol. 11, R2 (2010).

39. Donitz, J. et al. TrOn: an anatomical ontology for the beetle Tribolium castaneum. PLoS One 8, e70695 (2013).

40. Dnitz, J. et al. iBeetle-Base: a database for RNAi phenotypes in the red our beetle Tribolium castaneum. Nucleic. Acids. Res. 43, D721D726 (2015) doi:http://dx.doi.org/10.1093/nar/gku1054

Web End =10.1093/nar/gku1054 .

41. Kitzmann, P., Schwirz, J., Schmitt-Engel, C. & Bucher, G. RNAi phenotypes are inuenced by the genetic background of the injected strain. BMC Genomics 14, 5 (2013).

42. Chandler, C. H., Chari, S., Tack, D. & Dworkin, I. Causes and consequences of genetic background effects illuminated by integrative genomic analysis. Genetics 196, 13211336 (2014).

43. Arakane, Y. et al. The Tribolium chitin synthase genes TcCHS1 and TcCHS2 are specialized for synthesis of epidermal cuticle and midgut peritrophic matrix. Insect Mol. Biol. 14, 453463 (2005).

44. Ni, J. Q. et al. A genome-scale shRNA resource for transgenic RNAi in Drosophila. Nat. Methods 8, 405407 (2011).

45. Ni, J. Q. et al. A Drosophila resource of transgenic RNAi lines for neurogenetics. Genetics 182, 10891100 (2009).

46. Fu, J. et al. Asymmetrically expressed axin required for anterior development in Tribolium. Proc. Natl Acad. Sci. USA 109, 77827786 (2012).

47. Van der Zee, M., Stockhammer, O., von Levetzow, C., Nunes da Fonseca, R. & Roth, S. Sog/Chordin is required for ventral-to-dorsal Dpp/BMP transport and head formation in a short germ insect. Proc. Natl Acad. Sci. USA 103, 1630716312 (2006).

48. Kotkamp, K., Klingler, M. & Schoppmeier, M. Apparent role of Tribolium orthodenticle in anteroposterior blastoderm patterning largely reects novel functions in dorsoventral axis formation and cell survival. Development 137, 18531862 (2010).

49. Mullins, M. C., Hammerschmidt, M., Haffter, P. & Nsslein-Volhard, C. Large-scale mutagenesis in the zebrash: in search of genes controlling development in a vertebrate. Curr. Biol. 4, 189202 (1994).

50. Nusslein-Volhard, C. Of ies and shes. Science 266, 572574 (1994).51. Wieschaus, E., Nusslein-Volhard, C. & Jurgens, G. Mutations affecting the pattern of the larval cuticle in Drosophila melanogaster. 3. Zygotic loci on the X-chromosome and 4th chromosome. Wilhelm Rouxs Arch. Dev. Biol. 193, 296307 (1984).

52. Schpbach, T. & Wieschaus, E. Female sterile mutations on the second chromosome of Drosophila melanogaster. I. Maternal effect mutations. Genetics 121, 101117 (1989).

53. Brown, S. et al. A strategy for mapping bicoid on the phylogenetic tree. Curr. Biol. 11, R43R44 (2001).

54. Schoppmeier, M., Fischer, S., Schmitt-Engel, C., Lhr, U. & Klingler, M. An ancient anterior patterning system promotes caudal repression and head formation in Ecdysozoa. Curr. Biol. 19, 18111815 (2009).

55. Buescher, M., Hing, F. S. & Chia, W. Formation of neuroblasts in the embryonic central nervous system of Drosophila melanogaster is controlled by SoxNeuro. Development 129, 41934203 (2002).

56. Overton, P. M., Meadows, L. A., Urban, J. & Russell, S. Evidence for differential and redundant function of the Sox genes Dichaete and SoxN during CNS development in Drosophila. Development 129, 42194228 (2002).

57. Schmucker, D. et al. Drosophila Dscam is an axon guidance receptor exhibiting extraordinary molecular diversity. Cell 101, 671684 (2000).

58. Watson, F. L. et al. Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science 309, 18741878 (2005).

59. Eisner, T., Mchenry, F. & Salpeter, M. M. Defense mechanisms of arthropods.XV. Morphology of the quinone-producing glands of a tenebrionid beetle (ELEODES longicollis lec.). J. Morphol. 115, 355399 (1964).60. Li, J. et al. Odoriferous defensive stink gland transcriptome to identify novel genes necessary for quinone synthesis in the red our beetle, Tribolium castaneum. PLoS Genet. 9, e1003596 (2013).

61. Beutel, R. G., Leschen, R. A. B., Lawrence, J. F. &lipiski, S. A. Coleoptera, beetlesVol. 2 (De Gruyter, 2010).

62. Prout, M., Damania, Z., Soong, J., Fristrom, D. & Fristrom, J. W. Autosomal mutations affecting adhesion between wing surfaces in Drosophila melanogaster. Genetics 146, 275285 (1997).

NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8822

63. Walsh, E. P. & Brown, N. H. A screen to identify Drosophila genes required for integrin-mediated adhesion. Genetics 150, 791805 (1998).

64. Gelbart, W. M. et al. FlyBase: a Drosophila database. The FlyBase consortium. Nucleic Acids Res. 25, 6366 (1997).

65. Yang, J. et al. RBM24 is a major regulator of muscle-specic alternative splicing. Dev. Cell 31, 8799 (2014).

66. Poon, K. L. et al. RNA-binding protein RBM24 is required for sarcomere assembly and heart contractility. Cardiovasc. Res. 94, 418427 (2012).

67. Berghammer, A., Bucher, G., Maderspacher, F. & Klingler, M. A system to efciently maintain embryonic lethal mutations in the our beetle Tribolium castaneum. Dev. Genes Evol. 209, 382389 (1999).

Acknowledgements

We thank Deutsche Forschungsgemeinschaft for funding this project (FOR1234 iBeetle) and the Georg-August-University Gttingen and the Friedrich-Alexander-University Erlangen-Nrnberg for nancial and logistical support. We acknowledge support by the Open Access Publication Funds of the Georg-August-University Gttingen.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: M.T. and F.B. are afliated with Eupheria Biotech GmbH. All other authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Schmitt-Engel, C. et al. The iBeetle large-scale RNAi screen reveals gene functions for insect development and physiology. Nat. Commun. 6:7822 doi: 10.1038/ncomms8822 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

Author contributions

C.S.-E., U.M., J.S., D.S., N.S., N.T. and V.A.D. have contributed the vast bulk of the data: they performed the screening procedure, analysed and documented the data, annotated the phenotypes and contributed data to the NOF analysis. D.G. and T.R. performed the screening procedure, analysed and documented the data, annotated the phenotypes. C.S.-E., N.D.B.K., M.S., M.K. and G.B. developed the screening procedure. M.T., J.D. and B.M. developed the database collaborating with C.S.-E. and G.B. L.G. and M.S. dened the transcripts for template production. M.T. and F.B. dened and produced the iB

templates and produced the dsRNAs. I.S., E.K., J.T., S.K. helped with different aspects of the screen. S.L., J.U., J.S., Y.H., K.A.P. and R.S. contributed to the NOF analysis. K.A.P., R.S., M.F., S.R., E.A.W. were involved in the design of the study and the analysis of the data. C.S.-E., M.K., M.S. and G.B. designed the study, led the screen, headed the analysis and wrote the paper

10 NATURE COMMUNICATIONS | 6:7822 | DOI: 10.1038/ncomms8822 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.