1
|
Moore B, Herrera M, Gairin E, Li C, Miura S, Jolly J, Mercader M, Izumiyama M, Kawai E, Ravasi T, Laudet V, Ryu T. The chromosome-scale genome assembly of the yellowtail clownfish Amphiprion clarkii provides insights into the melanic pigmentation of anemonefish. G3 (BETHESDA, MD.) 2023; 13:6982751. [PMID: 36626199 PMCID: PMC9997566 DOI: 10.1093/g3journal/jkad002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/25/2022] [Revised: 11/25/2022] [Accepted: 12/12/2022] [Indexed: 01/11/2023]
Abstract
Anemonefish are an emerging group of model organisms for studying genetic, ecological, evolutionary, and developmental traits of coral reef fish. The yellowtail clownfish Amphiprion clarkii possesses species-specific characteristics such as inter-species co-habitation, high intra-species color variation, no anemone specificity, and a broad geographic distribution, that can increase our understanding of anemonefish evolutionary history, behavioral strategies, fish-anemone symbiosis, and color pattern evolution. Despite its position as an emerging model species, the genome of A. clarkii is yet to be published. Using PacBio long-read sequencing and Hi-C chromatin capture technology, we generated a high-quality chromosome-scale genome assembly initially comprised of 1,840 contigs with an N50 of 1,203,211 bp. These contigs were successfully anchored into 24 chromosomes of 843,582,782 bp and annotated with 25,050 protein-coding genes encompassing 97.0% of conserved actinopterygian genes, making the quality and completeness of this genome the highest among all published anemonefish genomes to date. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further analyses revealed higher copy numbers of erbb3b (a gene involved in melanocyte development) in A. clarkii compared with other anemonefish, thus suggesting a possible link between erbb3b and the natural melanism polymorphism observed in A. clarkii. The publication of this high-quality genome, along with A. clarkii's many unique traits, position this species as an ideal model organism for addressing scientific questions across a range of disciplines.
Collapse
Affiliation(s)
- Billy Moore
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Marcela Herrera
- Marine Eco-Evo-Devo Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Emma Gairin
- Marine Eco-Evo-Devo Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Chengze Li
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Saori Miura
- Marine Eco-Evo-Devo Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Jeffrey Jolly
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Manon Mercader
- Marine Eco-Evo-Devo Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Michael Izumiyama
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Erina Kawai
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Timothy Ravasi
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan.,Australian Research Council Centre of Excellence for Coral Reef Studies, James Cook University, Townsville, QLD 4811, Australia
| | - Vincent Laudet
- Marine Eco-Evo-Devo Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan.,Marine Research Station, Institute of Cellular and Organismic Biology, Academia Sinica, I-Lan 262, Taiwan
| | - Taewoo Ryu
- Marine Climate Change Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| |
Collapse
|
2
|
Spatial transcriptome analysis of long non-coding RNAs reveals tissue specificity and functional roles in cancer. J Zhejiang Univ Sci B 2023; 24:15-31. [PMID: 36632748 PMCID: PMC9837373 DOI: 10.1631/jzus.b2200206] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Long non-coding RNAs (lncRNAs) play a significant role in maintaining tissue morphology and functions, and their precise regulatory effectiveness is closely related to expression patterns. However, the spatial expression patterns of lncRNAs in humans are poorly characterized. Here, we constructed five comprehensive transcriptomic atlases of human lncRNAs covering thousands of major tissue samples in normal and disease states. The lncRNA transcriptomes exhibited high consistency within the same tissues across resources, and even higher complexity in specialized tissues. Tissue-elevated (TE) lncRNAs were identified in each resource and robust TE lncRNAs were refined by integrative analysis. We detected 1 to 4684 robust TE lncRNAs across tissues; the highest number was in testis tissue, followed by brain tissue. Functional analyses of TE lncRNAs indicated important roles in corresponding tissue-related pathways. Moreover, we found that the expression features of robust TE lncRNAs made them be effective biomarkers to distinguish tissues; TE lncRNAs also tended to be associated with cancer, and exhibited differential expression or were correlated with patient survival. In summary, spatial classification of lncRNAs is the starting point for elucidating the function of lncRNAs in both maintenance of tissue morphology and progress of tissue-constricted diseases.
Collapse
|
3
|
Roh H, Kim DH. Identification, classification and functional characterization of HSP70s in rainbow trout (Oncorhynchus mykiss) through multi-omics approaches. FISH & SHELLFISH IMMUNOLOGY 2022; 121:205-214. [PMID: 34990808 DOI: 10.1016/j.fsi.2021.12.059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Revised: 12/28/2021] [Accepted: 12/30/2021] [Indexed: 06/14/2023]
Abstract
Heat shock protein 70s (HSP70s) are known to play vital biological processes in rainbow trout. However, information on the numerous roles and classification of many different HSP70s is insufficient. The purpose of this study was to investigate the characteristics of all HSP70s in rainbow trout using multi-dimensional genomic and transcriptomic analyses for inspecting HSP70 homologs, phylogenetic characteristics, DNA motifs, and transcription factor binding sites (TFBSs). Also, the transcriptomic results in conditions of acute thermal stress and Ichthyophthirius multifiliis infection were used to characterize the expression of all HSP70 homologs, and the isoforms of the most sensitive HSP70 were predicted in silico. A total of 23 HSP70s were identified, and they were divided into seven evolutionary groups (groups 1-7). Groups 1 and 2 had relatively longer phylogenetic distances compared to the other groups, which can speculate origin of groups 1 and 2 HSP70s would be different compared to others. With transcriptomic profiling, most HSPs belonging to group 3 showed highly sensitive responses to I. multifiliis infection, not thermal stress, but the group 6 HSP70s had the opposite expression tendencies. Likewise, the composition of the TFBS in each HSP70 was consistent with its group classification. Since TFBSs are widely known to influence transcriptomic expression, they could be one of the major reasons for the different patterns of expression within the HSP70 groups. Moreover, this study demonstrated several isoforms of HSP70a, by far the most sensitive HSP70s, under several stress environments such as hypoxia, thermal, and overcrowding stress. This is an important fundamental study to expand the understanding of HSP70s in rainbow trout as well as for selecting the most sensitive biomarkers for types of stress.
Collapse
Affiliation(s)
- HyeongJin Roh
- Department of Aquatic Life Medicine, College of Fisheries Sciences, Pukyong National University, 45, Yongso-ro, Nam-Gu, Busan, South Korea
| | - Do-Hyung Kim
- Department of Aquatic Life Medicine, College of Fisheries Sciences, Pukyong National University, 45, Yongso-ro, Nam-Gu, Busan, South Korea.
| |
Collapse
|
4
|
Comprehensive Characterization of Multitissue Expression Landscape, Co-Expression Networks and Positive Selection in Pikeperch. Cells 2021; 10:cells10092289. [PMID: 34571938 PMCID: PMC8471114 DOI: 10.3390/cells10092289] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Revised: 08/27/2021] [Accepted: 08/29/2021] [Indexed: 11/19/2022] Open
Abstract
Promising efforts are ongoing to extend genomics resources for pikeperch (Sander lucioperca), a species of high interest for the sustainable European aquaculture sector. Although previous work, including reference genome assembly, transcriptome sequence, and single-nucleotide polymorphism genotyping, added a great wealth of genomic tools, a comprehensive characterization of gene expression across major tissues in pikeperch still remains an unmet research need. Here, we used deep RNA-Sequencing of ten vital tissues collected in eight animals to build a high-confident and annotated trancriptome atlas, to detect the tissue-specificity of gene expression and co-expression network modules, and to investigate genome-wide selective signatures in the Percidae fish family. Pathway enrichment and protein–protein interaction network analyses were performed to characterize the unique biological functions of tissue-specific genes and co-expression modules. We detected strong functional correlations and similarities of tissues with respect to their expression patterns—but also significant differences in the complexity and composition of their transcriptomes. Moreover, functional analyses revealed that tissue-specific genes essentially play key roles in the specific physiological functions of the respective tissues. Identified network modules were also functionally coherent with tissues’ main physiological functions. Although tissue specificity was not associated with positive selection, several genes under selection were found to be involved in hypoxia, immunity, and gene regulation processes, that are crucial for fish adaption and welfare. Overall, these new resources and insights will not only enhance the understanding of mechanisms of organ biology in pikeperch, but also complement the amount of genomic resources for this commercial species.
Collapse
|
5
|
Histological and biochemical evaluation of skeletal muscle in the two salmonid species Coregonus maraena and Oncorhynchus mykiss. PLoS One 2021; 16:e0255062. [PMID: 34383783 PMCID: PMC8360549 DOI: 10.1371/journal.pone.0255062] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Accepted: 07/08/2021] [Indexed: 12/03/2022] Open
Abstract
The growth of fishes and their metabolism is highly variable in fish species and is an indicator for fish fitness. Therefore, somatic growth, as a main biological process, is ecologically and economically significant. The growth differences of two closely related salmonids, rainbow trout (Oncorhynchus mykiss) and maraena whitefsh (Coregonus maraena), have not been adequately studied as a comparative study and are therefore insufficiently understood. For this reason, our aim was to examine muscle growth in more detail and provide a first complex insight into the growth and muscle metabolism of these two fish species at slaughter size. In addition to skeletal muscle composition (including nuclear counting and staining of stem and progenitor cells), biochemical characteristics, and enzyme activity (creatine kinase, lactate dehydrogenase, isocitrate dehydrogenase) of rainbow trout and maraena whitefish were determined. Our results indicate that red muscle contains cells with a smaller diameter compared to white muscle and those fibres had more stem and progenitor cells as a proportion of total nuclei. Interestingly, numerous interspecies differences were identified; in rainbow trout muscle RNA content, intermediate fibres and fibre diameter and in whitefish red muscle cross-sectional area, creatine kinase activity were higher compared to the other species at slaughter weight. The proportional reduction in red muscle area, accompanied by an increase in DNA content and a lower activity of creatine kinase, exhibited a higher degree of hypertrophic growth in rainbow trout compared to maraena whitefish, which makes this species particularly successful as an aquaculture species.
Collapse
|
6
|
Ali A, Thorgaard GH, Salem M. PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes. Front Genet 2021; 12:683408. [PMID: 34335690 PMCID: PMC8321248 DOI: 10.3389/fgene.2021.683408] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Accepted: 06/14/2021] [Indexed: 01/04/2023] Open
Abstract
Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.
Collapse
Affiliation(s)
- Ali Ali
- Department of Animal and Avian Sciences, University of Maryland, College Park, College Park, MD, United States
| | - Gary H. Thorgaard
- School of Biological Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA, United States
| | - Mohamed Salem
- Department of Animal and Avian Sciences, University of Maryland, College Park, College Park, MD, United States
| |
Collapse
|
7
|
Li X, Wu J, Xiao X, Rong Y, Yang H, Li J, Zhou Q, Zhou W, Shi J, Qi H, Du H. Characterization and complexity of transcriptome in Gymnocypris przewalskii using single-molecule long-read sequencing and RNA-seq. DNA Res 2021; 28:6275749. [PMID: 33989386 PMCID: PMC8320875 DOI: 10.1093/dnares/dsab005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Accepted: 05/11/2021] [Indexed: 11/30/2022] Open
Abstract
The Tibetan Schizothoracinae fish Gymnocypris przewalskii has the ability to adapt to the extreme plateau environment, making it an ideal biological material for evolutionary biology research. However, the lack of well-annotated reference genomes has limited the study of the molecular genetics of G. przewalskii. To characterize its transcriptome features, we first used long-read sequencing technology in combination with RNA-seq for transcriptomic analysis. A total of 159,053 full-length (FL) transcripts were captured by Iso-Seq, having a mean length of 3,445 bp with N50 value of 4,348. Of all FL transcripts, 145,169 were well-annotated in the public database and 134,537 contained complete open reading frames. There were 4,149 pairs of alternative splicing events, of which three randomly selected were defined by RT–PCR and sequencing, and 13,293 long non-coding RNAs detected, based on all-vs.-all BLAST. A total of 118,185 perfect simple sequence repeats were identified from FL transcripts. The FL transcriptome might provide basis for further research of G. przewalskii.
Collapse
Affiliation(s)
- Xindan Li
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China.,College of Marine Science, Shanghai Ocean University, Shanghai 201306, China
| | - Jinming Wu
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China
| | - Xinping Xiao
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China
| | - Yifeng Rong
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China.,College of Marine Science, Shanghai Ocean University, Shanghai 201306, China
| | - Haile Yang
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China
| | - Junyi Li
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China
| | - Qiong Zhou
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China
| | - Weiguo Zhou
- The Rescue and Rehabilitation Center of Naked Carps in Lake Qinghai, Xining, Qinghai 810016, China
| | - Jianquan Shi
- The Rescue and Rehabilitation Center of Naked Carps in Lake Qinghai, Xining, Qinghai 810016, China
| | - Hongfang Qi
- The Rescue and Rehabilitation Center of Naked Carps in Lake Qinghai, Xining, Qinghai 810016, China
| | - Hao Du
- Key Laboratory of freshwater biodiversity conservation, Ministry of Agriculture, Chinese Academy of Fishery Sciences, Wuhan, Hubei 430223, China.,College of Marine Science, Shanghai Ocean University, Shanghai 201306, China.,The Rescue and Rehabilitation Center of Naked Carps in Lake Qinghai, Xining, Qinghai 810016, China
| |
Collapse
|
8
|
Ma H, Han YC, Palti Y, Gao G, Liu S, Palmquist DE, Wiens GD, Shepherd BS. Structure and regulation of the NK-lysin (1-4) and NK-lysin like (a and b) antimicrobial genes in rainbow trout (Oncorhynchus mykiss). DEVELOPMENTAL AND COMPARATIVE IMMUNOLOGY 2021; 116:103961. [PMID: 33301795 DOI: 10.1016/j.dci.2020.103961] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 12/03/2020] [Accepted: 12/03/2020] [Indexed: 06/12/2023]
Abstract
Nk-lysin (Nkl), an antimicrobial peptide (AMP) product of natural killer cells and cytotoxic T cells in mammals, has recently been characterized in a number of finfish species. In this study, we identified six genes with sequence homology to Nkl and characterized their patterns of mRNA expression and abundances in rainbow trout (Oncorhynchus mykiss). The cDNA sequences for the six Nkls encoded precursor peptides of 128-133 aa in length, and mature peptides of 109-111 aa in length. Genomic DNA of the nkl1-4 genes consisted of five exons and four introns, whereas the nkl-like a & b genes consisted of four exons and three introns. Chromosomal locations of these peptides show that nkl1 was located on chromosome arm 25q, whereas the other five nkl genes were clustered on chromosome arm 19q. Phylogenetic analysis revealed a conserved structure of Nkls among the teleosts and further protein sequence analyses suggests that all six nkl genes fall within the Nkl sub-family of the Saposin family of proteins. Patterns of tissue-specific mRNA expression were asymmetric among the six trout Nkl homologues, with nkl1, nkl3, and nkl-like a & b occurring in immune competent organs such as spleen, gill, intestine and kidney, as well as pineal gland, brain and oocytes. However, nkl2 and nkl4, showed primary abundances in brain, pineal gland and oocyte tissues. Using mRNA sequencing, in whole-body pools of juvenile trout fry (1 g bw) exposed to Flavobacterium psychrophilum infection, we observed modest up-regulation (2-3 fold) of five (nkl 2-4 and nkl-like a & b) of the six nkl mRNAs over the five-day post-challenge time-course. However, no upregulation could be recorded in spleen tissue measured by qPCR in juvenile trout (270 g bw). Using mRNA sequencing again, mRNA abundances were determined in gill of juvenile trout (~57.7 g bw) exposed to various aquaculture stressors. The results indicated that all six nkls (nkl1-4 and nkl-like a and nkl-like b) were downregulated when exposed to high temperature, and that nkl1 was significantly downregulated following salinity challenge. Overall, these newly characterized AMPs may contribute to host innate immunity as they are modulated following pathogen challenge and by physiological stressors.
Collapse
Affiliation(s)
- Hao Ma
- USDA-ARS-NADC-Ruminant Diseases and Immunology Research Unit, 1920 Dayton Ave, Ames, IA, 50010, USA; USDA-ARS-National Center for Cool and Cold Water Aquaculture, 11861 Leetown Rd., Leetown, WV, 25430, USA
| | - Yueh-Chiang Han
- USDA-ARS-School of Freshwater Sciences, 600 E. Greenfield Ave., Milwaukee, WI, 53204, USA
| | - Yniv Palti
- USDA-ARS-National Center for Cool and Cold Water Aquaculture, 11861 Leetown Rd., Leetown, WV, 25430, USA
| | - Guangtu Gao
- USDA-ARS-National Center for Cool and Cold Water Aquaculture, 11861 Leetown Rd., Leetown, WV, 25430, USA
| | - Sixin Liu
- USDA-ARS-National Center for Cool and Cold Water Aquaculture, 11861 Leetown Rd., Leetown, WV, 25430, USA
| | - Debra E Palmquist
- USDA/ARS-Midwest Area Statistics Unit, 1815 N. Street, Peoria, IL, 61604, USA
| | - Gregory D Wiens
- USDA-ARS-National Center for Cool and Cold Water Aquaculture, 11861 Leetown Rd., Leetown, WV, 25430, USA
| | - Brian S Shepherd
- USDA-ARS-School of Freshwater Sciences, 600 E. Greenfield Ave., Milwaukee, WI, 53204, USA.
| |
Collapse
|
9
|
Gan W, Chung-Davidson YW, Chen Z, Song S, Cui W, He W, Zhang Q, Li W, Li M, Ren J. Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish. Sci Rep 2021; 11:1815. [PMID: 33469041 PMCID: PMC7815744 DOI: 10.1038/s41598-020-80168-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Accepted: 12/14/2020] [Indexed: 02/06/2023] Open
Abstract
Goldfish is an ornamental fish with diverse phenotypes. However, the limited genomic resources of goldfish hamper our understanding of the genetic basis for its phenotypic diversity. To provide enriched genomic resources and infer possible mechanisms underlying skin pigmentation, we performed a large-scale transcriptomic sequencing on 13 adult goldfish tissues, larvae at one- and three-days post hatch, and skin tissues with four different color pigmentation. A total of 25.52 Gb and 149.80 Gb clean data were obtained using the PacBio and Illumina platforms, respectively. Onto the goldfish reference genome, we mapped 137,674 non-redundant transcripts, of which 5.54% was known isoforms and 78.53% was novel isoforms of the reference genes, and the remaining 21,926 isoforms are novel isoforms of additional new genes. Both skin-specific and color-specific transcriptomic analyses showed that several significantly enriched genes were known to be involved in melanogenesis, tyrosine metabolism, PPAR signaling pathway, folate biosynthesis metabolism and so on. Thirteen differentially expressed genes across different color skins were associated with melanogenesis and pteridine synthesis including mitf, ednrb, mc1r, tyr, mlph and gch1, and xanthophore differentiation such as pax7, slc2a11 and slc2a15. These transcriptomic data revealed pathways involved in goldfish pigmentation and improved the gene annotation of the reference genome.
Collapse
Affiliation(s)
- Wu Gan
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China
| | - Yu-Wen Chung-Davidson
- Department of Fisheries and Wildlife, Michigan State University, East Lansing, MI, 48824, USA
| | - Zelin Chen
- South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China
| | - Shiying Song
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China
| | - Wenyao Cui
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China
| | - Wei He
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China
| | - Qinghua Zhang
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China
- International Research Center for Marine Biosciences, Ministry of Science and Technology, Shanghai Ocean University, Shanghai, 201306, China
| | - Weiming Li
- Department of Fisheries and Wildlife, Michigan State University, East Lansing, MI, 48824, USA
| | - Mingyou Li
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China.
| | - Jianfeng Ren
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, 201306, China.
- International Research Center for Marine Biosciences, Ministry of Science and Technology, Shanghai Ocean University, Shanghai, 201306, China.
| |
Collapse
|
10
|
Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper ( Epinephelus tukula). BIOLOGY 2021; 10:biology10010036. [PMID: 33430356 PMCID: PMC7825770 DOI: 10.3390/biology10010036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Revised: 12/25/2020] [Accepted: 01/05/2021] [Indexed: 12/16/2022]
Abstract
Simple Summary The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits. Abstract The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F1 juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.
Collapse
|
11
|
Abstract
Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.
Collapse
Affiliation(s)
- M. Lisette Delgado
- Department of Biology, Dalhousie University, Halifax, Nova Scotia, B3H 4R2, Canada
| | - Daniel E. Ruzzante
- Department of Biology, Dalhousie University, Halifax, Nova Scotia, B3H 4R2, Canada
| |
Collapse
|
12
|
Yamaguchi T, Chang CJ, Karger A, Keller M, Pfaff F, Wangkahart E, Wang T, Secombes CJ, Kimoto A, Furihata M, Hashimoto K, Fischer U, Dijkstra JM. Ancient Cytokine Interleukin 15-Like (IL-15L) Induces a Type 2 Immune Response. Front Immunol 2020; 11:549319. [PMID: 33193315 PMCID: PMC7658486 DOI: 10.3389/fimmu.2020.549319] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2020] [Accepted: 08/28/2020] [Indexed: 12/20/2022] Open
Abstract
Related interleukin-2, -15, and -15-like (IL-2, -15, and -15L) are ancient cytokines, with all three genes surviving in extant fish and some mammals. The present study is the first to identify IL-15L functions, namely in rainbow trout. In isolated trout splenocytes, and in vivo, purified recombinant IL-15L+IL-15Rα molecules induced expression of IL-4 and IL-13 homologs, which are markers of type 2 immunity. In contrast, trout IL-15 stimulated type 1 immunity markers, thus IL-15 and IL-15L can have opposing functions. Trout IL-15L was more dependent on "in trans" presentation by the receptor chain IL-15Rα than IL-15, and stimulated CD4-CD8-(IgM-) lymphocytes from thymus and spleen. We propose an important role for IL-15L early in the type 2 immunity cytokine cascade. Trout IL-2 and IL-15 exhibited features reminiscent of their mechanistic and functional dichotomy observed in mammals; for example, IL-15 but not IL-2 required a receptor alpha chain (only IL-15Rα in the case of fish) for its stability, and only IL-15 was efficient in stimulating lymphocytes from mucosal tissues. Data suggest that IL-15L and IL-15 may be particularly effective in stimulating innate lymphocyte type 2 cells (ILC2) and natural killer (NK) cells, respectively, but further identification of the cell types is needed. An interesting finding different from in mammals was the efficient stimulation of CD4+CD8+ thymocytes by IL-2. In short, this study presents fundamental information on the evolution of the IL-2/15/15L cytokine family.
Collapse
Affiliation(s)
- Takuya Yamaguchi
- Institute of Infectology, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Chia Jung Chang
- Institute of Infectology, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Axel Karger
- Institute of Molecular Virology and Cell Biology, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Markus Keller
- Institute of Novel and Emerging Infectious Diseases, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Florian Pfaff
- Institute of Diagnostic Virology, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Eakapol Wangkahart
- Scottish Fish Immunology Research Centre, School of Biological Sciences, University of Aberdeen, Aberdeen, United Kingdom
| | - Tiehui Wang
- Scottish Fish Immunology Research Centre, School of Biological Sciences, University of Aberdeen, Aberdeen, United Kingdom
| | - Christopher J Secombes
- Scottish Fish Immunology Research Centre, School of Biological Sciences, University of Aberdeen, Aberdeen, United Kingdom
| | - Azusa Kimoto
- Institute for Comprehensive Medical Science, Fujita Health University, Toyoake, Japan
| | - Mitsuru Furihata
- Nagano Prefectural Fisheries Experimental Station, Nagano, Japan
| | - Keiichiro Hashimoto
- Institute for Comprehensive Medical Science, Fujita Health University, Toyoake, Japan
| | - Uwe Fischer
- Institute of Infectology, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Johannes M Dijkstra
- Institute for Comprehensive Medical Science, Fujita Health University, Toyoake, Japan
| |
Collapse
|
13
|
Ocalewicz K, Gurgul A, Polonis M, Dobosz S. Preliminary Identification of Candidate Genes Related to Survival of Gynogenetic Rainbow Trout ( Oncorhynchus mykiss) Based on Comparative Transcriptome Analysis. Animals (Basel) 2020; 10:ani10081326. [PMID: 32751994 PMCID: PMC7459965 DOI: 10.3390/ani10081326] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Revised: 07/28/2020] [Accepted: 07/29/2020] [Indexed: 11/16/2022] Open
Abstract
In the present research, the eggs from four rainbow trout females were used to provide four groups of gynogenetic doubled haploids (DHs). The quality of the eggs from different clutches was comparable, however, interclutch differences were observed in the gynogenetic variants of the experiment and the survival of DH specimens from different groups varied from 3% to 57% during embryogenesis. Transcriptome analysis of the eggs from different females exhibited inter-individual differences in the maternal genes' expression. Eggs originating from females whose gynogenetic offspring had the highest survival showed an increased expression of 46 genes when compared to the eggs from three other females. Eggs with the highest survival of gynogenetic embryos showed an up-regulation of genes that are associated with cell survival, migration and differentiation (tyrosine-protein kinase receptor TYRO3-like gene), triglyceride metabolism (carnitine O-palmitoyltransferase 1 gene), biosynthesis of polyunsaturated fat (3-oxoacyl-acyl-carrier-protein reductase gene), early embryogenic development (protein argonaute-3 gene, leucine-rich repeat-containing protein 3-like gene), 5S RNA binding (ribosome biogenesis regulatory protein homolog) as well as senescence and aging (telomerase reverse transcriptase, TERT gene), among others. Positive correlation between the genotypic efficiency and egg transcriptome profiles indicated that at least some of the differentially expressed genes should be considered as potential candidate genes for the efficiency of gynogenesis in rainbow trout.
Collapse
Affiliation(s)
- Konrad Ocalewicz
- Department of Marine Biology and Ecology, Institute of Oceanography, Faculty of Oceanography and Geography, University of Gdansk, M. Piłsudskiego 46 Av, 81-378 Gdynia, Poland;
- Correspondence:
| | - Artur Gurgul
- Centre for Experimental and Innovative Medicine, University of Agriculture in Kraków, Rędzina 1c, 30-248 Kraków, Poland;
| | - Marcin Polonis
- Department of Marine Biology and Ecology, Institute of Oceanography, Faculty of Oceanography and Geography, University of Gdansk, M. Piłsudskiego 46 Av, 81-378 Gdynia, Poland;
| | - Stefan Dobosz
- Department of Salmonid Research, Inland Fisheries Institute in Olsztyn, Rutki, 83-330 Żukowo, Poland;
| |
Collapse
|
14
|
Feng X, Jia Y, Zhu R, Chen K, Chen Y. Characterization and analysis of the transcriptome in Gymnocypris selincuoensis on the Qinghai-Tibetan Plateau using single-molecule long-read sequencing and RNA-seq. DNA Res 2020; 26:353-363. [PMID: 31274170 PMCID: PMC6704404 DOI: 10.1093/dnares/dsz014] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2019] [Accepted: 06/11/2019] [Indexed: 12/26/2022] Open
Abstract
The lakes on the Qinghai-Tibet Plateau (QTP) are the largest and highest lake group in the world. Gymnocypris selincuoensis is the only cyprinid fish living in lake Selincuo, the largest lake on QTP. However, its genetic resource is still blank, limiting studies on molecular and genetic analysis. In this study, the transcriptome of G. selincuoensis was first generated by using PacBio Iso-Seq and Illumina RNA-seq. A full-length (FL) transcriptome with 75,435 transcripts was obtained by Iso-Seq with N50 length of 3,870 bp. Among all transcripts, 75,016 were annotated to public databases, 64,710 contain complete open reading frames and 2,811 were long non-coding RNAs. Based on all- vs.-all BLAST, 2,069 alternative splicing events were detected, and 80% of them were validated by reverse transcription polymerase chain reaction (RT-PCR). Tissue gene expression atlas showed that the number of detected expressed transcripts ranged from 37,397 in brain to 19,914 in muscle, with 10,488 transcripts detected in all seven tissues. Comparative genomic analysis with other cyprinid fishes identified 77 orthologous genes with potential positive selection (Ka/Ks > 0.3). A total of 56,696 perfect simple sequence repeats were identified from FL transcripts. Our results provide valuable genetic resources for further studies on adaptive evolution, gene expression and population genetics in G. selincuoensis and other congeneric fishes.
Collapse
Affiliation(s)
- Xiu Feng
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| | - Yintao Jia
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| | - Ren Zhu
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| | - Kang Chen
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Yifeng Chen
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| |
Collapse
|
15
|
Mehra R, Jasrotia RS, Mahajan A, Sharma D, Iquebal MA, Kaul S, Dhar MK. Transcriptome analysis of Snow Mountain Garlic for unraveling the organosulfur metabolic pathway. Genomics 2020; 112:99-107. [DOI: 10.1016/j.ygeno.2019.07.014] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2019] [Revised: 07/05/2019] [Accepted: 07/25/2019] [Indexed: 12/20/2022]
|
16
|
Syahputra K, Kania PW, Al-Jubury A, Jafaar RM, Dirks RP, Buchmann K. Transcriptomic analysis of immunity in rainbow trout (Oncorhynchus mykiss) gills infected by Ichthyophthirius multifiliis. FISH & SHELLFISH IMMUNOLOGY 2019; 86:486-496. [PMID: 30513380 DOI: 10.1016/j.fsi.2018.11.075] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2018] [Revised: 11/22/2018] [Accepted: 11/30/2018] [Indexed: 06/09/2023]
Abstract
The parasite Ichthyophthirius multifiliis infecting skin, fins and gills of a wide range of freshwater fish species, including rainbow trout, is known to induce a protective immune response in the host. Although a number of studies have reported activation of several immune genes in infected fish host, the immune response picture is still considered incomplete. In order to address this issue, a comparative transcriptomic analysis was performed on infected versus uninfected rainbow trout gills and it showed that a total of 3352 (7.2%) out of 46,585 identified gene sequences were significantly regulated after parasite infection. Of differentially expressed gene sequences, 1796 genes were up-regulated and 1556 genes were down-regulated. These were classified into 61 Gene Ontology (GO) terms and mapped to 282 reference canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Infection of I. multifiliis induced a clear differential expression of immune genes, related to both innate and adaptive immunity. A total of 268 (6.86%) regulated gene sequences were known to take part in 16 immune-related pathways. These involved pathways related to the innate immunity such as the Chemokine signaling pathway, Platelet activation, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, and Leukocyte transendothelial migration. Elevated transcription of genes encoding the TLR 8 gene and chemokines (CCL4, CCL19, CCL28, CXCL8, CXCL11, CXCL13, CXCL14) was recorded indicating their roles in recognition of I. multifiliis and subsequent induction of the inflammatory response, respectively. A number of upregulated genes in infected gills were associated with antigen processing/presentation and T and B cell receptor signaling (including B cell marker CD22 involved in B cell development). Overall the analysis supports the notion that I. multifiliis induces a massive and varied innate response upon which a range of adaptive immune responses are established which may contribute to the long lasting protection of immunized rainbow trout.
Collapse
Affiliation(s)
- Khairul Syahputra
- Department of Veterinary and Animal Science, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg C, Denmark.
| | - Per W Kania
- Department of Veterinary and Animal Science, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg C, Denmark
| | - Azmi Al-Jubury
- Department of Veterinary and Animal Science, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg C, Denmark
| | - Rzgar M Jafaar
- Department of Veterinary and Animal Science, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg C, Denmark
| | - Ron P Dirks
- Future Genomics Technologies B.V., Leiden, the Netherlands
| | - Kurt Buchmann
- Department of Veterinary and Animal Science, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg C, Denmark
| |
Collapse
|
17
|
Watanabe L, Gomes F, Vianez J, Nunes M, Cardoso J, Lima C, Schneider H, Sampaio I. De novo transcriptome based on next-generation sequencing reveals candidate genes with sex-specific expression in Arapaima gigas (Schinz, 1822), an ancient Amazonian freshwater fish. PLoS One 2018; 13:e0206379. [PMID: 30372461 PMCID: PMC6205615 DOI: 10.1371/journal.pone.0206379] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 10/11/2018] [Indexed: 12/24/2022] Open
Abstract
Background The Arapaima (Arapaima gigas) is one of the world's largest freshwater bony fish, and is found in the rivers of the Amazon basin. This species is a potential aquaculture resource, although reproductive management in captivity is limited in particular due to the lack of external sexual dimorphism. In this study, using the 454 Roche platform (pyrosequencing) techniques, we evaluated a major portion of the transcriptome of this important Amazonian species. Results Four libraries obtained from the liver and skin tissue of juvenile specimens (representing males and females separately) were sequenced, yielding 5,453,919 high-quality reads. The de novo transcriptome assembly resulted in 175,792 contigs, with 51,057 significant blast hits. A total of 38,586 transcripts were mapped by Gene Ontology using Blast2GO. We identified 20,219 genes in the total transcriptome (9,551 in the liver and 16,818 in the skin). The gene expression analyses indicated 105 genes in the liver and 204 in the skin with differentiated expression profiles, with 95 being over-expressed in the females and 214 in the males. The log2 Fold Change and heatmap based on Reads Per Kilobase per Million mapped reads (RPKM) revealed that the gene expression in the skin is highly differentiated between male and female arapaima, while the levels of expression in the liver are similar between the sexes. Conclusion Transcriptome analysis based on pyrosequencing proved to be a reliable tool for the identification of genes with differentiated expression profiles between male and female arapaima. These results provide useful insights into the molecular pathways of sexual dimorphism in this important Amazonian species, and for comparative analyses with other teleosts.
Collapse
Affiliation(s)
- Luciana Watanabe
- Laboratório de Genética e Biologia Molecular, Instituto de Estudos Costeiros (IECOS), Universidade Federal do Pará, Campus de Bragança, Pará, Brazil
- * E-mail:
| | - Fátima Gomes
- Laboratório de Genética e Biologia Molecular, Instituto de Estudos Costeiros (IECOS), Universidade Federal do Pará, Campus de Bragança, Pará, Brazil
| | - João Vianez
- Centro de Inovações Tecnológicas (CIT), Instituto Evandro Chagas (IEC), Ananindeua, Pará, Brazil
| | - Márcio Nunes
- Centro de Inovações Tecnológicas (CIT), Instituto Evandro Chagas (IEC), Ananindeua, Pará, Brazil
| | - Jedson Cardoso
- Centro de Inovações Tecnológicas (CIT), Instituto Evandro Chagas (IEC), Ananindeua, Pará, Brazil
| | - Clayton Lima
- Centro de Inovações Tecnológicas (CIT), Instituto Evandro Chagas (IEC), Ananindeua, Pará, Brazil
| | - Horacio Schneider
- Laboratório de Genética e Biologia Molecular, Instituto de Estudos Costeiros (IECOS), Universidade Federal do Pará, Campus de Bragança, Pará, Brazil
| | - Iracilda Sampaio
- Laboratório de Genética e Biologia Molecular, Instituto de Estudos Costeiros (IECOS), Universidade Federal do Pará, Campus de Bragança, Pará, Brazil
| |
Collapse
|
18
|
Shaw CH, Gao G, Wiens GD. Differential expression and evolution of three tandem, interleukin-1 receptor-like 1 genes in rainbow trout (Oncorhynchus mykiss). DEVELOPMENTAL AND COMPARATIVE IMMUNOLOGY 2018; 87:193-203. [PMID: 29886053 DOI: 10.1016/j.dci.2018.06.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Revised: 06/06/2018] [Accepted: 06/06/2018] [Indexed: 06/08/2023]
Abstract
Interleukin-1 receptor-like 1 (Il1rl1 or ST2), a member of the interleukin-1 receptor family, has pleiotropic roles including tissue homeostasis, inflammation, immune polarization, and disease resistance in mammals. A single orthologue was previously described in salmonid fish; however, a recently improved genome assembly of rainbow trout (Oncorhynchus mykiss) revealed three adjacent, tandem il1rl1 orthologues on chromosome Omy 03. Here, we report the genomic organization and evolution of the three il1rl1 genes (il1rl1α, il1rl1β, il1rl1γ), and use both RNA-seq and gene-specific qPCR methods to quantify expression patterns. Nucleotide sequence homology between the three genes is >95% and each predicted protein contains three IG/IG-like domains, a transmembrane region and a TIR domain. The amino acid sequence homology of the rainbow trout il1rl1 genes are highly related to two functional copies in Atlantic and Coho salmon (∼94%) but relatively low (22-26%) with avian and mammalian species. Transcript abundance measured by RNA-seq in 15 tissues of healthy adult rainbow trout indicate constitutive expression of each gene. In whole body lysates, il1rl1α was shown to have >20 fold mRNA expression compared to il1rl1β and il1rl1γ as measured by qPCR assays specific to il1rl1α or il1rl1γ, as well as a multi-gene qPCR assay (il1rl1α,β,γ). Unrooted phylogenetic trees grouped the rainbow trout il1rl1 genes apart from other interleukin-1 receptor family genes and genomic comparisons identify preserved synteny between mammals, birds and salmonids albeit a pseudogene is present in both Atlantic salmon and Coho salmon. Phylogenetic analyses suggest that the three genes arose by tandem duplication but are inconclusive whether these events occurred prior-to or after salmonid speciation. These findings further the understanding of interleukin receptor family evolution and their contribution to teleost immune function.
Collapse
Affiliation(s)
- Cassidy H Shaw
- National Center for Cool and Cold Water Aquaculture, Agricultural Research Service, USDA, Kearneysville, WV 25430, USA.
| | - Guangtu Gao
- National Center for Cool and Cold Water Aquaculture, Agricultural Research Service, USDA, Kearneysville, WV 25430, USA
| | - Gregory D Wiens
- National Center for Cool and Cold Water Aquaculture, Agricultural Research Service, USDA, Kearneysville, WV 25430, USA.
| |
Collapse
|
19
|
Transcriptome Analysis of Rainbow Trout ( Oncorhynchus mykiss) Eggs Subjected to the High Hydrostatic Pressure Treatment. Int J Genomics 2018; 2018:5197126. [PMID: 30225244 PMCID: PMC6129359 DOI: 10.1155/2018/5197126] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2017] [Revised: 05/08/2018] [Accepted: 07/15/2018] [Indexed: 11/18/2022] Open
Abstract
High hydrostatic pressure (HHP) causes depolymerization of the spindle microtubules. HHP applied to fish eggs prevents extrusion of the second polar body and inhibits the first cell cleavage, and it is used to produce triploids and diploid gynogenetic and androgenetic individuals. HHP has been also found to affect biomolecules including nucleic acids, and it may be presumed that HHP administered to the rainbow trout (Oncorhynchus mykiss) eggs disturbs cytoplasmic maternal RNA indispensable for the early embryogenesis. To verify this assumption, quality and quantity of RNA extracted from the rainbow trout eggs subjected to the high hydrostatic pressure shock were analyzed. Provided results exhibited that maternal transcriptome was resistant to a three-minute exposure to 65.5 MPa of HHP treatment. Some trend showing increase of the RNA integrity was observed in the HHP-treated eggs; however, the difference was not statistically significant. Alterations in the expression profiles in the rainbow trout eggs subjected to HHP were also negligible. Greater differences in the maternal gene expression were observed between eggs from different clutches than between HHP-treated and untreated eggs from the same clutch. It may be assumed that exposure to HHP shock was too short to modify significantly maternal transcripts in the rainbow trout eggs.
Collapse
|
20
|
Mohamed AR, King H, Evans B, Reverter A, Kijas JW. Multi-Tissue Transcriptome Profiling of North American Derived Atlantic Salmon. Front Genet 2018; 9:369. [PMID: 30271423 PMCID: PMC6146974 DOI: 10.3389/fgene.2018.00369] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2018] [Accepted: 08/22/2018] [Indexed: 12/23/2022] Open
Abstract
The availability of a reference genome assembly for Atlantic salmon, Salmo salar, SNP genotyping platforms and low cost sequencing are enhancing the understanding of both life history and production-related traits in this important commercial species. We collected and analyzed transcriptomes from selected tissues of Atlantic salmon to inform future functional and comparative genomics studies. Messenger RNA (mRNA) was isolated from pituitary gland, brain, ovary, and liver before Illumina sequencing produced a total of 640 million 150-bp paired-end reads. Following read mapping, feature counting, and normalization, cluster analysis identified genes highly expressed in a tissue-specific manner. We identified a cluster of 508 tissue specific genes for pituitary gland, 3395 for brain, 2939 for ovary, and 539 for liver. Functional profiling identified gene clusters describing the unique functions of each tissue. Moreover, highly-expressed transcription factors (TFs) present in each tissue-specific gene cluster were identified. TFs belonging to homeobox and bhlh families were identified for pituitary gland, pou and zf-c2h2 families for brain, arid, and zf-c2h2 for ovary and rxr-like family for liver. The data and analysis presented are relevant to the emerging Functional Annotation of All Salmonid Genomes (FAASG) initiative that is seeking to develop a detailed understanding of both salmonid evolution and the genomic elements that drive gene expression and regulation.
Collapse
Affiliation(s)
- Amin R Mohamed
- Commonwealth Scientific and Industrial Research Organisation Agriculture and Food, Queensland Bioscience Precinct, St Lucia, QLD, Australia.,Zoology Department, Faculty of Science, Benha University, Benha, Egypt
| | - Harry King
- Commonwealth Scientific and Industrial Research Organisation Agriculture, Hobart, TAS, Australia
| | | | - Antonio Reverter
- Commonwealth Scientific and Industrial Research Organisation Agriculture and Food, Queensland Bioscience Precinct, St Lucia, QLD, Australia
| | - James W Kijas
- Commonwealth Scientific and Industrial Research Organisation Agriculture and Food, Queensland Bioscience Precinct, St Lucia, QLD, Australia
| |
Collapse
|
21
|
Ali A, Al-Tobasei R, Kenney B, Leeds TD, Salem M. Integrated analysis of lncRNA and mRNA expression in rainbow trout families showing variation in muscle growth and fillet quality traits. Sci Rep 2018; 8:12111. [PMID: 30108261 PMCID: PMC6092380 DOI: 10.1038/s41598-018-30655-8] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2018] [Accepted: 08/04/2018] [Indexed: 02/07/2023] Open
Abstract
Muscle yield and quality traits are important for the aquaculture industry and consumers. Genetic selection for these traits is difficult because they are polygenic and result from multifactorial interactions. To study the genetic architecture of these traits, phenotypic characterization of whole body weight (WBW), muscle yield, fat content, shear force and whiteness were measured in ~500 fish representing 98 families from a growth-selected line. RNA-Seq was used to sequence the muscle transcriptome of different families exhibiting divergent phenotypes for each trait. We have identified 240 and 1,280 differentially expressed (DE) protein-coding genes and long noncoding RNAs (lncRNAs), respectively, in fish families exhibiting contrasting phenotypes. Expression of many DE lncRNAs (n = 229) was positively correlated with overlapping, neighboring or distantly located protein-coding genes (n = 1,030), resulting in 3,392 interactions. Three DE antisense lncRNAs were co-expressed with sense genes known to impact muscle quality traits. Forty-four DE lncRNAs had potential sponge functions to miRNAs that affect muscle quality traits. This study (1) defines muscle quality associated protein-coding and noncoding genes and (2) provides insight into non-coding RNAs involvement in regulating growth and fillet quality traits in rainbow trout.
Collapse
Affiliation(s)
- Ali Ali
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA
| | - Rafet Al-Tobasei
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA.,Department of Biostatistics, University of Alabama at Birmingham, Birmingham, AL, 35294-0022, USA
| | - Brett Kenney
- Division of Animal and Nutritional Science, West Virginia University, Morgantown, WV, 26506-6108, USA
| | - Timothy D Leeds
- The National Center for Cool and Cold Water Aquaculture, USDA Agricultural Research Service, Kearneysville, WV, 25430, USA
| | - Mohamed Salem
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA. .,Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA.
| |
Collapse
|
22
|
Paneru B, Ali A, Al-Tobasei R, Kenney B, Salem M. Crosstalk among lncRNAs, microRNAs and mRNAs in the muscle 'degradome' of rainbow trout. Sci Rep 2018; 8:8416. [PMID: 29849185 PMCID: PMC5976669 DOI: 10.1038/s41598-018-26753-2] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2018] [Accepted: 05/18/2018] [Indexed: 01/17/2023] Open
Abstract
In fish, protein-coding and noncoding genes involved in muscle atrophy are not fully characterized. In this study, we characterized coding and noncoding genes involved in gonadogenesis-associated muscle atrophy, and investigated the potential functional interplay between these genes. Using RNA-Seq, we compared expression pattern of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs of atrophying skeletal muscle from gravid females and control skeletal muscle from age-matched sterile individuals. A total of 852 mRNAs, 1,160 lncRNAs and 28 microRNAs were differentially expressed (DE) between the two groups. Muscle atrophy appears to be mediated by many genes encoding ubiquitin-proteasome system, autophagy related proteases, lysosomal proteases and transcription factors. Transcripts encoding atrogin-1 and mir-29 showed exceptional high expression in atrophying muscle, suggesting an important role in bulk muscle proteolysis. DE genes were co-localized in the genome with strong expression correlation, and they exhibited extensive 'lncRNA-mRNA', 'lncRNA-microRNA', 'mRNA-microRNA' and 'lncRNA-protein' physical interactions. DE genes exhibiting potential functional interactions comprised the highly correlated 'lncRNA-mRNA-microRNA' gene network described as 'degradome'. This study pinpoints extensive coding and noncoding RNA interactions during muscle atrophy in fish, and provides valuable resources for future mechanistic studies.
Collapse
Affiliation(s)
- Bam Paneru
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA
| | - Ali Ali
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA
| | - Rafet Al-Tobasei
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA.,Department of Biostatistics, University of Alabama at Birmingham, Birmingham, AL, 35294-0022, USA
| | - Brett Kenney
- Division of Animal and Nutritional Science, West Virginia University, Morgantown, 26506-6108, West Virginia, USA
| | - Mohamed Salem
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA. .,Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA.
| |
Collapse
|
23
|
Abernathy J, Overturf K. Expression of Antisense Long Noncoding RNAs as Potential Regulators in Rainbow Trout with Different Tolerance to Plant-Based Diets. Anim Biotechnol 2018; 30:87-94. [PMID: 29300121 DOI: 10.1080/10495398.2017.1401546] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Reformulation of aquafeeds in salmonid diets to include more plant proteins is critical for sustainable aquaculture. However, increasing plant proteins can lead to stunted growth and enteritis. Toward an understanding of the regulatory mechanisms behind plant protein utilization, directional RNA sequencing of liver tissues from a rainbow trout strain selected for growth on an all plant-protein diet and a control strain, both fed a plant diet for 12 weeks, were utilized to construct long noncoding RNAs. Antisense long noncoding RNAs were selected for differential expression and functional analyses since they have been shown to have regulatory actions within a genome. A total of 142 unique antisense long noncoding RNAs were differentially expressed between strains, 60 of which could be mapped to a gene. Genes underlying these noncoding RNAs are indicated in lipid metabolism and immunity. Six noncoding transcripts were also found to overlap with differentially expressed protein-coding genes, all of which were co-expressed. Associating variation in regulatory elements between rainbow trout strains with differing tolerance to plant-protein diets will assist in future studies toward increased gains throughout carnivorous aquaculture.
Collapse
Affiliation(s)
- Jason Abernathy
- a USDA, Agricultural Research Service , Harry K. Dupree Stuttgart National Aquaculture Research Center , Stuttgart , AR , USA
| | - Ken Overturf
- b USDA, Agricultural Research Service , Hagerman Fish Culture Experiment Station , Hagerman , ID , USA
| |
Collapse
|
24
|
Rawal HC, Kumar S, Mithra S V A, Solanke AU, Nigam D, Saxena S, Tyagi A, V S, Yadav NR, Kalia P, Singh NP, Singh NK, Sharma TR, Gaikwad K. High Quality Unigenes and Microsatellite Markers from Tissue Specific Transcriptome and Development of a Database in Clusterbean (Cyamopsis tetragonoloba, L. Taub). Genes (Basel) 2017; 8:genes8110313. [PMID: 29120386 PMCID: PMC5704226 DOI: 10.3390/genes8110313] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Revised: 10/23/2017] [Accepted: 11/06/2017] [Indexed: 12/23/2022] Open
Abstract
Clusterbean (Cyamopsis tetragonoloba L. Taub), is an important industrial, vegetable and forage crop. This crop owes its commercial importance to the presence of guar gum (galactomannans) in its endosperm which is used as a lubricant in a range of industries. Despite its relevance to agriculture and industry, genomic resources available in this crop are limited. Therefore, the present study was undertaken to generate RNA-Seq based transcriptome from leaf, shoot, and flower tissues. A total of 145 million high quality Illumina reads were assembled using Trinity into 127,706 transcripts and 48,007 non-redundant high quality (HQ) unigenes. We annotated 79% unigenes against Plant Genes from the National Center for Biotechnology Information (NCBI), Swiss-Prot, Pfam, gene ontology (GO) and KEGG databases. Among the annotated unigenes, 30,020 were assigned with 116,964 GO terms, 9984 with EC and 6111 with 137 KEGG pathways. At different fragments per kilobase of transcript per millions fragments sequenced (FPKM) levels, genes were found expressed higher in flower tissue followed by shoot and leaf. Additionally, we identified 8687 potential simple sequence repeats (SSRs) with an average frequency of one SSR per 8.75 kb. A total of 28 amplified SSRs in 21 clusterbean genotypes resulted in polymorphism in 13 markers with average polymorphic information content (PIC) of 0.21. We also constructed a database named ‘ClustergeneDB’ for easy retrieval of unigenes and the microsatellite markers. The tissue specific genes identified and the molecular marker resources developed in this study is expected to aid in genetic improvement of clusterbean for its end use.
Collapse
Affiliation(s)
- Hukam C Rawal
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Shrawan Kumar
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Amitha Mithra S V
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Amolkumar U Solanke
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Deepti Nigam
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Swati Saxena
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Anshika Tyagi
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Sureshkumar V
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Neelam R Yadav
- Department of Biotechnology and Molecular Biology, CCS Haryana Agricultural University, Hisar 125004, India.
| | - Pritam Kalia
- ICAR-Indian Agricultural Research Institute, New Delhi 110012, India.
| | | | | | - Tilak Raj Sharma
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| | - Kishor Gaikwad
- ICAR-National Research Centre on Plant Biotechnology, New Delhi 110012, India.
| |
Collapse
|
25
|
Oomen RA, Hutchings JA. Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing. Facets (Ott) 2017. [DOI: 10.1139/facets-2017-0015] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.
Collapse
Affiliation(s)
- Rebekah A. Oomen
- Department of Biology, Dalhousie University, Halifax, NS B3H 4J1, Canada
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, 0371 Oslo, Norway
- Institute of Marine Research, Flødevigen Research Station, 4817 His, Norway
| | - Jeffrey A. Hutchings
- Department of Biology, Dalhousie University, Halifax, NS B3H 4J1, Canada
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, 0371 Oslo, Norway
- Institute of Marine Research, Flødevigen Research Station, 4817 His, Norway
- Department of Natural Sciences, University of Agder, 4604 Kristiansand, Norway
| |
Collapse
|
26
|
Paneru BD, Al-Tobasei R, Kenney B, Leeds TD, Salem M. RNA-Seq reveals MicroRNA expression signature and genetic polymorphism associated with growth and muscle quality traits in rainbow trout. Sci Rep 2017; 7:9078. [PMID: 28831113 PMCID: PMC5567286 DOI: 10.1038/s41598-017-09515-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Accepted: 07/25/2017] [Indexed: 02/01/2023] Open
Abstract
The role of microRNA expression and genetic variation in microRNA-binding sites of target genes on growth and muscle quality traits is poorly characterized. We used RNA-Seq approach to investigate their importance on 5 growth and muscle quality traits: whole body weight (WBW), muscle yield, muscle crude-fat content, muscle shear force and whiteness. Phenotypic data were collected from 471 fish, representing 98 families (~5 fish/family) from a growth-selected line. Muscle microRNAs and mRNAs were sequenced from 22 families showing divergent phenotypes. Ninety microRNAs showed differential expression between families with divergent phenotypes, and their expression was strongly associated with variation in phenotypes. A total of 204 single nucleotide polymorphisms (SNPs) present in 3′ UTR of target genes either destroyed or created novel illegitimate microRNA target sites; of them, 78 SNPs explained significant variation in the aforementioned 5 muscle traits. Majority of the phenotype-associated SNPs were present in microRNA-binding sites of genes involved in energy metabolism and muscle structure. These findings suggest that variation in microRNA expression and/or sequence variation in microRNA binding sites in target genes play an important role in mediating differences in fish growth and muscle quality phenotypes.
Collapse
Affiliation(s)
- Bam Dev Paneru
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States
| | - Rafet Al-Tobasei
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States
| | - Brett Kenney
- Division of Animal and Nutritional Science, West Virginia University, Morgantown, 26506-6108, West Virginia, United States
| | - Timothy D Leeds
- The National Center for Cool and Cold Water Aquaculture, USDA Agricultural Research Service, Kearneysville, WV, 25430, United States
| | - Mohamed Salem
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States. .,Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States.
| |
Collapse
|
27
|
Al-Tobasei R, Ali A, Leeds TD, Liu S, Palti Y, Kenney B, Salem M. Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout. BMC Genomics 2017; 18:582. [PMID: 28784089 PMCID: PMC5547479 DOI: 10.1186/s12864-017-3992-z] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2016] [Accepted: 08/01/2017] [Indexed: 12/22/2022] Open
Abstract
Background Coding/functional SNPs change the biological function of a gene and, therefore, could serve as “large-effect” genetic markers. In this study, we used two bioinformatics pipelines, GATK and SAMtools, for discovering coding/functional SNPs with allelic-imbalances associated with total body weight, muscle yield, muscle fat content, shear force, and whiteness. Phenotypic data were collected for approximately 500 fish, representing 98 families (5 fish/family), from a growth-selected line, and the muscle transcriptome was sequenced from 22 families with divergent phenotypes (4 low- versus 4 high-ranked families per trait). Results GATK detected 59,112 putative SNPs; of these SNPs, 4798 showed allelic imbalances (>2.0 as an amplification and <0.5 as loss of heterozygosity). SAMtools detected 87,066 putative SNPs; and of them, 4962 had allelic imbalances between the low- and high-ranked families. Only 1829 SNPs with allelic imbalances were common between the two datasets, indicating significant differences in algorithms. The two datasets contained 7930 non-redundant SNPs of which 4439 mapped to 1498 protein-coding genes (with 6.4% non-synonymous SNPs) and 684 mapped to 295 lncRNAs. Validation of a subset of 92 SNPs revealed 1) 86.7–93.8% success rate in calling polymorphic SNPs and 2) 95.4% consistent matching between DNA and cDNA genotypes indicating a high rate of identifying SNPs with allelic imbalances. In addition, 4.64% SNPs revealed random monoallelic expression. Genome distribution of the SNPs with allelic imbalances exhibited high density for all five traits in several chromosomes, especially chromosome 9, 20 and 28. Most of the SNP-harboring genes were assigned to important growth-related metabolic pathways. Conclusion These results demonstrate utility of RNA-Seq in assessing phenotype-associated allelic imbalances in pooled RNA-Seq samples. The SNPs identified in this study were included in a new SNP-Chip design (available from Affymetrix) for genomic and genetic analyses in rainbow trout. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3992-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Rafet Al-Tobasei
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA
| | - Ali Ali
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA
| | - Timothy D Leeds
- National Center for Cool and Cold Water Aquaculture, ARS-USDA, Kearneysville, WV, 25430, USA
| | - Sixin Liu
- National Center for Cool and Cold Water Aquaculture, ARS-USDA, Kearneysville, WV, 25430, USA
| | - Yniv Palti
- National Center for Cool and Cold Water Aquaculture, ARS-USDA, Kearneysville, WV, 25430, USA
| | - Brett Kenney
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506, USA
| | - Mohamed Salem
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA. .,Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, USA.
| |
Collapse
|
28
|
Abernathy J, Brezas A, Snekvik KR, Hardy RW, Overturf K. Integrative functional analyses using rainbow trout selected for tolerance to plant diets reveal nutrigenomic signatures for soy utilization without the concurrence of enteritis. PLoS One 2017; 12:e0180972. [PMID: 28723948 PMCID: PMC5517010 DOI: 10.1371/journal.pone.0180972] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2016] [Accepted: 06/24/2017] [Indexed: 12/13/2022] Open
Abstract
Finding suitable alternative protein sources for diets of carnivorous fish species remains a major concern for sustainable aquaculture. Through genetic selection, we created a strain of rainbow trout that outperforms parental lines in utilizing an all-plant protein diet and does not develop enteritis in the distal intestine, as is typical with salmonids on long-term plant protein-based feeds. By incorporating this strain into functional analyses, we set out to determine which genes are critical to plant protein utilization in the absence of gut inflammation. After a 12-week feeding trial with our selected strain and a control trout strain fed either a fishmeal-based diet or an all-plant protein diet, high-throughput RNA sequencing was completed on both liver and muscle tissues. Differential gene expression analyses, weighted correlation network analyses and further functional characterization were performed. A strain-by-diet design revealed differential expression ranging from a few dozen to over one thousand genes among the various comparisons and tissues. Major gene ontology groups identified between comparisons included those encompassing central, intermediary and foreign molecule metabolism, associated biosynthetic pathways as well as immunity. A systems approach indicated that genes involved in purine metabolism were highly perturbed. Systems analysis among the tissues tested further suggests the interplay between selection for growth, dietary utilization and protein tolerance may also have implications for nonspecific immunity. By combining data from differential gene expression and co-expression networks using selected trout, along with ontology and pathway analyses, a set of 63 candidate genes for plant diet tolerance was found. Risk loci in human inflammatory bowel diseases were also found in our datasets, indicating rainbow trout selected for plant-diet tolerance may have added utility as a potential biomedical model.
Collapse
Affiliation(s)
- Jason Abernathy
- Hagerman Fish Culture Experiment Station, USDA-ARS, Hagerman, Idaho, United States of America
| | - Andreas Brezas
- Aquaculture Research Institute, University of Idaho, Hagerman, Idaho, United States of America
| | - Kevin R. Snekvik
- Department of Veterinary Microbiology and Pathology, Washington State University, Pullman, Washington, United States of America
| | - Ronald W. Hardy
- Aquaculture Research Institute, University of Idaho, Hagerman, Idaho, United States of America
| | - Ken Overturf
- Hagerman Fish Culture Experiment Station, USDA-ARS, Hagerman, Idaho, United States of America
| |
Collapse
|
29
|
Seasonal differences in the testicular transcriptome profile of free-living European beavers (Castor fiber L.) determined by the RNA-Seq method. PLoS One 2017; 12:e0180323. [PMID: 28678806 PMCID: PMC5498055 DOI: 10.1371/journal.pone.0180323] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2017] [Accepted: 06/14/2017] [Indexed: 11/20/2022] Open
Abstract
The European beaver (Castor fiber L.) is an important free-living rodent that inhabits Eurasian temperate forests. Beavers are often referred to as ecosystem engineers because they create or change existing habitats, enhance biodiversity and prepare the environment for diverse plant and animal species. Beavers are protected in most European Union countries, but their genomic background remains unknown. In this study, gene expression patterns in beaver testes and the variations in genetic expression in breeding and non-breeding seasons were determined by high-throughput transcriptome sequencing. Paired-end sequencing in the Illumina HiSeq 2000 sequencer produced a total of 373.06 million of high-quality reads. De novo assembly of contigs yielded 130,741 unigenes with an average length of 1,369.3 nt, N50 value of 1,734, and average GC content of 46.51%. A comprehensive analysis of the testicular transcriptome revealed more than 26,000 highly expressed unigenes which exhibited the highest homology with Rattus norvegicus and Ictidomys tridecemlineatus genomes. More than 8,000 highly expressed genes were found to be involved in fundamental biological processes, cellular components or molecular pathways. The study also revealed 42 genes whose regulation differed between breeding and non-breeding seasons. During the non-breeding period, the expression of 37 genes was up-regulated, and the expression of 5 genes was down-regulated relative to the breeding season. The identified genes encode molecules which are involved in signaling transduction, DNA repair, stress responses, inflammatory processes, metabolism and steroidogenesis. Our results pave the way for further research into season-dependent variations in beaver testes.
Collapse
|
30
|
Zhu C, Pan Z, Wang H, Chang G, Wu N, Ding H. De novo assembly, characterization and annotation for the transcriptome of Sarcocheilichthys sinensis. PLoS One 2017; 12:e0171966. [PMID: 28196101 PMCID: PMC5308828 DOI: 10.1371/journal.pone.0171966] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2016] [Accepted: 01/28/2017] [Indexed: 11/30/2022] Open
Abstract
The Chinese lake gudgeon Sarcocheilichthys sinensis is a small cyprinid fish with great aquaculture potential both for its edible and ornamental values. Nevertheless, available genomic and transcriptomic information for this fish is extremely deficient. In this study, a normalized cDNA library was constructed using 13 mixed tissues of an adult male S. sinensis, and was sequenced by the Illumina HiSeq2500 platform. De novo assembly was performed using 38,911,511 obtained clean reads, and a total of 147,282 unigenes with an average length of 900 bp were finally achieved. 96.2% of these unigenes were annotated in 9 public databases, and 16 segments of growth-related genes were identified for future studies. In addition, 28,493 unigenes were assigned to 61 subcategories of Gene Ontology (GO), and 10,483 unigenes were assigned to 25 categories of Cluster of Orthologous Group (COG). Moreover, 14,943 unigenes were classified into 225 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 30,666 microsatellites were detected from 17,627 unigenes with an average distribution density of 1:2405 bp. This transcriptome data set will be valuable for researches on discovery, expression and evolution on genes of interest. Meanwhile, the identified microsatellites would be useful tools for genetic and genomic studies in S. sinensis.
Collapse
Affiliation(s)
- Chuankun Zhu
- Jiangsu Engineering Laboratory for Breeding of Special Aquatic Organisms, Huaiyin Normal University, Huai’an, China
- Jiangsu Collaborative Innovation Center of Regional Modern Agriculture & Environmental Protection, Huaiyin Normal University, Huai’an, China
- * E-mail:
| | - Zhengjun Pan
- Jiangsu Engineering Laboratory for Breeding of Special Aquatic Organisms, Huaiyin Normal University, Huai’an, China
- Jiangsu Collaborative Innovation Center of Regional Modern Agriculture & Environmental Protection, Huaiyin Normal University, Huai’an, China
| | - Hui Wang
- Jiangsu Engineering Laboratory for Breeding of Special Aquatic Organisms, Huaiyin Normal University, Huai’an, China
- Jiangsu Collaborative Innovation Center of Regional Modern Agriculture & Environmental Protection, Huaiyin Normal University, Huai’an, China
| | - Guoliang Chang
- Jiangsu Engineering Laboratory for Breeding of Special Aquatic Organisms, Huaiyin Normal University, Huai’an, China
- Jiangsu Collaborative Innovation Center of Regional Modern Agriculture & Environmental Protection, Huaiyin Normal University, Huai’an, China
| | - Nan Wu
- Jiangsu Engineering Laboratory for Breeding of Special Aquatic Organisms, Huaiyin Normal University, Huai’an, China
- Jiangsu Collaborative Innovation Center of Regional Modern Agriculture & Environmental Protection, Huaiyin Normal University, Huai’an, China
| | - Huaiyu Ding
- Jiangsu Engineering Laboratory for Breeding of Special Aquatic Organisms, Huaiyin Normal University, Huai’an, China
- Jiangsu Collaborative Innovation Center of Regional Modern Agriculture & Environmental Protection, Huaiyin Normal University, Huai’an, China
| |
Collapse
|
31
|
Paneru B, Al-Tobasei R, Palti Y, Wiens GD, Salem M. Differential expression of long non-coding RNAs in three genetic lines of rainbow trout in response to infection with Flavobacterium psychrophilum. Sci Rep 2016; 6:36032. [PMID: 27786264 PMCID: PMC5081542 DOI: 10.1038/srep36032] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2016] [Accepted: 10/10/2016] [Indexed: 12/15/2022] Open
Abstract
Bacterial cold-water disease caused by Flavobacterium psychrophilum is one of the major causes of mortality of salmonids. Three genetic lines of rainbow trout designated as ARS-Fp-R (resistant), ARS-Fp-C (control) and ARS-Fp-S (susceptible) have significant differences in survival rate following F. psychrophilum infection. Previous study identified transcriptome differences of immune-relevant protein-coding genes at basal and post infection levels among these genetic lines. Using RNA-Seq approach, we quantified differentially expressed (DE) long non-coding RNAs (lncRNAs) in response to F. psychrophilum challenge in these genetic lines. Pairwise comparison between genetic lines and different infection statuses identified 556 DE lncRNAs. A positive correlation existed between the number of the differentially regulated lncRNAs and that of the protein-coding genes. Several lncRNAs showed strong positive and negative expression correlation with their overlapped, neighboring and distant immune related protein-coding genes including complement components, cytokines, chemokines and several signaling molecules involved in immunity. The correlated expressions and genome-wide co-localization suggested that some lncRNAs may be involved in regulating immune-relevant protein-coding genes. This study provides the first evidence of lncRNA-mediated regulation of the anti-bacterial immune response in a commercially important aquaculture species and will likely help developing new genetic markers for rainbow trout disease resistance.
Collapse
Affiliation(s)
- Bam Paneru
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, U.S
| | - Rafet Al-Tobasei
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, U.S
| | - Yniv Palti
- The National Center for Cool and Cold Water Aquaculture, USDA Agricultural Research Service, Kearneysville, WV 25430, U.S
| | - Gregory D Wiens
- The National Center for Cool and Cold Water Aquaculture, USDA Agricultural Research Service, Kearneysville, WV 25430, U.S
| | - Mohamed Salem
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, U.S.,Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, U.S
| |
Collapse
|
32
|
Osachoff HL, Brown LL, Tirrul L, van Aggelen GC, Brinkman FS, Kennedy CJ. Time course of hepatic gene expression and plasma vitellogenin protein concentrations in estrone-exposed juvenile rainbow trout (Oncorhynchus mykiss). COMPARATIVE BIOCHEMISTRY AND PHYSIOLOGY D-GENOMICS & PROTEOMICS 2016; 19:112-119. [DOI: 10.1016/j.cbd.2016.02.002] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2015] [Revised: 02/08/2016] [Accepted: 02/21/2016] [Indexed: 12/31/2022]
|
33
|
Hébert FO, Grambauer S, Barber I, Landry CR, Aubin-Horth N. Transcriptome sequences spanning key developmental states as a resource for the study of the cestode Schistocephalus solidus, a threespine stickleback parasite. Gigascience 2016; 5:24. [PMID: 27259971 PMCID: PMC4891850 DOI: 10.1186/s13742-016-0128-3] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Accepted: 05/11/2016] [Indexed: 01/06/2023] Open
Abstract
Background Schistocephalus solidus is a well-established model organism for studying the complex life cycle of cestodes and the mechanisms underlying host-parasite interactions. However, very few large-scale genetic resources for this species are available. We have sequenced and de novo-assembled the transcriptome of S. solidus using tissues from whole worms at three key developmental states - non-infective plerocercoid, infective plerocercoid and adult plerocercoid - to provide a resource for studying the evolution of complex life cycles and, more specifically, how parasites modulate their interactions with their hosts during development. Findings The de novo transcriptome assembly reconstructed the coding sequence of 10,285 high-confidence unigenes from which 24,765 non-redundant transcripts were derived. 7,920 (77 %) of these unigenes were annotated with a protein name and 7,323 (71 %) were assigned at least one Gene Ontology term. Our raw transcriptome assembly (unfiltered transcripts) covers 92 % of the predicted transcriptome derived from the S. solidus draft genome assembly currently available on WormBase. It also provides new ecological information and orthology relationships to further annotate the current WormBase transcriptome and genome. Conclusion This large-scale transcriptomic dataset provides a foundation for studies on how parasitic species with complex life cycles modulate their response to changes in biotic and abiotic conditions experienced inside their various hosts, which is a fundamental objective of parasitology. Furthermore, this resource will help in the validation of the S solidus gene features that have been predicted based on genomic sequence. Electronic supplementary material The online version of this article (doi:10.1186/s13742-016-0128-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- François Olivier Hébert
- Institut de Biologie Intégrative et des Systèmes (IBIS), Département de Biologie, Université Laval, Pavillon Charles-Eugène-Marchand, Québec, G1V 0A6, Canada.
| | - Stephan Grambauer
- Department of Neuroscience, Psychology and Behaviour, Adrian Building, University of Leicester, University Road, Leicester, LE1 7RH, UK
| | - Iain Barber
- Department of Neuroscience, Psychology and Behaviour, Adrian Building, University of Leicester, University Road, Leicester, LE1 7RH, UK
| | - Christian R Landry
- Institut de Biologie Intégrative et des Systèmes (IBIS), Département de Biologie, Université Laval, Pavillon Charles-Eugène-Marchand, Québec, G1V 0A6, Canada
| | - Nadia Aubin-Horth
- Institut de Biologie Intégrative et des Systèmes (IBIS), Département de Biologie, Université Laval, Pavillon Charles-Eugène-Marchand, Québec, G1V 0A6, Canada
| |
Collapse
|
34
|
Wang J, Fu L, Koganti PP, Wang L, Hand JM, Ma H, Yao J. Identification and Functional Prediction of Large Intergenic Noncoding RNAs (lincRNAs) in Rainbow Trout (Oncorhynchus mykiss). MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2016; 18:271-282. [PMID: 26864089 DOI: 10.1007/s10126-016-9689-5] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/24/2015] [Accepted: 11/30/2015] [Indexed: 06/05/2023]
Abstract
Long noncoding RNAs (lncRNAs) have been recognized in recent years as key regulators of diverse cellular processes. Genome-wide large-scale projects have uncovered thousands of lncRNAs in many model organisms. Large intergenic noncoding RNAs (lincRNAs) are lncRNAs that are transcribed from intergenic regions of genomes. To date, no lincRNAs in non-model teleost fish have been reported. In this report, we present the first reference catalog of 9674 rainbow trout lincRNAs based on analysis of RNA-Seq data from 15 tissues. Systematic analysis revealed that lincRNAs in rainbow trout share many characteristics with those in other mammalian species. They are shorter and lower in exon number and expression level compared with protein-coding genes. They show tissue-specific expression pattern and are typically co-expressed with their neighboring genes. Co-expression network analysis suggested that many lincRNAs are associated with immune response, muscle differentiation, and neural development. The study provides an opportunity for future experimental and computational studies to uncover the functions of lincRNAs in rainbow trout.
Collapse
Affiliation(s)
- Jian Wang
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506-6108, USA
| | - Liyuan Fu
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506-6108, USA
| | - Prasanthi P Koganti
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506-6108, USA
| | - Lei Wang
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506-6108, USA
| | - Jacqelyn M Hand
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506-6108, USA
| | - Hao Ma
- USDA/ARS National Center for Cool and Cold Water Aquaculture, Kearneysville, WV, 25430, USA
| | - Jianbo Yao
- Division of Animal and Nutritional Sciences, West Virginia University, Morgantown, WV, 26506-6108, USA.
| |
Collapse
|
35
|
Al-Tobasei R, Paneru B, Salem M. Genome-Wide Discovery of Long Non-Coding RNAs in Rainbow Trout. PLoS One 2016; 11:e0148940. [PMID: 26895175 PMCID: PMC4764514 DOI: 10.1371/journal.pone.0148940] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2015] [Accepted: 01/23/2016] [Indexed: 02/07/2023] Open
Abstract
The ENCODE project revealed that ~70% of the human genome is transcribed. While only 1–2% of the RNAs encode for proteins, the rest are non-coding RNAs. Long non-coding RNAs (lncRNAs) form a diverse class of non-coding RNAs that are longer than 200nt. Emerging evidence indicates that lncRNAs play critical roles in various cellular processes including regulation of gene expression. LncRNAs show low levels of gene expression and sequence conservation, which make their computational identification in genomes difficult. In this study, more than two billion Illumina sequence reads were mapped to the genome reference using the TopHat and Cufflinks software. Transcripts shorter than 200nt, with more than 83–100 amino acids ORF, or with significant homologies to the NCBI nr-protein database were removed. In addition, a computational pipeline was used to filter the remaining transcripts based on a protein-coding-score test. Depending on the filtering stringency conditions, between 31,195 and 54,503 lncRNAs were identified, with only 421 matching known lncRNAs in other species. A digital gene expression atlas revealed 2,935 tissue-specific and 3,269 ubiquitously-expressed lncRNAs. This study annotates the lncRNA rainbow trout genome and provides a valuable resource for functional genomics research in salmonids.
Collapse
Affiliation(s)
- Rafet Al-Tobasei
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States of America
| | - Bam Paneru
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States of America
| | - Mohamed Salem
- Computational Science Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States of America
- Department of Biology and Molecular Biosciences Program, Middle Tennessee State University, Murfreesboro, TN, 37132, United States of America
- * E-mail:
| |
Collapse
|
36
|
Deng T, Pang C, Lu X, Zhu P, Duan A, Tan Z, Huang J, Li H, Chen M, Liang X. De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery. PLoS One 2016; 11:e0147132. [PMID: 26766209 PMCID: PMC4713091 DOI: 10.1371/journal.pone.0147132] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2015] [Accepted: 12/29/2015] [Indexed: 01/11/2023] Open
Abstract
The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.
Collapse
Affiliation(s)
- Tingxian Deng
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Chunying Pang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Xingrong Lu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Peng Zhu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Anqin Duan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Zhengzhun Tan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Jian Huang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Hui Li
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Mingtan Chen
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Xianwei Liang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
- * E-mail:
| |
Collapse
|
37
|
|