Thyasirid species composition (Bivalvia: Thyasiridae) and genetic connectivity of Parathyasira equalis (A. E. Verrill & K. J. Bush, 1898) in deep basins of sub-Arctic fjords.

BACKGROUND: Thyasirid bivalves are often recorded as a dominant component of macrobenthic infaunal communities in depositional environments such as fjord basins. Fjord basins comprise patchy soft-bottom habitats bounded by steep walls and sills; however, little is known how this semi-isolated nature of fjords affects benthic populations. Accordingly, data on the composition and population connectivity of thyasirids can provide valuable information on the ecology of these ecosystems. RESULTS: The species composition of thyasirid bivalves has been studied in the basins of three sub-Arctic fjords (Nordland, Northern Norway). Overall, six thyasirid species were recorded: Parathyasira equalis, Parathyasira dunbari, Mendicula ferruginosa, Genaxinus eumyarius, Thyasira sarsii, and Thyasira obsoleta. The species composition remained stable within the basins during the sampling period (2013-2020) and suggested the importance of local reproduction over advection of individuals for population dynamics. Only one species, Parathyasira equalis, was common in all fjords. We have further investigated the population genetics of this species by combining two types of genetic markers: a 579 bp fragment of the cytochrome c oxidase subunit I (COI) gene and 4043 single-nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing. The latter provided a more in-depth resolution on the population genetics of this species and revealed a weak but significant differentiation of populations within fjords, further indicating limited connectivity between basins. CONCLUSION: Based on our findings, we conclude that limited dispersal between the basin communities results in weakly connected populations and might be an important structuring factor for macrobenthic communities.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts.

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Host habitat rather than evolutionary history explains gut microbiome diversity in sympatric stickleback species.

Host-associated microbiota can influence host phenotypic variation, fitness and potential to adapt to local environmental conditions. In turn, both host evolutionary history and the abiotic and biotic environment can influence the diversity and composition of microbiota. Yet, to what extent environmental and host-specific factors drive microbial diversity remains largely unknown, limiting our understanding of host-microbiome interactions in natural populations. Here, we compared the intestinal microbiota between two phylogenetically related fishes, the three-spined stickleback (Gasterosteus aculeatus) and the nine-spined stickleback (Pungitius pungitius) in a common landscape. Using amplicon sequencing of the V3-V4 region of the bacterial 16S rRNA gene, we characterised the α and ß diversity of the microbial communities in these two fish species from both brackish water and freshwater habitats. Across eight locations, α diversity was higher in the nine-spined stickleback, suggesting a broader niche use in this host species. Habitat was a strong determinant of ß diversity in both host species, while host species only explained a small fraction of the variation in gut microbial composition. Strong habitat-specific effects overruled effects of geographic distance and historical freshwater colonisation, suggesting that the gut microbiome correlates primarily with local environmental conditions. Interestingly, the effect of habitat divergence on gut microbial communities was stronger in three-spined stickleback than in nine-spined stickleback, possibly mirroring the stronger level of adaptive divergence in this host species. Overall, our results show that microbial communities reflect habitat divergence rather than colonisation history or dispersal limitation of host species.

The fitness landscape of a community of Darwin's finches.

Divergent natural selection should lead to adaptive radiation-that is, the rapid evolution of phenotypic and ecological diversity originating from a single clade. The drivers of adaptive radiation have often been conceptualized through the concept of "adaptive landscapes," yet formal empirical estimates of adaptive landscapes for natural adaptive radiations have proven elusive. Here, we use a 17-year dataset of Darwin's ground finches (Geospiza spp.) at an intensively studied site on Santa Cruz (Galápagos) to estimate individual apparent lifespan in relation to beak traits. We use these estimates to model a multi-species fitness landscape, which we also convert to a formal adaptive landscape. We then assess the correspondence between estimated fitness peaks and observed phenotypes for each of five phenotypic modes (G. fuliginosa, G. fortis [small and large morphotypes], G. magnirostris, and G. scandens). The fitness and adaptive landscapes show 5 and 4 peaks, respectively, and, as expected, the adaptive landscape was smoother than the fitness landscape. Each of the five phenotypic modes appeared reasonably close to the corresponding fitness peak, yet interesting deviations were also documented and examined. By estimating adaptive landscapes in an ongoing adaptive radiation, our study demonstrates their utility as a quantitative tool for exploring and predicting adaptive radiation.

CircPrime: a web-based platform for design of specific circular RNA primers.

BACKGROUND: Circular RNAs (circRNAs) are covalently closed-loop RNAs with critical regulatory roles in cells. Tens of thousands of circRNAs have been unveiled due to the recent advances in high throughput RNA sequencing technologies and bioinformatic tools development. At the same time, polymerase chain reaction (PCR) cross-validation for circRNAs predicted by bioinformatic tools remains an essential part of any circRNA study before publication. RESULTS: Here, we present the CircPrime web-based platform, providing a user-friendly solution for DNA primer design and thermocycling conditions for circRNA identification with routine PCR methods. CONCLUSIONS: User-friendly CircPrime web platform ( http://circprime.elgene.net/ ) works with outputs of the most popular bioinformatic predictors of circRNAs to design specific circular RNA primers. CircPrime works with circRNA coordinates and any reference genome from the National Center for Biotechnology Information database).

The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia.

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.

Contrasting population genetic responses to migration barriers in two native and an invasive freshwater fish.

Habitat fragmentation impacts the distribution of genetic diversity and population genetic structure. Therefore, protecting the evolutionary potential of species, especially in the context of the current rate of human-induced environmental change, is an important goal. In riverine ecosystems, migration barriers affect the genetic structure of native species, while also influencing the spread of invasive species. In this study, we compare genetic patterns of two native and one highly invasive riverine fish species in a Belgian river basin, namely the native three-spined stickleback (Gasterosteus aculeatus) and stone loach (Barbatula barbatula), and the non-native and invasive topmouth gudgeon (Pseudorasbora parva). We aimed to characterize both natural and anthropogenic determinants of genetic diversity and population genetic connectivity. Genetic diversity was highest in topmouth gudgeon, followed by stone loach and three-spined stickleback. The correlation between downstream distance and genetic diversity, a pattern often observed in riverine systems, was only marginally significant in stone loach and three-spined stickleback, while genetic diversity strongly declined with increasing number of barriers in topmouth gudgeon. An Isolation-By-Distance pattern characterizes the population genetic structure of each species. Population differentiation was only associated with migration barriers in the invasive topmouth gudgeon, while genetic composition of all species seemed at least partially determined by the presence of migration barriers. Among the six barrier types considered (watermills, sluices, tunnels, weirs, riverbed obstructions, and others), the presence of watermills was the strongest driver of genetic structure and composition. Our results indicate that conservation and restoration actions, focusing on conserving genetic patterns, cannot be generalized across species. Moreover, measures might target either on restoring connectivity, while risking a rapid spread of the invasive topmouth gudgeon, or not restoring connectivity, while risking native species extinction in upstream populations.

Complete mitochondrial genomes and updated divergence time of the two freshwater clupeids endemic to Lake Tanganyika (Africa) suggest intralacustrine speciation.

BACKGROUND: The hydrogeological history of Lake Tanganyika paints a complex image of several colonization and adaptive radiation events. The initial basin was formed around 9-12 million years ago (MYA) from the predecessor of the Malagarasi-Congo River and only 5-6 MYA, its sub-basins fused to produce the clear, deep waters of today. Next to the well-known radiations of cichlid fishes, the lake also harbours a modest clade of only two clupeid species, Stolothrissa tanganicae and Limnothrissa miodon. They are members of Pellonulini, a tribe of clupeid fishes that mostly occur in freshwater and that colonized West and Central-Africa during a period of high sea levels during the Cenozoic. There is no consensus on the phylogenetic relationships between members of Pellonulini and the timing of the colonization of Lake Tanganyika by clupeids. RESULTS: We use short-read next generation sequencing of 10X Chromium libraries to sequence and assemble the full mitochondrial genomes of S. tanganicae and L. miodon. We then use Maximum likelihood and Bayesian inference to place them into the phylogeny of Pellonulini and other clupeiforms, taking advantage of all available full mitochondrial clupeiform genomes. We identify Potamothrissa obtusirostris as the closest living relative of the Tanganyika sardines and confirm paraphyly for Microthrissa. We estimate the divergence of the Tanganyika sardines around 3.64 MYA [95% CI: 0.99, 6.29], and from P. obtusirostris around 10.92 MYA [95% CI: 6.37-15.48]. CONCLUSIONS: These estimates imply that the ancestor of the Tanganyika sardines diverged from a riverine ancestor and entered the proto-lake Tanganyika around the time of its formation from the Malagarasi-Congo River, and diverged into the two extant species at the onset of deep clearwater conditions. Our results prompt a more thorough examination of the relationships within Pellonulini, and the new mitochondrial genomes provide an important resource for the future study of this tribe, e.g. as a reference for species identification, genetic diversity, and macroevolutionary studies.

The terroir of the finch: How spatial and temporal variation shapes phenotypic traits in DARWIN'S finches.

The term terroir is used in viticulture to emphasize how the biotic and abiotic characteristics of a local site influence grape physiology and thus the properties of wine. In ecology and evolution, such terroir (i.e., the effect of space or "site") is expected to play an important role in shaping phenotypic traits. Just how important is the pure spatial effect of terroir (e.g., differences between sites that persist across years) in comparison to temporal variation (e.g., differences between years that persist across sites), and the interaction between space and time (e.g., differences between sites change across years)? We answer this question by analyzing beak and body traits of 4388 medium ground finches (Geospiza fortis) collected across 10 years at three locations in Galápagos. Analyses of variance indicated that phenotypic variation was mostly explained by site for beak size (η 2 = 0.42) and body size (η 2 = 0.43), with a smaller contribution for beak shape (η 2 = 0.05) and body shape (η 2 = 0.12), but still higher compared to year and site-by-year effects. As such, the effect of terroir seems to be very strong in Darwin's finches, notwithstanding the oft-emphasized interannual variation. However, these results changed dramatically when we excluded data from Daphne Major, indicating that the strong effect of terroir was mostly driven by that particular population. These phenotypic results were largely paralleled in analyses of environmental variables (rainfall and vegetation indices) expected to shape terroir in this system. These findings affirm the evolutionary importance of terroir, while also revealing its dependence on other factors, such as geographical isolation.

Context-dependent parasite infection affects trophic niche in populations of sympatric stickleback species.

How parasites alter host feeding ecology remains elusive in natural populations. A powerful approach to investigate the link between infection and feeding ecology is quantifying unique and shared responses to parasite infection in related host species within a common environment. Here, 9 pairs of sympatric populations of the three-spined and nine-spined stickleback fishes were sampled across a range of freshwater and brackish habitats to investigate how parasites alter host feeding ecology: (i) biotic and abiotic determinants of parasite community composition, and (ii) to what extent parasite infection correlates with trophic niche specialization of the 2 species, using stable isotope analyses (δ15N and δ13C). It was determined that parasite community composition and host parasite load varied among sites and species and were correlated with dissolved oxygen. It was also observed that the digenean Cyathocotyle sp.'s abundance, a common directly infecting parasite with a complex life cycle, correlated with host δ13C in a fish species-specific manner. In 6 sites, correlations were found between parasite abundance and their hosts' feeding ecology. These effects were location-specific and occasionally host species or host size-specific. Overall, the results suggest a relationship between parasite infection and host trophic niche which may be an important and largely overlooked ecological factor. The population specificity and variation in parasite communities also suggest this effect is multifarious and context-dependent.

Association between Chromosome 4 and mercury accumulation in muscle of the three-spined stickleback (Gasterosteus aculeatus).

Anthropogenic stressors, such as pollutants, act as selective factors that can leave measurable changes in allele frequencies in the genome. Metals are of particular concern among pollutants, because of interference with vital biological pathways. We use the three-spined stickleback as a model for adaptation to mercury pollution in natural populations. We collected sticklebacks from 21 locations in Flanders (Belgium), measured the accumulated levels of mercury in the skeletal muscle tissue, and genotyped the fish by sequencing (GBS). The spread of muscle mercury content across locations was considerable, ranging from 21.5 to 327 ng/g dry weight (DW). We then conducted a genome-wide association study (GWAS) between 28,450 single nucleotide polymorphisms (SNPs) and the accumulated levels of mercury, using different approaches. Based on a linear mixed model analysis, the GWAS yielded multiple hits with a single top hit on Chromosome 4, with eight more SNPs suggestive of association. A second approach, a latent factor mixed model analysis, highlighted one single SNP on Chromosome 11. Finally, an outlier test identified one additional SNP on Chromosome 4 that appeared under selection. Out of all ten SNPs we identified as associated with mercury in muscle, three SNPs all located on Chromosome 4 and positioned within a 2.5 kb distance of an annotated gene. Based on these results and the genome coverage of our SNPs, we conclude that the selective effect of mercury pollution in Flanders causes a significant association with at least one locus on Chromosome 4 in three-spined stickleback.

Extensive standing genetic variation from a small number of founders enables rapid adaptation in Daphnia.

We lack a thorough understanding of the origin and maintenance of standing genetic variation that enables rapid evolutionary responses of natural populations. Whole genome sequencing of a resurrected Daphnia population shows that standing genetic variation in over 500 genes follows an evolutionary trajectory that parallels the pronounced and rapid adaptive evolution of multiple traits in response to predator-driven natural selection and its subsequent relaxation. Genetic variation carried by only five founding individuals from the regional genotype pool is shown to suffice at enabling the observed evolution. Our results provide insight on how natural populations can acquire the genomic variation, through colonization by a few regional genotypes, that fuels rapid evolution in response to strong selection pressures. While these evolutionary responses in our study population involved hundreds of genes, we observed no evidence of genetic erosion.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species.

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Reconciling seascape genetics and fisheries science in three codistributed flatfishes.

Uncertainty hampers innovative mixed-fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species-specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population-specific information. We conclude that the management of fish stocks has to consider species-specific genetic structures and may benefit from the documentation of the genetic seascape and life-history traits.

Weak population structure and recent demographic expansion of the monogenean parasite Kapentagyrus spp. infecting clupeid fishes of Lake Tanganyika, East Africa.

Lake Tanganyika, East Africa, is the oldest and deepest African Great Lake and harbours one of the most diverse fish assemblages on earth. Two clupeid fishes, Limnothrissa miodon and Stolothrissa tanganicae, constitute a major part of the total fish catch, making them indispensable for local food security. Parasites have been proposed as indicators of stock structure in highly mobile pelagic hosts. We examined the monogeneans Kapentagyrus limnotrissae and Kapentagyrus tanganicanus (Dactylogyridae) infecting these clupeids to explore the parasites' lake-wide population structure and patterns of demographic history. Samples were collected at seven sites distributed across three sub-basins of the lake. Intraspecific morphological variation of the monogeneans (n = 380) was analysed using morphometrics and geomorphometrics of sclerotised structures. Genetic population structure of both parasite species (n = 246) was assessed based on a 415 bp fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Overall, we observed a lack of clear geographical morphological differentiation in both parasites along a north-south axis. This lack of geographical population structure was also reflected by a large proportion of shared haplotypes, and a pattern of seemingly unrestricted gene flow between populations. Significant morphological and genetic differentiation between some populations might reflect temporal differentiation rather than geographical isolation. Overall, the shallow population structure of both species of Kapentagyrus reflects the near-panmictic population structure of both host species as previously reported. Morphological differences related to host species identity of K. tanganicanus were consistent with incipient speciation at the genetic level. Both parasite species experienced a recent demographic expansion, which might be linked to paleohydrological events. Finally, interspecific hybridisation was found in Kapentagyrus, representing the first case in dactylogyrid monogeneans.

Sea ice loss increases genetic isolation in a high Arctic ungulate metapopulation.

Sea ice loss may have dramatic consequences for population connectivity, extinction-colonization dynamics, and even the persistence of Arctic species subject to climate change. This is of particular concern in face of additional anthropogenic stressors, such as overexploitation. In this study, we assess the population-genetic implications of diminishing sea ice cover in the endemic, high Arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) by analyzing the interactive effects of landscape barriers and reintroductions (following harvest-induced extirpations) on their metapopulation genetic structure. We genotyped 411 wild reindeer from 25 sampling sites throughout the entire subspecies' range at 19 microsatellite loci. Bayesian clustering analysis showed a genetic structure composed of eight populations, of which two were admixed. Overall population genetic differentiation was high (mean FST  = 0.21). Genetic diversity was low (allelic richness [AR] = 2.07-2.58; observed heterozygosity = 0.23-0.43) and declined toward the outer distribution range, where populations showed significant levels of inbreeding. Coalescent estimates of effective population sizes and migration rates revealed strong evolutionary source-sink dynamics with the central population as the main source. The population genetic structure was best explained by a landscape genetics model combining strong isolation by glaciers and open water, and high connectivity by dispersal across winter sea ice. However, the observed patterns of natural isolation were strongly modified by the signature of past harvest-induced extirpations, subsequent reintroductions, and recent lack of sea ice. These results suggest that past and current anthropogenic drivers of metapopulation dynamics may have interactive effects on large-scale ecological and evolutionary processes. Continued loss of sea ice as a dispersal corridor within and between island systems is expected to increase the genetic isolation of populations, and thus threaten the evolutionary potential and persistence of Arctic wildlife.

Temporally varying disruptive selection in the medium ground finch (Geospiza fortis).

Disruptive natural selection within populations exploiting different resources is considered to be a major driver of adaptive radiation and the production of biodiversity. Fitness functions, which describe the relationships between trait variation and fitness, can help to illuminate how this disruptive selection leads to population differentiation. However, a single fitness function represents only a particular selection regime over a single specified time period (often a single season or a year), and therefore might not capture longer-term dynamics. Here, we build a series of annual fitness functions that quantify the relationships between phenotype and apparent survival. These functions are based on a 9-year mark-recapture dataset of over 600 medium ground finches (Geospiza fortis) within a population bimodal for beak size. We then relate changes in the shape of these functions to climate variables. We find that disruptive selection between small and large beak morphotypes, as reported previously for 2 years, is present throughout the study period, but that the intensity of this selection varies in association with the harshness of environment. In particular, we find that disruptive selection was strongest when precipitation was high during the dry season of the previous year. Our results shed light on climatic factors associated with disruptive selection in Darwin's finches, and highlight the role of temporally varying fitness functions in modulating the extent of population differentiation.

Urbanization erodes niche segregation in Darwin's finches.

Urbanization is influencing patterns of biological evolution in ways that are only beginning to be explored. One potential effect of urbanization is in modifying ecological resource distributions that underlie niche differences and that thus promote and maintain species diversification. Few studies have assessed such modifications, or their potential evolutionary consequences, in the context of ongoing adaptive radiation. We study this effect in Darwin's finches on the Galápagos Islands, by quantifying feeding preferences and diet niche partitioning across sites with different degrees of urbanization. We found higher finch density in urban sites and that feeding preferences and diets at urban sites skew heavily toward human food items. Furthermore, we show that finches at urban sites appear to be accustomed to the presence of people, compared with birds at sites with few people. In addition, we found that human behavior via the tendency to feed birds at non-urban but tourist sites is likely an important driver of finch preferences for human foods. Site differences in diet and feeding behavior have resulted in larger niche breadth within finch species and wider niche overlap between species at the urban sites. Both factors effectively minimize niche differences that would otherwise facilitate interspecies coexistence. These findings suggest that both human behavior and ongoing urbanization in Galápagos are starting to erode ecological differences that promote and maintain adaptive radiation in Darwin's finches. Smoothing of adaptive landscapes underlying diversification represents a potentially important yet underappreciated consequence of urbanization. Overall, our findings accentuate the fragility of the initial stages of adaptive radiation in Darwin's finches and raise concerns about the fate of the Galápagos ecosystems in the face of increasing urbanization.

An exploration of the links between parasites, trophic ecology, morphology, and immunogenetics in the Lake Tanganyika cichlid radiation.

Differences in habitat and diet between species are often associated with morphological differences. Habitat and trophic adaptation have therefore been proposed as important drivers of speciation and adaptive radiation. Importantly, habitat and diet shifts likely impose changes in exposure to different parasites and infection risk. As strong selective agents influencing survival and mate choice, parasites might play an important role in host diversification. We explore this possibility for the adaptive radiation of Lake Tanganyika (LT) cichlids. We first compare metazoan macroparasites infection levels between cichlid tribes. We then describe the cichlids' genetic diversity at the major histocompatibility complex (MHC), which plays a key role in vertebrate immunity. Finally, we evaluate to what extent trophic ecology and morphology explain variation in infection levels and MHC, accounting for phylogenetic relationships. We show that different cichlid tribes in LT feature partially non-overlapping parasite communities and partially non-overlapping MHC diversity. While morphology explained 15% of the variation in mean parasite abundance, trophic ecology accounted for 16% and 22% of the MHC variation at the nucleotide and at the amino acid level, respectively. Parasitism and immunogenetic adaptation may thus add additional dimensions to the LT cichlid radiation.

First genomic study on Lake Tanganyika sprat Stolothrissa tanganicae: a lack of population structure calls for integrated management of this important fisheries target species.

BACKGROUND: Clupeid fisheries in Lake Tanganyika (East Africa) provide food for millions of people in one of the world's poorest regions. Due to climate change and overfishing, the clupeid stocks of Lake Tanganyika are declining. We investigate the population structure of the Lake Tanganyika sprat Stolothrissa tanganicae, using for the first time a genomic approach on this species. This is an important step towards knowing if the species should be managed separately or as a single stock. Population structure is important for fisheries management, yet understudied for many African freshwater species. We hypothesize that distinct stocks of S. tanganicae could be present due to the large size of the lake (isolation by distance), limnological variation (adaptive evolution), or past separation of the lake (historical subdivision). On the other hand, high mobility of the species and lack of obvious migration barriers might have resulted in a homogenous population. RESULTS: We performed a population genetic study on wild-caught S. tanganicae through a combination of mitochondrial genotyping (96 individuals) and RAD sequencing (83 individuals). Samples were collected at five locations along a north-south axis of Lake Tanganyika. The mtDNA data had low global FST and, visualised in a haplotype network, did not show phylogeographic structure. RAD sequencing yielded a panel of 3504 SNPs, with low genetic differentiation (FST = 0.0054; 95% CI: 0.0046-0.0066). PCoA, fineRADstructure and global FST suggest a near-panmictic population. Two distinct groups are apparent in these analyses (FST = 0.1338 95% CI: 0.1239,0.1445), which do not correspond to sampling locations. Autocorrelation analysis showed a slight increase in genetic difference with increasing distance. No outlier loci were detected in the RADseq data. CONCLUSION: Our results show at most very weak geographical structuring of the stock and do not provide evidence for genetic adaptation to historical or environmental differences over a north-south axis. Based on these results, we advise to manage the stock as one population, integrating one management strategy over the four riparian countries. These results are a first comprehensive study on the population structure of these important fisheries target species, and can guide fisheries management.