It's not all black and white: investigating colour polymorphism in manta rays across Indo-Pacific populations

Distribution of the reef manta ray Mobula alfredi and the oceanic manta ray Mobula birostris in the Philippines: a collaborative effort for conservation

Little is known about manta ray population size, structure and connectivity in the Philippines. In collaboration with dive operators, non-governmental organizations and authorities, sightings of manta rays were collated into a single national database. Using in-water photographs and videos gathered through citizen science and dedicated research efforts, this study compiled sightings between 2004 and 2020, showing 22 separate sites throughout the archipelago with manta rays present. A total of 392 individual reef manta rays (Mobula alfredi) and 107 oceanic manta rays (Mobula birostris) were identified from the collected footage. Four specific sites in the provinces of Masbate and Palawan together hosted 89% of all identified individuals and accounted for 95% of sightings, highlighting these areas are key aggregation sites. This study also reports the movements of M. birostris within the Philippines, based on photo-identification of three individuals moving 150 km between Cebu and Masbate. Despite the growing number of recreational divers in Daanbantayan and San Jacinto, an 80% decline in M. birostris sightings was observed at these sites. To ensure effective future conservation, it is recommended that efforts focus on the identification and protection of manta ray hotspots and migratory corridors, the creation of a sustainable tourism framework and, most important, the implementation of mitigation strategies to reduce fisheries interactions.

Gudgeon fish with and without genetically determined countershading coexist in heterogeneous littoral environments of an ancient lake

Countershading, characterized by a darker dorsal surface and lighter ventral surface, is common among many animals. This dorsoventral pigment polarity is often thought to be adaptive coloration for camouflage. By contrast, noncountershaded (melanistic) morphs often occur within a species due to genetic color polymorphism in terrestrial animals. However, the polymorphism with either countershaded or melanistic morphs is poorly known in wild aquatic animals. This study explored the genetic nature of diverged color morphs of a lineage of gudgeon fish (genus Sarcocheilichthys) in the ancient Lake Biwa and propose this system as a novel model for testing hypotheses of functional aspects of countershading and its loss in aquatic environments. This system harbors two color morphs that have been treated taxonomically as separate species; Sarcocheilichthys variegatus microoculus which occurs throughout the littoral zone and Sarcocheilichthys biwaensis which occurs in and around rocky areas. First, we confirmed that the divergence of dorsoventral color patterns between the two morphs is under strict genetic control at the levels of chromatophore distribution and melanin-related gene expression under common garden rearing. The former morph displayed sharp countershading coloration, whereas the latter morph exhibited a strong tendency toward its loss. The crossing results indicated that this divergence was likely controlled by a single locus in a two-allele Mendelian inheritance pattern. Furthermore, our population genomic and genome-wide association study analyses detected no genome-wide divergence between the two morphs, except for one region near a locus that may be associated with the color divergence. Thus, these morphs are either in a state of intraspecific color polymorphism or two incipient species. Evolutionary forces underlying this polymorphism appear to be associated with heterogeneous littoral environments in this lake. Future ecological genomic research will provide insight into adaptive functions of this widespread coloration, including the eco-evolutionary drivers of its loss, in the aquatic world.

High frequency of an otherwise rare phenotype in a small and isolated tiger population

Most endangered species exist today in small populations, many of which are isolated. Evolution in such populations is largely governed by genetic drift. Empirical evidence for drift affecting striking phenotypes based on substantial genetic data are rare. Approximately 37% of tigers (Panthera tigris) in the Similipal Tiger Reserve (in eastern India) are pseudomelanistic, characterized by wide, merged stripes. Camera trap data across the tiger range revealed the presence of pseudomelanistic tigers only in Similipal. We investigated the genetic basis for pseudomelanism and examined the role of drift in driving this phenotype's frequency. Whole-genome data and pedigree-based association analyses from captive tigers revealed that pseudomelanism cosegregates with a conserved and functionally important coding alteration in Transmembrane Aminopeptidase Q (Taqpep), a gene responsible for similar traits in other felid species. Noninvasive sampling of tigers revealed a high frequency of the Taqpep p.H454Y mutation in Similipal (12 individuals, allele frequency = 0.58) and absence from all other tiger populations (395 individuals). Population genetic analyses confirmed few (minimal number) tigers in Similipal, and its genetic isolation, with poor geneflow. Pairwise FST (0.33) at the mutation site was high but not an outlier. Similipal tigers had low diversity at 81 single nucleotide polymorphisms (mean heterozygosity = 0.28, SD = 0.27). Simulations were consistent with founding events and drift as possible drivers for the observed stark difference of allele frequency. Our results highlight the role of stochastic processes in the evolution of rare phenotypes. We highlight an unusual evolutionary trajectory in a small and isolated population of an endangered species.

Genome-wide SNPs detect no evidence of genetic population structure for reef manta rays (Mobula alfredi) in southern Mozambique

Little is known about the extent of genetic connectivity along continuous coastlines in manta rays, or whether site visitation is influenced by relatedness. Such information is pertinent to defining population boundaries and understanding localized dispersal patterns and behaviour. Here, we use 3057 genome-wide single-nucleotide polymorphisms (SNPs) to evaluate population genetic structure and assess the levels of relatedness at aggregation sites of reef manta rays (Mobula alfredi) in southern Mozambique (n = 114). Contrary to indications of limited dispersal along the southern Mozambican coastline inferred from photo-identification and telemetry studies, our results show no evidence of population structure (non-significant F_ST < 0.001) for M. alfredi along this coast. We also found no evidence that individuals sampled at the same site were more related than expected by chance for males, females or across both sexes, suggesting that kinship may not influence visitation patterns at these sites. We estimated the effective population size (N_e) of this population to be 375 (95% CI = 369-380). Comparison to a distant eastern Indian Ocean site (Western Australia, n = 15) revealed strong genetic differentiation between Mozambique and Western Australia (F_ST = 0.377), identifying the Indian Ocean basin as a barrier to dispersal. Our findings show that genetic connectivity in M. alfredi extends for several hundred kilometres along continuous coastlines. We therefore recommend that the population in Mozambique be considered a discrete management unit, and future conservation plans should prioritize integrated strategies along the entire southern coastline.

It's not all black and white: investigating colour polymorphism in manta rays across Indo-Pacific populations

Intraspecific colour polymorphisms have been the focus of numerous studies, yet processes affecting melanism in the marine environment remain poorly understood. Arguably, the most prominent example of melanism in marine species occurs in manta rays (Mobula birostris and Mobula alfredi). Here, we use long-term photo identification catalogues to document the frequency variation of melanism across Indo-Pacific manta ray populations and test for evidence of selection by predation acting on colour morph variants. We use mark-recapture modelling to compare survivorship of typical and melanistic colour morphs in three M. alfredi populations and assess the relationship between frequency variation and geographical distance. While there were large differences in melanism frequencies among populations of both species (0-40.70%), apparent survival estimates revealed no difference in survivorship between colour morphs. We found a significant association between phenotypic and geographical distance in M. birostris, but not in M. alfredi. Our results suggest that melanism is not under selection by predation in the tested M. alfredi populations, and that frequency differences across populations of both species are a consequence of neutral genetic processes. As genetic colour polymorphisms are often subjected to complex selection mechanisms, our findings only begin to elucidate the underlying evolutionary processes responsible for the maintenance and frequency variation of melanism in manta ray populations.