Ecological divergence and speciation between lemur (Eulemur) sister species in Madagascar

Understanding ecological niche evolution over evolutionary timescales is crucial to elucidating the biogeographic history of organisms. Here, we used, for the first time, climate‐based ecological niche models (ENMs) to test hypotheses about ecological divergence and speciation processes between sister species pairs of lemurs (genus Eulemur) in Madagascar. We produced ENMs for eight species, all of which had significant validation support. Among the four sister species pairs, we found nonequivalent niches between sisters, varying degrees of niche overlap in ecological and geographic space, and support for multiple divergence processes. Specifically, three sister‐pair comparisons supported the null model that niches are no more divergent than the available background region. These findings are consistent with an allopatric speciation model, and for two sister pairs (E. collaris–E. cinereiceps and E. rufus–E. rufifrons), a riverine barrier has been previously proposed for driving allopatric speciation. However, for the fourth sister pair E. flavifrons–E. macaco, we found support for significant niche divergence, and consistent with their parapatric distribution on an ecotone and the lack of obvious geographic barriers, these findings most strongly support a parapatric model of speciation. These analyses thus suggest that various speciation processes have led to diversification among closely related Eulemur species.     

The importance (or lack thereof) of niche divergence to the maintenance of a northern species complex: the case of the long‐toed salamander (Ambystoma macrodactylum Baird)

The relative importance of ecological vs. nonecological factors for the origin and maintenance of species is an open question in evolutionary biology. Young lineages – such as the distinct genetic groups that make up the ranges of many northern species – represent an opportunity to study the importance of ecological divergence during the early stages of diversification. Yet, few studies have examined the extent of niche divergence between lineages in previously glaciated regions and the role of ecology in maintaining the contact zones between them. In this study, we used tests of niche overlap in combination with ecological niche models to explore the extent of niche divergence between lineages of the long‐toed salamander (Ambystoma macrodactylum Baird) species complex and to determine whether contact zones correspond to (divergent) niche limits. We found limited evidence for niche divergence between the different long‐toed salamander lineages, substantial overlap in the predicted distribution of suitable climatic space for all lineages and range limits that are independent of niche limits. These results raise questions as to the importance of ecological divergence to the development of this widespread species complex and highlight the potential for non‐ecological factors to play a more important role in the maintenance of northern taxa.     

Insights from ecological niche modeling on the taxonomic distinction and niche differentiation between the black‐spotted and red‐spotted tokay geckoes (Gekko gecko)

The black‐spotted tokay and the red‐spotted tokay are morphologically distinct and have largely allopatric distributions. The black‐spotted tokay is characterized by a small body size and dark skin with sundry spots, while the red‐spotted tokay has a relatively large body size and red spots. Based on morphological, karyotypic, genetic, and distribution differences, recent studies suggested their species status; however, their classifications remain controversial, and additional data such as ecological niches are necessary to establish firm hypotheses regarding their taxonomic status. We reconstructed their ecological niches models using climatic and geographic data. We then performed niche similarity tests (niche identity and background tests) and point‐based analyses to explore whether ecological differentiation has occurred, and whether such differences are sufficient to explain the maintenance of their separate segments of environmental ranges. We found that both niche models of the black‐ and the red‐spotted tokay had a good fit and a robust performance, as indicated by the high area under the curve (AUC) values (“black” = 0.982, SD = ± 0.002, “red” = 0.966 ± 0.02). Significant ecological differentiation across the entire geographic range was found, indicating that the involvement of ecological differentiation is important for species differentiation. Divergence along the environmental axes is highly associated with climatic conditions, with isothermality being important for the “black” form, while temperature seasonality, precipitation of warmest quarter, and annual temperature range together being important for the “red” form. These factors are likely important factors in niche differentiation between the two forms, which result in morphological replacement. Overall, beside morphological and genetic differentiation information, our results contribute to additional insights into taxonomic distinction and niche differentiation between the black‐ and the red‐spotted tokay.     

Genetic,ecological, and behavioural divergence between two sibling snapping shrimp species (Crustacea: Decapoda: Alpheus)

Examination of genetic and ecological relationships within sibling species complexes can provide insights into species diversity and speciation processes. Alpheus angulatus and A. armillatus, two snapping shrimp species with overlapping ranges in the north-western Atlantic, are similar in morphology, exploit similar ecological niches and appear to represent recently diverged sibling species. We examined phylogenetic and ecological relationships between these two species with: (i) sequence data from two mitochondrial genes (16S rRNA and COI); (ii) data on potential differences in microhabitat distribution for A. armillatus and A. angulatus; and (iii) data from laboratory experiments on the level of reproductive isolation between the two species. DNA sequence data suggest A. armillatus and A. angulatus are sister species that diverged subsequent to the close of the Isthmus of Panama, and that haplotype diversity is lower in A. armillatus than in A. angulatus. Both species are distantly related to A. heterochaelis and A. estuariensis, two species with which A. angulatus shares some similarities in coloration. Ecological data on the distribution of A. angulatus and A. armillatus from two locations revealed differences in distribution of the two species between habitat patches, with each patch dominated by one or the other species. However, there was no apparent difference in distribution of the two species within habitat patches with respect to microhabitat location. Ecological data also revealed that heterospecific individuals often occur in close proximity (i.e. within metres or centimetres) where sympatric. Behavioural data indicated that these species are reproductively isolated, which is consistent with speciation in transient allopatry followed by post-divergence secondary contact. Our data further resolve taxonomic confusion between the sibling species, A. armillatus and A. angulatus, and suggest that sympatry in areas of range overlap and exploitation of similar ecological niches by these two recently diverged species have selected for high levels of behavioural incompatibility.     

Ecological divergence and habitat isolation between two migratory forms of Japanese threespine stickleback (Gasterosteus aculeatus)

When two closely related species migrate to divergent spawning sites, divergent use of spawning habitats can directly reduce heterospecific mating. Furthermore, adaptations to divergent spawning habitats can promote speciation as a by‐product of ecological divergence. Here, we investigated habitat isolation and ecological divergence between two anadromous forms of threespine stickleback (Gasterosteus aculeatus), the Japan Sea and Pacific Ocean forms. In several coastal regions of eastern Hokkaido, Japan, these forms migrate to the same watershed to spawn. Our field surveys in a single watershed revealed that segregation of distinct spawning sites between the two forms was maintained within the watershed across multiple years. These spawning sites diverged in salinity and predator composition. Morphological and physiological divergence between the forms also occurs in the direction predicted by ecological differences between the spawning sites. Our data indicate that migration into divergent spawning habitats can be an important mechanism contributing to speciation and phenotypic divergence in anadromous fishes.     

Host‐driven divergence in the parasitic plant Orobanche minor Sm. (Orobanchaceae)

Many parasitic angiosperms have a broad host range and are therefore considered to be host generalists. Orobanche minor is a nonphotosynthetic root parasite that attacks a range of hosts from taxonomically disparate families. In the present study, we show that O. minor sensu lato may comprise distinct, genetically divergent races isolated by the different ecologies of their hosts. Using a three‐pronged approach, we tested the hypothesis that intraspecific taxa O. minor var. minor and O. minor ssp. maritima parasitizing either clover (Trifolium pratense) or sea carrot (Daucus carota ssp. gummifer), respectively, are in allopatric isolation. Morphometric analysis revealed evidence of divergence but this was insufficient to define discrete, host‐specific taxa. Intersimple sequence repeat (ISSR) marker‐based data provided stronger evidence of divergence, suggesting that populations were isolated from gene flow. Phylogenetic analysis, using sequence‐characterized amplified region (SCAR) markers derived from ISSR loci, provided strong evidence for divergence by clearly differentiating sea carrot‐specific clades and mixed‐host clades. Low levels of intrapopulation SCAR marker sequence variation and floral morphology suggest that populations on different hosts are probably selfing and inbreeding. Morphologically cryptic Orobanche taxa may therefore be isolated from gene flow by host ecology. Together, these data suggest that host specificity may be an important driver of allopatric speciation in parasitic plants.     

Isotopic niche differs between seal and fish‐eating killer whales (Orcinus orca) in northern Norway

Ecological diversity has been reported for killer whales (Orcinus orca) throughout the North Atlantic but patterns of prey specialization have remained poorly understood. We quantify interindividual dietary variations in killer whales (n = 38) sampled throughout the year in 2017–2018 in northern Norway using stable isotopic nitrogen (δ¹⁵N: ¹⁵N/¹⁴N) and carbon (δ¹³C: ¹³C/¹²C) ratios. A Gaussian mixture model assigned sampled individuals to three differentiated clusters, characterized by disparate nonoverlapping isotopic niches, that were consistent with predatory field observations: seal‐eaters, herring‐eaters, and lumpfish‐eaters. Seal‐eaters showed higher δ¹⁵N values (mean ± SD: 12.6 ± 0.3‰, range = 12.3–13.2‰, n = 10) compared to herring‐eaters (mean ± SD: 11.7 ± 0.2‰, range = 11.4–11.9‰, n = 19) and lumpfish‐eaters (mean ± SD: 11.6 ± 0.2‰, range = 11.3–11.9, n = 9). Elevated δ¹⁵N values for seal‐eaters, regardless of sampling season, confirmed feeding at high trophic levels throughout the year. However, a wide isotopic niche and low measured δ¹⁵N values in the seal‐eaters, compared to that of whales that would eat solely seals (δ_N‐measured = 12.6 vs. δ_N‐expected = 15.5), indicated a diverse diet that includes both fish and mammal prey. A narrow niche for killer whales sampled at herring and lumpfish seasonal grounds supported seasonal prey specialization reflective of local peaks in prey abundance for the two fish‐eating groups. Our results, thus, show differences in prey specialization within this killer whale population in Norway and that the episodic observations of killer whales feeding on prey other than fish are a consistent behavior, as reflected in different isotopic niches between seal and fish‐eating individuals.     

Phylogenomics,biogeography, and evolution of morphology and ecological niche of the eastern Asian–eastern North American Nyssa (Nyssaceae)

Nyssa (Nyssaceae, Cornales) represents a classical example of the well‐known eastern Asian–eastern North American floristic disjunction. The genus consists of three species in eastern Asia, four species in eastern North America, and one species in Central America. Species of the genus are ecologically important trees in eastern North American and eastern Asian forests. The distribution of living species and a rich fossil record of the genus make it an excellent model for understanding the origin and evolution of the eastern Asian–eastern North American floristic disjunction. However, despite the small number of species, relationships within the genus have remained unclear and have not been elucidated using a molecular approach. Here, we integrate data from 48 nuclear genes, fossils, morphology, and ecological niche to resolve species relationships, elucidate its biogeographical history, and investigate the evolution of morphology and ecological niches, aiming at a better understanding of the well‐known EA–ENA floristic disjunction. Results showed that the Central American (CAM) Nyssa talamancana was sister to the remaining species, which were divided among three, rapidly diversified subclades. Estimated divergence times and biogeographical history suggested that Nyssa had an ancestral range in Eurasia and western North America in the late Paleocene. The rapid diversification occurred in the early Eocene, followed by multiple dispersals between and within the Erasian and North American continents. The genus experienced two major episodes of extinction in the early Oligocene and end of Neogene, respectively. The Central American N. talamancana represents a relic lineage of the boreotropical flora in the Paleocene/Eocene boundary that once diversified in western North America. The results supported the importance of both the North Atlantic land bridge and the Bering land bridge (BLB) for the Paleogene dispersals of Nyssa and the Neogene dispersals, respectively, as well as the role of Central America as refugia of the Paleogene flora. The total‐evidence‐based dated phylogeny suggested that the pattern of macroevolution of Nyssa coincided with paleoclimatic changes. We found a number of evolutionary changes in morphology (including wood anatomy and leaf traits) and ecological niches (precipitation and temperature) between the EA–ENA disjunct, supporting the ecological selection driving trait evolutions after geographic isolation. We also demonstrated challenges in phylogenomic studies of lineages with rapid diversification histories. The concatenation of gene data can lead to inference of strongly supported relationships incongruent with the species tree. However, conflicts in gene genealogies did not seem to impose a strong effect on divergence time dating in our case. Furthermore, we demonstrated that rapid diversification events may not be recovered in the divergence time dating analysis using BEAST if critical fossil constraints of the relevant nodes are not available. Our study provides an example of complex bidirectional exchanges of plants between Eurasia and North America in the Paleogene, but “out of Asia” migrations in the Neogene, to explain the present disjunct distribution of Nyssa in EA and ENA.     

Taxonomy and ecological niche modeling: Implications for the conservation of wood partridges (genus Dendrortyx)

Ecological niche models (ENMs) have a wide range of biological applications, particularly in conservation. To build these models, two sources of information are needed: occurrence records for the species of interest and environmental variables. However, taxonomic limits are often unclear, and the selection of occurrence data depends on the species concept being used. In this study we generated ENMs based on different taxonomic levels within the Dendrortyx group, which is comprised of three species and several subspecies; we analyzed the geographic and ecological distribution patterns and discuss the implications for the biogeography and conservation of this group. Our results suggest that the area with suitable climate depends on the taxonomic category used in the model, which in turn affects the interpretation of the importance of different biogeographic barriers and introduces variation into the potential differentiation of Dendrortyx. In terms of conservation, Dendrortyx macroura and Dendrortyx leucophrys are in a low risk category, that of “least concern,” although they may be amended to a higher category when their allopatric lineages are considered as the units for modeling. We suggest carrying out an a priori taxonomic analysis to facilitate the empirical identification of the units to be modeled in order to allow for a better ecological and biogeographic interpretation and more sound conservation policies.     

Ecological niche modelling of montane mammals in the Great Basin, North America: examining past and present connectivity of species across basins and ranges

Aim The goal of this study was to determine the extent of suitable habitats across the basins and ranges of the Great Basin for 13 montane mammals in the present and during the Last Glacial Maximum (LGM). For all these mammal species, we test whether: (1) more suitable habitat was available in basin areas during the LGM; (2) suitable habitat shifted upwards in elevation between the LGM and the present; (3) more ranges have suitable habitat than are currently occupied; and (4) these species are currently restricted to suitable habitats at higher‐elevation range areas. We also examine whether and how much distributional response varies among these montane mammal species. Location The Great Basin of western North America. Methods We re‐examine the past and present distributions of 13 Great Basin montane mammals using ecological niche modelling techniques that utilize now widely available species occurrence data and new, fine‐scale past climatological GIS layers in the present and at the LGM. These methods provide an objective, repeatable means for visual comparison of past and present modelled distributions for species examined in previous biogeographical studies. Results Our results indicate greater areal and lower elevational suitable habitat in the LGM than at present for nearly all montane mammals, and that there is more suitable habitat at present than is currently occupied. Our results also show that lowland areas provide suitable dispersal routes between ranges for most of the montane mammals both at the LGM and at present. However, three of the 13 species have little to no predicted suitable habitat in the LGM near currently occupied ranges, in contrast to the pattern for the other 10. For these species, the model results support more recent long‐distance colonization. Main conclusions Our finding of suitable lowland dispersal routes in the present for most species supports and greatly extends similar findings from single‐species studies. Our results also provide a visually striking confirmation that changes in species distribution and colonization histories of Great Basin montane mammals vary in a fashion related to the tolerances and requirements of each of these species; this has previously been hypothesized but not rigorously tested for multiple montane mammals in the region.     

Pliocene intraspecific divergence and Plio‐Pleistocene range expansions within Picea likiangensis (Lijiang spruce), a dominant forest tree of the Qinghai‐Tibet Plateau

A knowledge of intraspecific divergence and range dynamics of dominant forest trees in response to past geological and climate change is of major importance to an understanding of their recent evolution and demography. Such knowledge is informative of how forests were affected by environmental factors in the past and may provide pointers to their response to future environmental change. However, genetic signatures of such historical events are often weak at individual loci due to large effective population sizes and long generation times of forest trees. This problem can be overcome by analysing genetic variation across multiple loci. We used this approach to examine intraspecific divergence and past range dynamics in the conifer Picea likiangensis, a dominant tree of forests occurring in eastern and southern areas of the Qinghai‐Tibet Plateau (QTP). We sequenced 13 nuclear loci, two mitochondrial DNA regions and three plastid (chloroplast) DNA regions in 177 individuals sampled from 22 natural populations of this species, and tested the hypothesis that its evolutionary history was markedly affected by Pliocene QTP uplifts and Quaternary climatic oscillations. Consistent with the taxonomic delimitation of the three morphologically divergent varieties examined, all individuals clustered into three genetic groups with intervariety admixture detected in regions of geographical overlap. Divergence between varieties was estimated to have occurred within the Pliocene and ecological niche modelling based on 20 ecological variables suggested that niche differentiation was high. Furthermore, modelling of population‐genetic data indicated that two of the varieties (var. rubescens and var. linzhiensis) expanded their population sizes after the largest Quaternary glaciation in the QTP, while expansion of the third variety (var. likiangensis) began prior to this, probably following the Pliocene QTP uplift. These findings point to the importance of geological and climatic changes during the Pliocene and Pleistocene as causes of intraspecific diversification and range shifts of dominant tree species in the QTP biodiversity hot spot region.     

Phylogenomics,biogeography, and evolution of the blue‐ or white‐fruited dogwoods (Cornus)—Insights into morphological and ecological niche divergence following intercontinental geographic isolation

The eastern Asian (EA)–eastern North American (ENA) floristic disjunction represents a major pattern of phytogeography of the Northern Hemisphere. Despite 20 years of studies dedicated to identification of taxa that display this disjunct pattern, its origin and evolution remain an open question, especially regarding post‐isolation evolution. The blue‐ or white‐fruited dogwoods (BW) are the most species‐rich among the four major clades of Cornus L., consisting of ~35 species divided into three subgenera (subg. Yinquania, subg. Mesomora, and subg. Kraniopsis). The BW group provides an excellent example of the EA–ENA floristic disjunction for biogeographic study due to its diversity distribution centered in eastern Asia and eastern North America, yet its species relationships and delineation have remained poorly understood. In this study, we combined genome‐wide markers from RAD‐seq, morphology, fossils, and climate data to understand species relationships, biogeographic history, and ecological niche and morphological evolution. Our phylogenomic analyses with RAxML and MrBayes recovered a strongly supported and well‐resolved phylogeny of the BW group with three intercontinental disjunct clades in EA and ENA or Eurasia and North America, of which two are newly identified within subg. Kraniopsis. These analyses also recovered a potential new species but failed to resolve relationships within the C. hemsleyi–C. schindleri complex. In an effort to develop an approach to reduce computation time, analysis of different nodal age settings in treePL suggests setting a node's minimum age constraint to the lower bound of a fossil's age range to obtain similar ages to that of BEAST. Divergence time analyses with BEAST and treePL dated the BW stem back to the very Late Cretaceous and the divergence of the three subgenera in the Paleogene. By integrating fossil ages and morphology, a total evidence‐based dating approach was used in conjunction with time‐slice probabilities of dispersal under a DEC model to resolve ancestral ranges of each disjunct in the Miocene: Eurasia and ENA (disjunct 1), EA and western North America (disjunct 2), and EA (disjunct 3). The dated biogeographic history supports dispersal via the North Atlantic Land Bridge in the late Paleogene in disjunct 1 and dispersal via the Bering Land Bridge in the Miocene for disjuncts 2 and 3. Character mapping with a stochastic model in phytools and comparison of ecological niche, morphospace, and rate of evolution indicated differential divergence patterns in morphology, ecological niche, and molecules between disjunct sisters. Although morphological stasis was observed in most of the characters, evolutionary changes in growth habit and some features of leaf, flower, and fruit morphology occurred in one or both sister clades. A significant differentiation of ecological habitats in temperature, precipitation, and elevation between disjunct sisters was observed, suggesting a role of niche divergence in morphological evolution post‐isolation. The patterns of evolutionary rate between morphology and molecules varied among disjunct clades and were not always congruent between morphology and molecules, suggesting cases of non‐neutral morphological evolution driven by ecological selection. Our phylogenetic evidence and comparisons of evolutionary rate among disjunct lineages lend new insights into the formation of the diversity anomaly between EA and ENA, with particular support of an early diversification in EA. These findings, in conjunction with previous studies, again suggest that the EA–ENA disjunct floras are an assembly of lineages descended from the Mesophytic Forests that evolved from the early Paleogene “boreotropical flora” through varied evolutionary pathways across lineages.     

Environmentally induced and (epi‐)genetically based physiological trait differentiation between Heliosperma pusillum and its polytopically evolved ecologically divergent descendent,H. veselskyi (Caryophyllaceae: Sileneae)

Evolution is driven by natural selection, favouring individuals adapted in morphology and physiology to the environmental conditions at their growing site. Here, we studied environmentally induced and (epi‐)genetically based components of divergence of photosynthesis in response to irradiance and temperature and cellular characteristics reflecting water availability in two reciprocally non‐monophyletic, parapatric and interfertile, but morphologically and ecologically strongly divergent species in their natural habitat and in a common garden. Heliosperma pusillum occurs in alpine wet rocky habitats, whereas H. veselskyi is restricted to montane shady habitats with dry soils below overhanging rocks. Microclimatic divergence of natural habitats was evident in photosynthetically active photon flux density and leaf temperature. Photosynthetic light response was in agreement with differences in irradiance climate, suggesting that different photosynthetic parameters are favoured at particular growing sites. Photosynthesis was adapted to lower daily minimum temperatures at the alpine site, but did not differ at high temperatures. In the common garden, both species adjusted their functional traits in response to the shift in environmental conditions. The different light response of photosynthesis suggests an (epi‐)genetically based component affecting photosynthetic parameters. Early‐stage speciation between H. pusillum and H. veselskyi is probably environmentally induced, as the close connection between trait values and microclimatic conditions suggests that functional differentiation is adaptive. This is in line with the independent evolution of similar phenotypes under similar ecological conditions in different regions.     

Incomplete lineage sorting or secondary admixture: disentangling historical divergence from recent gene flow in the Vinous‐throated parrotbill (Paradoxornis webbianus)

Although Pleistocene glaciations had a major impact on the population genetic patterns of many species in North America and Europe, it remains unclear how these climatic fluctuations contributed to species diversification in East Asia. One reason for this is the difficulty of distinguishing genetic admixture following secondary contact from incomplete lineage sorting, both of which can generate similar patterns of genetic variation. Using a combination of multilocus analyses and coalescent simulation, we explore how these two processes occurred in the Pleistocene evolutionary history of a widespread East Asian bird, the Vinous‐throated parrotbill, Paradoxornis webbianus. Maximum likelihood (ML) tree identified two major mitochondrial lineages, which are geographically separated in most parts of its range, but are sympatric at a few sampling sites. NJ tree and Structure analysis of microsatellite data set revealed an extensive level of admixture and little population structure, suggesting recent admixture between two formerly separated groups. Networks from nuclear DNA data sets, however, did not indicate any geographically isolated groups but rather a panmictic population, thus support incomplete lineage sorting. By using coalescent simulation approaches, we show that both processes did occur, although at different temporal scales. During the Pleistocene glaciations, probably around 0.1–0.5 Ma (the Marine Isotope Stage 6, MIS6), P. webbianus contracted into two separate refugia, and subsequently accumulated genetic divergence. During the interglacial MIS5, the species expanded into previously glaciated areas allowing the once separated groups to come into contact and become admixed. Taken together, our results indicate the current genetic variation within P. webbianus is a combination pattern of widespread distribution in pre‐Pleistocene, then contraction and fragmentation into separated refugia during glacial advance, followed by recently postglacial expansion and admixture.     

r/K‐like trade‐off and voltinism discreteness: The implication to allochronic speciation in the fall webworm,Hyphantria cunea complex (Arctiidae)

North America has distinct types of Hyphantria moths (Arctiidae) characterized by red (RD)‐ and black (BL)‐headed larvae, of which the taxonomic status is unresolved. Genetic divergence of 26 populations, based on 710 bp of the mtCOI sequence, showed two phylogenetic lineages, which could not be connected in the haplotype network with 95% confidence. The two lineages are separated by 3.1% sequence divergence and should be considered for full species status. The estimated split occurred 1.2–1.6 million years ago. The range of the RD type covered most of the continent, whereas that of the BL type was limited to eastern deciduous forests. Several biological characteristics were differentiated in the zone of cohabitation where BL had more annual generations than RD. Spring emergence of BL precedes that of RD in the field by at least 1 month, because the diapause in BL was shallow, whereas it was deep in RD. Voltinism requires discreteness of numbers, which functions as a sink of hybrids between the two parental lines that have distinct but equally adaptive reproductive strategies; BL may be more r‐strategist‐like and RD more K‐strategist‐like, because fast‐developing BL has multivoltine life cycle, investing less silk proteins as the round‐the‐clock feeder, and slow‐developing RD univoltine one investing more silk as the nocturnal feeder. Also, intensity of diapause, deep in RD and weak in BL, was grossly different, which may enforce segregation of spring adults. Allochronic speciation avoiding coincidental occurrence of adult stages is therefore the most likely scenario. Because the adults never meet in nature, large morphological differentiation is not required.