Population genetic structure of Indian oil sardine, Sardinella longiceps along Indian coast

Multiple merger coalescent inference of effective population size

Variation in a sample of molecular sequence data informs about the past evolutionary history of the sample's population. Traditionally, Bayesian modelling coupled with the standard coalescent is used to infer the sample's bifurcating genealogy and demographic and evolutionary parameters such as effective population size and mutation rates. However, there are many situations where binary coalescent models do not accurately reflect the true underlying ancestral processes. Here, we propose a Bayesian non-parametric method for inferring effective population size trajectories from a multifurcating genealogy under the [Formula: see text]-coalescent. In particular, we jointly estimate the effective population size and the model parameter for the Beta-coalescent model, a special type of [Formula: see text]-coalescent. Finally, we test our methods on simulations and apply them to study various viral dynamics as well as Japanese sardine population size changes over time. The code and vignettes can be found in the phylodyn package.This article is part of the theme issue '"A mathematical theory of evolution": phylogenetic models dating back 100 years'.

Genetic analyses reveal panmixia in Indian waters and population subdivision across Indian Ocean and Indo-Malay Archipelago for Decapterus russelli

The Indian Scad, Decapterus russelli is an important pelagic carangid widely distributed throughout the Indian Ocean and the Indo-West Pacific. Despite being widely distributed in the Indian Ocean, the information regarding genetic structuring and diversity has been lacking compared to its Indo Malay Archipelago counterparts. The present study was conducted to investigate the genetic stock structure of D. russelli based on mitochondrial (Cyt b) and nuclear (DrAldoB1) markers along Indian waters. The results indicated the presence of a single panmictic stock across the Indian Ocean region. High haplotype diversity associated with low nucleotide diversity suggested a population bottleneck followed by rapid population growth. Phylogenetic analysis revealed the absence of geographical clustering of lineages with the most common haplotype distributed globally. The pelagic life style, migratory capabilities, and larval dispersal may be the contributing factors to the observed spatial homogeneity of D. russelli. However, significant genetic differentiation was observed between the populations from Indian Ocean and Indo-Malay Archipelago. Hierarchical molecular variance analysis (AMOVA), pairwise FST comparisons and SAMOVA showed existence of two distinct genetic stocks of D. russelli in the Indian Ocean and IMA. The observed interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the Indo-Malay Archipelago and the Indian Ocean regions suggest the influence of historic isolation, ocean surface currents and biotic and abiotic features of the ocean. Also, there was a significant relationship between genetic distance and geographical distance between population pairs in a manner consistent with isolation-by-distance. These resulted in the evolution of a phylogeographic break for this species between these regions. The findings of these results suggest that D. russelli from the Indian Ocean shall be managed in its entire area of distribution as a single stock. Further, the Indian Ocean and Indo-Malayan stocks can be managed separately.

The sequence and de novo assembly of the genome of the Indian oil sardine, Sardinella longiceps

The Indian oil sardine, Sardinella longiceps, is a widely distributed and commercially important small pelagic fish of the Northern Indian Ocean. The genome of the Indian oil sardine has been characterized using Illumina and Nanopore platforms. The assembly is 1.077 Gb (31.86 Mb Scaffold N50) in size with a repeat content of 23.24%. The BUSCO (Benchmarking Universal Single Copy Orthologues) completeness of the assembly is 93.5% when compared with Actinopterygii (ray finned fishes) data set. A total of 46316 protein coding genes were predicted. Sardinella longiceps is nutritionally rich with high levels of omega-3 polyunsaturated fatty acids (PUFA). The core genes for omega-3 PUFA biosynthesis, such as Elovl 1a and 1b,Elovl 2, Elovl 4a and 4b,Elovl 8a and 8b,and Fads 2, were observed in Sardinella longiceps. The presence of these genes may indicate the PUFA biosynthetic capability of Indian oil sardine, which needs to be confirmed functionally.

Genomic investigations provide insights into the mechanisms of resilience to heterogeneous habitats of the Indian Ocean in a pelagic fish

The adaptive genetic variation in response to heterogeneous habitats of the Indian Ocean was investigated in the Indian oil sardine using ddRAD sequencing to understand the subpopulation structure, stock complexity, mechanisms of resilience, and vulnerability in the face of climate change. Samples were collected from different ecoregions of the Indian ocean and ddRAD sequencing was carried out. Population genetic analyses revealed that samples from the Gulf of Oman significantly diverged from other Indian Ocean samples. SNP allele-environment correlation revealed the presence of candidate loci correlated with the environmental variables like annual sea surface temperature, chlorophyll-a, and dissolved oxygen concentration which might represent genomic regions allegedly diverging as a result of local adaptation. Larval dispersal modelling along the southwest coast of India indicated a high dispersal rate. The two major subpopulations (Gulf of Oman and Indian) need to be managed regionally to ensure the preservation of genetic diversity, which is crucial for climatic resilience.

Population Genetic Structure of Ornate Threadfin Bream (Nemipterus hexodon) in Thailand

Ornate threadfin bream (Nemipterus hexodon) is an economically important fishery species in Southeast Asia. In Thailand, N. hexodon decreased dramatically due to overexploitation for commercial purposes. To construct an effective sustainable management plan, genetic information is necessary. Thus, in our study, the population genetic structure and demographic history of N. hexodon were investigated using 419 bp of the mitochondrial DNA sequence in cytochrome oxidase subunit I gene (mtDNA COI). A total of 142 samples was collected from nine localities in the Gulf of Thailand (Chonburi, Samut Songkhram, Surat Thani, Nakhon Si Thammarat, Songkhla), and the Andaman Sea (Satun, Trang, Krabi, Phang Nga). Fourteen polymorphic sites defined 18 haplotypes, revealing a high haplotype diversity and low nucleotide diversity among nine localities. The analysis of molecular variance (AMOVA) analysis, pairwise F _ST , and minimum spanning network result revealed that the genetic structure of N. hexodon was separated into two populations: the Gulf of Thailand and the Andaman Sea population. The genetic structure of N. hexodon can be explained by a disruption of gene flow from the geographic barrier and the Pleistocene isolation of the marine basin hypothesis. Neutrality tests, Bayesian skyline analysis, mismatch distribution, and the estimated values of population growth suggested that N. hexodon had experienced a population expansion. The genetic information would certainly help us gain insight into the population genetic structure of N. hexodon living on the coast of Thailand.

Signals of selection in the mitogenome provide insights into adaptation mechanisms in heterogeneous habitats in a widely distributed pelagic fish

Oceans are vast, dynamic, and complex ecosystems characterized by fluctuations in environmental parameters like sea surface temperature (SST), salinity, oxygen availability, and productivity. Environmental variability acts as the driver of organismal evolution and speciation as organisms strive to cope with the challenges. We investigated the evolutionary consequences of heterogeneous environmental conditions on the mitogenome of a widely distributed small pelagic fish of Indian ocean, Indian oil sardine, Sardinella longiceps. Sardines were collected from different eco-regions of the Indian Ocean and selection patterns analyzed in coding and non-coding regions. Signals of diversifying selection were observed in key functional regions involved in OXPHOS indicating OXPHOS gene regulation as the critical factor to meet enhanced energetic demands. A characteristic control region with 38–40 bp tandem repeat units under strong selective pressure as evidenced by sequence conservation and low free energy values was also observed. These changes were prevalent in fishes from the South Eastern Arabian Sea (SEAS) followed by the Northern Arabian Sea (NAS) and rare in Bay of Bengal (BoB) populations. Fishes belonging to SEAS exhibited accelerated substitution rate mainly due to the selective pressures to survive in a highly variable oceanic environment characterized by seasonal hypoxia, variable SST, and food availability.

Contemporary and historic patterns of intraspecific diversity in Indian anchovy, Stolephorus indicus, from Indian peninsular waters

We analyzed intraspecific diversity of Indian anchovy, Stolephorus indicus, a commercially and ecologically important species, using mitochondrial DNA markers so as to derive insights into population structuring and historical demography. Analyses were carried out on 128 and 138 individuals collected from 5 locations along the range of distribution using mitochondrial ATPase (843 bp) and COI (663 bp) sequences respectively. Significant connectivity and gene flow was detected among fishes collected from all the geographic locations as indicated by lack of structuring in Bayesian clustering analysis along with insignificant Φ_ST values. Oceanographic features of the Bay of Bengal, Arabian Sea and Andaman Sea may be favorable for the dispersal of anchovy larvae and subsequent gene flow. Historical demographic analyses indicated a demographic and spatial expansion taken place approximately during 125,000 years before present, the Pleistocene epoch. Indian Ocean witnessed emergence of upwelling events and consequent increase in productivity during the Pleistocene epoch causing a demographic and spatial expansion of anchovies. Management measures for this species should be devised considering it as a single stock along its entire range of distribution.

Population Structure of Blue Marlin, Makaira nigricans, in the Pacific and Eastern Indian Oceans

Hui Chen, Chia-Hao Chang, Chi-Lu Sun, Kwang-Tsao Shao, Su-Zan Yeh, and Gerard DiNardo (2016) Blue marlin Makaira nigricans is economically important for fisheries worldwide. However, overfishing has substantially reduced the stock size. Better knowledge of blue marlin population genetics will help improve management and conservation. Previous genetic studies concluded that the Pacific blue marlin should be considered a single stock. This study investigated the population genetic structure of blue marlin inhabiting the Pacific and eastern Indian oceans based on mtDNA cytochrome b (cyt b) and control region (CR) sequence variation. We collected tissue samples (n = 183) from three Pacific and one Indian Ocean, and determined the sequences of 1140 bp of cyt b and 905 bp of CR. Phylogenetic analysis revealed that blue marlin contain two clades, the Atlantic clade and the ubiquitous clade, and that all the eastern Indian and Pacific individuals collected for this study belonged to the ubiquitous clade. All eastern Indian and Pacific blue marlin possess extremely high haplotype diversity (h) and low nucleotide diversity (π). The results of pairwise ΦST, hierarchical analysis of molecular variance (AMOVA) and spatial analysis of molecular variance (SAMOVA) all support that there is no population differentiation among eastern Indian and Pacific blue marlin. Neutrality tests and pairwise mismatch distribution analysis both indicated that eastern Indian and Pacific blue marlin have undergone a rapid population expansion on the order of 0.30 to 0.65 million years ago (mya). This study demonstrates that blue marlin in the Pacific and eastern Indian oceans constitute a single stock. International cooperation will be required to preserve blue marlin as a resource; moreover, the high genetic variation of blue marlin in this region suggests that unique haplotypes in the population are sensitive to high harvesting levels and could disappear.