Development of resident and migratory three-spined stickleback, Gasterosteus aculeatus.

The three-spined stickleback (Gasterosteus aculeatus) is a teleost fish and a model organism in evolutionary ecology, useful for both laboratory and natural experiments. It is especially valued for the substantial intraspecific variation in morphology, behaviour and genetics. Classic work of Swarup (1958) has described the development in the laboratory of embryos from a single freshwater population, but this was carried out at higher temperature than many stickleback would encounter in the wild and variation between populations was not addressed. Here we describe the development of embryos from two sympatric, saltwater ecotypes of stickleback from North Uist, Scotland raised at 14°C, the approximate temperature of North Uist lochs in the breeding season. The two ecotypes were (a) a large, migratory form in which the adults are completely plated with bony armour and (b) a smaller, low-plated form that is resident year-round in saltwater lagoons. By monitoring embryos every 24-hours post fertilisation, important characteristics of development were observed and photographed to provide a reference for North Uist ecotypes at this temperature. Hatching success was greater than 85% and did not differ between resident and migratory stickleback, but migratory eggs hatched significantly earlier than the resident ecotype. Our work provides a framework that can now be used to compare stickleback populations that may also grow in distinct environmental conditions, to help understand the breadth of normal developmental features and to characterise abnormal development.

Physiological responses of Atlantic cod to climate change indicate that coastal ecotypes may be better adapted to tolerate ocean stressors.

Healthy ecosystems and species have some degree of resilience to changing conditions, however as the frequency and severity of environmental changes increase, resilience may be diminished or lost. In Sweden, one example of a species with reduced resilience is the Atlantic cod (Gadus morhua). This species has been subjected to overfishing, and with additional pressures such as habitat degradation and changing environmental conditions there has been little to no recovery, despite more than a decade of management actions. Given the historical ecological, economical, and cultural significance of cod, it is important to understand how Atlantic cod respond to global climate change to recover and sustainably manage this species in the future. A multi-stressor experiment was conducted to evaluate physiological responses of juvenile cod exposed to warming, ocean acidification, and freshening, changes expected to occur in their nursery habitat. The response to single drivers showed variable effects related to fish biometrics and increased levels of oxidative stress dependent parameters. Importantly, two separate responses were seen within a single treatment for the multi-stressor and freshening groups. These within-treatment differences were correlated to genotype, with the offshore ecotype having a heightened stress response compared to the coastal ecotype, which may be better adapted to tolerate future changes. These results demonstrate that, while Atlantic cod have some tolerance for future changes, ecotypes respond differently, and cumulative effects of multiple stressors may lead to deleterious effects for this important species.

A promoter polymorphism defines distinct roles in anther development for Col-0 and Ler-0 alleles of Arabidopsis ACYL-COA BINDING PROTEIN3.

Acyl-CoA-Binding Proteins (ACBPs) bind acyl-CoA esters and function in lipid metabolism. Although acbp3-1, the ACBP3 mutant in Arabidopsis thaliana ecotype Col-0, displays normal floral development, the acbp3-2 mutant from ecotype Ler-0 characterized herein exhibits defective adaxial anther lobes and improper sporocyte formation. To understand these differences and identify the role of ERECTA in ACBP3 function, the acbp3 mutants and acbp3-erecta (er) lines were analyzed by microscopy for anther morphology and high-performance liquid chromatography for lipid composition. Defects in Landsberg anther development were related to the ERECTA-mediated pathway because the progenies of acbp3-2 × La-0 and acbp3-1 × er-1 in Col-0 showed normal anthers, contrasting to that of acbp3-2 in Ler-0. Polymorphism in the regulatory region of ACBP3 enabled its function in anther development in Ler-0 but not Col-0 which harbored an AT-repeat insertion. ACBP3 expression and anther development in acbp3-2 were restored using ACBP3pro (Ler)::ACBP3 not ACBP3pro (Col)::ACBP3. SPOROCYTELESS (SPL), a sporocyte formation regulator activated ACBP3 transcription in Ler-0 but not Col-0. For anther development, the ERECTA-related role of ACBP3 is required in Ler-0, but not Col-0. The disrupted promoter regulatory region for SPL binding in Col-0 eliminates the role of ACBP3 in anther development.

Envirotyping within a multi-environment trial allowed identifying genetic determinants of winter oilseed rape yield stability.

KEY MESSAGE: A comprehensive environmental characterization allowed identifying stable and interactive QTL for seed yield: QA09 and QC09a were detected across environments; whereas QA07a was specifically detected on the most stressed environments. A main challenge for rapeseed consists in maintaining seed yield while adapting to climate changes and contributing to environmental-friendly cropping systems. Breeding for cultivar adaptation is one of the keys to meet this challenge. Therefore, we propose to identify the genetic determinant of seed yield stability for winter oilseed rape using GWAS coupled with a multi-environmental trial and to interpret them in the light of environmental characteristics. Due to a comprehensive characterization of a multi-environmental trial using 79 indicators, four contrasting envirotypes were defined and used to identify interactive and stable seed yield QTL. A total of four QTLs were detected, among which, QA09 and QC09a, were stable (detected at the multi-environmental trial scale or for different envirotypes and environments); and one, QA07a, was specifically detected into the most stressed envirotype. The analysis of the molecular diversity at QA07a showed a lack of genetic diversity within modern lines compared to older cultivars bred before the selection for low glucosinolate content. The results were discussed in comparison with other studies and methods as well as in the context of breeding programs.

Ecotype variation in the endemic tree Callicarpa subpubescens on small oceanic islands: genetic, phenotypic, and environmental insights.

Callicarpa subpubescens, endemic to the Ogasawara Islands, is suggested to have multiple ecotypes in the Hahajima Islands, specifically in the central part of the Ogasawara Islands. In this study, associations between genetic groups and spatial distribution, habitat, leaf morphology, size structure, and flowering time of each genetic group were investigated on Hahajima and the satellite Imoutojima Islands. Genetic groups were identified using EST-SSR markers, revealing four ecotypes named based on morphological features: Dwarf (D), Glabrescent (G), Tall (T), and Middle (M), with M being a result of the hybridization of G and T. Ecotype D, adapted to dry environments, is characterized by small tree size, dense thick leaves with abundant hairs, and is distributed in dry scrub. Ecotype G, adapted to understory of mesic forests, lacks leaf hairs. Ecotype T, adapted to the canopy of mesic forests, has hairy leaves and is tall in tree height. Ecotype M, adapted to the canopy of mesic scrub or edges of mesic forests, has hairy leaves but with a shorter tree height than ecotype T. Flowering peaks differed among all ecotype pairs except G and M, but the flowering times more or less overlapped among all ecotypes, suggesting that pre-mating isolation among ecotypes is not perfect. Post-mating isolation is considered absent, as there were no differences in the results, germination, and survival rates of one-year seedlings among inter- and intra-ecotype crossings. The existence of such ecotypes provides valuable insights into the ongoing speciation processes adapting to the oceanic island environments.

Skeletonema marinoi ecotypes show specific habitat-related responses to fluctuating light supporting high potential for growth under photobioreactor light regime.

Diatoms are a diverse group of phytoplankton usually dominating areas characterized by rapidly shifting light conditions. Because of their high growth rates and interesting biochemical profile, their biomass is considered for various commercial applications. This study aimed at identifying strains with superior growth in a photobioreactor (PBR) by screening the natural intraspecific diversity of ecotypes isolated from different habitats. We investigated the effect of PBR light fluctuating on a millisecond scale (FL, simulating the light in a PBR) on 19 ecotypes of the diatom Skeletonema marinoi isolated from the North Sea-Baltic Sea area. We compare growth, pigment ratios, phylogeny, photo-physiological variables and photoacclimation strategies between all strains and perform qPCR and absorption spectra analysis on a subset of strains. Our results show that the ecotypes responded differently to FL, and have contrasting photo-physiological and photoprotective strategies. The strains from Kattegat performed better in FL, and shared common photoacclimation and photoprotection strategies that are the results of adaptation to the specific light climate of the Kattegat area. The strains that performed better with FL conditions had a high light (HL)-acclimated phenotype coupled with unique nonphotochemical quenching features. Based on their characteristics, three strains were identified as good candidates for growth in PBRs.

Phenotypic characterization and analysis of genetic diversity between commercial crossbred and indigenous chickens from three different agro-ecological zones using DArT-Seq technology.

Indigenous and were used to study genetic diversity and population structure analyses. Polymorphism information content (PIC) values ranged from 0.0 to 0.5, with 21,285 SNP markers (35%) being in the lowest PIC value range (0 to 0.15) while 13,511 (commercial chickens have developed unique adaptations to their environments, which may include nutrition, pathogens, and thermal stress. Besides, environmental pressures and artificial selection have generated significant genome-wide divergence in chickens, as those selection pressures contribute a considerable evolutionary force to phenotypic and genotypic differentiation. Herein, we determined genomic diversity of indigenous chickens from semi-deciduous rainforest (SDR), coastal savannah (CS) and Guinea savannah (GS) agro-ecological zones (AEZs) in Ghana and commercial crossbreds (CC) reared at the Kwame Nkrumah University of Science and Technology (KNUST). We generated SNP markers from 82 chickens (62 indigenous chicken ecotypes and 26 commercial crossbred ecotype) using DArT-Seq technology. A total of 85,396 SNP markers were generated and after filtering the data, 58,353 markers 21%) were in the highest PIC value range (0.45 to 0.50). The CC were more genetically diverse than the indigenous birds, with the highest expected heterozygosity value of 0.220. Between the commercial crossbreds population and the indigenous ecotypes, pairwise FST values were estimated to be 0.105 between CS, 0.096 between SDF, and 0.133 between GS. Furthermore, PCA analysis showed that the CC, SDF and GS chickens clustered together and are genetically distant from the commercial crossbred. We herein show that chickens from the AEZs studied can be considered as one population. However, due the abundance of agro-byproducts in the SDR compared to the CS and GS, chickens from the SDR AEZ had better growth compared to their counterparts. It is suggested that the genetic diversity within the local ecotypes could form the basis for genetic improvement.

Identifying potential provenances for climate-change adaptation using spatially variable coefficient models.

BACKGROUND: Selection of climate-change adapted ecotypes of commercially valuable species to date relies on DNA-assisted screening followed by growth trials. For trees, such trials can take decades, hence any approach that supports focussing on a likely set of candidates may save time and money. We use a non-stationary statistical analysis with spatially varying coefficients to identify ecotypes that indicate first regions of similarly adapted varieties of Douglas-fir (Pseudotsuga menziesii (Mirbel) Franco) in North America. For over 70,000 plot-level presence-absences, spatial differences in the survival response to climatic conditions are identified. RESULTS: The spatially-variable coefficient model fits the data substantially better than a stationary, i.e. constant-effect analysis (as measured by AIC to account for differences in model complexity). Also, clustering the model terms identifies several potential ecotypes that could not be derived from clustering climatic conditions itself. Comparing these six identified ecotypes to known genetically diverging regions shows some congruence, as well as some mismatches. However, comparing ecotypes among each other, we find clear differences in their climate niches. CONCLUSION: While our approach is data-demanding and computationally expensive, with the increasing availability of data on species distributions this may be a useful first screening step during the search for climate-change adapted varieties. With our unsupervised learning approach being explorative, finely resolved genotypic data would be helpful to improve its quantitative validation.

Sustainable Agrochemical Exploitation of Ligurian Aromatic Plant Ecotypes.

The characterization of local improved varieties as well as the reduction of synthetic chemical fertilizers are sustainable approaches in the vision of a new precision Farming. Aim of our study was to improve the geographical characterization of local ecotypes and to identify peculiar features of new crops in terms of bioactive compounds. NMR and LC-MS metabolite profiling approaches followed by multivariate data analysis were applied to characterize local rosemary and garlic ecotypes. With the aim of applying for a protected designation of origin, orthogonal partial least squares discriminant analysis (OPLS-DA) was used to identify representative sensory quality indicators for Vessalico garlic and rosemary "Eretto Liguria" local ecotypes, Variable Influence on Projections (VIP) values of OPLS-DA indicated six metabolites as quality indicators for Vessalico garlic and sixteen metabolites as quality indicators for rosemary "Eretto Liguria". Finally, to discover and utilize new ecotypes in a sustainable way, Vessalico garlic extracts antiviral activity, previously evaluated against Tomato brown rugose fruit virus (ToBRFV), a Tobamovirus affecting tomato crops, was extended to Pepino mosaic virus (PepMV) with positive results.

Chromosomal inversions from an initial ecotypic divergence drive a gradual repeated radiation of Galápagos beetles.

Island faunas exhibit some of the most iconic examples where similar forms repeatedly evolve within different islands. Yet, whether these deterministic evolutionary trajectories within islands are driven by an initial, singular divergence and the subsequent exchange of individuals and adaptive genetic variation between islands remains unclear. Here, we study a gradual, repeated evolution of low-dispersive highland ecotypes from a dispersive lowland ecotype of Calosoma beetles along the island progression of the Galápagos. We show that repeated highland adaptation involved selection on multiple shared alleles within extensive chromosomal inversions that originated from an initial adaptation event on the oldest island. These highland inversions first spread through dispersal of highland individuals. Subsequent admixture with the lowland ecotype resulted in polymorphic dispersive populations from which the highland populations evolved on the youngest islands. Our findings emphasize the significance of an ancient divergence in driving repeated evolution and highlight how a mixed contribution of inter-island colonization and within-island evolution can shape parallel species communities.

Chromosome-scale Genome Assembly of the Rough Periwinkle Littorina saxatilis.

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation.

Metabolomic profiling of wild rooibos (Aspalathus linearis) ecotypes and their antioxidant-derived phytopharmaceutical potential.

INTRODUCTION: Aspalathus linearis (commonly known as rooibos) is endemic to the Cape Floristic Region of South Africa and is a popular herbal drink and skin phytotherapeutic ingredient, with health benefits derived primarily from its unique phenolic content. Several, seemingly habitat-specific ecotypes from the Cederberg (Western Cape) and Northern Cape have morphological, ecological, genetic and biochemical differences. OBJECTIVES AND METHODS: Despite the commercial popularity of the cultivated variety, the uncultivated ecotypes are largely understudied. To address gaps in knowledge about the biochemical constituency, ultra-performance liquid chromatography-mass spectrometry analysis of fifteen populations was performed, enabling high-throughput metabolomic fingerprinting of 50% (v/v) methanolic extracts. Antioxidant screening of selected populations was performed via three assays and antimicrobial activity on two microbial species was assessed. The metabolomic results were corroborated with total phenolic and flavonoid screening of the extracts. RESULTS AND DISCUSSION: Site-specific chemical lineages of rooibos ecotypes were confirmed via multivariate data analyses. Important features identified via PLS-DA disclosed higher relative abundances of certain tentative metabolites (e.g., rutin, aspalathin and apiin) present in the Dobbelaarskop, Blomfontein, Welbedacht and Eselbank sites, in comparison to other locations. Several unknown novel metabolites (e.g., m/z 155.0369, 231.0513, 443.1197, 695.2883) are responsible for metabolomic separation of the populations, four of which showed higher amounts of key metabolites and were thus selected for bioactivity analysis. The Welbedacht and Eselbank site 2 populations consistently displayed higher antioxidant activities, with 2,2-azino-bis (3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) radical scavenging activities of 679.894 ± 3.427 µmol Trolox/g dry matter and 635.066 ± 5.140 µmol Trolox/g dry matter, respectively, in correlation with a high number of phenolic and flavonoid compounds. The contribution of the individual metabolites to the pharmacological effectiveness of rooibos remains unknown and as such, further structural elucidation and phytopharmacological testing is thus urgently needed.

Phytochemical Composition and Antimicrobial Properties of New Lavandula angustifolia Ecotypes.

The purpose of this study was to characterize ethanol extracts from leaves and flowers of two ecotypes (PL-intended for industrial plantations and KC-intended for cut flowers) of Lavandula angustifolia Mill. The plant was cultivated in 2019 in southern Poland as part of a long-term research plan to develop new varieties resistant to difficult environmental conditions. The collected leaves and flowers were used to prepare ethanol extracts, which were then analyzed in terms of phytochemical composition and antioxidant, bactericidal, and fungicidal properties. Using UPLC techniques, 22 compounds belonging to phenolic acids and flavonoids were identified. UPLC test results indicated that ethanol extracts from leaves and flowers differ in phytochemical composition. Lower amounts of phenolic acids and flavonoids were identified in leaf extracts than in flower extracts. The predominant substances in the flower extracts were rosmarinic acid (829.68-1229.33 µg/g), ferulic acid glucoside III (810.97-980.55 µg/g), and ferulic acid glucoside II (789.30-885.06 µg/g). Ferulic acid glucoside II (3981.95-6561.19 µg/g), ferulic acid glucoside I (2349.46-5503.81 µg/g), and ferulic acid glucoside III (1303.84-2774.17 µg/g) contained the highest amounts in the ethanol extracts of the leaves. The following substances were present in the extracts in trace amounts or at low levels: apigenin, kaempferol, and caftaric acid. Leaf extracts of the PL ecotype quantitatively (µg/g) contained more phytochemicals than leaf extracts of the KC ecotype. The results obtained in this study indicate that antioxidant activity depends on the ecotype. Extracts from the PL ecotype have a better ability to eliminate free radicals than extracts from the KC ecotype. At the same time, it was found that the antioxidant activity (total phenolic content, ABTS•+, DPPH•, and FRAP) of PL ecotype leaf extracts was higher (24.49, 177.75, 164.88, and 89.10 µmol (TE)/g) than that determined in flower extracts (15.84, 125.05, 82.35, and 54.64 µmol (TE)/g). The test results confirmed that leaf and flower extracts, even at low concentrations (0.313-0.63%), significantly inhibit the growth of selected Gram-negative and Gram-positive bacteria and Candida yeasts. Inhibition of mold growth was observed at a dose extract of at least 1 mL/100 mL.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization.

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

A single Prochlorococcus ecotype dominates the tropical Bay of Bengal with ultradian growth.

The Bay of Bengal (BoB) spans >2.2 million km2 in the northeastern Indian Ocean and is bordered by dense populations that depend upon its resources. Over recent decades, a shift from larger phytoplankton to picoplankton has been reported, yet the abundance, activity, and composition of primary producer communities are not well-characterized. We analysed the BoB regions during the summer monsoon. Prochlorococcus ranged up to 3.14 × 105 cells mL-1 in the surface mixed layer, averaging 1.74 ± 0.46 × 105 in the upper 10 m and consistently higher than Synechococcus and eukaryotic phytoplankton. V1-V2 rRNA gene amplicon analyses showed the High Light II (HLII) ecotype formed 98 ± 1% of Prochlorococcus amplicons in surface waters, comprising six oligotypes, with the dominant oligotype accounting for 65 ± 4% of HLII. Diel sampling of a coherent water mass demonstrated evening onset of cell division and rapid Prochlorococcus growth between 1.5 and 3.1 div day-1, based on cell cycle analysis, as confirmed by abundance-based estimates of 2.1 div day-1. Accumulation of Prochlorococcus produced by ultradian growth was restricted by high loss rates. Alongside prior Arabian Sea and tropical Atlantic rates, our results indicate Prochlorococcus growth rates should be reevaluated with greater attention to latitudinal zones and influences on contributions to global primary production.

Multicellular ecotypes shape progression of lung adenocarcinoma from ground-glass opacity toward advanced stages.

Lung adenocarcinoma is a type of cancer that exhibits a wide range of clinical radiological manifestations, from ground-glass opacity (GGO) to pure solid nodules, which vary greatly in terms of their biological characteristics. Our current understanding of this heterogeneity is limited. To address this gap, we analyze 58 lung adenocarcinoma patients via machine learning, single-cell RNA sequencing (scRNA-seq), and whole-exome sequencing, and we identify six lung multicellular ecotypes (LMEs) correlating with distinct radiological patterns and cancer cell states. Notably, GGO-associated neoantigens in early-stage cancers are recognized by CD8+ T cells, indicating an immune-active environment, while solid nodules feature an immune-suppressive LME with exhausted CD8+ T cells, driven by specific stromal cells such as CTHCR1+ fibroblasts. This study also highlights EGFR(L858R) neoantigens in GGO samples, suggesting potential CD8+ T cell activation. Our findings offer valuable insights into lung adenocarcinoma heterogeneity, suggesting avenues for targeted therapies in early-stage disease.

Morphometric differentiation of three chicken ecotypes of Ethiopia using multivariate analysis.

Twenty-one morphometric traits were measured on 770 extensively managed indigenous chickens in the western zone of Tigray, comprising 412 hens and 358 cocks in three agro-ecologies. The quantitative traits for male and female chicken ecotypes were separately analyzed using multivariate analysis with SAS 2008. Four and seven principal components accounted for about 74.26% and 69.77% of the total variability in morphometric traits for males and females, respectively. Earlobe length, wingspan, skull length, and shank length werethe most important traits for discriminating among female chicken ecotypes, while wingspan, neck length, earlobe length, spur length, body length, and shank length were the most important discriminatory traits among male chicken ecotypes. The discriminant analysis accurately classified 97.3% of female and 100% of male chicken ecotypes. Cluster analysis revealed the genetic heterogeneity of indigenous chicken populations in both sexes. This finding suggests the presence of morphological variations among the indigenous chicken populations in the different agro-ecological zones, classified as distinct indigenous chicken ecotypes (Lowland, Midland, and Highland). Further DNA-based studies are needed to confirm and complement these morphological variations for effective conservation and the development of sustainable genetic improvement strategies for indigenous chicken populations in the region.

Rediversification following ecotype isolation reveals hidden adaptive potential.

Microbial communities play a critical role in ecological processes, and their diversity is key to their functioning. However, little is known about whether communities can regenerate ecological diversity following ecotype removal or extinction and how the rediversified communities would compare to the original ones. Here, we show that simple two-ecotype communities from the E. coli long-term evolution experiment (LTEE) consistently rediversified into two ecotypes following the isolation of one of the ecotypes, coexisting via negative frequency-dependent selection. Communities separated by more than 30,000 generations of evolutionary time rediversify in similar ways. The rediversified ecotype appears to share a number of growth traits with the ecotype it replaces. However, the rediversified community is also different from the original community in ways relevant to the mechanism of ecotype coexistence-for example, in stationary phase response and survival. We found substantial variation in the transcriptional states between the two original ecotypes, whereas the differences within the rediversified community were comparatively smaller, although the rediversified community showed unique patterns of differential expression. Our results suggest that evolution may leave room for alternative diversification processes even in a maximally reduced community of only two strains. We hypothesize that the presence of alternative evolutionary pathways may be even more pronounced in communities of many species where there are even more potential niches, highlighting an important role for perturbations, such as species removal, in evolving ecological communities.

The genomic consequences of selection across development.

Understanding how natural selection drives diversification in nature has been at the forefront of biological research for over a century. The main idea is simple: natural selection favours individuals best suited to pass on their genes. However, the journey from birth to reproduction is complex as organisms experience multiple developmental stages, each influenced by genetic and environmental factors (Orr, 2009). These complexities compound even further as each stage of development might be governed by a unique underlying set of alleles and genes. In this issue of Molecular Ecology, Goebl et al. (2022) examine the role of natural selection in driving ecotypic divergence across different life history stages of the prairie sunflower Helianthus petiolaris. The authors used reciprocal transplant experiments, demographic models, and genomic sequencing to explore fitness variation across developmental stages. They show how natural selection impacts population divergence across multiple life history stages and evaluate the resulting allele frequency changes. Goebl et al. link these results to the role of chromosomal inversions, thus furthering our understanding of how ecological divergence proceeds in the face of gene flow. Below, we explore these results in detail and complement their interpretation by considering the evolution of genetic correlations amongst traits governing fitness.

Division time (td) for in situ growth measurements demonstrates thermal ecotypes of Karlodinium veneficum.

The toxic dinoflagellate Karlodinium veneficum forms fish killing blooms in temperate estuaries worldwide. These blooms have variable toxicity which may be related to bloom stage and in situ growth rates of the constituent K. veneficum cells. Measurement of in situ growth rates is challenging and methods such as the mitotic index technique require knowledge of the dynamics of cell division. In order to better understand these dynamics, we determined the duration of cell division (td) in four geographically distinct laboratory strains of K. veneficum at three different environmentally relevant temperatures. The results demonstrated that the td value for each strain, growing at strain-specific optimal temperatures, was 1.6 ± 0.1 h. This value corresponded to a range of growth rates from 0.17 ± 0.08 d-1 to 0.62 ± 0.07 d-1. Equivalent values of td spread across four geographically distinct laboratory strains and a nearly fourfold range of growth rates implies that 1.6 h represents the td value of K. veneficum. Additionally, temperature conditions yielding this value for td and the highest growth rates varied among strains, indicating cold-adapted (Norway), warm-adapted (Florida, USA), and eurythermally-adapted (Maryland, USA) strains. These differences have been apparently retained in culture over many years, indicating a conserved genetic basis that suggests distinct thermal ecotypes of the morphospecies K. veneficum. This knowledge together with the first estimate of td for K. veneficum will be useful in future field studies aimed at correlating bloom toxicity with in situ growth rate using the mitotic index technique.