1
|
Poikela N, Laetsch DR, Hoikkala V, Lohse K, Kankare M. Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana. Genome Biol Evol 2024; 16:evae024. [PMID: 38482698 PMCID: PMC10972691 DOI: 10.1093/gbe/evae024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/30/2024] [Indexed: 04/01/2024] Open
Abstract
Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.
Collapse
Affiliation(s)
- Noora Poikela
- Department of Biological and Environmental Science, University of Jyväskylä, FI-40014, Jyväskylä, Finland
| | - Dominik R Laetsch
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Ville Hoikkala
- Department of Biological and Environmental Science, University of Jyväskylä, FI-40014, Jyväskylä, Finland
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Maaria Kankare
- Department of Biological and Environmental Science, University of Jyväskylä, FI-40014, Jyväskylä, Finland
| |
Collapse
|
2
|
Poikela N, Laetsch DR, Kankare M, Hoikkala A, Lohse K. Experimental introgression in Drosophila: Asymmetric postzygotic isolation associated with chromosomal inversions and an incompatibility locus on the X chromosome. Mol Ecol 2023; 32:854-866. [PMID: 36461113 PMCID: PMC10107139 DOI: 10.1111/mec.16803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 11/16/2022] [Accepted: 11/23/2022] [Indexed: 12/04/2022]
Abstract
Interspecific gene flow (introgression) is an important source of new genetic variation, but selection against it can reinforce reproductive barriers between interbreeding species. We used an experimental approach to trace the role of chromosomal inversions and incompatibility genes in preventing introgression between two partly sympatric Drosophila virilis group species, D. flavomontana and D. montana. We backcrossed F1 hybrid females from a cross between D. flavomontana female and D. montana male with the males of the parental species for two generations and sequenced pools of parental strains and their reciprocal second generation backcross (BC2 mon and BC2 fla) females. Contrasting the observed amount of introgression (mean hybrid index, HI) in BC2 female pools along the genome to simulations under different scenarios allowed us to identify chromosomal regions of restricted and increased introgression. We found no deviation from the HI expected under a neutral null model for any chromosome for the BC2 mon pool, suggesting no evidence for genetic incompatibilities in backcrosses towards D. montana. In contrast, the BC2 fla pool showed high variation in the observed HI between different chromosomes, and massive reduction of introgression on the X chromosome (large X-effect). This observation is compatible with reduced recombination combined with at least one dominant incompatibility locus residing within the X inversion(s). Overall, our study suggests that genetic incompatibilities arising within chromosomal inversions can play an important role in speciation.
Collapse
Affiliation(s)
- Noora Poikela
- Department of Biological and Environmental Science, University of Jyväskylä, Jyvaskyla, Finland
| | - Dominik R Laetsch
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Maaria Kankare
- Department of Biological and Environmental Science, University of Jyväskylä, Jyvaskyla, Finland
| | - Anneli Hoikkala
- Department of Biological and Environmental Science, University of Jyväskylä, Jyvaskyla, Finland
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
3
|
Xia T, Nishimura T, Nagata N, Kubota K, Sota T, Takami Y. Reproductive isolation via divergent genital morphology due to cascade reinforcement in Ohomopterus ground beetles. J Evol Biol 2023; 36:169-182. [PMID: 36357996 DOI: 10.1111/jeb.14116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Revised: 08/29/2022] [Accepted: 09/26/2022] [Indexed: 11/12/2022]
Abstract
Secondary contact between incipient species and selection against maladaptive hybridization can drive reinforcement between populations in contact and result in reproductive character displacement (RCD). Resultant divergence in mating traits within a species may generate downstream reproductive isolation between populations with displaced and non-displaced traits, referred to as the cascade reinforcement hypothesis. We examined this hypothesis using three allopatric populations of the ground beetle Carabus maiyasanus with a genital lock-and-key system. This species shows RCD in male and female genital morphologies in populations in contact with the sister species C. iwawakianus. In a reciprocal mating experiment using three allopatric populations with differences in male and female genital sizes, insemination failure increased as the difference in genital size increased. Based on the reproductive isolation index, insemination failure was the major postmating-prezygotic isolation barrier, at least in one population pair with comparable total isolation to those of other species pairs. By contrast, there was only incomplete premating isolation among populations. These results suggest that RCD in genital morphologies drives incipient allopatric speciation, supporting the cascade reinforcement hypothesis. These findings provide insight into the roles of interspecific interactions and subsequent trait diversification in speciation processes.
Collapse
Affiliation(s)
- Tian Xia
- Graduate School of Human Development and Environment, Kobe University, Nada, Kobe, Japan
| | - Taira Nishimura
- Graduate School of Human Development and Environment, Kobe University, Nada, Kobe, Japan
| | - Nobuaki Nagata
- Department of Zoology, Graduate School of Science, Kyoto University, Sakyo, Kyoto, Japan.,National Museum of Nature and Science, Tsukuba, Ibaraki, Japan
| | - Kohei Kubota
- Department of Forest Science, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo, Tokyo, Japan
| | - Teiji Sota
- Department of Zoology, Graduate School of Science, Kyoto University, Sakyo, Kyoto, Japan
| | - Yasuoki Takami
- Graduate School of Human Development and Environment, Kobe University, Nada, Kobe, Japan
| |
Collapse
|
4
|
Yusuf LH, Tyukmaeva V, Hoikkala A, Ritchie MG. Divergence and introgression among the virilis group of Drosophila. Evol Lett 2022; 6:537-551. [PMID: 36579165 PMCID: PMC9783487 DOI: 10.1002/evl3.301] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 09/23/2022] [Accepted: 10/12/2022] [Indexed: 12/03/2022] Open
Abstract
Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.
Collapse
Affiliation(s)
- Leeban H. Yusuf
- Centre for Biological Diversity, School of BiologyUniversity of St AndrewsSt AndrewsKY16 9THUnited Kingdom
| | - Venera Tyukmaeva
- Centre for Biological Diversity, School of BiologyUniversity of St AndrewsSt AndrewsKY16 9THUnited Kingdom,Department of Evolution, Ecology and BehaviourUniversity of LiverpoolLiverpoolL69 7ZBUnited Kingdom
| | - Anneli Hoikkala
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskylä40014Finland
| | - Michael G. Ritchie
- Centre for Biological Diversity, School of BiologyUniversity of St AndrewsSt AndrewsKY16 9THUnited Kingdom
| |
Collapse
|
5
|
Rushworth CA, Wardlaw AM, Ross-Ibarra J, Brandvain Y. Conflict over fertilization underlies the transient evolution of reinforcement. PLoS Biol 2022; 20:e3001814. [PMID: 36228022 PMCID: PMC9560609 DOI: 10.1371/journal.pbio.3001814] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Accepted: 09/02/2022] [Indexed: 11/16/2022] Open
Abstract
When two species meet in secondary contact, the production of low fitness hybrids may be prevented by the adaptive evolution of increased prezygotic isolation, a process known as reinforcement. Theoretical challenges to the evolution of reinforcement are generally cast as a coordination problem, i.e., "how can statistical associations between traits and preferences be maintained in the face of recombination?" However, the evolution of reinforcement also poses a potential conflict between mates. For example, the opportunity costs to hybridization may differ between the sexes or species. This is particularly likely for reinforcement based on postmating prezygotic (PMPZ) incompatibilities, as the ability to fertilize both conspecific and heterospecific eggs is beneficial to male gametes, but heterospecific mating may incur a cost for female gametes. We develop a population genetic model of interspecific conflict over reinforcement inspired by "gametophytic factors", which act as PMPZ barriers among Zea mays subspecies. We demonstrate that this conflict results in the transient evolution of reinforcement-after females adaptively evolve to reject gametes lacking a signal common in conspecific gametes, this gamete signal adaptively introgresses into the other population. Ultimately, the male gamete signal fixes in both species, and isolation returns to pre-reinforcement levels. We interpret geographic patterns of isolation among Z. mays subspecies considering these findings and suggest when and how this conflict can be resolved. Our results suggest that sexual conflict over fertilization may pose an understudied obstacle to the evolution of reinforcement.
Collapse
Affiliation(s)
- Catherine A. Rushworth
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, Minnesota, United States of America
- Department of Evolution and Ecology and Center for Population Biology, University of California, Davis, California, United States of America
- Department of Biology, Utah State University, Logan, Utah, United States of America
| | - Alison M. Wardlaw
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, Minnesota, United States of America
- Canada Revenue Agency—Agence du revenu du Canada, Ottawa, Ontario, Canada
| | - Jeffrey Ross-Ibarra
- Department of Evolution and Ecology and Center for Population Biology, University of California, Davis, California, United States of America
- Genome Center, University of California, Davis, California, United States of America
| | - Yaniv Brandvain
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, Minnesota, United States of America
- * E-mail:
| |
Collapse
|
6
|
Hoikkala A, Poikela N. Adaptation and ecological speciation in seasonally varying environments at high latitudes: Drosophila virilis group. Fly (Austin) 2022; 16:85-104. [PMID: 35060806 PMCID: PMC8786326 DOI: 10.1080/19336934.2021.2016327] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Living in high latitudes and altitudes sets specific requirements on species’ ability to forecast seasonal changes and to respond to them in an appropriate way. Adaptation into diverse environmental conditions can also lead to ecological speciation through habitat isolation or by inducing changes in traits that influence assortative mating. In this review, we explain how the unique time-measuring systems of Drosophila virilis group species have enabled the species to occupy high latitudes and how the traits involved in species reproduction and survival exhibit strong linkage with latitudinally varying photoperiodic and climatic conditions. We also describe variation in reproductive barriers between the populations of two species with overlapping distributions and show how local adaptation and the reinforcement of prezygotic barriers have created partial reproductive isolation between conspecific populations. Finally, we consider the role of species-specific chromosomal inversions and the X chromosome in the development of reproductive barriers between diverging lineages.
Collapse
Affiliation(s)
- Anneli Hoikkala
- Department of Biological and Environmental Sciences, University of Jyväskylä, Jyväskylä, Finland
| | - Noora Poikela
- Department of Biological and Environmental Sciences, University of Jyväskylä, Jyväskylä, Finland
| |
Collapse
|
7
|
Hickmann F, Cordeiro EG, Soares PL, Aurélio MSL, Schwertner CF, Corrêa AS. Reproductive Patterns Drive the Gene Flow and Spatial Dispersal of Euschistus heros (Hemiptera: Pentatomidae). JOURNAL OF ECONOMIC ENTOMOLOGY 2021; 114:2346-2354. [PMID: 34657956 DOI: 10.1093/jee/toab190] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Indexed: 06/13/2023]
Abstract
Euschistus heros (Hemiptera: Pentatomidae) has two allopatric strains with a hybrid zone in central Brazil. Asymmetric dispersal and gene flow between these strains of E. heros have been observed, where the South strain (SS) moves more quickly to the northern regions of the country than the North strain (NS) to the southern areas. In addition, SS generally has a bigger body size and presents dark brown coloration, and NS is usually smaller in size and presents light brown coloration. Here, we studied the reproductive behavior and tested for the presence of assortative mating and reproductive barriers between the two allopatric strains of E. heros. Nonrandom mating was observed in the SS strain based on mating choice trials and the reproductive isolation indexes. SS females and males prefer to mate with their co-specific (same strain) partner, while NS insects showed no mating preference. The insect's pronotum width was positively associated with the mating choice suggesting size-assortative mating in E. heros. Reciprocal crosses between strains yielded similar reproductive outputs when compared with pure strain crosses, suggesting similar fitness of hybrid pure strains. The asymmetric gene flow in the hybridization zone that favors SS seems to be associated with the reproductive behavior of the species, which favors the typical phenotype found in the SS populations.
Collapse
Affiliation(s)
- Frederico Hickmann
- Department of Entomology and Acarology, University of Sao Paulo, Luiz de Queiroz College of Agriculture (USP/ESALQ), Piracicaba, SP, 13418-900, Brazil
| | - Erick Goes Cordeiro
- Department of Entomology and Acarology, University of Sao Paulo, Luiz de Queiroz College of Agriculture (USP/ESALQ), Piracicaba, SP, 13418-900, Brazil
| | - Patrícia Lima Soares
- Department of Entomology and Acarology, University of Sao Paulo, Luiz de Queiroz College of Agriculture (USP/ESALQ), Piracicaba, SP, 13418-900, Brazil
| | - Mateus Souza L Aurélio
- Department of Entomology and Acarology, University of Sao Paulo, Luiz de Queiroz College of Agriculture (USP/ESALQ), Piracicaba, SP, 13418-900, Brazil
| | - Cristiano Feldens Schwertner
- Department of Ecology and Evolutionary Biology, Graduate Program in Ecology and Evolution, Federal University of São Paulo, Diadema, SP, 09972-270Brazil
| | - Alberto Soares Corrêa
- Department of Entomology and Acarology, University of Sao Paulo, Luiz de Queiroz College of Agriculture (USP/ESALQ), Piracicaba, SP, 13418-900, Brazil
| |
Collapse
|
8
|
Poikela N, Tyukmaeva V, Hoikkala A, Kankare M. Multiple paths to cold tolerance: the role of environmental cues, morphological traits and the circadian clock gene vrille. BMC Ecol Evol 2021; 21:117. [PMID: 34112109 PMCID: PMC8191109 DOI: 10.1186/s12862-021-01849-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Accepted: 06/03/2021] [Indexed: 11/10/2022] Open
Abstract
Background Tracing the association between insect cold tolerance and latitudinally and locally varying environmental conditions, as well as key morphological traits and molecular mechanisms, is essential for understanding the processes involved in adaptation. We explored these issues in two closely-related species, Drosophila montana and Drosophila flavomontana, originating from diverse climatic locations across several latitudes on the coastal and mountainous regions of North America. We also investigated the association between sequence variation in one of the key circadian clock genes, vrille, and cold tolerance in both species. Finally, we studied the impact of vrille on fly cold tolerance and cold acclimation ability by silencing it with RNA interference in D. montana. Results We performed a principal component analysis (PCA) on variables representing bioclimatic conditions on the study sites and used latitude as a proxy of photoperiod. PC1 separated the mountainous continental sites from the coastal ones based on temperature variability and precipitation, while PC2 arranged the sites based on summer and annual mean temperatures. Cold tolerance tests showed D. montana to be more cold-tolerant than D. flavomontana and chill coma resistance (CTmin) of this species showed an association with PC2. Chill coma recovery time (CCRT) of both species improved towards northern latitudes, and in D. flavomontana this trait was also associated with PC1. D. flavomontana flies were darkest in the coast and in the northern mountainous populations, but coloration showed no linkage with cold tolerance. Body size decreased towards cold environments in both species, but only within D. montana populations largest flies showed fastest recovery from cold. Finally, both the sequence analysis and RNAi study on vrille suggested this gene to play an essential role in D. montana cold resistance and acclimation, but not in recovery time. Conclusions Our study demonstrates the complexity of insect cold tolerance and emphasizes the need to trace its association with multiple environmental variables and morphological traits to identify potential agents of natural selection. It also shows that a circadian clock gene vrille is essential both for short- and long-term cold acclimation, potentially elucidating the connection between circadian clock system and cold tolerance. Supplementary Information The online version contains supplementary material available at 10.1186/s12862-021-01849-y.
Collapse
Affiliation(s)
- Noora Poikela
- Department of Biological and Environmental Science, University of Jyväskylä, P.O. Box 35, 40014, Jyväskylä, Finland.
| | - Venera Tyukmaeva
- Department of Biological and Environmental Science, University of Jyväskylä, P.O. Box 35, 40014, Jyväskylä, Finland.,Centre d'Ecologie Fonctionelle et Evolutive, CNRS, Montpellier, France
| | - Anneli Hoikkala
- Department of Biological and Environmental Science, University of Jyväskylä, P.O. Box 35, 40014, Jyväskylä, Finland
| | - Maaria Kankare
- Department of Biological and Environmental Science, University of Jyväskylä, P.O. Box 35, 40014, Jyväskylä, Finland
| |
Collapse
|
9
|
Yukilevich R. Reproductive Character Displacement Drives Diversification of Male Courtship Songs in Drosophila. Am Nat 2021; 197:690-707. [PMID: 33989143 DOI: 10.1086/714046] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
Abstract
AbstractMale secondary sexual traits are one of the most striking and diverse features of the animal kingdom. While these traits are often thought to evolve via sexual selection, many questions remain about their patterns of diversification and their role in speciation. To address these questions, I performed a comparative study of precopulatory male courtship songs of 119 Drosophila species across 10 distinct species groups. I related song divergence to genetic distances, geographic relationships, and sexual isolation between species. On the basis of pairwise Euclidean song distances, species groups typically retained their phylogenetic signal while species within groups diverged five times more in sympatry relative to allopatry, producing a pattern of reproductive character displacement. This occurred despite similar genetic distances in allopatry and sympatry, was exaggerated among younger species pairs, and was driven primarily by the parameter interpulse interval. While sexual isolation in sympatry was high even with low song divergence, these variables were correlated with each other and with increased divergence of female mating preferences in sympatry. The widespread pattern of character displacement implies that allopatric divergence due to processes like sexual selection are very slow relative to sympatric processes such as reinforcement and reproductive interference in driving song diversification across Drosophila.
Collapse
|
10
|
Larson EL, Brassil MM, Maslan J, Juárez D, Lilagan F, Tipton H, Schweitzer A, Skillman J, Monsen-Collar KJ, Peterson MA. The effects of heterospecific mating frequency on the strength of cryptic reproductive barriers. J Evol Biol 2019; 32:900-912. [PMID: 31162735 DOI: 10.1111/jeb.13495] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2019] [Revised: 05/22/2019] [Accepted: 05/28/2019] [Indexed: 01/02/2023]
Abstract
Heterospecific mating frequency is critical to hybrid zone dynamics and can directly impact the strength of reproductive barriers and patterns of introgression. The effectiveness of post-mating prezygotic (PMPZ) reproductive barriers, which include reduced fecundity via heterospecific matings and conspecific sperm precedence, may depend on the number, identity and order of mates. Studies of PMPZ barriers suggest that they may be important in many systems, but whether these barriers are effective at realistic heterospecific mating frequencies has not been tested. Here, we evaluate the strength of cryptic reproductive isolation in two leaf beetles (Chrysochus auratus and C. cobaltinus) in the context of a range of heterospecific mating frequencies observed in natural populations. We found both species benefited from multiple matings, but the benefits were greater in C. cobaltinus and extended to heterospecific matings. We found that PMPZ barriers greatly limited hybrid production by C. auratus females with moderate heterospecific mating frequencies, but that their effectiveness diminished at higher heterospecific mating frequencies. In contrast, there was no evidence for PMPZ barriers in C. cobaltinus females at any heterospecific mating frequency. We show that integrating realistic estimates of cryptic isolation with information on relative abundance and heterospecific mating frequency in the field substantially improves our understanding of the strong directional bias in F1 production previously documented in the Chrysochus hybrid zone. Our results demonstrate that heterospecific mating frequency is critical to understanding the impact of cryptic post-copulatory barriers on hybrid zone structure and dynamics, and that future studies of such barriers should incorporate field-relevant heterospecific mating frequencies.
Collapse
Affiliation(s)
- Erica L Larson
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA.,Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Margaret M Brassil
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Jonathan Maslan
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Danielle Juárez
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Flordeliza Lilagan
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Hallie Tipton
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Andrew Schweitzer
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | - Joe Skillman
- Biology Department, Western Washington University, Bellingham, Washington, USA
| | | | - Merrill A Peterson
- Biology Department, Western Washington University, Bellingham, Washington, USA
| |
Collapse
|