The origin and integration of novel traits are fundamental processes during the developmental evolution of complex organisms. Yet how novel traits integrate into pre-existing contexts remains poorly understood. Beetle horns represent a spectacular evolutionary novelty integrated within the context of the adult dorsal head, a highly conserved trait complex present since the origin of insects. We investigated whether otd1/2 and six3, members of a highly conserved gene network that instructs the formation of the anterior end of most bilaterians, also play roles in patterning more recently evolved traits. Using ablation-based fate-mapping, comparative larval RNA interference (RNAi) and transcript sequencing, we found that otd1/2, but not six3, play a fundamental role in the post-embryonic formation of the adult dorsal head and head horns of Onthophagus beetles. By contrast, neither gene appears to pattern the adult head of Tribolium flour beetles even though all are expressed in the dorsal head epidermis of both Onthophagus and Tribolium. We propose that, at least in beetles, the roles of otd genes during post-embryonic development are decoupled from their embryonic functions, and that potentially non-functional post-embryonic expression in the dorsal head facilitated their co-option into a novel horn-patterning network during Onthophagus evolution.
Onthophagus spp six3 and otd gene sequences
Sequence and associated annotations for O. taurus six3 and otd genes and O. sagittarius six3 and otd transcripts. Flat file Genbank format.
Onthophagus six3 and otd genes and primers.gb
Insect SIX protein family sequences and tree (data for Figure S1)
Aligned sequences and RaxML generated phylogram of selected insect SIX proteins.
InsectSIXtree.nex
Insect OTD protein sequences and tree (data for Figure S2)
Aligned sequence and FastTree generated phylogram of selected insect orthodenticle proteins.
InsectOTDtree.nex
R script to reproduce expression plots (Figures S3 and S8)
Annotated R script with code to read and plot raw data from qPCR and RNAseq estimations of stage specific and tissue specific levels of six3 and otd transcripts. Requires files "tcas_qPCR_stages.csv", "otau_rnaseq_stages.csv", "osag_rnaseq_stages.csv", "tcas_qpcr_tissues.csv" and "otau_rnaseq_tissues.csv".
Expression_plots.R
Tribolium castaneum qPCR by stage
Raw data table for T. castaneum stage-specific qPCR estimation of six3 and otd transcript levels. Use R script in "Expression_plots.R" for analysis and plotting.
tcas_qPCR_stages.csv
Onthophagus taurus RNAseq by stage
Raw data table for O. taurus stage-specific RNAseq estimation of six3 and otd transcript levels. Use R script in "Expression_plots.R" for analysis and plotting.
otau_rnaseq_stages.csv
Onthophagus sagittarius RNAseq by stage
Raw data table for O. sagittarius stage-specific RNAseq estimation of six3 and otd transcript levels. Use R script in "Expression_plots.R" for analysis and plotting.
osag_rnaseq_stages.csv
Tribolium castaneum qPCR by pupal tissue
Raw data table for T. castaneum tissue-specific qPCR estimation of six3 and otd transcript levels. Use R script in "Expression_plots.R" for analysis and plotting.
tcas_qpcr_tissues.csv
Onthophagus taurus RNAseq by pupal tissue
Raw data table for O. taurus tissue-specific RNAseq estimation of six3 and otd transcript levels. Use R script in "Expression_plots.R" for analysis and plotting.
otau_rnaseq_tissues.csv
R script to reproduce RNAi phenotype plots (Figures S4 and S5))
Annotated R script with code to read and plot raw data from phenotypic measurements of dsRNA treated beetles . Requires files "tcas_RNAi.csv", "otau_RNAi.csv" and "osag_RNAi.csv".
RNAiAnalyses.R
Tribolium castaneum RNAi phenotype data
Phenotypic measurements from dsRNA injected T. castaneum. For analysis and plotting using the R code in file "RNAiAnalyses.R".
tcas_RNAi.csv
Onthophagus taurus RNAi phenotype data
Phenotypic measurements from dsRNA injected O. taurus. For analysis and plotting using the R code in file "RNAiAnalyses.R".
otau_RNAi.csv
Onthophagus sagitarius RNAi phenotype data
Phenotypic measurements from dsRNA injected O. sagittarius. For analysis and plotting using the R code in file "RNAiAnalyses.R".
osag_RNAi.csv
Tribolium castaneum qPCR validation of RNAi KD efficiency (Figure S9)
Raw data from qPCR experiments to validate the knockdown efficiency of dsRNA injections in Tribolium castaneum.
RNAiknockdownefficiency_qPCR_rawdata.xlsx