1
|
Rahi ML, Mather PB, de Bello Cioffi M, Ezaz T, Hurwood DA. Genomic Basis of Freshwater Adaptation in the Palaemonid Prawn Genus Macrobrachium: Convergent Evolution Following Multiple Independent Colonization Events. J Mol Evol 2023; 91:976-989. [PMID: 38010517 DOI: 10.1007/s00239-023-10149-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Accepted: 11/14/2023] [Indexed: 11/29/2023]
Abstract
Adaptation to different salinity environments can enhance morphological and genomic divergence between related aquatic taxa. Species of prawns in the genus Macrobrachium naturally inhabit different osmotic niches and possess distinctive lifecycle traits associated with salinity tolerance. This study was conducted to investigate the patterns of adaptive genomic divergence during freshwater colonization in 34 Macrobrachium species collected from four continents; Australia, Asia, North and South America. Genotyping-by-sequencing (GBS) technique identified 5018 loci containing 82,636 single nucleotide polymorphisms (SNPs) that were used to reconstruct a phylogenomic tree. An additional phylogeny was reconstructed based on 43 candidate genes, previously identified as being potentially associated with freshwater adaptation. Comparison of the two phylogenetic trees revealed contrasting topologies. The GBS tree indicated multiple independent continent-specific invasions into freshwater by Macrobrachium lineages following common marine ancestry, as species with abbreviated larval development (ALD), i.e., species having a full freshwater life history, appeared reciprocally monophyletic within each continent. In contrast, the candidate gene tree showed convergent evolution for all ALD species worldwide, forming a single, well-supported clade. This latter pattern is likely the result of common evolutionary pressures selecting key mutations favored in continental freshwater habitats Results suggest that following multiple independent invasions into continental freshwaters at different evolutionary timescales, Macrobrachium taxa experienced adaptive genomic divergence, and in particular, convergence in the same genomic regions with parallel shifts in specific conserved phenotypic traits, such as evolution of larger eggs with abbreviated larval developmental.
Collapse
Affiliation(s)
- Md Lifat Rahi
- Fisheries and Marine Resource Technology Discipline, Khulna University, Khulna, Bangladesh.
| | - Peter B Mather
- Faculty of Science, Queensland University of Technology (QUT), Brisbane, QLD, 4001, Australia
| | - Marcelo de Bello Cioffi
- Department of Genetics and Evolution, Federal University of Sao Carlos, São Carlos, SP, Brazil
| | - Tariq Ezaz
- Institute for Applied Ecology (IAE), University of Canberra (UC), Canberra, ACT, 2617, Australia
| | - David A Hurwood
- Faculty of Science, Queensland University of Technology (QUT), Brisbane, QLD, 4001, Australia
| |
Collapse
|
2
|
Chowdhury MAA, Islam MR, Amin A, Mou SN, Ullah KN, Baten A, Shoyaib M, Ali AA, Chowdhury FT, Rahi ML, Khan H, Amin MA, Islam MR. Integrated transcriptome catalog of Tenualosa ilisha as a resource for gene discovery and expression profiling. Sci Data 2023; 10:214. [PMID: 37062771 PMCID: PMC10106452 DOI: 10.1038/s41597-023-02132-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 04/03/2023] [Indexed: 04/18/2023] Open
Abstract
The silver pride of Bangladesh, migratory shad, Tenualosa ilisha (Hilsa), makes the highest contribution to the total fish production of Bangladesh. Despite its noteworthy contribution, a well-annotated transcriptome data is not available. Here we report a transcriptomic catalog of Hilsa, constructed by assembling RNA-Seq reads from different tissues of the fish including brain, gill, kidney, liver, and muscle. Hilsa fish were collected from different aquatic habitats (fresh, brackish, and sea water) and the sequencing was performed in the next generation sequencing (NGS) platform. De novo assembly of the sequences obtained from 46 cDNA libraries revealed 462,085 transcript isoforms that were subsequently annotated using the Universal Protein Resource Knowledgebase (UniPortKB) as a reference. Starting from the sampling to final annotation, all the steps along with the workflow are reported here. This study will provide a significant resource for ongoing and future research on Hilsa for transcriptome based expression profiling and identification of candidate genes.
Collapse
Affiliation(s)
- Md Arko Ayon Chowdhury
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh
| | - Md Rakibul Islam
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh
| | - Al Amin
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh
| | - Sadia Noor Mou
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh
| | - Kazi Newaz Ullah
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh
- Department of Zoology, Jagannath University, Dhaka, 1100, Bangladesh
| | - Abdul Baten
- Institute of Precision Medicine and Bioinformatics, Sydney Local Health District, Royal Prince Alfred Hospital, Camperdown, Australia
| | - Mohammad Shoyaib
- Institute of Information Technology (IIT), University of Dhaka, Dhaka, 1000, Bangladesh
| | - Amin Ahsan Ali
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh
| | - Farhana Tasnim Chowdhury
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh
| | - Md Lifat Rahi
- Fisheries and Marine Resource Technology (FMRT) Discipline, Khulna University, Khulna, 9208, Bangladesh
| | - Haseena Khan
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh
| | - M Ashraful Amin
- Center for Computational and Data Sciences (CCDS), Independent University, Bangladesh (IUB), Dhaka, Bangladesh.
| | - Mohammad Riazul Islam
- Molecular Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka, 1000, Bangladesh.
| |
Collapse
|
3
|
Rahi ML, Mather PB, Ezaz T, Hurwood DA. The Molecular Basis of Freshwater Adaptation in Prawns: Insights from Comparative Transcriptomics of Three Macrobrachium Species. Genome Biol Evol 2019; 11:1002-1018. [PMID: 30840062 PMCID: PMC6450038 DOI: 10.1093/gbe/evz045] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/12/2019] [Indexed: 12/17/2022] Open
Abstract
Elucidating the molecular basis of adaptation to different environmental conditions is important because adaptive ability of a species can shape its distribution, influence speciation, and also drive a variety of evolutionary processes. For crustaceans, colonization of freshwater habitats has significantly impacted diversity, but the molecular basis of this process is poorly understood. In the current study, we examined three prawn species from the genus Macrobrachium (M. australiense, M. tolmerum, and M. novaehollandiae) to better understand the molecular basis of freshwater adaptation using a comparative transcriptomics approach. Each of these species naturally inhabit environments with different salinity levels; here, we exposed them to the same experimental salinity conditions (0‰ and 15‰), to compare expression patterns of candidate genes that previously have been shown to influence phenotypic traits associated with freshwater adaptation (e.g., genes associated with osmoregulation). Differential gene expression analysis revealed 876, 861, and 925 differentially expressed transcripts under the two salinities for M. australiense, M. tolmerum, and M. novaehollandiae, respectively. Of these, 16 were found to be unannotated novel transcripts and may be taxonomically restricted or orphan genes. Functional enrichment and molecular pathway mapping revealed 13 functionally enriched categories and 11 enriched molecular pathways that were common to the three Macrobrachium species. Pattern of selection analysis revealed 26 genes with signatures of positive selection among pairwise species comparisons. Overall, our results indicate that the same key genes and similar molecular pathways are likely to be involved with freshwater adaptation widely across this decapod group; with nonoverlapping sets of genes showing differential expression (mainly osmoregulatory genes) and signatures of positive selection (genes involved with different life history traits).
Collapse
Affiliation(s)
- Md Lifat Rahi
- Science and Engineering Faculty, School of Earth Environment and Biological Sciences (EEBS), Queensland University of Technology (QUT), Brisbane, Queensland, Australia
| | - Peter B Mather
- Science and Engineering Faculty, School of Earth Environment and Biological Sciences (EEBS), Queensland University of Technology (QUT), Brisbane, Queensland, Australia
| | - Tariq Ezaz
- Wildlife Genetics Laboratory, Institute for Applied Ecology, University of Canberra, Australian Capital Territory, Australia
| | - David A Hurwood
- Science and Engineering Faculty, School of Earth Environment and Biological Sciences (EEBS), Queensland University of Technology (QUT), Brisbane, Queensland, Australia
| |
Collapse
|
4
|
Rahi ML, Amin S, Mather PB, Hurwood DA. Candidate genes that have facilitated freshwater adaptation by palaemonid prawns in the genus Macrobrachium: identification and expression validation in a model species ( M. koombooloomba). PeerJ 2017; 5:e2977. [PMID: 28194319 PMCID: PMC5301973 DOI: 10.7717/peerj.2977] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2016] [Accepted: 01/10/2017] [Indexed: 11/20/2022] Open
Abstract
Background The endemic Australian freshwater prawn, Macrobrachium koombooloomba, provides a model for exploring genes involved with freshwater adaptation because it is one of the relatively few Macrobrachium species that can complete its entire life cycle in freshwater. Methods The present study was conducted to identify potential candidate genes that are likely to contribute to effective freshwater adaptation by M. koombooloomba using a transcriptomics approach. De novo assembly of 75 bp paired end 227,564,643 high quality Illumina raw reads from 6 different cDNA libraries revealed 125,917 contigs of variable lengths (200–18,050 bp) with an N50 value of 1597. Results In total, 31,272 (24.83%) of the assembled contigs received significant blast hits, of which 27,686 and 22,560 contigs were mapped and functionally annotated, respectively. CEGMA (Core Eukaryotic Genes Mapping Approach) based transcriptome quality assessment revealed 96.37% completeness. We identified 43 different potential genes that are likely to be involved with freshwater adaptation in M. koombooloomba. Identified candidate genes included: 25 genes for osmoregulation, five for cell volume regulation, seven for stress tolerance, three for body fluid (haemolymph) maintenance, eight for epithelial permeability and water channel regulation, nine for egg size control and three for larval development. RSEM (RNA-Seq Expectation Maximization) based abundance estimation revealed that 6,253, 5,753 and 3,795 transcripts were expressed (at TPM value ≥10) in post larvae, juveniles and adults, respectively. Differential gene expression (DGE) analysis showed that 15 genes were expressed differentially in different individuals but these genes apparently were not involved with freshwater adaptation but rather were involved in growth, development and reproductive maturation. Discussion The genomic resources developed here will be useful for better understanding the molecular basis of freshwater adaptation in Macrobrachium prawns and other crustaceans more broadly.
Collapse
Affiliation(s)
- Md Lifat Rahi
- Science and Engineering Faculty, School of Earth Environment and Biological Sciences, Queensland University of Technology (QUT) , Brisbane , Queensland , Australia
| | - Shorash Amin
- Science and Engineering Faculty, School of Biomedical Sciences, Queensland University of Technology , Brisbane , Queensland , Australia
| | - Peter B Mather
- Science and Engineering Faculty, School of Earth Environment and Biological Sciences, Queensland University of Technology (QUT) , Brisbane , Queensland , Australia
| | - David A Hurwood
- Science and Engineering Faculty, School of Earth Environment and Biological Sciences, Queensland University of Technology (QUT) , Brisbane , Queensland , Australia
| |
Collapse
|
5
|
Moshtaghi A, Rahi ML, Nguyen VT, Mather PB, Hurwood DA. A transcriptomic scan for potential candidate genes involved in osmoregulation in an obligate freshwater palaemonid prawn ( Macrobrachium australiense). PeerJ 2016; 4:e2520. [PMID: 27761323 PMCID: PMC5068373 DOI: 10.7717/peerj.2520] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Accepted: 09/02/2016] [Indexed: 01/11/2023] Open
Abstract
Background Understanding the genomic basis of osmoregulation (candidate genes and/or molecular mechanisms controlling the phenotype) addresses one of the fundamental questions in evolutionary ecology. Species distributions and adaptive radiations are thought to be controlled by environmental salinity levels, and efficient osmoregulatory (ionic balance) ability is the main mechanism to overcome the problems related to environmental salinity gradients. Methods To better understand how osmoregulatory performance in freshwater (FW) crustaceans allow individuals to acclimate and adapt to raised salinity conditions, here we (i), reviewed the literature on genes that have been identified to be associated with osmoregulation in FW crustaceans, and (ii), performed a transcriptomic analysis using cDNA libraries developed from mRNA isolated from three important osmoregulatory tissues (gill, antennal gland, hepatopancreas) and total mRNA from post larvae taken from the freshwater prawn, Macrobrachium australiense using Illumina deep sequencing technology. This species was targeted because it can complete its life cycle totally in freshwater but, like many Macrobrachium sp., can also tolerate brackish water conditions and hence should have genes associated with tolerance of both FW and saline conditions. Results We obtained between 55.4 and 65.2 million Illumina read pairs from four cDNA libraries. Overall, paired end sequences assembled into a total of 125,196 non-redundant contigs (≥200 bp) with an N50 length of 2,282 bp and an average contig length of 968 bp. Transcriptomic analysis of M. australiense identified 32 different gene families that were potentially involved with osmoregulatory capacity. A total of 32,597 transcripts were specified with gene ontology (GO) terms identified on the basis of GO categories. Abundance estimation of expressed genes based on TPM (transcript per million) ≥20 showed 1625 transcripts commonly expressed in all four libraries. Among the top 10 genes expressed in four tissue libraries associated with osmoregulation, arginine kinase and Na+/K+- ATPase showed the highest transcript copy number with 7098 and 660, respectively in gill which is considered to be the most important organ involved in osmoregulation. Discussion The current study provides the first broad transcriptome from M. australiense using next generation sequencing and identifies potential candidate genes involved in salinity tolerance and osmoregulation that can provide a foundation for investigating osmoregulatory capacity in a wide variety of freshwater crustaceans.
Collapse
Affiliation(s)
- Azam Moshtaghi
- Science and Engineering Faculty, Queensland University of Technology , Brisbane , Queensland , Australia
| | - Md Lifat Rahi
- Science and Engineering Faculty, Queensland University of Technology , Brisbane , Queensland , Australia
| | - Viet Tuan Nguyen
- School of Science and Engineering, University of the Sunshine Coast , Sippy Downs , Queensland , Australia
| | - Peter B Mather
- Science and Engineering Faculty, Queensland University of Technology , Brisbane , Queensland , Australia
| | - David A Hurwood
- Science and Engineering Faculty, Queensland University of Technology , Brisbane , Queensland , Australia
| |
Collapse
|