Received 30 November 2015; accepted 26 January 2016; published 29 January 2016
Many beetles in the family Coccinellidae are identified by the number and pattern of spots appearing on the forewings, or elytra. Both the common and scientific names of many species describe the numbers of spots; for example, Adalia bipunctata (Linneaus) is commonly called the two-spotted lady beetle, Coccinella septempunctata Linnaeus is the seven-spotted lady beetle, and Coccinella novemnotata Herbst is the nine-spotted lady beetle. The species Coleomegilla maculata (DeGeer) (Coleoptera: Coccinellidae) is commonly called the pink lady beetle or the twelve-spotted lady beetle. It is a beneficial omnivore found in US agroecosystems and occurs throughout much of the North and South American continents. The name “twelve-spotted lady beetle” describes the typical spot pattern on C. maculata: dark spots on a lighter background that is pink to red or orange. While numerous color and pattern variations in other species of coccinellids have been described, and the inheritance of those patterns have been analyzed  , C. maculata is not a species that exhibits dramatic polymorphism in wild populations, and therefore it has not been used as a model for phenotype inheritance.
To facilitate genetics and biological control research, colonies of C. maculata were kept in continuous culture and inbred for over 64 generations since 2009. One result of inbreeding was the discovery and selection of novel phenotypes unique to laboratory strains. A stable homozygous strain of pale eyed beetles with pale yellow coloration in the cuticle, ye, was described earlier  . After further laboratory inbreeding, an additional phenotypic trait was identified and selected and was described here. This trait involved the pattern of spots on the elytra of the adult beetles. The pattern was an expanded dark region of the two anterior spots on each elytron such that the spots merged. This trait was labeled “ten spotted” and the strain abbreviated 10sp. The pattern was observed in both the wild type strain and in the ye strain. Phenotypically distinct individuals exhibiting both the ye and 10sp traits combined were selected, bred to stability, and crossed back to the parental twelve spotted phenotype to determine the heritability of the trait. To evaluate the effects of the trait on overall fitness of the insects, selected biometric measurements of the ye/10sp strain were compared to those of the wild type and the ye laboratory strains. Two sequences that were predicted as potential candidate genes contributing to the phenotype were identified, and compared with genes from model insects.
The ye/10sp strain was analyzed using classical Mendelian breeding and documented by digital image collection. Insects were maintained as previously described   . Individual insects were isolated for reciprocal straincrosses, and putative heterozygous first generation offspring (F1) crosses were pooled. Observed phenotype data were collected and entered into a spreadsheet, and phenotypic ratios were compared to expected Mendelian ratios for a single locus incomplete dominant allele using chi-squared distribution test (Excel®, Microsoft Corporation, Redmond, WA, USA). The expected phenotype of all offspring from a homozygous 10sp parent mated to an individual with wild type twelve spotted pattern was all 10sp, but with a heterozygous moderate expression pattern. Offspring from heterozygous parents were expected to result in 75% 10sp and 25% twelve spotted wild type phenotypes; offspring from a heterozygous parent mated to a wild type individual were expected to result in half 10sp and half wild type phenotypes. Images were collected using a Nikon digital camera, DMX 1200, with factory supplied ACT-1 software. The camera was mounted on a Nikon Stereomicroscope SMZ1500 (Nikon Corporation, Tokyo, Japan) with aperture fully closed to provide maximum depth of field. Objective lens was WD54 1x, oculars C-W10xA/22. Zoom was set at 0.75, 1, and 2 for various images. Shutter speeds were 1/50 to 1/75 sec. A Nikon NI-150 high intensity double gooseneck illuminator set at 75% or higher intensity illuminated the subjects from two opposing sides.
To estimate strain fecundity, the number of eggs in egg masses were counted as harvested from stable strain colony cages. A grouping of more than four eggs clustered together was considered a mass. To estimate fecundity and fertility on a finer scale, mature and apparently gravid individual females were isolated from each strain and eggs were collected and counted on a daily basis over ten consecutive days. Eggs were observed daily and the number of hatched neonates was counted. Pupae were weighed individually using a Sartorius CP2P-F analytical balance. Data were analyzed by one way analysis of variance using SigmaPlot, version 12software (SSI, San Jose, CA USA).
As research on coccinellids advances, high-throughput sequencing is expected to play a greater role in both gene expression study and defining molecular markers  . A high-throughput sequencing project was undertaken utilizing total RNA extracted from highly inbred (six isofemale selection steps) specimens of C. maculata, resulting in two adult transcriptomes that assembled into over 33,000 sequences each  . Two predicted complete protein coding sequences identified from the C. maculata transcriptomes  were compared to sequences in NCBI GenBank using BLAST alignment  . Characterized protein sequences from model insects were aligned with the predicted C. maculata sequences using DNASTARLasergene 8 MegAlign software (Madison, WI USA).
Inheritance of the 10sp trait is autosomal and exhibits incomplete dominance. Reciprocal crosses of insects with the trait and wild type or ye insects having the characteristic twelve spot phenotype produced offspring with expansion and merging of the anterior spots as heterozygotes, following expected Mendelian inheritance ratios (Table 1). All crosses resulted in ratios of offspring that did not differ from the expected ratios for an incomplete dominant autosomal allele (p > 0.05). The offspring of heterozygotes, when mated, produced offspring with a range of phenotypes including the normal spotting pattern, with twelve spots, and enlarged spots that were partially merged, and fully merged spots. Strain types are shown in Figure 1, with the homozygous characteristic pattern of the 10sp strain shown in Figure 1(B), and an example of a heterozygous ten spotted pattern shown in Figure 1(C). The differences between the homozygous and heterozygous patterns were subtle and for counts of the inherited trait, any variation of the trait was counted as ten spotted. Digital images of individual insects in homozygous and heterozygous form were compared by measuring the smallest region of the constriction between merged spots, the “waist”, and the widest region of the anterior spot, and using the measurements to estimate a ratio (Figure 2). The spot ratio of homozygous beetles was 0.886 ± 0.0347 (n = 13) while the ratio of heterozygous beetles was 0.542 ± 0.0819 (n = 8). These ratios were significantly different (student t-test, p < 0.001).
As represented in Figure 3, the strain ye/10sp was similar in fecundity and fertility to the laboratory wild type strain and the ye strain. Egg masses collected from the ovipositing colony cages did not differ significantly in size (F = 0.613, df = 2, p = 0.543). Because the egg masses collected from the colony cages might represent interrupted oviposition events, a second fecundity estimate was measured. Fecundity was measured by the number of eggs produced in ten days by individual gravid females from each of three strains. Individual gravid females produced quantities of eggs over a ten day period that were not significantly different (F = 3.052, df = 2, p = 0.067). However, wild type fertile egg masses had a higher hatch rate compared with the ye strain (F = 3.882, df = 2, p = 0.046), and the hatch rate was significant among the treatment groups (p = 0.039). However, the
Table 1. Segregation of adult pattern in the progeny of crosses between 10sp and 12 spotted parents.
Figure 1. The ten spotted (10sp) phenotype of Coleomegilla maculata. (A) Two specimens of the yellow eyes and elytra (ye) phenotype exhibiting the typical twelve spotted pattern (right) and the 10sp form (left) with merged spots on elytra. (B) A specimen exhibiting the homozygous form of the trait; the merged spots have very little constriction between them. (C) A specimen exhibiting the heterozygous trait form. The dark area between the spots is narrowed.
Figure 2. Diagram of measurement of the spot shape of individuals to identify incomplete dominance characteristics, with images of insect elytra flanking diagrammatic representations. Left: strong homozygous phenotype, with a wide waist region of melanization. Right, moderate heterozygous phenotype, with a narrow waist region of melanization.
difference between the ye/10sp and ye rate was not statistically significant. The mass of pupae from the three strains were compared and while the mean mass of the individuals from the ye strain were found to be significantly smaller than those of the wild type strain (F = 10.172, df = 2, p < 0.001), the difference was small (<10%) and not visually obvious.
Sequences involved in melanin biosynthesis (predicted N-b alanyldopamine synthase, ebony) and patterning (predicted transcription factor bric a brac) were identified from C. maculata transcriptomes  . Predicted translations of the sequences were very similar to those from the red flour beetle, Tribolium castaneum, and somewhat similar to genes from the fly Drosophila melanogaster, both genetically robust model organisms. Table 2 summarizes the similarities by alignment of the C. maculata sequences to genetics model insects. Translated sequences and ClustalW alignments are included as Supplementary Figures S1-S4.
The dark colors of insects are often based on pigments in the cuticle. One of the best known dark pigment
Figure 3. Bionomic measurements of colony strain fitness. Columns depict means with ±standard error bars. Means labeled with the same letter, or unlabeled columns are not significantly different. Upper left: mean number of eggs per egg mass collected from colonies of a stable homozygous strain of 10 spotted yellow phenotype of Coleomegilla maculata (ye/10sp), a stable homozygous strain of twelve (wild type) spotted yellow phenotype of C. maculata (ye), and a wild type strain (pink with twelve spots) of C. maculata (wt). Upper right, mean number of eggs per individual isolated gravid female from a stable homozygous strain of 10 spotted yellow phenotype of C. maculata (ye/10sp), a stable homozygous strain of twelve (wild type) spotted yellow phenotype of C. maculata (ye), and a wild type strain (pink with twelve spots) of C. maculata (wt). Lower left: mean ratio of hatched eggs to unhatched eggs from masses collected from the three strains: ye/10sp, ye, and wt. Lower right: mean mass of individual pupae (mg) collected from the three strains: ye/10sp, ye, and wt.
compounds is melanin, a derivative of the amino acid tyrosine. Melanization in insects is influenced by multiple genes at multiple loci. In Drosophila melanogaster pigment metabolism involves tyrosine hydroxylase, encoded by the gene pale; phenol oxidases and dopa decarboxylase; genes in the Yellow family; the genes black and tan; and N-b alanyldopamine synthase, encoded by the gene ebony. In Bicyclus anynana butterflies, different mutations implicating an entirely different enzyme, cysteine sulfinic acid decarboxylase, produce melanic mutants in the larva and adult forms of the insect  . The gene ebony has been described in the beetle Tribolium castaneum, and when it is disrupted the resulting phenotype is uniformly dark  . Melanic forms of insects, or “dark morphs” may be adapted to resist diseases. For example, a dark form of the wax moth, Galleria mellonella, had a thicker cuticle and was resistant to infection by the entomopathogen Beauveria bassiana, a fungus used commercially for biological control of insect pests  . Melanin in insects is also associated with temperature modulation, as demonstrated in the lady beetle Adalia bipunctata. Forms with dark elytra may be better adapted for climate tolerance, particularly in cold conditions, and an associated increase in activity may influences mate choice in some populations  . Recent research on the model organism Drosophila melanogaster indicates that dark pigmentation patterning is associated with resistance to ultraviolet radiation exposure, and melanic forms are most closely correlated with high levels of exposure  . Toxic defensive compounds and elytra coloration are correlated in Harmonia axiridis  and Coccinella septempunctata  ; whether this is also true for C. maculata is a topic for further research. Both the ye and ye/10sp insects produce reflex hemolymph leakage accompanied by a distinctively unpleasant odor when disturbed (empirical observation), but variation between the phenotypic strains has not yet been examined or analyzed.
Patterning in insects may be regulated by transcription factors such as optomotor-blind (omb) or bric a brac (bab)  . Color patterns in wings of butterflies are complex, and may be the product of the co-option of developmental pathways, as exemplified by the eye development gene optix which is correlated with wing patterns of Heliconius butterflies  . The genes involved in color patterning in beetles have not yet been discovered.
The ladybird beetle C. maculata has appealing characteristics for use as a genetic model organism. It is relatively easy to find and maintain in culture, it is visually appealing and unthreatening, and has a relatively rapid reproductive rate. For molecular genetics, it has a small genome, and a pair of transcriptomes of the adult life stage have been sequenced  . The species exhibits high genetic polymorphism across its wide geographic distribution  , and laboratory cultures kept for gene sequencing are more suitable for sample submission when inbred through multiple isofemale selections to increase homozygosity. Naturally occurring mutant strains of other insect species, and of course in domesticated vertebrates, have proven extremely valuable in genetic discovery. The ye and 10sp strains of C. maculata do not appear to arise in natural populations, but arose spontaneously during laboratory inbreeding. Recently, a visible pigmentation marker was constructed based on the ebony gene from the silk moth Bombyx mori for use in insect transgenesis  . The gene sequences presented herein could be used for similar biotechnology applications research. Coloration in lady beetles is indicative of a wide range physiological and ecological traits of interest. The combination of distinctive phenotypic laboratory strains and increasing genetic sequence availability provide valuable scientific resources for studying complex and fascinating interactive relationships of organisms with the environment, with each other, and in trophic relationships with other organisms.
Comparison of two sequences that are similar to the genes ebony and bric-a-brac from the C. maculata transcriptomes (Supplementary Figures S1-S4) with sequences from the model insects D. melanogaster, B. mori, and T. castaneum (Table 2) indicate that similar genes involved in melanization and patterning are present in C.
Table 2. Comparisons of predicted ebony and bric-a-brac sequences from transcriptomes to model insects in GenBank.
maculata. While these particular genes may not play a role in the phenotypes described here, they provide an initiation point for further studies on pigmentation biochemistry and pattern formation in a novel and important beneficial lady beetle. Predictably, the C. maculata sequences are more similar to T. castaneum sequences than to those of D. melanogaster or B. mori, because the genes of two beetles, albeit distantly related, should logically be more similar to each other than to those from other insect orders, either Diptera or Lepidoptera (respectively). As more genes from non-model organisms are annotated and curated, the conservation and divergence of gene evolution and function will be better understood.
The cost of producing genetic strains is inbreeding depression. The fitness assessments presented here, indicating relatively robust fecundity, fertility, and size, bode well for continuing availability of the strains described here, and other strains. Laboratory inbreeding of biological control agents may change predation characteristics  , therefore it is important to maintain and improve awareness and understanding of those traits that could lead to changes in the effectiveness of beneficial insects.
Additional studies of coccinellid pigments will require chemical or physical extraction and isolation methods and molecular biological methods such as gene disruption or over expression.
Further studies of C. maculata and related coccinellids and the genetic and biochemical processes responsible for elytra coloration will be facilitated by this unique strain of beetles. Molecular genetic markers will facilitate mark and recapture studies to evaluate immigration and emigration to and from natural and managed ecosystems  , and even from plant type or species may be possible. This strain is stable in the homozygous form. The strain will be useful for molecular genetic and biochemical studies of insect pigments, immunity, and evolution of gene regulatory networks.
An earlier version of this manuscript was reviewed by Jian Chen and Jonathan G. Lundgren, and anonymous reviewers. The comments and suggestions provided helpful improvement. The author thanks Jeff Gore and Brenda Yant for sharing field collected wild insects to found the original C. maculata colony from Stoneville, MS. The author thanks Mary Elizabeth Huddleston, Joseph Grey Ballenger, and Morgan Holmes for technical support and assistance with insect maintenance. The United States Government has the right to retain a non-exclusive, royalty-free license in and to any copyright of this article. This article reports the results of research only. Mention of a commercial or proprietary product does not constitute an endorsement of the product by the United States Department of Agriculture. USDA is an equal opportunity provider and employer.
The author claims no competing interests.
Margaret L. Allen designed and analyzed all experiments and data.
List of Abbreviations Used
10sp = Ten spotted pattern (strain)
wt = wild type
f = female
m = male
G = generation
R = reciprocal
ye = yellow eyes and elytra (strain)
F = filial
avg = average
Figure S1. Nucleotide sequence and translated amino acid sequence for predicted ebony gene transcript from transcriptomes of two individual adult C. maculata specimens. Illumina sequences were assembled individually and combined, and all assemblies resulted in the same predicted sequence. Lower case letters indicate untranslated sequence.
Figure S2. ClustalW alignment (Gonnet method) of predicted translated ebony sequence with predicted sequence from Tribolium castaneum genome project, and validated sequence from Drosophila melanogaster.
Figure S3. Nucleotide sequence and translated amino acid sequence for predicted bric-a-brac gene transcript from transcriptomes of two individual adult C. maculata specimens. Illumina sequences were assembled individually and combined, and all assemblies resulted in the same predicted sequence. Lower case letters indicate untranslated sequence.
Figure S4. ClustalW alignment (Gonnet method) of predicted translated bric-a-brac sequence with predicted sequence from Tribolium castaneum genome project, and validated sequence from Drosophila melanogaster.