1
|
Nielsen ME, Nylin S, Wiklund C, Gotthard K. Evolution of butterfly seasonal plasticity driven by climate change varies across life stages. Ecol Lett 2023; 26:1548-1558. [PMID: 37366181 DOI: 10.1111/ele.14280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Revised: 05/17/2023] [Accepted: 05/23/2023] [Indexed: 06/28/2023]
Abstract
Photoperiod is a common cue for seasonal plasticity and phenology, but climate change can create cue-environment mismatches for organisms that rely on it. Evolution could potentially correct these mismatches, but phenology often depends on multiple plastic decisions made during different life stages and seasons that may evolve separately. For example, Pararge aegeria (Speckled wood butterfly) has photoperiod-cued seasonal life history plasticity in two different life stages: larval development time and pupal diapause. We tested for climate change-associated evolution of this plasticity by replicating common garden experiments conducted on two Swedish populations 30 years ago. We found evidence for evolutionary change in the contemporary larval reaction norm-although these changes differed between populations-but no evidence for evolution of the pupal reaction norm. This variation in evolution across life stages demonstrates the need to consider how climate change affects the whole life cycle to understand its impacts on phenology.
Collapse
Affiliation(s)
| | - Sören Nylin
- Zoology Department, Stockholm University, Stockholm, Sweden
| | | | - Karl Gotthard
- Zoology Department, Stockholm University, Stockholm, Sweden
| |
Collapse
|
2
|
Kimmitt AA, Pegan TM, Jones AW, Wacker KS, Brennan CL, Hudon J, Kirchman JJ, Ruegg K, Benz BW, Herman R, Winger BM. Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum. Proc Biol Sci 2023; 290:20221334. [PMID: 36695033 PMCID: PMC9874272 DOI: 10.1098/rspb.2022.1334] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Accepted: 12/19/2022] [Indexed: 01/26/2023] Open
Abstract
Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.
Collapse
Affiliation(s)
- Abigail A. Kimmitt
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI 48109, USA
| | - Teresa M. Pegan
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI 48109, USA
| | - Andrew W. Jones
- Department of Ornithology, Cleveland Museum of Natural History, Cleveland, OH 44106, USA
| | - Kristen S. Wacker
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI 48109, USA
| | - Courtney L. Brennan
- Department of Ornithology, Cleveland Museum of Natural History, Cleveland, OH 44106, USA
| | - Jocelyn Hudon
- Royal Alberta Museum, Edmonton, Alberta Canada, T5J 0G2
| | | | - Kristen Ruegg
- Biology Department, Colorado State University, Fort Collins, CO 80521, USA
| | - Brett W. Benz
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI 48109, USA
| | - Rachael Herman
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI 48109, USA
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY 11794, USA
| | - Benjamin M. Winger
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan, Ann Arbor, MI 48109, USA
| |
Collapse
|
3
|
Lindestad O, Nylin S, Wheat CW, Gotthard K. Local adaptation of life cycles in a butterfly is associated with variation in several circadian clock genes. Mol Ecol 2021; 31:1461-1475. [PMID: 34931388 DOI: 10.1111/mec.16331] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Revised: 11/30/2021] [Accepted: 12/07/2021] [Indexed: 12/25/2022]
Abstract
Many insects exhibit geographical variation in voltinism, the number of generations produced per year. This includes high-latitude species in previously glaciated areas, meaning that divergent selection on life cycle traits has taken place during or shortly after recent colonization. Here, we use a population genomics approach to compare a set of nine Scandinavian populations of the butterfly Pararge aegeria that differ in life cycle traits (diapause thresholds and voltinism) along both north-south and east-west clines. Using a de novo-assembled genome, we reconstruct colonization histories and demographic relationships. Based on the inferred population structure, we then scan the genome for candidate loci showing signs of divergent selection potentially associated with population differences in life cycle traits. The identified candidate genes include a number of components of the insect circadian clock (timeless, timeless2, period, cryptochrome and clockwork orange). Most notably, the gene timeless, which has previously been experimentally linked to life cycle regulation in P. aegeria, is here found to contain a novel 97-amino acid deletion unique to, and fixed in, a single population. These results add to a growing body of research framing circadian gene variation as a potential mechanism for generating local adaptation of life cycles.
Collapse
Affiliation(s)
- Olle Lindestad
- Department of Zoology, Stockholm University, Stockholm, Sweden
| | - Sören Nylin
- Department of Zoology, Stockholm University, Stockholm, Sweden
| | | | - Karl Gotthard
- Department of Zoology, Stockholm University, Stockholm, Sweden
| |
Collapse
|
4
|
Walton W, Stone GN, Lohse K. Discordant Pleistocene population size histories in a guild of hymenopteran parasitoids. Mol Ecol 2021; 30:4538-4550. [PMID: 34252238 DOI: 10.1111/mec.16074] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Revised: 06/23/2021] [Accepted: 07/06/2021] [Indexed: 01/03/2023]
Abstract
Signatures of past changes in population size have been detected in genome-wide variation in many species. However, the causes of such demographic changes and the extent to which they are shared across co-distributed species remain poorly understood. During Pleistocene glacial maxima, many temperate European species were confined to southern refugia. While vicariance and range expansion processes associated with glacial cycles have been widely documented, it is unclear whether refugial populations of co-distributed species have experienced shared histories of population size change. We analyse whole-genome sequence data to reconstruct and compare demographic histories during the Quaternary for Iberian refuge populations in a single ecological guild (seven species of chalcid parasitoid wasps associated with oak cynipid galls). For four of these species, we find support for large changes in effective population size (Ne ) through the Pleistocene that coincide with major climate events. However, there is little evidence that the timing, direction and magnitude of demographic change are shared across species, suggesting that demographic histories in this guild are largely idiosyncratic, even at the scale of a single glacial refugium.
Collapse
Affiliation(s)
- William Walton
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Graham N Stone
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
5
|
Sarabia C, vonHoldt B, Larrasoaña JC, Uríos V, Leonard JA. Pleistocene climate fluctuations drove demographic history of African golden wolves (Canis lupaster). Mol Ecol 2020; 30:6101-6120. [PMID: 33372365 DOI: 10.1111/mec.15784] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Revised: 12/03/2020] [Accepted: 12/14/2020] [Indexed: 12/31/2022]
Abstract
Pleistocene climate change impacted entire ecosystems throughout the world. In the northern hemisphere, the distribution of Arctic species expanded during glacial periods, while more temperate and mesic species contracted into climatic refugia, where isolation drove genetic divergence. Cycles of local cooling and warming in the Sahara region of northern Africa caused repeated contractions and expansions of savannah-like environments which connected mesic species isolated in refugia during interglacial times, possibly driving population expansions and contractions; divergence and geneflow in the associated fauna. Here, we use whole genome sequences of African golden wolves (Canis lupaster), a generalist mesopredator with a wide distribution in northern Africa to estimate their demographic history and past episodes of geneflow. We detect a correlation between divergence times and cycles of increased aridity-associated Pleistocene glacial cycles. A complex demographic history with responses to local climate change in different lineages was found, including a relict lineage north of the High Atlas Mountains of Morocco that has been isolated for more than 18,000 years, possibly a distinct ecotype.
Collapse
Affiliation(s)
- Carlos Sarabia
- Conservation and Evolutionary Genetics Group, Estación Biológica de Doñana (EBD-CSIC, Seville, Spain
| | - Bridgett vonHoldt
- Faculty of Ecology and Evolutionary Biology, University of Princeton, Princeton, NJ, USA
| | | | - Vicente Uríos
- Vertebrate Zoology Research Group, University of Alicante, Alicante, Spain
| | - Jennifer A Leonard
- Conservation and Evolutionary Genetics Group, Estación Biológica de Doñana (EBD-CSIC, Seville, Spain
| |
Collapse
|
6
|
Pruisscher P, Nylin S, Gotthard K, Wheat CW. Genetic variation underlying local adaptation of diapause induction along a cline in a butterfly. Mol Ecol 2018; 27:3613-3626. [PMID: 30105798 DOI: 10.1111/mec.14829] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2017] [Revised: 07/17/2018] [Accepted: 07/19/2018] [Indexed: 01/23/2023]
Abstract
Diapause is a life history strategy allowing individuals to arrest development until favourable conditions return, and it is commonly induced by shortened day length that is latitude specific for local populations. Although understanding the evolutionary dynamics of a threshold trait like diapause induction provides insights into the adaptive process and adaptive potential of populations, the genetic mechanism of variation in photoperiodic induction of diapause is not well understood. Here, we investigate genetic variation underlying latitudinal variation in diapause induction and the selection dynamics acting upon it. Using a genomewide scan for divergent regions between two populations of the butterfly Pararge aegeria that differ strongly in their induction thresholds, we identified and investigated the patterns of variation in those regions. We then tested the association of these regions with diapause induction using between-population crosses, finding significant SNP associations in four genes present in two chromosomal regions, one with the gene period, and the other with the genes kinesin, carnitine O-acetyltransferase and timeless. Patterns of allele frequencies in these two regions in population samples along a latitudinal cline suggest strong selection against heterozygotes at two genes within these loci (period, timeless). Evidence for additional loci modifying the diapause decision was found in patterns of allelic change in relation to induction thresholds over the cline, as well as in backcross analyses. Taken together, population-specific adaptations of diapause induction appear to be due to a combination of alleles of larger and smaller effect size, consistent with an exponential distribution of effect sizes involved in local adaption.
Collapse
Affiliation(s)
| | - Sören Nylin
- Department of Zoology, Stockholm University, Stockholm, Sweden
| | - Karl Gotthard
- Department of Zoology, Stockholm University, Stockholm, Sweden
| | | |
Collapse
|
7
|
Andersen JC, Havill NP, Caccone A, Elkinton JS. Postglacial recolonization shaped the genetic diversity of the winter moth ( Operophtera brumata) in Europe. Ecol Evol 2017; 7:3312-3323. [PMID: 28515868 PMCID: PMC5433974 DOI: 10.1002/ece3.2860] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Revised: 01/20/2017] [Accepted: 01/28/2017] [Indexed: 12/05/2022] Open
Abstract
Changes in climate conditions, particularly during the Quaternary climatic oscillations, have long been recognized to be important for shaping patterns of species diversity. For species residing in the western Palearctic, two commonly observed genetic patterns resulting from these cycles are as follows: (1) that the numbers and distributions of genetic lineages correspond with the use of geographically distinct glacial refugia and (2) that southern populations are generally more diverse than northern populations (the “southern richness, northern purity” paradigm). To determine whether these patterns hold true for the widespread pest species the winter moth (Operophtera brumata), we genotyped 699 individual winter moths collected from 15 Eurasian countries with 24 polymorphic microsatellite loci. We find strong evidence for the presence of two major genetic clusters that diverged ~18 to ~22 ka, with evidence that secondary contact (i.e., hybridization) resumed ~ 5 ka along a well‐established hybrid zone in Central Europe. This pattern supports the hypothesis that contemporary populations descend from populations that resided in distinct glacial refugia. However, unlike many previous studies of postglacial recolonization, we found no evidence for the “southern richness, northern purity” paradigm. We also find evidence for ongoing gene flow between populations in adjacent Eurasian countries, suggesting that long‐distance dispersal plays an important part in shaping winter moth genetic diversity. In addition, we find that this gene flow is predominantly in a west‐to‐east direction, suggesting that recently debated reports of cyclical outbreaks of winter moth spreading from east to west across Europe are not the result of dispersal.
Collapse
Affiliation(s)
- Jeremy C Andersen
- Department of Environmental Conservation University of Massachusetts Amherst Amherst MA USA.,Present address: Jeremy C. Andersen, Department of Environmental Science Policy and Management University of California Berkeley Berkeley CA USA
| | | | - Adalgisa Caccone
- Department of Ecology & Evolutionary Biology Yale University New Haven CT USA
| | - Joseph S Elkinton
- Department of Environmental Conservation University of Massachusetts Amherst Amherst MA USA
| |
Collapse
|
8
|
Longdon B, Day JP, Schulz N, Leftwich PT, de Jong MA, Breuker CJ, Gibbs M, Obbard DJ, Wilfert L, Smith SCL, McGonigle JE, Houslay TM, Wright LI, Livraghi L, Evans LC, Friend LA, Chapman T, Vontas J, Kambouraki N, Jiggins FM. Vertically transmitted rhabdoviruses are found across three insect families and have dynamic interactions with their hosts. Proc Biol Sci 2017; 284:20162381. [PMID: 28100819 PMCID: PMC5310039 DOI: 10.1098/rspb.2016.2381] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2016] [Accepted: 12/20/2016] [Indexed: 11/17/2022] Open
Abstract
A small number of free-living viruses have been found to be obligately vertically transmitted, but it remains uncertain how widespread vertically transmitted viruses are and how quickly they can spread through host populations. Recent metagenomic studies have found several insects to be infected with sigma viruses (Rhabdoviridae). Here, we report that sigma viruses that infect Mediterranean fruit flies (Ceratitis capitata), Drosophila immigrans, and speckled wood butterflies (Pararge aegeria) are all vertically transmitted. We find patterns of vertical transmission that are consistent with those seen in Drosophila sigma viruses, with high rates of maternal transmission, and lower rates of paternal transmission. This mode of transmission allows them to spread rapidly in populations, and using viral sequence data we found the viruses in D. immigrans and C. capitata had both recently swept through host populations. The viruses were common in nature, with mean prevalences of 12% in C. capitata, 38% in D. immigrans and 74% in P. aegeria We conclude that vertically transmitted rhabdoviruses may be widespread in a broad range of insect taxa, and that these viruses can have dynamic interactions with their hosts.
Collapse
Affiliation(s)
- Ben Longdon
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
| | - Jonathan P Day
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
| | - Nora Schulz
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany
| | - Philip T Leftwich
- School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich NR4 7TJ, UK
| | - Maaike A de Jong
- School of Biological Sciences, University of Bristol, Bristol Life Sciences Building, 24 Tyndall Avenue, Bristol BS8 1TQ, UK
| | - Casper J Breuker
- Evolutionary Developmental Biology Research Group, Department of Biological and Medical Sciences, Faculty of Health and Life Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford OX3 0BP, UK
| | - Melanie Gibbs
- NERC Centre for Ecology and Hydrology, Crowmarsh Gifford, Maclean Building, Wallingford, Oxfordshire OX10 8BB, UK
| | - Darren J Obbard
- Institute of Evolutionary Biology, University of Edinburgh, Ashworth Laboratories, Charlotte Auerbach Road, Edinburgh EH9 3FL, UK
| | - Lena Wilfert
- Centre for Ecology and Conservation, Biosciences, College of Life and Environmental Sciences, University of Exeter, Penryn Campus TR10 9FE, UK
| | - Sophia C L Smith
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
| | - John E McGonigle
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
| | - Thomas M Houslay
- Centre for Ecology and Conservation, Biosciences, College of Life and Environmental Sciences, University of Exeter, Penryn Campus TR10 9FE, UK
| | - Lucy I Wright
- Centre for Ecology and Conservation, Biosciences, College of Life and Environmental Sciences, University of Exeter, Penryn Campus TR10 9FE, UK
- Zoological Society of London, Regent's Park, London NW1 4RY, UK
| | - Luca Livraghi
- Evolutionary Developmental Biology Research Group, Department of Biological and Medical Sciences, Faculty of Health and Life Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford OX3 0BP, UK
| | - Luke C Evans
- Evolutionary Developmental Biology Research Group, Department of Biological and Medical Sciences, Faculty of Health and Life Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford OX3 0BP, UK
- Ecology and Evolutionary Biology Research Division, School of Biological Sciences, University of Reading, Whiteknights, Reading RG6 6AS, UK
| | - Lucy A Friend
- School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich NR4 7TJ, UK
| | - Tracey Chapman
- School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich NR4 7TJ, UK
| | - John Vontas
- Lab Pesticide Science, Agricultural University of Athens, Iera Odos 75, 11855, Athens, Greece
- Molecular Entomology, Institute Molecular Biology and Biotechnology/Foundation for Research and Technology, Voutes, 70013, Heraklio, Crete, Greece
| | - Natasa Kambouraki
- Lab Pesticide Science, Agricultural University of Athens, Iera Odos 75, 11855, Athens, Greece
- Molecular Entomology, Institute Molecular Biology and Biotechnology/Foundation for Research and Technology, Voutes, 70013, Heraklio, Crete, Greece
| | - Francis M Jiggins
- Department of Genetics, University of Cambridge, Cambridge CB2 3EH, UK
| |
Collapse
|
9
|
Drag L, Hauck D, Bérces S, Michalcewicz J, Šerić Jelaska L, Aurenhammer S, Cizek L. Genetic differentiation of populations of the threatened saproxylic beetle Rosalia longicorn,Rosalia alpina(Coleoptera: Cerambycidae) in Central and South-east Europe. Biol J Linn Soc Lond 2015. [DOI: 10.1111/bij.12624] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Lukas Drag
- Faculty of Science; University of South Bohemia; Branisovska 31 37005 Ceske Budejovice Czech Republic
- Biology Centre CAS, v. v. i.; Institute of Entomology; Branisovska 31 37005 Ceske Budejovice Czech Republic
| | - David Hauck
- Biology Centre CAS, v. v. i.; Institute of Entomology; Branisovska 31 37005 Ceske Budejovice Czech Republic
| | - Sándor Bérces
- Duna-Ipoly National Park Directorate; Kolto utca 21 1121 Budapest Hungary
| | - Jakub Michalcewicz
- Institute of Forest Ecosystem Protection; Faculty of Forestry; University of Agriculture; Al. 29 Listopada 46 31-425 Kraków Poland
| | - Lucija Šerić Jelaska
- Faculty of Science; University of Zagreb; Rooseveltov trg 6 10000 Zagreb Croatia
| | - Sandra Aurenhammer
- Institute for Animal Ecology and Landscape Planning; Bergmanngasse 22 8010 Graz Austria
| | - Lukas Cizek
- Faculty of Science; University of South Bohemia; Branisovska 31 37005 Ceske Budejovice Czech Republic
- Biology Centre CAS, v. v. i.; Institute of Entomology; Branisovska 31 37005 Ceske Budejovice Czech Republic
| |
Collapse
|
10
|
Aalberg Haugen IM, Gotthard K. Diapause induction and relaxed selection on alternative developmental pathways in a butterfly. J Anim Ecol 2014; 84:464-72. [PMID: 25267557 DOI: 10.1111/1365-2656.12291] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2014] [Accepted: 09/21/2014] [Indexed: 11/29/2022]
Abstract
Seasonal phenotypic plasticity entails differential trait expression depending on the time of season. The facultative induction of winter diapause in temperate insects is a developmental switch mechanism often leading to differential expression in life-history traits. However, when there is a latitudinal shift from a bivoltine to univoltine life cycle, selection for pathway-specific expression is disrupted, which may allow drift towards less optimal trait values within the non-selected pathway. We use field- and experimental data from five Swedish populations of Pararge aegeria to investigate latitudinal variation in voltinism, local adaptation in the diapause switch and footprints of selection on pathway-specific regulation of life-history traits and sexual dimorphism in larval development. Field data clearly illustrated how natural populations gradually shift from bivoltinism to univoltinism as latitude increases. This was supported experimentally as the decrease in direct development at higher latitudes was accompanied by increasing critical daylengths, suggesting local adaptation in the diapause switch. The differential expression among developmental pathways in development time and growth rate was significantly less pronounced in univoltine populations. Univoltine populations showed no significant signs of protandry during larval development, suggesting that erosion of the direct development pathway under relaxed selection has led to the loss of its sex-specific modifications.
Collapse
Affiliation(s)
| | - Karl Gotthard
- Department of Zoology, Stockholm University, Stockholm, SE-10691, Sweden
| |
Collapse
|