Mitogenomic recognition of incognito lineages in the mud spiny lobster Panulirus polyphagus (Herbst, 1793): A tale of unique genetic structuring and diversification.

This study provides the first documentation of three deep conspecific lineages within Panulirus polyphagus in the Indian Ocean, bridging the gap in genetic research. Comparative mitogenomics between lineages (L) at both species and family levels, evolutionary relationships and heterogeneity of sequence divergence within Decapoda, and divergence time estimation were performed. The characterized mitogenomes ranged from 15,685-15,705 bp in size and exhibited a typical pancrustacean pattern. Among the three lineages, L1 predominated the Bay of Bengal, L2 the Arabian Sea, and L2.a, a less common lineage genetically closer to L2, was restricted to the latter region. A minor lineage L1.a, was observed in the Coral Triangle area. All PCGs displayed evidence of purifying selection across species and family levels. The largest genetic distance (K2P) between lineages was 9 %, notably between L1.a and L2.a. The phylogenetic tree subdivided the Achelates into Palinuridae and Scyllaridae, and the topology demonstrated a distinct pattern of lineage diversification within P. polyphagus. AliGROOVE analysis revealed no discernible divergence in Decapoda. The diversification of P. polyphagus appears to have occurred during Miocene, with further diversification in Pliocene. Furthermore, genetic stocks and population connectivity recognized here will provide valuable insight for spatial management planning of this dwindling resource.

The chromosome level genome assembly of the Asian green mussel, Perna viridis.

The Asian green mussel, Perna viridis is an important aquaculture species in the family Mytilidae contributing substantially to molluscan aquaculture. We generated a high-quality chromosome level assembly of this species by combining PacBio single molecule sequencing technique (SMRT), Illumina paired-end sequencing, high-throughput chromosome conformation capture technique (Hi-C) and Bionano mapping. The final assembly resulted in a genome of 723.49 Mb in size with a scaffold N50 of 49.74 Mb with 99% anchored into 15 chromosomes. A total of 49654 protein-coding genes were predicted from the genome. The presence of 634 genes associated with the cancer pathway and 408 genes associated with viral carcinogenesis indicates the potential of this species to be used as a model for cancer studies. The chromosome-level assembly of this species is also a valuable resource for further genomic selection and selective breeding for improving economically important aquaculture traits and augmenting aquaculture productivity.

Heat Stress Effects on Physiological and Milk Yield Traits of Lactating Holstein Friesian Crossbreds Reared in Tanga Region, Tanzania.

Global warming caused by climate change is a challenge for dairy farming, especially in sub-Saharan countries. Under high temperatures and relative humidity, lactating dairy cows suffer from heat stress. The objective of this study was to investigate the effects and relationship of heat stress (HS) measured by the temperature-humidity index (THI) regarding the physiological parameters and milk yield and composition of lactating Holstein Friesian crossbred dairy cows reared in the humid coastal region of Tanzania. A total of 29 lactating Holstein Friesian x Zebu crossbred dairy cows with 50% (HF50) and 75% (HF75) Holstein Friesian gene levels in the second and third months of lactation were used. The breed composition of Holstein Friesians was determined based on the animal recording system used at the Tanzania Livestock Research Institute (TALIRI), Tanga. The data collected included the daily temperature, relative humidity, daily milk yield, and physiological parameters (core body temperature, rectal temperature, respiratory rate, and panting score). THI was calculated using the equation of the National Research Council. The THI values were categorized into three classes, i.e., low THI (76-78), moderate THI (79-81), and high THI (82-84). The effects of THI on the physiological parameters and milk yield and composition were assessed. The effects of the genotype, the parity, the lactation month, and the interaction of these parameters with THI on the milk yield, milk composition, and physiological parameters were also investigated. The results show that THI and its interaction with genotypes, parity, and the lactation month had a highly significant effect on all parameters. THI influenced (p ˂ 0.05) the average daily milk yield and milk fat %, protein %, lactose %, and solids-not-fat %. As the THI increased from moderate to high levels, the average daily milk yield declined from 3.49 ± 0.04 to 3.43 ± 0.05 L/day, while the fat % increased from 2.66 ± 0.05% to 3.04 ± 0.06% and the protein decreased from 3.15 ± 0.02% to 3.13 ± 0.03%. No decline in lactose % was observed, while the solid-not-fat % declined from 8.56 ± 0.08% to 8.55 ± 0.10% as the THI values increased from moderate to high. Also, the THI influenced physiological parameters (p ˂ 0.05). The core body temperature (CBT), rectal temperature (RT), respiratory rate (RR) and panting score (PS) increased from 35.60 ± 0.01 to 36.00 ± 0.01 °C, 38.03 ± 0.02 to 38.30 ± 0.02 °C, 62.53 ± 0.29 to 72.35 ± 0.28 breaths/min, and 1.35 ± 0.01 to 1.47 ± 0.09, respectively, as the THI increased from low to high. The THI showed a weak positive correlation with the average daily milk yield and fat percentage, whereas the protein, lactose, and solids-not-fat percentages showed negative relationships with THI (p ≤ 0.05). CBT, RT, RR, and PS showed positive relationships (p ≤ 0.05) with THI. These negative relationships indicate that there is an antagonistic correlation between sensitivity to HS and the level of production. It is concluded that the THI, the genotype, the parity, and the lactation month, along with their interactions with THI, significantly influenced the milk yield, milk composition, and physiological parameters of lactating Holstein Friesian dairy crosses at THI thresholds ranging from 77 to 84.

Tidings from the Tides-De novo transcriptome assembly of the endemic estuarine bivalve Villorita cyprinoides.

The Indian black clam Villorita cyprinoides Gray, 1825, is an economically valuable estuarine bivalve that faces challenges from multiple stressors and anthropogenic pressures. However, limited genomic resources have hindered molecular investigations into the impact of these stressors on clam populations. Here, we have generated the first transcriptomic reference datasets for V. cyprinoides to address this knowledge gap. A total of 25,040,592 and 22,486,217 million Illumina paired-end reads generated from two individuals were assembled using Trinity and rnaSPAdes. From the 47,607 transcripts identified as Coding Domain Sequences, 37,487 returned positive BLAST hits against six different databases. Additionally, a total of 14,063 Single Sequence Repeats were identified using GMATA. This study significantly enhances the genetic understanding of V. cyprinoides, a potential candidate for aquaculture that supports the livelihoods of many people dependent on small-scale fisheries. The data generated provides insights into broader genealogical connections within the family Cyrenidae through comparative transcriptomics. Furthermore, this transcriptional profile serves as baseline data for future studies in toxicological and conservation genetics.

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.

Heat stress effects on milk yield traits and metabolites and mitigation strategies for dairy cattle breeds reared in tropical and sub-tropical countries.

Heat stress is an important problem for dairy industry in many parts of the world owing to its adverse effects on productivity and profitability. Heat stress in dairy cattle is caused by an increase in core body temperature, which affects the fat production in the mammary gland. It reduces milk yield, dry matter intake, and alters the milk composition, such as fat, protein, lactose, and solids-not-fats percentages among others. Understanding the biological mechanisms of climatic adaptation, identifying and exploring signatures of selection, genomic diversity and identification of candidate genes for heat tolerance within indicine and taurine dairy breeds is an important progression toward breeding better dairy cattle adapted to changing climatic conditions of the tropics. Identifying breeds that are heat tolerant and their use in genetic improvement programs is crucial for improving dairy cattle productivity and profitability in the tropics. Genetic improvement for heat tolerance requires availability of genetic parameters, but these genetic parameters are currently missing in many tropical countries. In this article, we reviewed the HS effects on dairy cattle with regard to (1) physiological parameters; (2) milk yield and composition traits; and (3) milk and blood metabolites for dairy cattle reared in tropical countries. In addition, mitigation strategies such as physical modification of environment, nutritional, and genetic development of heat tolerant dairy cattle to prevent the adverse effects of HS on dairy cattle are discussed. In tropical climates, a more and cost-effective strategy to overcome HS effects is to genetically select more adaptable and heat tolerant breeds, use of crossbred animals for milk production, i.e., crosses between indicine breeds such as Gir, white fulani, N'Dama, Sahiwal or Boran to taurine breeds such as Holstein-Friesian, Jersey or Brown Swiss. The results of this review will contribute to policy formulations with regard to strategies for mitigating the effects of HS on dairy cattle in tropical countries.

Genes and models for estimating genetic parameters for heat tolerance in dairy cattle.

Dairy cattle are highly susceptible to heat stress. Heat stress causes a decline in milk yield, reduced dry matter intake, reduced fertility rates, and alteration of physiological traits (e.g., respiration rate, rectal temperature, heart rates, pulse rates, panting score, sweating rates, and drooling score) and other biomarkers (oxidative heat stress biomarkers and stress response genes). Considering the significant effect of global warming on dairy cattle farming, coupled with the aim to reduce income losses of dairy cattle farmers and improve production under hot environment, there is a need to develop heat tolerant dairy cattle that can grow, reproduce and produce milk reasonably under the changing global climate and increasing temperature. The identification of heat tolerant dairy cattle is an alternative strategy for breeding thermotolerant dairy cattle for changing climatic conditions. This review synthesizes information pertaining to quantitative genetic models that have been applied to estimate genetic parameters for heat tolerance and relationship between measures of heat tolerance and production and reproductive performance traits in dairy cattle. Moreover, the review identified the genes that have been shown to influence heat tolerance in dairy cattle and evaluated the possibility of using them in genomic selection programmes. Combining genomics information with environmental, physiological, and production parameters information is a crucial strategy to understand the mechanisms of heat tolerance while breeding heat tolerant dairy cattle adapted to future climatic conditions. Thus, selection for thermotolerant dairy cattle is feasible.

Full-length transcriptome from different life stages of cobia (Rachycentron canadum, Rachycentridae).

Cobia (Rachycentron canadum, Rachycentridae) is one of the prospective species for mariculture. The transcriptome-based study on cobia was hampered by an inadequate reference genome and a lack of full-length cDNAs. We used a long-read based sequencing technology (PacBio Sequel II Iso-Seq3 SMRT) to obtain complete transcriptome sequences from larvae, juveniles, and various tissues of adult cobia, and a single SMRTcell generated 99 gigabytes of data and 51,205,946,694 bases. A total of 8609435, 7441673 and 9140164 subreads were generated from the larval, juvenile, and adult sample pools, with mean sub-read lengths of 2109.9, 1988.2 and 1996.2 bp, respectively. All samples were combined to increase transcript recovery and clustered into 35661 high-quality reads. This is the first report on a full-length transcriptome from R. canadum. Our results illustrate a significant increase in the identified amount of cobia LncRNAs and alternatively spliced transcripts, which will help improve genome annotation. Furthermore, this information will be beneficial for nutrigenomics and functional studies on cobia and other commercially important mariculture species.

Comparative mitogenomics of Clupeoid fish provides insights into the adaptive evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes and codon usage in the heterogeneous habitats.

Clupeoid fish can be considered excellent candidates to understand the role of mitochondrial DNA in adaptive evolution, as they have colonized different habitats (marine, brackish, freshwater, tropical and temperate regions) over millions of years. Here, we investigate patterns of tRNA location, codon usage bias, and lineage-specific diversifying selection signals to provide novel insights into how evolutionary improvements of mitochondrial metabolic efficiency have allowed clupeids to adapt to different habitats. Based on whole mitogenome data of 70 Clupeoids with a global distribution we find that purifying selection was the dominant force acting and that the mutational deamination pressure in mtDNA was stronger than the codon/amino acid constraints. The codon usage pattern appears evolved to achieve high translational efficiency (codon/amino acid-related constraints), as indicated by the complementarity of most codons to the GT-saturated tRNA anticodon sites (retained by deamination-induced pressure) and usage of the codons of the tRNA genes situated near to the control region (fixed by deamination pressure) where transcription efficiency was high. The observed shift in codon preference patterns between marine and euryhaline/freshwater Clupeoids indicates possible selection for improved translational efficiency in mitochondrial genes while adapting to low-salinity habitats. This mitogenomic plasticity and enhanced efficiency of the metabolic machinery may have contributed to the evolutionary success and abundance of Clupeoid fish.

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.

The signals of selective constraints on the mitochondrial non-coding control region: insights from comparative mitogenomics of Clupeoid fishes.

The vertebrate mitochondrial genome is characterized by an exceptional organization evolving towards a reduced size. However, the persistence of a non-coding and highly variable control region is against this evolutionary trend that is explained by the presence of conserved sequence motifs or binding sites for nuclear-organized proteins that regulate mtDNA maintenance and expression. We performed a comparative mitogenomic investigation of the non-coding control region to understand its evolutionary patterns in Clupeoid fishes which are widely distributed across oceans of the world, exhibiting exemplary evolutionary potential. We confirmed the ability of sequence flanking the conserved sequence motifs in the control region to form stable secondary structures. The existence of evolutionarily conserved secondary structures without primary structure conservation suggested the action of selective constraints towards maintaining the secondary structure. The functional secondary structure is maintained by retaining the frequency of discontinuous AT and TG repeats along with compensatory base substitutions in the stem forming regions which can be considered as a selective constraint. The nucleotide polymorphism along the flanking regions of conserved sequence motifs can be explained as errors during the enzymatic replication of secondary structure-forming repeat elements. The evidence for selective constraints on secondary structures emphasizes the role of the control region in mitogenome function. Maintenance of high frequency of discontinuous repeats can be proposed as a model of adaptive evolution against the mutations that break the secondary structure involved in the efficient regulation of mtDNA functions substantiating the efficient functioning of the control region even in a high nucleotide polymorphism environment.

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.

Characterization of the Whole Mitogenome of Largehead Hairtail Trichiurus lepturus (Trichiuridae): Insights into Special Characteristics.

Mitochondrial genomes are valuable for making inferences regarding population genetics, phylogenetics and recently in understanding adaptive variation in OXPHOS genes. Ribbon fish, Trichiurus lepturus is distributed in the Indian Ocean and contributes to a major share of pelagic fish landings. Considering the importance, we characterized the complete mitogenome of large head hairtail T. lepturus from Indian waters by PCR amplification and sequencing. The entire mitogenome is 16,840 bp in length, with 13 protein coding genes, 2 subunit ribosomal RNA genes, and 1 control region similar to other vertebrates. Genome contains only 21 tRNAs (tRNAPro is absent) as against 22 tRNAs reported in other vertebrates. The putative origin of L-strand replication (OL) was located in a cluster of five tRNA genes (WANCY region) and may not be functional due to the variation of conserved motif (5'-GCCGG-3') and loop. These two unusual features are common for Trichiurus mitogenomes. Phylogenetic analysis supported the monophyly of Trichiuridae. This study will provide information for further studies on taxonomy, phylogeny and patterns of adaptive evolution in the mitogenome of Trichiurus species.

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.

Isolation and characterization of Aquaporin 1 (AQP1), sodium/potassium-transporting ATPase subunit alpha-1 (Na/K-ATPase α1), Heat Shock Protein 90 (HSP90), Heat Shock Cognate 71 (HSC71), Osmotic Stress Transcription Factor 1 (OSTF1) and Transcription Factor II B (TFIIB) genes from a euryhaline fish, Etroplus suratensis.

The present study reports the complete sequences of Aquaporin 1 (AQP1) gene and partial sequences of genes, Sodium/potassium-transporting ATPase subunit alpha-1 (Na/K-ATPase α1 subunit), Osmotic Stress Transcription Factor 1 (OSTF1), Transcription Factor II B (TFIIB), Heat Shock Cognate 71 (HSC71) and Heat Shock Protein 90 (HSP90) obtained from mRNA and genomic DNA of Etroplus suratensis. They are candidate genes involved in stress responses of fishes. AQP1 gene was 2163 bp long. Its mRNA sequence has 55 bp 5' UTR, 783 bp open reading frame (ORF), 119 bp 3' UTR, three intronic regions and 90% identity with AQP1 of Oreochromis niloticus. The partial Na/K-ATPase α1subunit gene obtained 5998 bp length with an ORF of 2213 bp and 12 intronic regions. The partial OSTF1, TF IIB, HSC71 and HSP90 mRNA sequences obtained were 1473 bp, 587 bp, 1708 bp and 151 bp in length respectively. All the genes showed a high sequence similarity with respective genes reported from fishes. Comparison of AQP1 and Na/K-ATPase α1 genomic DNA sequence of E. suratensis collected from different water system showed two type of AQP1 with one synonymous mutation in exon-1 and higher sequence difference in intronic regions (including addition, deletion, transition and transversion mutations) with few synonymous and non-synonymous mutations in the exons of Na/K-ATPase α1. The sequence information of these major candidate genes involved in stress responses will help in further studies on population genetics, adaptive variations and genetic improvement programs of this cichlid species having aquaculture, ornamental and evolutionary importance.

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

Genetic stock structure and historical demography of Indian oil sardine, a commercially and ecologically important small pelagic fish, was studied using mitochondrial control region and Cytochrome C Oxidase I (COI) sequences. A 758 bp portion of the control region in 287 individuals and a 576 bp portion of the COI gene in 291 individuals from 10 locations along the distribution range were amplified resulting in 236 and 84 haplotypes, respectively. The high haplotype and low nucleotide diversity values (0.99 and 0.19 for control region and 0.85 and 0.004 for COI, respectively) are characteristic of populations having undergone a demographic expansion. Genetic differentiation, ΦST, was low and insignificant between populations using both control region and COI gene markers. Mismatch analysis showed a recent demographic and spatial expansion coinciding with the late Pleistocene epoch. Mantel tests revealed the lack of isolation by distance which is attributable either to high levels of migration overriding the effect of genetic drift or to insufficient time for accomplishing a balance between migration and drift after a recent range expansion.

Tick burden and prevalence of Theileria parva infection in Tarime zebu cattle in the lake zone of Tanzania.

This study was carried out to assess the distribution, abundance of different tick genera and prevalence of Theileria parva infection in Tarime zebu cattle kept in selected wards of Serengeti and Tarime districts in Mara region. Adult ticks were identified and counted from half body parts of 360 animals which were extensively managed in communal land with natural pastures. Concurrently, blood samples were collected and thereafter DNA extracted and a nested polymerase chain reaction (nPCR) was done using primers specific for p104 gene to detect the presence of T. parva DNA. Ticks were identified into four groups: Amblyomma genus, Boophilus sub-genus of Rhipicephalus genus, other species of Rhipicephalus, and Hyalomma genus. Rhipicephalus genus accounted for 71.8 % of the total ticks, whereas Amblyomma, Boophilus sub-genus of Rhipicephalus genus and Hyalomma constituted 14.1, 14.0 and 0.1 %, respectively. There were more animals (p < 0.05) infested with ticks in Tarime district (96.1 %) than in Serengeti (61.7 %). The average counts of ticks were higher in adult animals (p < 0.05) than in young animals. The overall prevalence of T. parva was 27.7 % and was higher (p < 0.05) in Serengeti (38.3 %) than in Tarime district (16.7 %). However, all animals tested positive for T. parva did not show any clinical signs of East Coast fever (ECF), suggesting the existence of subclinical infection in Tarime zebu. These results suggest that Tarime cattle can tolerate ECF infection and are likely to serve as potential carriers of T. parva to other less-tolerant cattle breeds in mixed herds. Since Tarime cattle are preferred by most farmers with mixed herds, routine screening for T. parva is highly recommended to minimize introduction of infected cattle into an immunologically naive population.