Find research datasets worth reusing

Data from: The impact of rate heterogeneity on inference of phylogenetic models of trait evolution

Rates of trait evolution are known to vary across phylogenies; however, standard evolutionary models assume a homogeneous process of trait change. These simple methods are widely applied in small-scale phylogenetic studies, whereas models of rate heterogeneity are not, so the prevalence and patterns of potential rate variation in groups up to hundreds of species remain unclear. The extent to which trait evolution is modelled accurately on a given phylogeny is also largely unknown because studies typically lack absolute model fit tests. We investigated these issues by applying both rate-static and variable-rates methods on (i) body mass data for 88 avian clades of 10–318 species, and (ii) data simulated under a range of rate-heterogeneity scenarios. Our results show that rate heterogeneity is present across small-scaled avian clades, and consequently applying only standard single-process models prompts inaccurate inferences about the generating evolutionary process. Specifically, these approaches underestimate rate variation, and systematically mislabel temporal trends in trait evolution. Conversely, variable-rates approaches have superior relative fit (they are the best model) and absolute fit (they describe the data well). We show that rate changes such as single internal branch variations, rate decreases and early bursts are hard to detect, even by variable-rates models. We also use recently developed absolute adequacy tests to highlight misleading conclusions based on relative fit alone (e.g. a consistent preference for constrained evolution when isolated terminal branch rate increases are present). This work highlights the potential for robust inferences about trait evolution when fitting flexible models in conjunction with tests for absolute model fit.

Data from: Plant community evenness responds to spatial plant-soil feedback heterogeneity primarily through the diversity of soil conditioning

1.Plant-soil feedback (PSF) has been identified as a key driver of local plant diversity and evenness in competitive communities. However, while it has been shown that spatial PSF heterogeneity can alter plant performance and competitive interactions, there is no proof of principle that spatial PSF heterogeneity enhances community diversity. 2.Using a grassland model system we separated two aspects of spatial heterogeneity: the number of species conditioning the soil and spatial distribution of the PSFs. 3.Our data show that PSFs promoted a higher plant evenness when the soil was conditioned by multiple species (mixed-conditioned), then when the soil was conditioned by a single species (mono-conditioned). On mono-conditioned soils, heterospecifics typically outperformed the focal species. In addition, there was a trend for increasing community evenness from uniform, via fine-grained to coarse-grained mixed-conditioned soils, but this was not significant. 4.On mixed-conditioned soils, performance of all competing species was intermediate to the best and the worst mono-conditioned soils, leading to higher community evenness. 5.Our data demonstrate that PSFs play a role in promoting plant evenness. Across mono-conditioned soils, PSF led to altered competitive hierarchies. However, on soils conditioned by multiple species, competitive ability among species was more similar and this led to higher plant evenness. The spatial distribution of the heterogeneity, on the other hand, did not significantly affect plant evenness. Our data therefore show that community evenness was more strongly related to the number of plant species that conditioned the soil than the spatial distribution of the PSF heterogeneity. Future studies need to investigate the importance of PSFs in the field across plant life-stages and multiple generations.

opencc-zeroDec 2018View details →

Data from: Do genetic structure and landscape heterogeneity impact color morph frequency in a polymorphic salamander?

Landscape heterogeneity plays an important role in population structure and divergence, particularly for species with limited vagility. Here, we used a landscape genetic approach to identify how landscape and environmental variables affect genetic structure and color morph frequency in a polymorphic salamander. The Eastern Red-backed Salamander, Plethodon cinereus, is widely distributed in northeastern North America and contains two common color morphs, striped and unstriped, that are divergent in ecology, behavior, and physiology. To quantify population structure, rates of gene flow, and genetic drift, we amplified 10 microsatellite loci from 648 individuals across 28 sampling localities. This study was conducted in northern Ohio, where populations of P. cinereus exhibit an unusually wide range of morph frequency variation. To test whether genetic distance was more correlated with morph frequency, elevation, canopy cover, waterways, ecological niche, or geographic distance, we used resistance distance and least cost path analyses. We then examined whether landscape and environmental variables, genetic distance, or geographic distance were correlated with variation in morph frequency. Tests for population structure revealed three genetic clusters across our sampling range, with one cluster monomorphic for the striped morph. Rates of gene flow and genetic drift were low to moderate across sites. Genetic distance was most correlated with ecological niche, elevation, and a combination of landscape and environmental variables. In contrast, morph frequency variation was correlated with waterways and geographic distance. Thus, our results suggest that selection is also an important evolutionary force across our sites, and a balance between gene flow, genetic drift, and selection interact to maintain the two color morphs.

Data from: Linking resource availability and heterogeneity to understorey species diversity through succession in boreal forest of Canada

Understorey vegetation hosts the majority of species diversity and contributes greatly to ecosystem functioning in natural systems. Although patterns of understorey abundance, species diversity and composition associated with forest stand development are well researched, mechanisms driving these patterns remain largely speculative. We sampled fire origin stands of varying stand ages and overstorey compositions on mesic sites of the boreal forest of Canada and used structural equation modeling (SEM) to link time since fire (stand age), light availability and heterogeneity, substrate heterogeneity and soil nitrogen to understorey vegetation cover and species diversity. The most parsimonious model for total understorey cover showed a positive direct effect of stand age (r = 0.43) and an indirect effect via mean light level (0.18) and shrub cover (-0.11), with a positive total effect (0.50); the percent broadleaf canopy had a direct negative effect (-0.22) and an indirect effect via shrub cover (-0.11). The model for total understorey species richness showed an indirect effect of stand age via mean light (0.24), light heterogeneity (0.10), and substrate heterogeneity (0.07), with a positive total effect (0.52); percent broadleaf canopy had an indirect effect via light heterogeneity (0.09), and substrate heterogeneity (-0.10). Soil nitrogen did not significantly influence either understorey cover or species richness. The models for vascular plants followed similar trends to those for total understorey cover and species richness; however, there was an opposite indirect effect of light heterogeneity for both cover and species richness of non-vascular plants. Shrub cover had positive direct and negative direct and indirect effects on both vascular and non-vascular cover and species richness. Our findings indicate that understorey cover and species diversity are driven by time since disturbance, light availability as influenced by overstorey and shrub layers, but with important additional effects mediated by light and substrate heterogeneity. Non-vascular understorey vegetation is more strongly determined by time since disturbance than vascular vegetation, and negatively affected by broadleaf tree abundance. The overall results highlight the importance of colonization, light availability and heterogeneity, substrate specialization and growth dynamics in determining successional patterns of boreal forest understorey vegetation.

opencc-zeroJun 2019View details →

Data from: Aβ42 fibril formation from predominantly oligomeric samples suggests a link between oligomer heterogeneity and fibril polymorphism

Aβ oligomers play a central role in the pathogenesis of Alzheimer's disease. Oligomers of different sizes, morphology, and structures have been reported in both in vivo and in vitro studies, but there is a general lack of understanding about where to place these oligomers in the overall process of Aβ aggregation and fibrillization. Here we show that Aβ42 spontaneously forms oligomers with a wide range of sizes in the same sample. These Aβ42 samples contain predominantly oligomers, and they quickly form fibrils upon incubation at 37°C. When fractionated using ultrafiltration filters, the samples enriched with smaller oligomers form fibrils at a faster rate than the samples enriched with larger oligomers, with both a shorter lag time and faster fibril growth rate. This observation is independent of Aβ42 batches and HFIP treatment. Furthermore, the fibrils formed by the samples enriched with larger oligomers are more readily solubilized by EGCG, a main catechin component of green tea. These results suggest that the fibrils formed by larger oligomers may adopt a different structure from fibrils formed by smaller oligomers, pointing to a link between oligomer heterogeneity and fibril polymorphism.

Data from: Selective breeding and selection mapping using a novel wild-derived heterogeneous stock mice revealed two closely-linked loci for tameness

Tameness is a major behavioral factor for domestication, and can be divided into two potential components: motivation to approach humans (active tameness) and reluctance to avoid humans (passive tameness). We identified genetic loci for active tameness through selective breeding, selection mapping, and association analysis. In previous work using laboratory and wild mouse strains, we found that laboratory strains were predominantly selected for passive tameness but not active tameness during their domestication. To identify genetic regions associated with active tameness, we applied selective breeding over 9 generations for contacting, a behavioural parameter strongly associated with active tameness. The prerequisite for successful selective breeding is high genetic variation in the target population, so we established and used a novel resource, wild-derived heterogeneous stock (WHS) mice from eight wild strains. The mice had genetic variations not present in other outbred mouse populations. Selective breeding of the WHS mice increased the contacting level through the generations. Selection mapping was applied to the selected population using a simulation based on a non-selection model and inferred haplotype data derived from single-nucleotide polymorphisms. We found a genomic signature for selection on chromosome 11 containing two closely linked loci.

Data from: Bayesian inference of selection in a heterogeneous environment from genetic time-series data

Evolutionary geneticists have sought to characterize the causes and molecular targets of selection in natural populations for many years. Although this research program has been somewhat successful, most statistical methods employed were designed to detect consistent, weak to moderate selection. In contrast, phenotypic studies in nature show that selection varies in time and that individual bouts of selection can be strong. Measurements of the genomic consequences of such fluctuating selection could help test and refine hypotheses concerning the causes of ecological specialization and the maintenance of genetic variation in populations. Herein, I proposed a Bayesian non-homogenous hidden Markov model to estimate effective population sizes and quantify variable selection in heterogeneous environments from genetic time-series data. The model is described and then evaluated using a series of simulated data, including cases where selection occurs on a trait with a simple or polygenic molecular basis. The proposed method accurately distinguished neutral loci from non-neutral loci under strong selection, but not from those under weak selection. Selection coefficients were accurately estimated when selection was constant or when the fitness values of genotypes varied linearly with the environment, but these estimates were less accurate when fitness was polygenic or the relationship between the environment and the fitness of genotypes was non-linear. Past studies of temporal evolutionary dynamics in lab populations have been remarkably successful. The proposed method makes similar analyses of genetic time-series data from natural populations more feasible, and thereby could help answer fun damental questions about the causes and consequences of evolution in the wild.

Data from: Differential temperature sensitivity of synaptic and firing processes in a neural mass model of epileptic discharges explains heterogeneous response of experimental epilepsy to focal brain cooling

Experiments with drug-induced epilepsy in rat brains and epileptic human brain region reveal that focal cooling can suppress epileptic discharges without affecting the brain's normal neurological function. Findings suggest a viable treatment for intractable epilepsy cases via an implantable cooling device. However, precise mechanisms by which cooling suppresses epileptic discharges are still not clearly understood. Cooling experiments in vitro presented evidence of reduction in neurotransmitter release from presynaptic terminals and loss of dendritic spines at post-synaptic terminals offering a possible synaptic mechanism. We show that termination of epileptic discharges is possible by introducing a homogeneous temperature factor in a neural mass model which attenuates the post-synaptic impulse responses of the neuronal populations. This result however may be expected since such attenuation leads to reduced post-synaptic potential and when the effect on inhibitory interneurons is less than on excitatory interneurons, frequency of firing of pyramidal cells is consequently reduced. While this is observed in cooling experiments in vitro, experiments in vivo exhibit persistent discharges during cooling but suppressed in magnitude. This leads us to conjecture that reduction in the frequency of discharges may be compensated through intrinsic excitability mechanisms. Such compensatory mechanism is modelled using a reciprocal temperature factor in the firing response function in the neural mass model. We demonstrate that the complete model can reproduce attenuation of both magnitude and frequency of epileptic discharges during cooling. The compensatory mechanism suggests that cooling lowers the average and the variance of the distribution of threshold potential of firing across the population. Bifurcation study with respect to the temperature parameters of the model reveals how heterogeneous response of epileptic discharges to cooling (termination or suppression only) is exhibited. Possibility of differential temperature effects on post-synaptic potential generation of different populations is also explored.

opencc-zeroDec 2017View details →

Data from: Heterogeneity in efflux pump expression predisposes antibiotic-resistant cells to mutation

<p>Antibiotic resistance is often the result of mutations that block drug activity; however, bacteria also evade antibiotics by transiently expressing genes such as multidrug efflux pumps. A crucial question is whether transient resistance can promote permanent genetic changes. Previous studies have established that antibiotic treatment can select tolerant cells that then mutate to achieve permanent resistance. Whether these mutations result from antibiotic stress or preexist within the population is unclear. To address this question, we focused on the multidrug pump AcrAB-TolC. Using time-lapse microscopy, we found that cells with higher acrAB expression have lower expression of the DNA mismatch repair gene mutS, lower growth rates, and higher mutation frequencies. Thus, transient antibiotic resistance from elevated acrAB expression can promote spontaneous mutations within single cells.</p>

opencc-zeroJun 2019View details →

Data from: Bayesian tip dating reveals heterogeneous morphological clocks in Mesozoic birds

Recently, comprehensive morphological datasets including nearly all the well-recognized Mesozoic birds became available, making it feasible for statistically rigorous methods to unveil finer evolutionary patterns during early avian evolution. Here, we exploited the advantage of Bayesian tip dating under relaxed morphological clocks to estimate both the divergence times and evolutionary rates while accounting for their uncertainties. We further subdivided the characters into six body regions (i.e., skull, axial skeleton, pectoral girdle and sternum, forelimb, pelvic girdle, and hindlimb) to assess evolutionary rate heterogeneity both along the lineages and across partitions. We observed extremely high rates of morphological character changes during early avian evolution and the clock rates are quite heterogeneous among the six regions. The branch subtending Pygostylia shows extremely high rate in the axial skeleton, while the branches subtending Ornithothoraces and Enantiornithes show notably high rates in the pectoral girdle and sternum, and moderately high rates in the forelimb. The extensive modifications in these body regions largely correspond to refinement of the flight capability. This study reveals the power and flexibility of Bayesian tip dating implemented in MrBayes to investigate evolutionary dynamics in deep time.

opencc-zeroDec 2017View details →

Data from: Heterogeneity in local density allows a positive evolutionary relationship between self-fertilisation and dispersal

Despite empirical evidence for a positive relationship between dispersal and self-fertilisation (selfing), theoretical work predicts that these traits should always be negatively correlated, and the Good Coloniser Syndrome of high dispersal and selfing (Cf. Baker's Law) should not evolve. Critically, previous work assumes that adult density is spatiotemporally homogeneous, so selfing results in identical offspring production for all patches, eliminating the benefit of dispersal for escaping from local resource competition. We investigate the joint evolution of dispersal and selfing in a demographically structured metapopulation model where local density is spatiotemporally heterogeneous due to extinction-recolonisation dynamics. Selfing alleviates outcrossing failure due to low local density (an Allee Effect) while dispersal alleviates competition through dispersal of propagules from high- to low-density patches. Because local density is spatiotemporally heterogenous in our model, selfing does not eliminate heterogeneity in competition, so dispersal remains beneficial even under full selfing. Hence the Good Coloniser Syndrome is evolutionarily stable under a broad range of conditions, and both negative and positive relationships between dispersal and selfing are possible, depending on the environment. Our model thus accommodates positive empirical relationships between dispersal and selfing not predicted by previous theoretical work and provides additional explanations for negative relationships.

Data from: Geographical variation in the heterogeneity of mutualistic networks

Plant–animal mutualistic networks are characterized by highly heterogeneous degree distributions. The majority of species interact with few partner species, while a small number are highly connected to form network hubs that are proposed to play an important role in community stability. It has not been investigated, however, if or how the degree distributions vary among types of mutualisms or communities, or between plants and animals in the same network. Here, we evaluate the degree distributions of pollination and seed-dispersal networks, which are two major types of mutualistic networks that have often been discussed in parallel, using an index based on Pielou's evenness. Among 56 pollination networks we found strong negative correlation of the heterogeneity between plants and animals, and geographical shifts of network hubs from plants in temperate regions to animals in the tropics. For 28 seed-dispersal networks, by contrast, the correlation was positive, and there is no comparable geographical pattern. These results may be explained by evolution towards specialization in the presence of context-dependent costs that occur if plants share the animal species as interaction partner. How the identity of network hubs affects the stability and resilience of the community is an important question for future studies.

Data from: Genomic heterogeneity of historical gene flow between two species of newts inferred from transcriptome data

The role of gene flow in species formation is a major unresolved issue in speciation biology. Progress in this area requires information on the long-term patterns of gene flow between diverging species. Here, we used thousands of single-nucleotide polymorphisms derived from transcriptome resequencing and a method modeling the joint frequency spectrum of these polymorphisms to reconstruct patterns of historical gene flow between two Lissotriton newts: L. vulgaris (Lv) and L. montandoni (Lm). We tested several models of divergence including complete isolation and various scenarios of historical gene flow. The model of secondary contact received the highest support. According to this model, the species split from their common ancestor ca. 5.5 million years (MY) ago, evolved in isolation for ca. 2 MY, and have been exchanging genes for the last 3.5 MY Demographic changes have been inferred in both species, with the current effective population size of ca. 0.7 million in Lv and 0.2 million in Lm. The postdivergence gene flow resulted in two-directional introgression which affected the genomes of both species, but was more pronounced from Lv to Lm. Interestingly, we found evidence for genomic heterogeneity of interspecific gene flow. This study demonstrates the complexity of long-term gene flow between distinct but incompletely reproductively isolated taxa which divergence was initiated millions of years ago.

Data from: Food web heterogeneity and succession in created saltmarshes

Ecological restoration must achieve functional as well as structural recovery. Functional metrics for re-establishment of trophic interactions can be used to complement traditional monitoring of structural attributes. In addition, topographic effects on food web structure provide added information within a restoration context; often, created sites may require spatial heterogeneity to effectively match structure and function of natural habitats. We addressed both of these issues in our study of successional development of benthic food web structure, with focus on bottom-up-driven changes in macroinvertebrate consumer assemblages in the saltmarshes of the Venice Lagoon, Italy. We combined quantified estimates of the changing community composition with stable isotope data (13C:12C and 15N:14N) to compare the general trophic structure between created (2–14 years) marshes and reference sites and along topographic elevation gradients within saltmarshes. Macrofaunal invertebrate consumers exhibited local, habitat-specific trophic patterns. Stable isotope-based trophic structure changed with increasing marsh age, in particular with regard to mid-elevation (Salicornia) habitats. In young marshes, the mid-elevation consumer signatures resembled those of unvegetated ponds. The mid-elevation of older and natural marshes had a more distinct Salicornia zone food web, occasionally resembling that of the highest (Sarcocornia-dominated) elevation. In summary, this indicates that primary producers and availability of vascular plant detritus structure consumer trophic interactions and the flow of carbon. Functionally different consumers, subsurface-feeding detritivores (Oligochaeta) and surface grazers (Hydrobia sp.), showed distinct but converging trajectories of isotopic change over time, indicating that successional development may be asymmetric between 'brown' (detrital) guilds and 'green' (grazing) guilds in the food web. Synthesis and applications. Created marsh food webs converged into a natural state over about a decade, with successional shifts seen in both consumer community composition and stable isotope space. Strong spatial effects were noted, highlighting the utility of stable isotopes to evaluate functional equivalence in spatially heterogeneous systems. Understanding the recovery of functional properties such as food web support, and their inherent spatial variability, is key to planning and managing successful habitat restoration.

Data from: Spatial heterogeneity lowers rather than increases host-parasite specialization

Abiotic environmental heterogeneity can promote the evolution of diverse resource specialists, which in turn may increase the degree of host-parasite specialization. We coevolved Pseudomonas fluorescens and lytic phage ϕ2 in spatially structured populations, each consisting of two interconnected subpopulations evolving in the same or different nutrient media (homogeneous and heterogeneous environments, respectively). Counter to the normal expectation, host-parasite specialization was significantly lower in heterogeneous compared with homogeneous environments. This result could not be explained by dispersal homogenizing populations, as this would have resulted in the heterogeneous treatments having levels of specialization equal to or greater than that of the homogeneous environments. We argue that selection for costly generalists is greatest when the coevolving species are exposed to diverse environmental conditions and that this can provide an explanation for our results. A simple coevolutionary model of this process suggests that this can be a general mechanism by which environmental heterogeneity can reduce rather than increase host-parasite specialization.

Data from: Bulk genotyping of biopsies can create spurious evidence for heterogeneity in mutation content

When multiple samples are taken from the neoplastic tissues of a single patient, it is natural to compare their mutation content. This is often done by bulk genotyping of whole biopsies, but the chance that a mutation will be detected in bulk genotyping depends on its local frequency in the sample. When the underlying mutation count per cell is equal, homogenous biopsies will have more high-frequency mutations, and thus more detectable mutations, than heterogeneous ones. Using simulations, we show that bulk genotyping of data simulated under a neutral model of somatic evolution generates strong spurious evidence for non-neutrality, because the pattern of tissue growth systematically generates differences in biopsy heterogeneity. Any experiment which compares mutation content across bulk-genotyped biopsies may therefore suggest mutation rate or selection intensity variation even when these forces are absent. We discuss computational and experimental approaches for resolving this problem.

opencc-zeroDec 2010View details →

Data from: Evolution of virulence in heterogeneous host communities under multiple trade-offs

Many pathogens and parasites are transmitted through hosts that differ in species, sex, genotype, or immune status. In addition, virulence (here defined as disease-induced mortality) and transmission can vary during the infectious period within hosts of different state. Most models of virulence evolution assume that transmission and virulence are constant over the infectious period and that the host population is homogenous. Here, we examine a multi-species SIR model where transmission occurs within and between species, and transmission and virulence varied during the infectious period. This allows us to understand virulence evolution in a broader range of situations that characterize many emerging diseases. Because emerging pathogens are by definition new to their host populations, they should be expected to rapidly adapt after emergence. We illustrate these evolutionary effects using the framework of adaptive dynamics to examine how virulence evolves after emergence in response to the relative strength of selection on pathogen fitness and mutational variance for virulence. We illustrate the role of evolution by simulating adaptive walks to an evolutionarily stable virulence. We found that the magnitude of between-species transmission and the relative timing of transmission and mortality across species were of primary importance for determining the evolutionarily stable virulence.

opencc-zeroDec 2010View details →

Data from: Abiotic heterogeneity drives parasite local adaptation in coevolving bacteria and phages

Spatial abiotic heterogenity can result in divergent selection, hence might increase the magnitude of host-parasite local adaptation (the mean difference in fitness of sympatric versus allopatric host-parasite combinations). We explicitly tested this hypothesis by measuring local adapation in experimentally coevolved populations of bacteria and viruses evolved in the same or different nutrient media. Consistent with previous work, we found that mean levels of evolved phage infectivity and bacteria resistance varied with nutrient concentration, with maximal levels at nutrient concentrations that supported the greatest densities of bacteria. Despite this variation in evolved mean infectivity and resistance between treatments, we found that parasite local adaptation was greatly increased when measured between populations evolved in different, compared to the same, media. This pattern is likely to have resulted from different media imposing divergent selection on bacterial hosts, and phages in turn adapting to their local hosts. These results demonstrate that the abiotic environment can play a strong and predictable role in driving patterns of local adaptation.

Data from: Quantifying individual heterogeneity and its influence on life-history trajectories: different methods for different questions and contexts

Heterogeneity among individuals influences the life-history trajectories we observe at the population level because viability selection, selective immigration and emigration processes, and ontogeny change the proportion of individuals with specific trait values with increasing age. Here, we review the two main approaches that have been proposed to account for these processes in life-history trajectories, contrasting how they quantify ontogeny and selection, and proposing ways to overcome some of their limitations. Nearly all existing approaches to model individual heterogeneity assume either a single normal distribution or a priori known groups of individuals. Ontogenetic processes, however, can vary across individuals through variation in life-history tactics. We show the usefulness of describing ontogenetic processes by modelling trajectories with a mixture model that focuses on heterogeneity in life-history tactics. Additionally, most methods examine individual heterogeneity in a single trait, ignoring potential correlations among multiple traits caused by latent common sources of individual heterogeneity. We illustrate the value of using a joint modelling approach to assess the presence of a shared latent correlation and its influence on life-history trajectories. We contrast the strengths and limitations of different methods for different research questions, and we exemplify the differences among methods using empirical data from long-term studies of ungulates.

Data from: Shear-induced orientational dynamics and spatial heterogeneity in suspensions of motile phytoplankton

Fluid flow, ubiquitous in natural and man-made environments, has the potential to profoundly impact the transport of microorganisms, including phytoplankton in aquatic habitats and bioreactors. Yet, the effect of ambient flow on the swimming behavior of phytoplankton has remained poorly understood, largely due to the difficulty of observing cell-flow interactions at the microscale. Here, we present microfluidic experiments where we tracked individual cells for four species of motile phytoplankton exposed to a spatially non-uniform fluid shear rate, characteristic of many flows in natural and artificial environments. We observed that medium to high mean shear rates (1 – 25 s^-1) produce heterogeneous cell concentrations in the form of regions of accumulation and regions of depletion. The location of these regions relative to the flow depends on the cells' propulsion mechanism, body shape, and flagellar arrangement, as captured by an effective aspect ratio. Species having a large effective aspect ratio accumulated in the high-shear regions, due to shear-induced alignment of the swimming orientation with the fluid streamlines. Species having an effective aspect ratio close to unity exhibited little preferential accumulation at low to moderate flow rates, but strongly accumulated in the low-shear regions under high flow conditions, potentially due to an active, behavioral response of cells to shear. These observations demonstrate that ambient fluid flow can strongly affect the motility and spatial distribution of phytoplankton and highlight the rich dynamics emerging from the interaction between motility, morphology, and flow.