Skip to main content
  • Research article
  • Open access
  • Published:

X-ray diffraction pattern from the flight muscle of Toxorhynchites towadensis reveals the specific phylogenic position of mosquito among Diptera

Abstract

Introduction

The Diptera are a group of insects with only a single pair of wings (forewings), and are considered monophyletic (originating from a common ancestor). The flight muscle in Diptera has features not observed in other insects, such as the long Pro-Ala-rich peptide associated with tropomyosin, not with troponin-I as in other insects, and the formation of a superlattice by myosin filaments analogous to that in vertebrate skeletal muscle.

Results

Here we describe X-ray diffraction patterns from the flight muscle of a mosquito, Toxorhynchites towadensis (Culicidae), belonging to a primitive group of Diptera. The diffraction pattern indicates that myosin filaments in the flight muscle of this species do not form a superlattice. X-ray diffraction also shows meridional reflections that are not observed in other dipterans, but are present in the patterns from bumblebee (Hymenoptera) flight muscle.

Conclusion

These observations suggest that the superlattice structure evolved after the common ancestor of Diptera had diverged from other insects. The flight muscle of mosquito may retain primitive structural features that are shared by Hymenoptera.

Introduction

Mosquitos belong to the insect order Diptera, which is characterized by having only one pair of wings. Diptera is one of the orders of insects that possess asynchronous flight muscles (asynchronous IFM), capable of driving their wings at high frequencies not achievable by the repetition of ordinary contraction-relaxation cycles. The reported frequencies of wing beat of mosquitos range between 320 and 800 Hz [17]. These previous studies have shown that the buzzing sound produced in this frequency range is also utilized for male-female communication/species recognition.

Although asynchronous IFM is believed to have occurred many times independently in the course of insect evolution, Diptera is considered to be monophyletic [8]. The asynchronous IFM of dipterans is known to have features not observed in other insects. One is the flight muscle-specific tropomyosin, which has a long proline-alanine-rich sequence at its C-terminus [911]. This sequence is observed in the flight muscles of all winged insects, but it is associated with troponin-I in other orders [9]. A second characteristic feature is the superlattice arrangement of the thick filaments [12]. In this arrangement, the position of the crown of myosin heads is axially shifted by one third of the 14.5-nm basic repeat with respect to its neighboring thick filaments (Fig. 1). This shift is not only inferred from X-ray diffraction, but is also confirmed by the analysis of electron micrographs (12). As a result, the thick filaments in a sarcomere form a lattice that has a larger unit cell size than that of a simple lattice in which all the thick filaments have the same axial position. Lattice structure is readily identified by recording X-ray diffraction patterns; in the case of the simple lattice, the intensity of the 1st myosin meridional reflection has a single peak on the meridian, whereas in the case of the superlattice, the intensity forms two peaks. This was first observed in Drosophila flight muscle [12], but was later confirmed in other species ([13] as well as in the present study).

Fig. 1
figure 1

Simple lattice and superlattice structures of myosin filaments in insect flight muscle. a, c simple lattice seen in insect orders other than Diptera; b, d superlattice seen in higher dipterans. a, b top view; c, d side view. Based on [12]

Mosquitos (Culicidae) are considered to be one of the most primitive groups of dipterans; it is thus of interest to know whether they also exhibit these features. However, their body sizes are generally small, and it is difficult to prepare flight muscle specimens suitable for recording high-quality X-ray diffraction patterns. Here, we describe the diffraction patterns from the IFM of the largest mosquito species in Japan, Toxorhynchites towadensis (elephant mosquito). The length of its dorsal longitudinal muscle fibers exceeds 2 mm, which makes it possible to use the mounting technique used for larger insects. The high-quality diffraction patterns recorded in this way revealed many features that are not observed in other (higher) dipteran species.

Materials and methods

Materials

Live mosquitoes (Toxorhynchites towadensis) were collected in a building close to the campus of the SPring-8 synchrotron facility in Harima, Japan (Fig. 2). This species was readily identified by its size, bluish metallic body color, and curved proboscis. Its IFM was glycerinated in a 50 % mixture of glycerol and a relaxing solution with protease inhibitors, and stored at –20 °C, as previously described [13]. On the day before the experiment, the dorsal longitudinal muscle fibers were mounted to the specimen chamber by clamping both ends with half-split gold meshes for electron microscopy [1315]. The composition of solutions (relaxing and rigor) was as described previously [16, 17]. For comparison, flight muscles from other dipterans and bumblebee (Bombus ignitus) were also used; these were prepared in the same manner. These dipterans include Phytomia zonata (hover fly), Promachus yesonicus (robber fly), Tipula coquilletti (crane fly) and Bibio sp. (march fly), and all these specimens were collected in or near the campus of SPring-8.

Fig. 2
figure 2

Picture of the giant mosquito (Toxorhynchites towadensis) that was used for the X-ray experiment. The size of the mesh is 5 mm

X-ray diffraction recording

X-ray diffraction patterns were recorded at the BL45XU beamline of SPring-8 [18]. The specimen-to-detector distance was 2.4 m, and the X-ray wavelength was 0.1 nm. As a detector, a cooled CCD camera (C4880-50, Hamamatsu Photonics, Hamamatsu, Japan; 1000 × 1018 pixels, pixel size, 150 × 150 μm) was used in combination with an image intensifier (VP5445-MOD, Hamamatsu Photonics). Some of the patterns from bumblebee muscle fibers were recorded on Imaging Plates (Fuji Film Co., Tokyo, Japan) with a pixel size of 50 × 50 μm. The diffraction patterns recorded under the same experimental condition were summed to improve the signal-to-noise ratio, and the background scattering was subtracted by the method described previously [15, 19].

Results

Diffraction patterns recorded from the mosquito IFM fibers are shown in Fig. 3. Like IFM diffraction patterns from other insects, the diffraction pattern from the mosquito IFM consists of a number of layer-line reflections (ladder-rung-like reflections originating from the helical arrangement of proteins along the filament axis) and meridional reflections (usually spot-like reflections that appear right on the meridian, originating from the axial repeat of monomers of myofilaments). Most of these reflections are indexable to the 38.7-nm long-pitch helical repeat of the actin filament with its associated regulatory proteins, and the 4-start helix of the myosin filament (116 nm basic repeat and 14.5-nm monomer pitch). These features have been described previously in detail for the IFM from the giant waterbug, Lethocerus [20]. A number of equatorial reflections, originating from the hexagonal arrangement of myofilaments within a sarcomere, are also observed. The intense 2,0 reflection and the much weaker 1,1 reflection indicate that the hexagonal lattice structure is such that a thin filament is positioned midway between the two neighboring thick filaments, like in other insects (for reasoning see ref. [21]). The unit cell size of the hexagonal lattice (the d-spacing for the 1,0 reflection) is ~45 nm, and the value is smaller in rigor (Fig. 3b) than in the relaxed state (Fig. 3a). Otherwise, the difference between the relaxed and rigor patterns is small, except for the slightly stronger actin-based layer line reflections in rigor.

Fig. 3
figure 3

Diffraction patterns recorded from the flight muscle of Toxorhynchites towadensis. a relaxed; b rigor. The low-angle area of the diffraction pattern in (a) is magnified in the box on the lower right. In the box, the innermost equatorial reflections (1,0, 1,1 and 2,0) and the two lowest-angle layer-line reflections (d-spacing, 58 and 38.7 nm) are indicated by green arrows. Note that the 1st and 2nd myosin meridional reflections (1MM and 2MM) have single peaks right on the meridian. The scale on the right edge of panel a indicates d-spacing in nm

In principle, layer line reflections that arise from the helical arrangement of objects have continuous intensity profiles. If the objects are arranged in a crystalline lattice, the layer lines have intensities only on the spots where the conditions of diffraction are met. This phenomenon is called lattice sampling.

The layer line reflections of the mosquito, especially those in the lower-angle region are finely sampled into discrete spots, indicating that the contractile proteins in the mosquito flight muscle are in crystalline order. This sampling is observed up to the myosin layer line at d = 7.2 nm, indicating that the 3-dimensional lattice structure of the mosquito IFM is very regular (in many other insects, the sampling is restricted to lower-order reflections). The pattern of the sampling is identical for all the layer-line reflections, and coincides with the equatorial reflections. This means that the contractile proteins in the mosquito IFM form a simple lattice (orientations of the proteins on the neighboring filaments are the same). This is in marked contrast to the myosin filament lattice of Drosophila IFM [12], which forms a superlattice with a greater lattice constant, similar to the superlattice of vertebrate skeletal muscle [22]. The normal insect lattice has a lattice constant (d1,0 spacing) of ~45 nm, and the superlattice has a lattice constant of ~78 nm \( \left(=45\times \sqrt{3}\right) \). The superlattice structure manifests itself most conspicuously in the myosin meridional reflection at d = 14.5 nm, whose intensity is split in the middle and has peaks at d ~ 78 nm when measured along the equator. In contrast, the myosin meridional reflection of the mosquito has a sharp peak right on the meridian (Fig. 3), as in non-dipteran insects, such as giant waterbug [20] and bumblebee [13, 15]. So far, all of the studied dipteran species belonging to suborder Brachycera, including Drosophila [12], Phytomia (hover fly, Fig. 4a), Tabanus [23] (horse fly), and Promachus (robber fly, Fig. 4b), and suborder Nematocera, including Ctenacroscelis [13], Tipula (these two are crane flies) and Bibio (march fly, Fig. 4d) are found to have the superlattice structure, and the mosquito as presented here is the first dipteran that does not have the superlattice structure.

Fig. 4
figure 4

Diffraction patterns from the flight muscles of other dipteran species. a Phytomia zonata (hover fly); b Promachus yesonicus (robber fly); c Tipula coquilletti (crane fly) and d Bibio sp. (march fly). Green arrowheads indicate peaks indexable to the thick-filament superlattice with a greater lattice constant. Blue arrowheads indicate that the intensity of the 1st myosin meridional reflection (at d = 14.5 nm) is weak on the meridian. The pattern in (a) was recorded in the relaxed state, and the others were recorded in the rigor state, in which the 1st myosin meridional reflection is weakened

The layer line reflections of the mosquito IFM, as shown in Fig. 3, are indexable to either the myosin helix with a basic repeat of 116 nm or the actin helix with a basic repeat of 38.7 nm (Fig. 5). Therefore, the mosquito myofilaments are believed to have the same helical symmetries as those of other insects. Usually, the layer line at d = 38.7 nm (combined actin 1st and myosin 3rd (based on the 116-nm basic repeat) is the strongest. The mosquito pattern is unusual in that the 2nd myosin layer line (at d = 58 nm) is very intense. This layer line is usually weak even if present, and is clearly recognized only in a few cases (e.g., in giant waterbug [24]).

Fig. 5
figure 5

Indexing of reflections in the diffraction pattern from the mosquito (left) and bumblebee (Bombus, right) flight muscles. The green lines indicate the positions of reflections expected from the myosin repeat, and the pink lines indicate those expected from the actin repeat. Note the resemblance of positions of meridional reflections that appear at angles lower than 14.5 nm−1

The meridional reflections of the mosquito IFM is very complex when compared with those of other insects (Fig. 5). In the giant waterbug and most of the dipterans (Fig. 4), the meridional reflections are ascribed to either the myosin axial repeat (14.5, 7.25, 4.83 nm, etc.) or the troponin axial repeat (12.9, 6.45, 4.3 nm, etc.). The actual axial troponin repeat is 38.7 nm, but because of its helical arrangement on the six thin filaments surrounding each thick filament, only reflections of the orders of multiples of three have intensities [20]. Therefore, the first meridional reflection that appears in the low-angle region is the 1st-order myosin reflection (8th of the 116-nm repeat) at d = 14.5 nm. However, many more reflections are observed on the meridian of the mosquito pattern. Many meridional reflections are observed in the region of d > 14.5 nm (Figs. 5 and 6), and they are not readily indexable to known structures. For example, the reflection at d = 16.7 nm may be indexed as the 7th of the 116-nm myosin repeat, but an intensity on the meridian is not expected from the 4-start helical symmetry of the myosin filament. The reflection at 24.4 nm is not indexable to any of the known repeats. This situation is similar to the bumblebee flight muscle, which also have many meridional reflections at d > 14.5 nm (Figs. 5 and 6). Originally, some of these lower-angle reflections in the bumblebee pattern were indexed to a 155 nm repeat (4 x actin’s 38.7-nm repeat), and were argued that they were of the thin-filament origin. However, higher-resolution recording using Imaging Plates (Fig. 6) reveal that the spacing of these reflections does not exactly coincide with that predicted from the actin repeat. The patterns of occurrence of meridional reflections are similar in mosquito and bumblebee, although they are not identical.

Fig. 6
figure 6

Comparison of intensity profiles of meridional reflections from the flight muscles from the mosquito and other insects. From the top, mosquito, bumblebee (Bombus) and hover fly (Phytomia). The green lines indicate the positions of reflections expected from the myosin repeat, and the pink lines indicate those expected from the actin repeat. Note that in both mosquito and bumblebee, there are a number of reflections that are not indexable to either myosin or actin repeat. The bumblebee data were recorded by using an Imaging Plate with a pixel size of 50 μm, while others were taken by using a cooled CCD camera with a pixel size of 150 μm. The intensity is normalized with respect to that of the 6th actin layer line at d = 5.9 nm

These observations suggest that the structure of the IFM of this primitive dipteran has a higher crystalline order, and is more complex than that of higher dipterans in that it contains more periodically arranged components.

Discussion

In this study, X-ray diffraction patterns from the flight muscle of a mosquito (Toxorhynchites towadensis), one of the most primitive dipterans, were described and compared with those of higher dipterans. It has been shown that the patterns from the mosquito are distinct from those from higher dipterans in that they have reflections that are not found or weak in other dipterans (especially in meridional reflections), and most importantly, they lack a myosin filament superlattice structure.

Dipterans have traditionally been classified into two major groups, the more primitive Nematocera (slender-bodied dipterans with long antennae) and Brachycera (stout flies with short antennae). A recent extensive study of the phylogenetic relationship of dipterans based on morphological and genomic traits [25] places a few groups of former Nematocera, including Culicomorpha (mosquitoes and allies), Psychodomorpha (sand flies and allies) and Tipulomorpha (crane flies and allies) as truly primitive groups of dipterans, while placing Bibionomorpha (march flies and allies) as a sister group of Brachycera. In support of this, the present study has shown that the diffraction patterns from a march fly, Bibio have traits similar to those of Brachycera species, i.e., the superlattice structure of the thick filaments and simple meridional reflections (Fig. 4). However, the crane flies (Ctenacroscelis and Tipula) are found to have these traits, suggesting a closer relationship between Tipulomorpha and Brachycera. To establish that the lack of the superlattice structure and the more complex set of meridional reflections are the primitive traits of dipteran IFM, diffraction patterns should be taken from other primitive families, such as Chilonomidae (midges) and Psychodidae (moth flies); however, these do not include species large enough to guarantee high-quality diffraction patterns. Finally, it is desirable to record diffraction patterns from dipterans belonging to the most primitive Deuterophlebiidae, located at the very base of the dipteran diversification. Again, there is little hope of obtaining diffraction patterns from these, because only one very tiny species (body length, ~2 mm) is known from Japan and it is endangered. In any event, it is clear from the present study that the superlattice structure of the thick filaments is not a trait present in the common ancestor, but it evolved after diversion of dipterans from other groups on insects.

The origin of the unusual meridional reflections from the mosquito IFM, not indexable to the repeats of either actin or myosin filaments, is currently unknown. In bumblebee IFM, similar meridional reflections are observed, and they are resistant to the digestion of the thin filaments by gelsolin (an F-actin-severing protein), suggesting that they are of the thick-filament origin (Iwamoto, unpublished). Therefore, at least in bumblebee, these may represent the repeat of unknown thick-filament proteins. Gelsolin treatment of mosquito IFM fibers is expected to clarify whether their unusual reflections are also of the thick-filament origin. A candidate for such an unknown protein is flightin, a myosin-binding protein found in Drosophila [26] and its periodicity is unknown. However, the unusual meridional reflections are not found in Drosophila and other Brachycera species, so that it may follow the basic repeat of myosin (116 nm) or the contrast of density along the filament may be so low so that they do not contribute to intensities of meridional reflections. It is also unclear whether the loss of the unknown thick-filament protein(s) and the acquisition of the superlattice structure are coupled.

Conclusions

The architecture of the myofilament lattice in the flight muscle of the mosquito, Toxorhynchites towadensis (Culicidae), retains features observed in other insect orders but not in higher dipterans (simple myosin-filament lattice and structures that give unusual meridional reflections). Examination of the IFM structures of other primitive dipterans is expected to give more insights about the early molecular evolution of dipterans.

Abbreviations

CCD:

Charge-coupled device

IFM:

Insect flight muscle

References

  1. Jackson JC, Robert D. Nonlinear auditory mechanism enhances female sounds for male mosquitoes. Proc Natl Acad Sci USA. 2006;103:16734–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Caprio MA, Huang JX, Faver MK, Moore A. Characterization of male and female wingbeat frequencies in the Anopheles quadrimaculatus complex in Mississippi. J Am Mosq Control Assoc. 2001;17:186–9.

    CAS  PubMed  Google Scholar 

  3. Charlwood JD, Thompson R, Madsen H. Observations on the swarming and mating behaviour of Anopheles funestus from southern Mozambique. Malar J. 2003;2:2.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Mankin RW. Acoustical detection of Aedes taeniorhynchus swarms and emergence exoduses in remote salt marshes. J Am Mosq Control Assoc. 1994;10:302–8.

    CAS  PubMed  Google Scholar 

  5. Gibson G, Warren B, Russell IJ. Humming in tune: sex and species recognition by mosquitoes on the wing. J Assoc Res Otolaryngol. 2010;11:527–40.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Pennetier C, Warren B, Dabire KR, Russell IJ, Gibson G. Singing on the wing as a mechanism for species recognition in the malarial mosquito Anopheles gambiae. Curr Biol. 2010;20:131–6.

    Article  CAS  PubMed  Google Scholar 

  7. Warren B, Gibson G, Russell IJ. Sex recognition through midflight mating duets in Culex mosquitoes is mediated by acoustic distortion. Curr Biol. 2009;19:485–91.

    Article  CAS  PubMed  Google Scholar 

  8. Beutel RG, Pohl H. Endopterygote systematics - where do we stand and what is the goal (Hexapoda, Arthropoda)? Syst Entom. 2006;31:202–19.

    Article  Google Scholar 

  9. Peckham M, Cripps R, White D, Bullard B. Mechanics and protein content of insect flight muscles. J Exp Biol. 1992;168:57–76.

    CAS  Google Scholar 

  10. Mogami K, Fujita SC, Hotta Y. Identification of Drosophila indirect flight muscle myofibrillar proteins by means of two-dimensional electrophoresis. J Biochem. 1982;91:643–50.

    CAS  PubMed  Google Scholar 

  11. Mateos J, Herranz R, Domingo A, Sparrow J, Marco R. The structural role of high molecular weight tropomyosins in dipteran indirect flight muscle and the effect of phosphorylation. J Muscle Res Cell Motil. 2006;27:189–201.

    Article  CAS  PubMed  Google Scholar 

  12. Squire JM, Bekyarova T, Farman G, Gore D, Rajkumar G, Knupp C, et al. The myosin filament superlattice in the flight muscles of flies: A-band lattice optimisation for stretch-activation? J Mol Biol. 2006;361:823–38.

    Article  CAS  PubMed  Google Scholar 

  13. Iwamoto H. Evidence for unique structural change of the thin filaments upon calcium-activation of insect flight muscle. J Mol Biol. 2009;390:99–111.

    Article  CAS  PubMed  Google Scholar 

  14. Iwamoto H, Oiwa K, Suzuki T, Fujisawa T. X-ray diffraction evidence for the lack of stereospecific protein interactions in highly activated actomyosin complex. J Mol Biol. 2001;305:863–74.

    Article  CAS  PubMed  Google Scholar 

  15. Iwamoto H, Inoue K, Yagi N. Fast X-ray recordings reveal dynamic action of contractile and regulatory proteins in stretch-activated insect flight muscle. Biophys J. 2010;99:184–92.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Iwamoto H. Strain sensitivity and turnover rate of low force cross-bridges in contracting skeletal muscle fibers in the presence of phosphate. Biophys J. 1995;68:243–50.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Iwamoto H. Influence of ionic strength on the actomyosin reaction steps in contracting skeletal muscle fibers. Biophys J. 2000;78:3138–49.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Fujisawa T, Inoue K, Oka T, Iwamoto H, Uruga T, Kumasaka T, et al. Small-angle X-ray scattering station at the SPring-8 RIKEN beamline. J Appl Cryst. 2000;33:797–800.

    Article  CAS  Google Scholar 

  19. Iwamoto H, Wakayama J, Fujisawa T, Yagi N. Static and dynamic X-ray diffraction recordings from living mammalian and amphibian skeletal muscles. Biophys J. 2003;85:2492–506.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Tregear RT, Edwards RJ, Irving TC, Poole KJV, Reedy MC, Schmitz H, et al. X-ray diffraction indicates that active cross-bridges bind to actin target zones in insect flight muscle. Biophys J. 1998;74:1439–51.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Iwamoto H. Flight muscle-specific Pro-Ala-rich extension of troponin is important for maintaining the insect-type myofilament lattice integrity. J Struct Biol. 2013;183:33–9.

    Article  CAS  PubMed  Google Scholar 

  22. Huxley HE, Brown W. The low-angle X-ray diagram of vertebrate striated muscle and its behaviour during contraction and rigor. J Mol Biol. 1967;30:383–434.

    Article  CAS  PubMed  Google Scholar 

  23. Iwamoto H, Trombitás K, Yagi N, Suggs JA, Bernstein SI. X-ray diffraction from flight muscle with a headless myosin mutation: implications for interpreting reflection patterns. Front Physiol Striated Muscle Physiol. 2014;5:416.

    Google Scholar 

  24. Bekyarova TI, Reedy MC, Baumann BAJ, Tregear RT, Ward A, Krzic U, et al. Reverse actin sliding triggers strong myosin binding that moves tropomyosin. Proc Natl Acad Sci USA. 2008;105:10372–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Wiegmann BM, Trautwein MD, Winkler IS, Barr NB, Kim J-W, Lambkin C, et al. Episodic radiations in the fly tree of life. Proc Natl Acad Sci USA. 2011;108:5690–5.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Vigoreaux JO, Saide JD, Valgeirsdottir K, Pardue ML. Flightin, a novel myofibrillar protein of Drosophila stretch-activated muscles. J Cell Biol. 1993;121:587–98.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The author expresses thanks to Dr. N. Yagi (JASRI) for his fruitful discussion, and Dr. T. Hikima (RIKEN) for his help at the beamline. The experiments were performed under approval of the SPring-8 Proposal Review Committee (Proposal Nos. 2009A1469, 2012A1272, 2012B1256, 2013B1333 and 2014A1298). Supported by Grant-in-Aid for Scientific Research, The Ministry of Education, Culture, Sports, Science and Technology, Japan, No. 23612009 and 26440185.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hiroyuki Iwamoto.

Additional information

Competing interests

The author declares that he has no competing interests.

Author’s contributions

HI designed the study and performed the experiments, analyses and the writing of the manuscript.

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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

Iwamoto, H. X-ray diffraction pattern from the flight muscle of Toxorhynchites towadensis reveals the specific phylogenic position of mosquito among Diptera. Zoological Lett 1, 24 (2015). https://doi.org/10.1186/s40851-015-0024-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40851-015-0024-1

Keywords