Skip to main content

Genetic basis of allochronic differentiation in the fall armyworm

Abstract

Background

Very little is known on how changes in circadian rhythms evolve. The noctuid moth Spodoptera frugiperda (Lepidoptera: Noctuidae) consists of two strains that exhibit allochronic differentiation in their mating time, which acts as a premating isolation barrier between the strains. We investigated the genetic basis of the strain-specific timing differences to identify the molecular mechanisms of differentiation in circadian rhythms.

Results

Through QTL analyses we identified one major Quantitative trait chromosome (QTC) underlying differentiation in circadian timing of mating activity. Using RADtags, we identified this QTC to be homologous to Bombyx mori C27, on which the clock gene vrille is located, which thus became the major candidate gene. In S. frugiperda, vrille showed strain-specific polymorphisms. Also, vrille expression differed significantly between the strains, with the rice-strain showing higher expression levels than the corn-strain. In addition, RT-qPCR experiments with the other main clock genes showed that pdp1, antagonist of vrille in the modulatory feedback loop of the circadian clock, showed higher expression levels in the rice-strain than in the corn-strain.

Conclusions

Together, our results indicate that the allochronic differentiation in the two strains of S. frugiperda is associated with differential transcription of vrille or a cis-acting gene close to vrille, which contributes to the evolution of prezygotic isolation in S. frugiperda.

Background

Virtually all life on earth experiences a similar day-night cycle, yet some species have evolved to be day-active, while others are night-active. Even closely related species living in the same habitat can differ from each other in their daily activity rhythms [1]. The molecular mechanisms of circadian rhythms have been revealed in a broad range of eukaryotic species from algae [2] to mammals (reviewed in SM Reppert and DR Weaver [3]), and also in insects with the increasing number of sequenced insect genomes [4–7]. Surprisingly, despite this growing knowledge on the molecular basis of circadian rhythms, the evolution of differentiation in daily activity patterns is largely unexplored. In insects, differences in diurnal mating times have been found to prevent gene flow between populations [8–14]. By determining the genetic basis of allochronic differentiation between closely related species, or even between divergent populations within species, the initial steps causing differentiation in daily activity rhythms can be discovered, which is important for an understanding of the evolution of circadian rhythms on a micro-evolutionary time-scale.

An ideal model organism for the study of the evolution of circadian rhythms is the noctuid moth Spodoptera frugiperda (Lepidoptera: Noctuidae), as it consists of two naturally occurring morphologically identical strains that exhibit strain-specific timing of mating in the night [15, 16]. The so-called corn- and rice-strains seem to be in the process of ecological speciation in sympatry [17]. Although the hybridization rate is up to 16% in the field [18], the two strains do not merge into one panmictic population, which is probably prevented by a combination of different isolation barriers [17]. So far, three possible prezygotic mating barriers have been described in this species: a) differential host plant choice [19–23], b) strain-specific timing of mating in the night [15, 16], and c) female sex pheromone differences [24–26]. Recent studies have shown that host preference in the field is not as specific as previously thought [27–29]. Therefore, habitat isolation seems to be a relatively weak prezygotic mating barrier. Differences in female sex pheromone composition are also likely to constitute a weak prezygotic mating barrier [26, 30]. As both strains consistently differ in their timing of reproductive activity at night [15, 16], allochronic divergence seems to be a major barrier separating the two S. frugiperda strains. The corn-strain calls, mates and oviposits approximately three hours earlier than the rice-strain, with only a small overlapping time-window between the strains [15, 16].

Allochronic speciation due to seasonal timing differences has been suggested for several insect species, e.g. crickets [31, 32], fruit flies [33, 34] and mosquitoes [13]. However, surprisingly little research has been conducted on the importance and exact genetic changes underlying temporal speciation (reviewed in AT Groot [1]). Recently, a study by Kaiser et al. [35] determined the genomic basis of circadian and circalunar timing adaptations in the midges Clunio marinus. Different naturally occurring strains of C. marinus emerge at different time points in the circadian as well as the circalunar rhythm and mate and oviposit shortly after. The abundance of different splice variants of the calcium/calmodulin-dependent kinase II.1 (CaMKII.1) is associated with these allochronic differences.

Also the two S. frugiperda strains differ in their diurnal mating patterns and we hypothesize that genetic and/or expression differences in one or more clock genes underlie their differences in timing of reproductive activity.

In general, biological clocks are a network of genes and gene products that enhance and suppress each other in a rhythmic manner, entrained by environmental factors such as light, temperature or tides [36, 37]. Within insects, the clock gene network is best described in the vinegar fly Drosophila melanogaster, where the network consists of two interlocked feedback loops [36, 38]: one involving the genes vrille (vri), PAR-domain protein 1 (pdp1), clock (clk) and cycle (cyc); and the other incorporating period (per), timeless (tim), clk and cyc. In addition, kinases phosphorylate clock proteins (e.g. phosphorylation of PER by DOUBLETIME (DBT) and CASEIN KINASE2α (CK2α)) and facilitate their accumulation [36], while cryptochrome 1 (cry1) functions as a circadian photoreceptor. Most of these genes are also present in Lepidoptera [5, 39–42], and are thus good candidate genes that may underlie the timing differences between the corn- and the rice-strain (see Fig. 1). Additionally, a second cryptochrome, cryptochrome 2 (cry2), is present in Lepidoptera and able to repress CLK: CYC mediated transcription [40, 41, Fig. 1]. Even though most of the clock genes and their organization in the circadian clock are well conserved across taxonomic groups [43, 44], differences among insects orders have been highlighted [41, 43, 45], especially in the daily oscillation patterns of gene expression [46–49]. Identifying the exact genetic changes and the mechanism(s) underlying allochronic differentiation in behavior will shed light on how temporal differentiation may evolve.

Fig. 1
figure 1

Two feedback loops that define the circadian rhythm in Danaus plexippus. Adapted from S Zhan et al. [5]. In the first feedback loop (green proteins), the CLK:CYC dimer promotes the transcription of Per, Tim and Cry2. PER, TIM and CRY2 proteins enter the nucleus where a PER:CRY2 dimer inhibits the binding of the CLK:CYC dimer to E-Boxes and thus inhibits the expression of genes with E-Box promoters. In the second, modulatory feedback loop (yellow proteins), the CLK:CYC dimer promotes the transcription of Vri and /or Pdp1. VRI and PDP1 proteins enter the nucleus, where PDP1 promotes Clk transcription while VRI inhibits Clk transcription

In this study, we determined the genetic basis of the major prezygotic isolation barrier, i.e. differentiation in the diurnal mating patterns of the two strains of S. frugiperda. To identify the genomic loci (chromosomes) that explain most of the variance in timing of mating activity in an unbiased way, we first conducted quantitative trait locus (QTL) analyses using AFLP markers and RADtags. To determine which of the candidate clock genes map to the major Quantitative trait chromosome (QTC), we indirectly mapped these genes onto our genetic map by homologizing the S. frugiperda linkage map to the B. mori chromosomes and identifying the location of the candidate genes on the homologous chromosomes. This approach is possible due to the high syntheny between S. frugiperda and B. mori [50]. Since vrille was located on the major QTC, we determined differences in expression levels, as well as sequence differences of this gene, between the two strains, and compared this to differences in the other main clock genes.

Methods

As S. frugiperda is a non-model organism, no assembled genome was available when we started investigating the genetic basis of the differences in mating time of the two strains. Hence, this manuscript presents results that were obtained with a variety of methods to overcome this lack of information. In the course of experiments, the genome of both the corn- and the rice-strain became available to the members of The Fall armyworm International Public Consortium, which facilitated genomic comparisons and expression analysis of all main clock genes. An overview of the chronology of the experiments is depicted in Additional file 1.

Insects

Individuals used for the QTL analysis descended from > 200 rice-strain larvae and > 100 corn-strain larvae, collected from different fields in Florida in 2003 and 2004, respectively (Additional file 2). These populations were reared for 10 (corn-strain) and 21 (rice-strain) generations in mass culture at the USDA-ARS in Gainesville, FL, before shipment to the Max Planck Institute for Chemical Ecology (MPICE) in 2007. These populations were also used by G Schöfl, DG Heckel and AT Groot [16]. We refer to these populations as CL1 and RL1 (Additional file 2). Unfortunately, these two populations died after six years of laboratory rearing. Therefore, we established new laboratory populations starting with ~ 300 larvae collected in Florida (rice-strain) and Puerto Rico (corn-strain) in 2010 (Additional file 2), which were shipped directly to MPICE, where all adults were screened for strain-specific COI markers [23], and separated accordingly into strain-specific colonies. We refer to these populations as CL2 and RL2 (Additional file 2). All populations were reared in climate chambers with reversed light:dark (L:D) cycle (photophase started at 10 pm, scotophase started at noon) and 14:10 L:D photoperiod at 26 °C and 70% RH. Adults were fed with a 10% honey-water solution and random single-pair-matings were set up to maintain the populations and minimize inbreeding.

Generation of backcrosses

For the QTL analysis, we generated female-informative backcrosses (Additional file 3). Single pair matings between pure corn- and rice-strain individuals were performed to obtain F1 hybrids (denoted CR from corn mothers and RC from rice mothers). Hybrid females were then backcrossed to pure rice-strain males to produce different backcross families (Additional file 3). Two backcross families (BCs) were used for the QTL analysis (BC_A: RCxR, BC_B: CRxR) (the first two letters of a backcross refer to the mother, the last letter to the father). The two rice-strain fathers used to generate both backcrosses were kin.

Phenotyping backcross families

To determine the phenotype for the QTL analysis, we observed the mating behavior of a) pure strain individuals in intra-strain (CxC, RxR) and inter-strain matings (CxR, RxC), b) hybrid females backcrossed to pure strain males (CRxC, CRxR, RCxC, RCxR), and c) female backcross offspring crossed to pure strain males (CR-RxC, CR-RxR, RC-RxC, RC-RxR). The observations of mating behavior were performed as described by [16] and summarized here. One to four day-old virgin females and males were set up in single pairs in clear plastic cups (16 oz.) and provided with 10% honey solution. All matings were set up simultaneously and placed in a walk-in climate chamber (26 °C, 70% RH, L:D 14:10) two hours before scotophase. In total, 320 to 400 couples were observed throughout the scotophase and one hour into photophase (in total 11 h), with a 30 min interval, i.e. each couple was observed once every 30 min. All pairs were observed for three consecutive nights starting at the first day of the mating. The onset time of the first mating was the phenotype used for the timing QTL analysis. After observation, all individuals were frozen at -80 °C for further genetic analysis.

Genetic map construction

DNA of 90 randomly chosen backcross females (44x RC-R, 46x CR-R) as well as of their parents and grandparents were used to generate AFLP markers. Female backcross individuals were chosen to construct the genetic map and conduct the QTL analysis, as Schöfl et al. 2009 [16] showed that the onset time of mating was mainly influenced by the female mating partner (i.e. mating time was significantly different between corn-strain and rice-strain females, irrespective of the strain-identity of their mating partner). After DNA extraction, AFLP markers were generated as described in Groot et al. [51], and summarized here: 200 ng DNA of each sample was digested with EcoRI and MseI (New England Biolabs, Ipswich, MA, USA), and EcoRI- and MseI-adapters were ligated to the fragments, preamplified and selectively amplified with different EcoRI- and MseI-primer combinations (Additional file 4). The generated AFLP fragments were analyzed on a 6.5% polyacrylamide gel using a LI-COR 4300 DNA analyzer (LI-COR Biosciences, Lincoln, NE, USA). AFLP gels were scored with AFLP-Quantar Pro 1.0 (KeyGene, Wageningen, The Netherlands). To identify corn-strain specific markers, we scored markers that were present in the corn-strain grandparent (C grandmother or grandfather), the hybrid mother (RC or CR), and half of the offspring females (heterozygote females), but absent in the rice-strain grandparent (R), the backcross male (R), and the homozygote backcross (CR-R and RC-R) females. For identification of rice-strain specific markers, we scored markers present in the rice-strain grandparent, the hybrid mother and the homozygote offspring females, but absent in the corn-strain grandparent, the father and the heterozygote backcross females. All markers were converted to the same phase by inverting the absence/presence patterns of all rice-strain specific markers.

After scoring at least 450 markers, we constructed a linkage map for each BC with MapMaker 3.0 (http://www.broadinstitute.org/ftp/distribution/software/mapmaker3/). Markers were clustered into linkage groups (LG) using a LOD of 4.5. In BC_A, 30 LGs were identified that refer to the 30 autosomes in a backcross family, as there is no crossing over in lepidopteran females [52]. In BC_B, 29 LGs were identified. The chromosome names (chromosome 1 to 30) were chosen arbitrarily for each linkage map, so that the same numbers in the different linkage maps are not necessarily homologous. Markers present in both backcrosses (Additional file 4) were used to homologize the chromosomes of these backcrosses.

QTL analysis

A QTL analysis based on AFLP markers was conducted with two female-informative backcross families. A total of 465 (in BC family A) and 514 (in BC family B) informative AFLP markers were used to identify the 30 S. frugiperda autosomes (Additional file 4). Each chromosome consisted of at least two AFLP markers from different primer combinations up to a maximum of 26 markers. Seventeen markers in each BC family did not map to any linkage group. Three chromosomes (corresponding to three linkage groups in BC_A and 2 linkage groups in BC_B) could not be homologized between the two linkage maps.

To identify candidate QTC, each chromosome was tested for a significant difference in the phenotype (i.e. onset time of first mating) between the homo- and heterozygote backcross females. The two backcrosses were also combined for this analysis, to increase the sample size and thus the possibility to detect QTC. Because of the absence of crossing-over in female Lepidoptera [52], each identified QTC corresponds to an individual chromosome, on average 1/31 of the genome. Statistical analysis was performed with R 2.5.0 (R-Development-Core-Team, 2007) and SAS® software (SAS institute, Cary, NC, USA, 2002-2008). To assess how much of the variance is explained by the different QTC (R2 value) we conducted a two-sided t-test and a GLM (The GLM procedure of the SAS software, with onset time of the first copulation as dependent variable). Chromosomes with a significant correlation (P < 0.05) were considered a QTC.

Homologizing linkage maps to Bombyx mori chromosomes, using RAD markers

To identify candidate genes on the QTC, the linkage map was homologized to the reference genome of B. mori, using restriction site associated DNA (RAD) analysis (see Baxter et al [53] and Groot et al. [54]). DNA of parents, female grandparents and 11 backcross individuals per backcross family was digested with the Sbf1 restriction enzyme (New England Biolabs, Ipswich, MA, USA), barcoded, pooled, sheared and amplified, following the procedure described in Groot et al [54]. The pool was paired-end sequenced (50 bp fw, 50 bp rev) by FASTERIS (Geneva, Switzerland) with a HiSeq Illumina sequencer, resulting in 76 million reads. The reads were separated by barcodes into pools per individual and filtered for quality (q10 = 99%). On average, there were 5-10 different paired-end reads per forward read (Additional file 5). The segregation patterns that were obtained with the AFLP markers for different linkage groups were utilized to identify RAD markers segregating in the same pattern with RAD tools [53]. The AFLP markers showed a specific presence (1) /absence (0) pattern in backcross individuals for each chromosome, thus RAD markers showing the same 1/0 pattern in the backcross individuals were identified as belonging to the same specific linkage group. All sequences matching an AFLP segregation pattern were pooled across the individuals, after which the paired-end sequences were retrieved, resulting in 30 FASTA files (one file per chromosome). Each group was assembled into RAD contigs using CLC Genomics Workbench (CLC bio version 5.0.1; www.clcbio.com). Sequences were trimmed for length and quality with standard settings (nucleotide mismatch cost = 2; in/del cost = 2; length fraction = 0.35; similarity = 0.9; when bases conflicted, the base with highest frequency was chosen) and assembled de novo. Contig sizes ranged from 89 – 598 bp, with the majority of contigs being 100 – 250 bp long.

Resulting contigs from the paired-end RAD sequences were BLASTed against the scaffolds of the corn-variant assembly 3.0 and the rice-variant assembly 1.0 of the S. frugiperda genome [55], http://bipaa.genouest.org/is/lepidodb/spodoptera_frugiperda/. The scaffolds with the best BLAST hits were BLASTed in SilkDB (http://silkworm.genomics.org.cn/) and KAIKObase [56], http://sgp.dna.affrc.go.jp/KAIKObase/]. We considered the S. frugiperda chromosomes and B. mori chromosomes to be homologous when a) multiple contigs of the same Sf chromosome produced significant BLAST hits (e-value < E-10) to the same Bm chromosome, or b) in cases where multiple Bm chromosomes hit contigs of one Sf chromosome, the hit with the lowest e-value was chosen (Additional files 6 and 7).

After homologizing the S. frugiperda chromosomes to the B. mori chromosomes, we assessed the location of candidate genes involved in the circadian rhythm (Fig. 1), using KAIKObase (http://sgp.dna.affrc.go.jp/KAIKObase) and the homology table (Additional file 6). The position of vri on the timing QTC Sf_C25 (Bm_C27) was verified by mapping it via single nucleotide polymorphisms (SNPs) to the combined linkage map, using the segregation patterns of the SNPs in the grandparents, parents, and 8 backcross females of both backcross families (see Additional file 8).

Expression analysis of the main candidate gene, vrille

To determine strain-specific expression differences and an overview of the daily oscillation in the candidate gene vri, as well as the main other circadian clock genes period, timeless, cryptochrome 2, PAR-domain protein 1, clock and cycle (Fig. 1), we conducted two reverse transcription-quantitative real-time PCR (RT-qPCR) experiments with mRNA from heads of female S. frugiperda of both strains. We chose female heads for this experiment, as the mating process is started by female calling, followed by male courtship and then copulation. Thus, the mating time is influenced significantly more by the female partner than by the male partner [16]. In the first experiment, 15 females of both strains were transferred from the rearing cups to a 10 ml Falcon tube, immediately frozen in liquid nitrogen and kept at -80 °C, which was repeated every hour for 24 h. RNA was isolated from three pools of five heads, providing three biological replicates per strain per time point. RNA extraction, cDNA synthesis and RT-qPCR reaction were conducted, as described in AT Groot et al. [54] and summarized here. Pools of five heads were ground with mortar and pestle in liquid nitrogen, RNA was isolated using Direct-zol™ RNA MiniPrep (Zymo Research corp.), and DNase was digested by adding 10 μl 10x Turbo DNase buffer and 1 μl Turbo DNase (Ambion, LIFE TECHNOLOGIES, Darmstadt, Germany). cDNA was synthesized from 1000 ng RNA using Verso cDNA synthesis kit (Thermo Fisher Scientific, Schwerte, Germany). RT-qPCR experiments were conducted with 5 ng cDNA per reaction, 2 technical replicates on each plate, using ABsolute Blue QPCR SYBR Green Low Rox Mix (Thermo Fisher Scientific, Schwerte, Germany) and Applied Biosystems 7500 Real-Time PCR (ThermoFisher Scientific) (see Additional file 8 for details). The Elongation Initiation Factor 1α (eIF1α) was used as the reference gene and amplified for all samples. eIF1α was chosen as the reference gene as it showed the most stable expression over the different time points in a pre-experiment. Relative expression levels were calculated as copy numbers per 1000 copies eIF1α.

To verify our first RT-qPCR results, we conducted a second qPCR experiment where we focused on the 12 most relevant time points in which the two strains showed differences in clock expression levels, i.e. from 5 h before (-5 h) until 6 h into scotophase (+6 h). For this second experiment, we again observed mating couples for one night, as described above, and extracted RNA extraction in the second night, for which we chose corn-strain females that showed reproductive behavior (female calling, copulation) early in the night and rice-strain females that showed reproductive behavior later in the night. Every hour, 6 - 10 females of both strains were transferred from the rearing cups to a 10 ml Falcon tube, immediately frozen in liquid nitrogen and kept at -80 °C. This time, RNA was isolated individually from six heads, providing six biological replicates per strain per time point. RNA extraction, cDNA synthesis and RT-qPCR reaction were conducted on the individual samples. Heads were ground with mortar and pestle in liquid nitrogen, RNA was isolated using innuPREP RNA Mini Kit (Analytik Jena, Germany). RT-qPCR experiments were conducted with 5 ng cDNA per reaction, this time 3 technical replicates on each plate, using 5X Hot firepol EvaGreen® qPCR mix plus (ROX) (Thermo Fisher Scientific, Schwerte, Germany) and a Bio-Rad CFX machine (Biotum) (see Additional file 8 for details). eIF1A was again used as the reference gene and amplified for all samples. Relative expression levels were calculated as above.

Statistical differences in strain-specific expression levels of the clock genes were tested in R (R Development Core Team 2010). For each time point, we tested the normality of the data with a Shapiro-Wilk test (shapiro.test function under R). To ensure that the variances were equal in both rice and corn dataset, we ran a Fisher test (var.test function under R). Finally, we used a Student test to compare the expression levels of two samples for each gene (t.test function in R) with a Bonferroni correction for multiple comparisons.

We also tested for significant difference between the whole expression pattern of the corn-strain and the rice-strain for individual clock genes. For this, we first used a Shapiro-Wilk normality test for each dataset (dataset = expression data 1 gene for one strain for each timepoint) followed by an F test to compare the variances of both strains per clock gene. Finally, we used a Two Sample t-test to test for a significant difference between the corn-strain expression and the rice-strain expression for each clock gene.

All graphs were made in R with the ggplot2 package [57].

Sequence differences of the main candidate gene, vrille

To assess strain-specific sequence differences in vri, the sequence of the gene was established stepwise. First degenerate primers based on insect ESTs and genomic sequences (gb|AY526608.1, gb|AY576272.1, gb|AADK01019845.1) were used to obtain partial sequences. After obtaining the sequences, primers were designed to sequence further. The DNA Walking SpeedUp™ Kit II (SEEGENE, Eschborn, Germany) was used to obtain the sequence upstream of the coding sequence (see Additional file 9 for all primers used). To determine exon/intron structure, the coding region was sequenced from cDNA. Subsequently, parts of the gene were sequenced in 88 different samples (including backcross individuals and corn- and rice-strain individuals from different regions; Additional file 2), using Sanger-sequencing and Sequencher 4.10.1 for analysis. All obtained sequences are available in GenBank (accession numbers KM675483 - KM675658). Subsequently, the S. frugiperda genome for both strains (http://www6.inra.fr/lepidodb/SfruDB) and an in-house RNAseq database of larval guts became available [55]. With the full length mRNA acquired from the RNAseq database and BLASTed against the genome, the full sequence of vri was obtained, including a large intron in the 5’UTR. The corn-strain genome was not complete in this region, thus two BAC clones (AUA0AAA25YL06FM1, AUA0AAA20YH15RM1) spanning the region were obtained from the Centre National de Ressources Génomiques Végétales (CNRGV, Toulouse, France) and shotgun sequenced using Sanger sequencing and Sequencher for analysis. Based on the alignment of the rice-strain genome from SfruDB and the BAC clone sequences, additional parts of the intron were sequenced in 12 corn-strain and 12 rice-strain individuals from the CL_1 and RL_1 populations, as well as the parental and F1 generations of the backcross families.

When the genome of S. frugiperda became available, the other main clock genes were annotated in both strains as part of the WGS project [55]. The coding sequences of each gene of both strains were aligned to identify polymorphisms between the strains. The protein domains were identified on http://prosite.expasy.org/ and using Protein BLAST on http://blast.ncbi.nlm.nih.gov/.

Results

QTL analysis

QTCs were identified by testing each linkage group for a significant association with the phenotypic trait, i.e. onset time of first mating. Resulting P-values and R 2-values (see Additional file 6) refers to the combined analysis of both BCs to use the biggest sample size possible. One QTC (Sf_C25, P < 0.0001, R 2 = 0.19) had a major effect on the variance in the strain-specific timing of mating. This QTC is homologous to Bombyx mori chromosome 27 (Bm_C27) and explained 19% of the variance of the onset time of first mating (Fig. 2). Bm_C27 is 14.5 Mb in size (52.8 cM) and represents 3.3% of the total B. mori genome [56, 58]. The difference in the onset time of mating between heterozygous and homozygous individuals (carrying only rice-strain copies of the chromosome) for Sf_C27 did not differ between the combined or individual analysis of backcrosses (Additional file 10). The LOD scores of all linkage groups are shown in Additional File 11.

Fig. 2
figure 2

a Location of clock-related genes on the homologous chromosomes of S. frugiperda and B. mori. Chromosomes that do not show a QTL are depicted light grey, minor QTL chromosomes dark grey and the major QTL chromosome black. Names of clock-related genes are positioned on the right of the chromosome they are situated on in B. mori. b Phenotype (onset time of first mating) for major QTL chromosome Sf_C25 (=Bm_C27) of pure corn-strain and rice-strain individuals vs. heterozygous CR and homozygous RR backcross individuals

Even though the major QTC Sf_C25 explained 19% of the variance between the strains, it only achieved a power of 0.55, i.e. only 55% of QTCs of this magnitude could be detected with our setup on average (Additional file 12). This is due to the small sample size and implies a chance of missing or underestimating minor QTCs. However, the finding of the same major QTC in two independent backcrosses strengthens our finding. We did detect three minor QTCs in the combined analysis that were not consistently present in both families when analyzed individually: Sf_C28 (Bm_C2, P = 0.014, R 2 = 0.08), Sf_C30 (Bm_C6, P = 0.0104, R 2 = 0.08) and Sf_C20 (Bm_C12, P = 0.023, R2 = 0.07).

Mapping the candidate genes

The main circadian rhythm genes are located on the following chromosomes (Additional file 6): per, clk, cyc and PdP1 on the sex chromosome (Bm_C01), jetlag on Bm_C3 (Sf_C10), tim on Bm_C4 (Sf_C31), CK2α on Bm_C5 (Sf_C09), cry2 on Bm_C15 (Sf_C23), CK2β and cry1 on Bm_C15 (Sf_C23), dbt on Bm_C17 (Sf_C17), shaggy on Bm_C18 (Sf_C13), clockwork orange on Bm_C22 (Sf_C6), slimb on Bm_C24 (Sf_C12, 32), vri on Bm_C27 (Sf_C25) and CK1α on Bm_scaf256 (which has not been mapped to a B. mori chromosome, and thus cannot be homologized). Thus, of all candidate genes, only vri mapped to the one major QTC Bm_C27 (Sf_C25). Since the reciprocal F1 hybrids (CR and RC) did not differ in their onset time of mating (Additional file 13; [16]), the involvement of the sex chromosome in the timing differentiation between the two strains can be excluded, which thus excludes per, clk, cyc and PdP1, which are located on the sex chromosome.

Expression differences of the main candidate gene vrille

The main candidate gene vrille showed a clear daily oscillation pattern in its expression in both strains and the overall expression of vrille differed between the strains (Fig. 3). We found significant differences in expression levels between the two strains in the first 24 h experiment (experiment 1, P = 0.00065), as well as in the second 12 h experiment (experiment 2, P <0.0001) (Fig. 3). These differences were due to overall significantly higher expression levels of vri in the rice-strain compared to the corn-strain.

Fig. 3
figure 3

Expression levels of vrille for the two strains of S. frugiperda (corn-strain in red and rice-strain in blue). The dots represent individual measurements; the line connects the means for each time point and each strain; the grey shaded area represents the confidence interval. Experiment 1 was conducted over 24 h, experiment 2 over 12 time points. Significant difference between the overall expression pattern of the corn-strain and the rice-strain is indicated by the asterisks

In comparison, the other clock genes showed rhythmic expression levels as well, except for clk and cyc which did not show a clear daily oscillation pattern (Fig. 4a-c). Three clock genes, per, tim and cry2 did not differ between the strains in either experiment (Fig. 4a1-6a3). While in the first 24 h experiment three additional clock genes did show significant differences between the two strains (pdp1 P = 0.00543, clk P = 0.00222, and cyc P = 0.00002; Fig. 4b, c), only pdp1 also showed significant differences in the second 12 h experiment (P = 0.0038). All differences between the two strains were due to significantly higher expression levels in the rice-strain than in the corn-strain (Figs. 3 and 4).

Fig. 4
figure 4

Expression levels of key clock genes for the two strains of S. frugiperda (corn-strain in red and rice-strain in blue). The dots represent individual measurements; the line connects the means for each time point and each strain; the grey shaded area represents the confidence interval. Experiment 1 was conducted over 24 h, experiment 2 over 12 time points. Significant difference between the overall expression pattern of the corn-strain and the rice-strain is indicated by the asterisks. A1-A3. Genes of the first, main feedback loop (period, timeless, cryptochrome 2); B1-B2: Core clock genes connecting both feedback loops (cycle, clock); C1-C2. Genes of the second, modulatory feedback loop (vrille, PAR-domain protein 1; Fig. 4 = C1 and was repeated to facilitate an overview of the key clock genes)

Sequence differences of the candidate gene vrille

Vrille is a short gene without introns in the protein coding region, coding for a 367 aa protein, followed by a 1234 bp 3’ UTR. The 5’ UTR is divided into a 45 bp segment and a 375 bp segment by an intron that contains regulatory elements, namely 11 Ebox elements (Ebox A – Ebox K) with the core sequence CACGTG (Fig. 3). Near Eboxes E, F and I, five strain-specific polymorphisms were identified in the investigated corn-strain and rice-strain populations and the maternal grandmothers of BC_A and BC_B (Table 1).

Table 1 Variation in the regulatory intron in the 5’UTR of vri

The availability of a corn-strain and a rice-strain variant of the S. frugiperda genome assembly enabled us to compare the sequences of the main clock genes in these assemblies (see Additional file 14). Sequence differences of the main clock genes were assessed by identifying non-synonymous as well as synonymous SNPs in the protein coding regions. We mostly found synonymous (syn) SNPs between the corn-strain and the rice-strain genome assembly, with the exception of clk, cyc and tim, which showed two, two and one non-synonymous (non-syn) SNPs, respectively (Table 2). All non-synonymous SNPs were located in non-conserved domains of these genes (Additional file 14).

Table 2 Number of SNPs between the corn-strain and rice-strain variant of the S. frugiperda genome assembly

Discussion

In this study, we aimed to identify the genetic basis of the main prezygotic isolation barrier between the two strains of S. frugiperda, i.e. allochronic differentiation. We found one consistent Quantitative trait chromosome (QTC) that significantly accounted for the difference in the onset time of mating in the two strains, Sf_C25, which is homologous to Bm_C27. Detecting only one major QTC is of note because the timing of behavior is a complex trait, depending on the complex network of the circadian clock and its interlocked feedback loops of transcription and translation (see Fig. 1). Thus, one might anticipate that the difference in mating time between the corn and rice strains would represent a polygenic trait affected by multiple loci of modest to small effect. Yet, differences in one or more genes on Sf_C25 explain 19% of the timing variance observed between the strains, i.e. most likely there are more loci affecting this phenotype that we did not detect in this study. QTC Sf_C25 is a single autosome and the homologous chromosome in B. mori (Bm_C27) is 14.5 Mb (52.8 cM). This size is comparable to QTL intervals found in other QTL studies [59–62]. In B. mori only one candidate clock gene is known to be located on this chromosome, namely vri. All other known clock genes map to different chromosomes in Lepidoptera.

The major QTC Sf_C25 explained 19% of the variance between the strains, while three minor QTCs in the combined analysis were not consistently present in both families when analyzed individually. In addition, no known genes involved in the circadian rhythm are located on the Bombyx homologs of these minor QTCs. Possibly, allochronic differentiation between the two strains is affected by an interaction between different factors involved in the circadian rhythm regulation.

A limitation of our indirect mapping approach is the different number of autosomes in B. mori (27) and S. frugiperda (30). When homologizing Helicoverpa armigera and B. mori, Sahara et al. [63] found Bm chromosomes 11, 23 and 24 to be merged from two chromosomes each. For Bm chromosome 23, we identified two homologous chromosomes in S. frugiperda: Sf_C12 and Sf_C29. Whether Bm_C11 and Bm_C24 are also represented by two chromosomes in S. frugiperda remains to be determined. The incomplete homology did not affect our result, because a) we confirmed the position of vrille on our major Sf QTC, b) all minor QTC have a confident homologue in B. mori (see Fig. 2, Additional file 6), none of which contain known clock genes, and c) we homologized all chromosomes with known clock genes. The high synteny level between B. mori and S. frugiperda [50] also supports our conclusion that vrille is the only clock gene located on a QTC in S. frugiperda.

Candidate gene vrille

Within the network of the circadian clock genes in insects, vri is a powerful player (e.g. in fire ants [64], pea aphids [65] and bean bugs [66]) and best described in Drosophila [36, 38, 67, 68]. VRI inhibits clk transcription, and since a dimer of CLK and CYC promotes many E-Box promoted genes, clk inhibition represses transcription of the core clock genes. Consequently, vri mutants have altered behavioral rhythms [68]. Hence, in S. frugiperda a strain-specific difference in vri structure or expression may cause a strain-specific expression difference in other clock genes, leading to a timing difference in behavior.

Our qPCR results show that vri is consistently higher expressed in the rice-strain compared to the corn-strain (Fig. 3). Since VRI inhibits the transcription of clk, which is needed for the expression of most circadian rhythm genes, differences in vri expression may impact the expression of the downstream players of the interlocked feedback loops. Interestingly, we only found strain-specific expression differences in the genes that constitute the second, modulatory feedback loop: vri, pdp1, clk and cyc with the rice-strain showing an overall higher expression of these genes in at least one of the experiments. Pdp1, clk and cyc are all located on the sex chromosome, which is not involved in the timing differentiation between the strains (see Additional file 13). Thus, the observed expression pattern differences reflect downstream temporal regulation dependent upon the initial upstream genetic difference in the expression of vrille. The fact that only one feedback loop shows expression differences between the strains is surprising, as the CLK:CYC dimer connects and promotes the transcription of the downstream genes in both feedback loops (Fig. 1). Yet, mRNA expression differences do not directly translate to differences at the protein level. For the generally low turnover of the CLK and CYC proteins, the total effect of mRNA expression differences on protein abundance can be subtle. Additionally, there is always competition for the CLK:CYC dimers between the E-Box promoter genes, and basic helix-loop-helix transcription factors like CLK and CYC exhibit differences in their binding specificities to E-Box binding sites [69]. Since per and tim showed higher expression peaks compared to pdp1 and vri, it is possible that per- and tim E-Boxes attract CLK:CYC dimers more easily and are thus less affected by a limited CLK:CYC dimer abundance.

Differences in expression levels of clock genes have been shown to be linked to different activity patterns. For example, the migratory songbird Emberiza melanocephala shows significant differences in clock gene expression levels, rather than phase shifts, between life history states that differ in their daily activity patterns (e.g. night activity in the migratory life state and day activity in the pre-migratory state) [70]. In insects, the cricket Gryllus bimaculatus shows lower expression of tim and per in the nocturnal nymphs compared to the diurnal adults [71]. Whether and how the differences in transcription levels of clock genes in S. frugiperda strains may affect behavioral variation in timing of sexual activities remains to be determined.

In our search for sequence differences in vri that might account for the timing difference, we found five strain-specific polymorphisms surrounding E-boxes between the corn-strain and a rice-strain population from Florida and in the parents of the backcross families used for the QTL analysis (Fig. 5). Since the binding specificity of basic helix-loop-helix transcription factors, such as CLK and CYC, is influenced by the genomic region surrounding the E-box binding site [69], a less efficient binding of a transcription factor to the active vri E-Box element (s) in e.g. the corn-strain could change the expression of vri. Alternatively, a cis-regulatory element regulating this gene could be situated on the same chromosome at a more distant region that we did not yet sequence. Mutations in cis-regulatory elements generally cause expression differences [72, 73] and are hypothesized to be key elements of evolutionary changes [74].

Fig. 5
figure 5

Structure of vri in the corn- and rice-strain of S. frugiperda and strain-specific polymorphisms in the intron in the 5’UTR

When aligning the protein coding regions of the other annotated clock genes in the corn-strain and rice-strain variant of the S. frugiperda genome, we identified a number of synonymous SNPs in every clock gene, and five non-synonymous SNP in clk, cyc and tim. However, none of these non-synonymous SNPs affected conserved protein domains of these genes. Also, these SNPs are based solely on the genome alignments and were not confirmed by testing multiple corn- and rice-strain individuals and may thus not be strain-specific. Tim is located on Sf_C31 (Bm_C04), which is not a QTC for the timing differences. Also, tim did not exhibit any strain-specific expression differences (Fig. 4a). Clk and cyc are situated on the sex-chromosome, which does not underlie the strain-specific timing differences, based on the observation of reciprocal backcrosses. The expression differences of these genes can be explained as a downstream-effect of differences in the expression of vri and pdp1.

Comparisons of clock gene expression patterns between species

To compare the expression patterns of the main clock genes in the night-active Lepidoptera S. frugiperda to expression data in other species is difficult, as expression studies addressing the clock genes are often conducted under Zeitgeber time (ZT) conditions (L:D 12:12), while we used a longer day (L:D 14:10). With this difference in mind, some expression levels can be compared.

In the day-active migratory Lepidoptera Danaus plexippus, and in the day-active Drosophila melanogaster, vri expression is highest early in the scotophase [38, 68]. In S. frugiperda, vri peaks earlier, namely 4 h before the scotophase in the rice-strain and 3 – 1 h before the scotophase in the corn-strain (24 h experiment). In S. frugiperda, pdp1 shows an expression peak 8 h after vri in the rice-strain and 4-7 h after vri in the corn-strain. In Drosophila expression levels of pdp1 peaked 6 h after vri [38] and thus more comparable to the rice-strain. These different times of peak expression levels are expected, as VRI and PDP1 proteins both bind to the promoter region of clock at different time points to inhibit or facilitate clock expression, respectively.

As for other clock genes, in both D. plexippus and D. melanogaster, tim expression peaks 2 h into scotophase (ZT 14) [75], while we found the highest expression levels of tim at 4 h into scotophase. In addition, in D. plexippus, per expression peaks 2 h into scotophase (ZT 14) [75], while we found the highest expression levels of per at 6-7 h (corn-strain) to 10 h into scotophase (rice-strain; 24 h experiment). Since TIM, PER and CRY2 proteins form a trimer to enter the nucleus, similar expression patterns are expected in these genes, so that the peak differences between per and tim in S. frugiperda is surprising.

Expression differences in clock genes are especially interesting in the context of diurnal and nocturnal activity patterns. Martin-Fairey et al. [76] investigated PER protein expression levels in diurnal grass rats Arvicanthis niloticus. When comparing individuals that exhibited the usual day-active behavior to individuals that had adopted a night-active behavior, they found PER protein expression levels comparable to nocturnal and diurnal rodents [76]: The expression of PER protein peaks in the early morning in diurnal grass rats [76, 77] while it peaks in the late night in nocturnal rodent species [78, 79]. To be able to make such comparisons in Lepidoptera, clock gene expression studies across species should be conducted under comparable conditions.

Conclusion

In summary, we identified one major QTC for the timing difference in mating between the two S. frugiperda strains. The clock gene vrille (vri) is located on this QTC and thus the major candidate for the strain-specific timing differences. Strain-specific expression differences, as well as strain-specific polymorphisms in the regulatory region of vri, support the hypothesis that vri plays a key role in the timing differentiation of these two strains. As allochronic, diurnal differentiation is likely the major isolation barrier driving divergence between S. frugiperda populations, the mechanism by which vri and the other circadian clock genes influence this differentiation should be elucidated, which is possible through fine-scale mapping and functional analyses. This will advance our understanding of the molecular basis of incipient speciation in sympatry.

References

  1. Groot AT: Circadian rhythms of sexual activities in moths: a review. Front Ecol Evol. 2014;2:43.

  2. Schulze T, Prager K, Dathe H, Kelm J, Kiessling P, Mittag M. How the green alga Chlamydomonas reinhardtii keeps time. Protoplasma. 2010;244(1-4):3–14.

    Article  CAS  PubMed  Google Scholar 

  3. Reppert SM, Weaver DR. Coordination of circadian timing in mammals. Nature. 2002;418(6901):935–41.

    Article  CAS  PubMed  Google Scholar 

  4. Hardin PE. Molecular genetic analysis of circadian timekeeping in Drosophila. Adv Genet. 2011;74:141–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Zhan S, Merlin C, Boore JL, Reppert SM. The monarch butterfly genome yields insights into long-distance migration. Cell. 2011;147(5):1171–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Weinstock GM, Robinson GE, Gibbs RA, Worley KC, Evans JD, Maleszka R, Robertson HM, Weaver DB, Beye M, Bork P, et al. Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006;443(7114):931–49.

    Article  CAS  Google Scholar 

  7. Cong Q, Borek D, Otwinowski Z, Grishin NV. Tiger swallowtail genome reveals mechanisms for speciation and caterpillar chemical defense. Cell Rep. 2015;10(6):910–9.

    Article  CAS  Google Scholar 

  8. Devries PJ, Austin GT, Martin NH. Diel activity and reproductive isolation in a diverse assemblage of Neotropical skippers (Lepidoptera : Hesperiidae). Biol J Linn Soc. 2008;94(4):723–36.

    Article  Google Scholar 

  9. Konno Y. Time-Lag between Sex-Pheromone Content and the Calling Behavior in the Yellow Peach Moth, Conogethes-Punctiferalis (Guenee) (Lepidoptera, Pyralidae). Appl Entomol Zool. 1986;21(4):622–4.

    CAS  Google Scholar 

  10. Miyatake T, Shimizu T. Genetic correlations between life-history and behavioral traits can cause reproductive isolation. Evolution. 1999;53(1):201–8.

    Article  Google Scholar 

  11. Miyatake T. Insect quality control: synchronized sex, mating system, and biological rhythm. Appl Entomol Zool. 2011;46(1):3–14.

    Article  Google Scholar 

  12. Monti L, Genermont J, Malosse C, LalanneCassou B. A genetic analysis of some components of reproductive isolation between two closely related species, Spodoptera latifascia (Walker) and S-descoinsi (Lalanne-Cassou and Silvain) (Lepidoptera: Noctuidae). J Evol Biol. 1997;10(1):121–34.

    Article  Google Scholar 

  13. Rund SS, Lee SJ, Bush BR, Duffield GE. Strain- and sex-specific differences in daily flight activity and the circadian clock of Anopheles gambiae mosquitoes. J Insect Physiol. 2012;58(12):1609–19.

    Article  CAS  PubMed  Google Scholar 

  14. Tauber E, Kyriacou CP, et al. Molecular evolution and population genetics of circadian clock genes. In: Circadian Rhythms, vol. 393. 2005. p. 797–817.

    Chapter  Google Scholar 

  15. Pashley DP, Hammond AM, Hardy TN. Reproductive isolating mechanisms in fall armyworm host strains (Lepidoptera, Noctuidae). Ann Entomol Soc Am. 1992;85(4):400–5.

    Article  Google Scholar 

  16. Schöfl G, Heckel DG, Groot AT. Time-shifted reproductive behaviours among fall armyworm (Noctuidae: Spodoptera frugiperda) host strains: evidence for differing modes of inheritance. J Evol Biol. 2009;22(7):1447–59.

    Article  PubMed  Google Scholar 

  17. Groot AT, Marr M, Heckel DG, Schöfl G. The roles and interactions of reproductive isolation mechanisms in fall armyworm (Lepidoptera: Noctuidae) host strains. Ecol Entomol. 2010;35:105–18.

    Article  Google Scholar 

  18. Prowell DP, McMichael M, Silvain JF. Multilocus genetic analysis of host use, introgression, and speciation in host strains of fall armyworm (Lepidoptera : Noctuidae). Ann Entomol Soc Am. 2004;97(5):1034–44.

    Article  CAS  Google Scholar 

  19. Pashley DP. Host-associated genetic differentiation in fall armyworm (Lepidoptera, Noctuidae) - a sibling species complex. Ann Entomol Soc Am. 1986;79(6):898–904.

    Article  Google Scholar 

  20. Meagher RL, Gallo-Meagher M. Identifying host strains of fall armyworm (Lepidoptera: Noctuidae) in Florida using mitochondrial markers. Fla Entomol. 2003;86(4):450–5.

    Article  Google Scholar 

  21. Machado V, Wunder M, Baldissera VD, Oliveira JV, Fiuza LM, Nagoshi RN. Molecular characterization of host strains of Spodoptera frugiperda (Lepidoptera : Noctuidae) in Southern Brazil. Ann Entomol Soc Am. 2008;101(3):619–26.

    Article  CAS  Google Scholar 

  22. Nagoshi RN, Silvie P, Meagher RL, Lopez J, Machados V. Identification and comparison of fall armyworm (Lepidoptera : Noctuidae) host strains in Brazil, Texas, and Florida. Ann Entomol Soc Am. 2007;100(3):394–402.

    Article  CAS  Google Scholar 

  23. Nagoshi RN, Meagher RL, Adamczyk JJ, Braman SK, Brandenburg RL, Nuessly G. New restriction fragment length polymorphisms in the cytochrome oxidase I gene facilitate host strain identification of fall armyworm (Lepidoptera: Noctuidae) populations in the southeastern United States. J Econ Entomol. 2006;99(3):671–7.

    Article  CAS  PubMed  Google Scholar 

  24. Groot AT, Marr M, Schöfl G, Lorenz S, Svatos A, Heckel DG: Host strain specific sex pheromone variation in Spodoptera frugiperda. Front Zool. 2008;5:21.

  25. Lima ER, McNeil JN. Female sex pheromones in the host races and hybrids of the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae). Chemoecology. 2009;19(1):29–36.

    Article  CAS  Google Scholar 

  26. Unbehend M, Hänniger S, Meagher RL, Heckel DG, Groot AT. Pheromonal divergence between two strains of Spodoptera frugiperda. J Chem Ecol. 2013;39(3):364–76.

    Article  CAS  PubMed  Google Scholar 

  27. Juárez ML, Murua MG, Garcia MG, Ontivero M, Vera MT, Vilardi JC, Groot AT, Castagnaro AP, Gastaminza G, Willink E. Host association of Spodoptera frugiperda (Lepidoptera: Noctuidae) corn and rice strains in Argentina, Brazil, and Paraguay. J Econ Entomol. 2012;105(2):573–82.

    Article  PubMed  Google Scholar 

  28. Juárez ML, Schöfl G, Vera M, Vilardi J, Murua M, Willink E, Hänniger S, Heckel DG, Groot AT: Population structure of Spodoptera frugiperda maize and rice host forms in South America: are they host strains? Entomol Exp Appl 2014.

  29. Groot AT, Unbehend M, Hänniger S, Juárez ML, Kost S, Heckel DG, et al. Evolution of reproductive isolation of Spodoptera frugiperda. In: Pheromone Communication in Moths: Evolution, Behavior and Application. Oakland, CA: UC Press; 2016.

    Google Scholar 

  30. Unbehend M, Hänniger S, Vásquez GM, Juárez ML, Reisig D, McNeil JN, Meagher RL, Jenkins DA, Heckel DG, Groot AT. 2013. Geographic variation in sexual attraction of S podoptera frugiperda corn- and rice-strain males to pheromone lures. Plos One. 9;e89255.

  31. Danley PD, Decarvalho TN, Fergus DJ, Shaw KL. Reproductive asynchrony and the divergence of Hawaiian crickets. Ethology. 2007;113(12):1125–32.

    Article  Google Scholar 

  32. Fergus DJ, Shaw KL. Circadian rhythms and period expression in the Hawaiian cricket genus Laupala. Behav Genet. 2013;43(3):241–53.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Prabhakaran PM, Sheeba V. Sympatric Drosophilid species melanogaster and ananassae differ in temporal patterns of activity. J Biol Rhythms. 2012;27(5):365–76.

    Article  PubMed  Google Scholar 

  34. Tauber E, Roe H, Costa R, Hennessy JM, Kyriacou CP. Temporal mating isolation driven by a behavioral gene in Drosophila. Curr Biol. 2003;13(2):140–5.

    Article  CAS  PubMed  Google Scholar 

  35. Kaiser TS, Poehn B, Szkiba D, Preussner M, Sedlazeck FJ, Zrim A, Neumann T, Nguyen LT, Betancourt AJ, Hummel T et al: The genomic basis of circadian and circalunar timing adaptations in a midge. Nature 2016, 540(7631):69- + .

  36. Hardin PE. The circadian timekeeping system of Drosophila. Curr Biol. 2005;15(17):R714–22.

    Article  CAS  PubMed  Google Scholar 

  37. Kaiser TS, Neumann D, Heckel DG: Timing the tides: Genetic control of diurnal and lunar emergence times is correlated in the marine midge Clunio marinus. BMC Genet. 2011;12:49.

  38. Cyran SA, Buchsbaum AM, Reddy KL, Lin MC, Glossop NRJ, Hardin PE, Young MW, Storti RV, Blau J. vrille, Pdp1, and dClock form a second feedback loop in the Drosophila circadian clock. Cell. 2003;112(3):329–41.

    Article  CAS  PubMed  Google Scholar 

  39. Zhu H, Gegear RJ, Casselman A, Kanginakudru S, Reppert SM. Defining behavioral and molecular differences between summer and migratory monarch butterflies. BMC Biol. 2009;7:14.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Zhu H, Yuan Q, Briscoe AD, Froy O, Casselman A, Reppert SM. The two CRYs of the butterfly. Curr Biol. 2005;15(23):R953–954.

    Article  CAS  PubMed  Google Scholar 

  41. Yuan Q, Metterville D, Briscoe AD, Reppert SM. Insect cryptochromes: gene duplication and loss define diverse ways to construct insect circadian clocks. Mol Biol Evol. 2007;24(4):948–55.

    Article  CAS  PubMed  Google Scholar 

  42. Trang LTD, Sehadova H, Ichihara N, Iwai S, Mita K, Takeda M. Casein kinases I of the silkworm, Bombyx mori: Their possible roles in circadian timing and developmental determination. J Biol Rhythms. 2006;21(5):335–49.

    Article  CAS  Google Scholar 

  43. Sandrelli F, Costa R, Kyriacou CP, Rosato E. Comparative analysis of circadian clock genes in insects. Insect Mol Biol. 2008;17(5):447–63.

    Article  CAS  PubMed  Google Scholar 

  44. Tomioka K, Matsumoto A. Circadian molecular clockworks in non-model insects. Curr Opin Insect Sci. 2015;7:58–64.

    Article  Google Scholar 

  45. Goto SG, Denlinger DL. Short-day and long-day expression patterns of genes involved in the flesh fly clock mechanism: period, timeless, cycle and cryptochrome. J Insect Physiol. 2002;48(8):803–16.

    Article  CAS  PubMed  Google Scholar 

  46. Reppert SM, Tsai T, Roca AL, Sauman I. Cloning of a Structural and Functional Homolog of the Circadian Clock Gene Period from the Giant Silkmoth Antheraea-Pernyi. Neuron. 1994;13(5):1167–76.

    Article  CAS  PubMed  Google Scholar 

  47. Froy O, Gotter AL, Casselman AL, Reppert SM. Illuminating the circadian clock in monarch butterfly migration. Science. 2003;300(5623):1303–5.

    Article  CAS  PubMed  Google Scholar 

  48. Iwai S, Fukui Y, Fujiwara Y, Takeda M. Structure and expressions of two circadian clock genes, period and timeless in the commercial silkmoth. Bombyx mori. J Insect Physiol. 2006;52(6):625–37.

    Article  CAS  PubMed  Google Scholar 

  49. Zhu H, Sauman I, Yuan Q, Casselman A, Emery-Le M, Emery P, Reppert SM. Cryptochromes define a novel circadian clock mechanism in monarch butterflies that may underlie sun compass navigation. PLoS Biol. 2008;6(1):e4.

    Article  PubMed  PubMed Central  Google Scholar 

  50. d'Alencon E, Sezutsu H, Legeai F, Permal E, Bernard-Samain S, Gimenez S, Gagneur C, Cousserans F, Shimomura M, Brun-Barale A, et al. Extensive synteny conservation of holocentric chromosomes in Lepidoptera despite high rates of local genome rearrangements. Proc Natl Acad Sci U S A. 2010;107(17):7680–5.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Groot AT, Estock ML, Horovitz JL, Hamilton J, Santangelo RG, Schal C, Gould F. QTL analysis of sex pheromone blend differences between two closely related moths: Insights into divergence in biosynthetic pathways. Insect Biochem Mol Biol. 2009;39(8):568–77.

    Article  CAS  PubMed  Google Scholar 

  52. Heckel DG. Comparative genetic linkage mapping in insects. Annu Rev Entomol. 1993;38:381–408.

    Article  Google Scholar 

  53. Baxter SW, Davey JW, Johnston JS, Shelton AM, Heckel DG, Jiggins CD, Blaxter ML. Linkage mapping and comparative genomics using next-generation RAD sequencing of a non-model organism. Plos One. 2011;6(4):11.

    Google Scholar 

  54. Groot AT, Staudacher H, Barthel A, Inglis O, Schöfl G, Santangelo RG, Gebauer-Jung S, Vogel H, Emerson J, Schal C, et al. One quantitative trait locus for intra- and interspecific variation in a sex pheromone. Mol Ecol. 2013;22(4):1065–80.

    Article  CAS  PubMed  Google Scholar 

  55. Gouin A, Bretaudeau A, Nam K, Gimenez S, Aury J-M, Duvic B, Hilliou F, Durand N, Montagné N, Darboux I et al: Extreme gene family expansion underlies adaptation for polyphagy in Spodoptera frugiperda. 2017; submitted.

  56. Shimomura M, Minami H, Suetsugu Y, Ohyanagi H, Satoh C, Antonio B, Nagamura Y, Kadono-Okuda K, Kajiwara H, Sezutsu H, et al. KAIKObase: An integrated silkworm genome database and data mining tool. BMC Genom. 2009;10:8.

    Article  Google Scholar 

  57. Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer; 2009.

    Book  Google Scholar 

  58. Xia QY, Wang J, Zhou ZY, Li RQ, Fan W, Cheng DJ, Cheng TC, Qin JJ, Duan J, Xu HF, et al. The genome of a lepidopteran model insect, the silkworm Bombyx mori. Insect Biochem Mol Biol. 2008;38(12):1036–45.

    Article  Google Scholar 

  59. Gleason JM, Ritchie MG. Do quantitative trait loci (QTL) for a courtship song difference between Drosophila simulans and D. sechellia coincide with candidate genes and intraspecific QTL? Genetics. 2004;166(3):1303–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Moehring AJ, Mackay TFC. The quantitative genetic basis of male mating behavior in Drosophila melanogaster. Genetics. 2004;167(3):1249–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Shaw KL, Parsons YM, Lesnick SC. QTL analysis of a rapidly evolving speciation phenotype in the Hawaiian cricket Laupala. Mol Ecol. 2007;16(14):2879–92.

    Article  CAS  PubMed  Google Scholar 

  62. Gleason JM, James RA, Wicker-Thomas C, Ritchie MG. Identification of quantitative trait loci function through analysis of multiple cuticular hydrocarbons differing between Drosophila simulans and Drosophila sechellia females. Heredity. 2009;103(5):416–24.

    Article  CAS  PubMed  Google Scholar 

  63. Sahara K, Yoshido A, Shibata F, Fujikawa-Kojima N, Okabe T, Tanaka-Okuyama M, Yasukochi Y. FISH identification of Helicoverpa armigera and Mamestra brassicae chromosomes by BAC and fosmid probes. Insect Biochem Mol Biol. 2013;43(8):644–53.

    Article  CAS  PubMed  Google Scholar 

  64. Ingram KK, Kutowoi A, Wurm Y, Shoemaker D, Meier R, Bloch G. The molecular clockwork of the fire ant Solenopsis invicta. Plos One. 2012;7(11):1–11.

    Article  Google Scholar 

  65. Cortes T, Ortiz-Rivas B, Martinez-Torres D. Identification and characterization of circadian clock genes in the pea aphid Acyrthosiphon pisum. Insect Mol Biol. 2010;19:123–39.

    Article  CAS  PubMed  Google Scholar 

  66. Ikeno T, Numata H, Goto SG. Molecular characterization of the circadian clock genes in the bean bug, Riptortus pedestris, and their expression patterns under long- and short-day conditions. Gene. 2008;419(1-2):56–61.

    Article  CAS  PubMed  Google Scholar 

  67. Glossop NRJ, Houl JH, Zheng H, Ng FS, Dudek SM, Hardin PE. VRILLE feeds back to control circadian transcription of Clock in the Drosophila circadian oscillator. Neuron. 2003;37(2):249–61.

    Article  CAS  PubMed  Google Scholar 

  68. Blau J, Young MW. Cycling vrille expression is required for a functional Drosophila clock. Cell. 1999;99(6):661–71.

    Article  CAS  PubMed  Google Scholar 

  69. Gordan R, Shen N, Dror I, Zhou T, Horton J, Rohs R, Bulyk ML. Genomic regions flanking E-Box binding sites influence DNA binding specificity of bHLH transcription factors through DNA shape. Cell Rep. 2013;3(4):1093–104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Singh D, Trivedi AK, Rani S, Panda S, Kumar V. Circadian timing in central and peripheral tissues in a migratory songbird: dependence on annual life-history states. FASEB J. 2015;29(10):4248–55.

    Article  CAS  PubMed  Google Scholar 

  71. Uryu O, Tomioka K. Post-Embryonic Development of Circadian Oscillations Within and Outside the Optic Lobe in the Cricket. Gryllus bimaculatus. Zool Sci. 2014;31(4):237–43.

    Article  CAS  PubMed  Google Scholar 

  72. Wittkopp PJ, Haerum BK, Clark AG. Independent effects of cis- and trans-regulatory variation on gene expression in Drosophila melanogaster. Genetics. 2008;178(3):1831–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Wittkopp PJ, Haerum BK, Clark AG. Regulatory changes underlying expression differences within and between Drosophila species. Nat Genet. 2008;40(3):346–50.

    Article  CAS  PubMed  Google Scholar 

  74. Wray GA. The evolutionary significance of cis-regulatory mutations. Nat Rev Genet. 2007;8(3):206–16.

    Article  CAS  PubMed  Google Scholar 

  75. Merlin C, Gegear RJ, Reppert SM. Antennal circadian clocks coordinate sun compass orientation in migratory monarch butterflies. Science. 2009;325(5948):1700–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Martin-Fairey CA, Ramanathan C, Stowie A, Walaszczyk E, Smale L, Nunez AA. Plastic oscillators and fixed rhythms: Changes in the phase of clock-gene rhythms in the PVN are not reflected in the phase of the melatonin rhythm of grass rats. Neuroscience. 2015;288:178–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Ramanathan C, Stowie A, Smale L, Nunez AA. Phase preference for the display of activity is associated with the phase of extra-suprachiasmatic nucleus oscillators within and between species. Neuroscience. 2010;170(3):758–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Angeles-Castellanos M, Mendoza J, Escobar C. Restricted feeding schedules phase shift daily rhythms of c-Fos and protein Per1 immunoreactivity in corticolimbic regions in rats. Neuroscience. 2007;144(1):344–55.

    Article  CAS  PubMed  Google Scholar 

  79. Feillet CA, Mendoza J, Albrecht U, Pevet P, Challet E. Forebrain oscillators ticking with different clock hands. Mol Cell Neurosci. 2008;37(2):209–21.

    Article  CAS  PubMed  Google Scholar 

  80. Hänniger S, Dumas P, Schöfl G, Gebauer-Jung S, Vogel H, Unbehend M, Heckel DG, Groot A. Data from: Genetic basis of allochronic differentiation in the fall armyworm. http://dx.doi.org/10.17617/3.k

Download references

Acknowledgements

We thank Domenica Schnabelrauch, Susanne Donnerhacke and Antje Schmaltz for their support in molecular analyses; Steffen Reifarth and Johannes Fleischmann for their support with insect observations, Simon Baxter for donating RAD P1 adapters and for help with RAD sequencing, Rob Meagher, Carlos Blanco and Laura Juarez for insects and DNA samples.

Funding

This research was funded by the Deutsche Forschungsgemeinschaft (P.S.GR362721), the National Science Foundation (awards IOS-1052238 and IOS-1456973), and the Max-Planck-Gesellschaft.

Availability of data and material

The Illumina RAD sequences are available from the Edmond respository under http://dx.doi.org/10.17617/3.k [80]. Sequence data is accessible via GenBank (KM675483 - KM675658).

Authors’ contributions

ATG, GS, SH, PD and DGH designed the experiments. SH, PD, ATG, GS conducted the experiments. SGJ, HV, DGH and SH conducted the bioinformatics analysis of the NGS data. PD, SH, DHG and ATG analyzed the data. SH, ATG, MU and PD drafted the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors have no financial or non-financial competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sabine Hänniger.

Additional files

Additional file 1:

Chronological outline of experiments. (PDF 86 kb)

Additional file 2:

Spodoptera frugiperda populations. (PDF 68 kb)

Additional file 3:

Generation of female-informative backcross families for QTL analysis. (PDF 61 kb)

Additional file 4:

Numbers of scored AFLP markers. (PDF 56 kb)

Additional file 5:

Coverage of RAD sequences. (PDF 224 kb)

Additional file 6:

Homologous chromosomes of S. frugiperda and B. mori. (PDF 20 kb)

Additional file 7:

Details of BLAST hits of Sf BC contigs to scaffolds of Sf genome and of BLAST hits of Sf genome scaffolds to Bm genome. (PDF 57 kb)

Additional file 8:

PCR conditions used in the mentioned experiments. (PDF 47 kb)

Additional file 9:

Primer combinations and annealing temperatures (Ta). (PDF 71 kb)

Additional file 10:

Mean onset time of mating for individuals that are homozygous or heterozygous for the QTC Sf_C25. (PDF 107 kb)

Additional file 11:

LOD scores for all linkage groups. (PDF 36 kb)

Additional file 12:

Power analysis for backcross families. (PDF 317 kb)

Additional file 13:

Mating time in S. frugiperda hybrids. (PDF 29 kb)

Additional file 14:

Position of SNPs between the corn-strain and rice-strain of S. frugiperda in the coding regions of the major clock genes. (PDF 76 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hänniger, S., Dumas, P., Schöfl, G. et al. Genetic basis of allochronic differentiation in the fall armyworm. BMC Evol Biol 17, 68 (2017). https://doi.org/10.1186/s12862-017-0911-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-017-0911-5

Keywords