Search

Viral nervous necrosis (VNN) disease caused by the nervous necrosis virus (NNV) is a major disease, leading to a huge economic loss in aquaculture. Previous GWAS and QTL mapping have identified a major QTL for NNV resistance in linkage group 20 in Asian seabass. However, no causative gene for NNV resistance has been identified. In this study, RNA-seq from brains of Asian seabass fingerlings challenged with NNV at four time points (5, 10, 15 and 20 days post-challenge) identified 1228, 245, 189 and 134 DEGs, respectively. Eight DEGs, including rrm1, were located in the major QTL for NNV resistance. An association study in 445 survived and 608 dead fingerlings after NNV challenge revealed that the SNP in rrm1 were significantly associated with NNV resistance. Therefore, rrm1 was selected for functional analysis, as a candidate gene for NNV resistance. The expression of rrm1 was significantly increased in the gill, liver, spleen and muscle, and was suppressed in the brain, gut and skin after NNV challenge. The rrm1 protein was localized in the nuclear membrane. Over-expression of rrm1 significantly decreased viral RNA and titer in NNV-infected Asian seabass cells, whereas knock-down of rrm1 significantly increased viral RNA and titer in NNV-infected Asian seabass cells. The rrm1 knockout heterozygous zebrafish was more susceptible to NNV infection. Our study suggests that rrm1 is one of the causative genes for NNV resistance and the SNP in the gene may be applied for accelerating genetic improvement for NNV resistance.

Understanding the functions of genes related to disease resistance and identifying polymorphisms in these genes are essential in molecular breeding for disease resistance. Viral nervous necrosis (VNN) is one of the major diseases in the Asian seabass, Lates calcarifer. Our previous works on QTL mapping, GWAS and cell-line transcriptome analysis of the Asian seabass after NNV challenge revealed that the gene GAB3 might be a candidate gene for VNN resistance. In this study, we cloned and characterized GAB3, and identified SNPs in the gene of the Asian seabass. The cDNA of the gene was 2165 bp, containing an ORF of 1674 bp encoding 557 amino acids. The gene consisted of 10 exons and nine introns. It was ubiquitously expressed in normal fish. An analysis of the association between two SNPs in the second intron and NNV resistance in 1035 fish descended from 43 families revealed that the two SNPs were significantly associated with VNN resistance. After NNV infection, the expression of GAB3 was significantly increased in the brain, spleen, muscle and gut, and was suppressed in the liver. The GAB3 protein was localized in the nucleus. Overexpression of GAB3 with specific GAB3-pcDNA was positively correlated to increased viral RNA and titer in NNV-infected Asian seabass cells. Our study provides new evidence to support that GAB3 may be an important gene related to NNV resistance. In addition, the SNPs provide DNA markers for the selection of candidate genes resistance to NNV at the juvenile stage of Asian seabass.

Tilapias are important aquaculture species. Male-biased sexual size dimorphism is very common and males are preferred for aquaculture in tilapia. However, the mechanisms underlying sexual dimorphism remain to be elucidated. One hundred and thirty-six sex-biased genes, of which 23 were male-biased and 113 were female-biased, were identified via reanalysis of a muscle transcriptome data using the latest reference genome assembly. These genes were mapped to KEGG pathways that are related to somatic cell metabolism, growth and differentiations, such as MAPK, FoxO4 and metabolism pathways as well as developmental processes responsible for skeletal muscle development. Pathways related to cell growth and proliferations such as MAPK signaling pathways are upregulated in males while pathways regulating cell division such as FoxO4 are upregulated in females. Sex-biased genes in tilapia skeletal muscle have higher evolution rates (dN/dS) compared to unbiased genes. Female-biased and male-biased genes showed 17.4% and 13.5% higher dN/dS, respectively, compared to unbiased genes. Our results suggest that some of the male and female sex-biased genes were under selection pressures. Three SNPs located in the promoter region of one sex-biased gene RASGRF1 on LG1 were associated with bodyweight differences in the hybrid tilapia. These sex-biased genes identified in this study may serve as candidate genes for future functional analysis on sexual size dimorphism and for developing DNA markers for selecting fast-growing saline tilapia in aquaculture. Keywords: Sexual dimorphism, Body size, RNA-Seq, Tilapia, Aquaculture, Evolution

14 pages, 4 figures, supplementary material https://doi.org/10.1093/molbev/msab110.-- Data Availability: The genome sequences and related annotations of fighting fish are hosted by the web server of Temasek Lifesciences Laboratory (https://genhua.tll.org.sg/) and archived in China National GeneBank (CNGB) and DDBJ with BioProject accession nos. CNP0001745 and PRJDB7253, respectively. Sequences used for whole-genome sequencing, RNA sequencing, and RAD sequencing are available with the DDBJ Sequencing Read Archive (SRA) through BioProject ID PRJDB7253–PRJDB7255, Resolving the genomic basis underlying phenotypic variations is a question of great importance in evolutionary biology. However, understanding how genotypes determine the phenotypes is still challenging. Centuries of artificial selective breeding for beauty and aggression resulted in a plethora of colors, long-fin varieties, and hyper-aggressive behavior in the air-breathing Siamese fighting fish (Betta splendens), supplying an excellent system for studying the genomic basis of phenotypic variations. Combining whole-genome sequencing, quantitative trait loci mapping, genome-wide association studies, and genome editing, we investigated the genomic basis of huge morphological variation in fins and striking differences in coloration in the fighting fish. Results revealed that the double tail, elephant ear, albino, and fin spot mutants each were determined by single major-effect loci. The elephant ear phenotype was likely related to differential expression of a potassium ion channel gene, kcnh8. The albinotic phenotype was likely linked to a cis-regulatory element acting on the mitfa gene and the double-tail mutant was suggested to be caused by a deletion in a zic1/zic4 coenhancer. Our data highlight that major loci and cis-regulatory elements play important roles in bringing about phenotypic innovations and establish Bettas as new powerful model to study the genomic basis of evolved changes, This project was supported by the Temasek Lifesciences Laboratory, Singapore, With the institutional support of the ‘Severo OchoaCentre of Excellence’ accreditation (CEX2019-000928-S)

Fish color is an important trait involved in many biological functions and related to commercial value of fish. However, our understanding of the molecular mechanisms underlying color variations in fish species is limited. In this study, we found 12 natural golden mutants in a Mozambique tilapia population. Analysis of the inheritance pattern of the golden and blackish skin colors in nine backcross families revealed that the golden phenotype is a recessive mutant determined by a single locus and follows simple Mendelian inheritance. Using linkage mapping, fine mapping and analyses of gene expressions, we uncovered that the Pmel17 gene was responsible for golden and black skin colors. Further sequence analysis of the gene discovered that an insertion in 3′ UTR of Pmel17 was fully associated with the golden skin color in 1835 tilapias. Expression analysis showed that Pmel17 with the insertion reduced the transcripts of Pmel17. These results indicate that the insertion in 3′ UTR of Pmel17 is responsible for the golden skin color in tilapia. The insertion can be used as a reliable DNA marker to accelerate the breeding for golden color.

Background: Palm oil is an important feedstock for biofuel. Palm oil yield is seriously affected by drought stress. However, not much is known about the molecular responses of oil palm to drought stress.Results: We studied the root transcriptomic responses of oil palm seedlings under normal culture and 14-day drought stress using RNA-seq and bioinformatics analysis. We identified 1293 differentially expressed genes (DEGs), involved in several molecular processes, such as cell wall biogenesis and functions, phenylpropanoid biosynthesis and metabolisms, ion transport and homeostasis and cellular ketone metabolic process, as well as small molecule biosynthetic process. We observed that DEGs were significantly enriched into the two categories: hormone regulation and metabolism, as well as ABC transporters. In addition, we identified three protein-protein interaction networks involved in the response to drought stress, including ion transport, reactive nitrogen species metabolic process and nitrate assimilation. Finally, 96 differentially expressed transcription factors were detected to be associated with drought stress response, which were classified into 28 families.Conclusions : The transcriptomic responses of oil palm seedlings to drought stress were systematically analysed, revealing important genes, pathways, networks and transcription factors involved in drought stress responses. These results provide new insights into the mechanisms of drought stress responses in economic crops. The genes and pathways identified in this study provide valuable genomic resources to improve drought tolerance of oil palm by both genetic modification and selective breeding.

Tilapia is one of most important foodfish species. The low omega-3 to omega-6 fatty acid ratio in freshwater tilapia meat is disadvantageous for human health. Increasing omega-3 content is an important task in breeding to increase the nutritional value of tilapia. However, conventional breeding to increase omega-3 content is difficult and slow. To accelerate the increase of omega-3 through marker-assisted selection (MAS), we conducted QTL mapping for fatty acid contents and profiles in a F2 family of saline tilapia generated by crossing red tilapia and Mozambique tilapia. The total omega-3 content in F2 hybrid tilapia was 2.5 ± 1.0 mg/g, higher than that (2.00 mg/g) in freshwater tilapia. Genotyping by sequencing (GBS) technology was used to discover and genotype SNP markers, and microsatellites were also genotyped. We constructed a linkage map with 784 markers (151 microsatellites and 633 SNPs). The linkage map was 2076.7 cM long and consisted of 22 linkage groups. Significant and suggestive QTL for total lipid content were mapped on six linkage groups (LG3, -4, -6, -8, -13, and -15) and explained 5.8–8.3% of the phenotypic variance. QTL for omega-3 fatty acids were located on four LGs (LG11, -18, -19, and -20) and explained 5.0 to 7.5% of the phenotypic variance. Our data suggest that the total lipid and omega-3 fatty acid content were determined by multiple genes in tilapia. The markers flanking the QTL for omega-3 fatty acids can be used in MAS to accelerate the genetic improvements of these traits in salt-tolerant tilapia. NRF (Natl Research Foundation, S’pore)

Breeding for disease resistance is an important task for genetic improvement in aquaculture. Identification of DNA markers associated with disease resistance is a critical step towards molecular breeding for disease resistance. Iridovirus causes mass mortality in hatchery of an important food fish Asian seabass (Lates calcarifer). In this study, an SNX2 (Sorting nexin 2) homolog (LcaSNX2) from L. calcarifer was cloned. Its roles in the infection of Singapore grouper iridovirus (SGIV) were characterized. The gene contained an ORF of 1542 bp encoding 513 amino acids, a 5’ UTR of 242 bp and a 3′ UTR of 2650 bp. Two SNPs in the 3’ UTR of the gene were identified, and they were significantly associated with SGIV resistance. Thus, they provide useful DNA markers for selecting Asian seabass juveniles for resistance against SGIV. The SNX2 gene was expressed ubiquitously in healthy fish, and highly expressed in the eye, gill, fin and brain. The expression of the gene was induced after SGIV infection in vitro, suggesting that LcaSNX2 plays a role in virus-host interaction during SGIV infection. Overexpression of LcaSNX2 increased the infection and replication of SGIV. Knockdown of this gene reduced viral gene transcriptional level, cytopathic effect and virus titer, indicating that suppression of the gene reduces SGIV infection. These data indicate that LcaSNX2 gene is essential for efficient replication of SGIV.

Asian seabass, an important food fish in Southeast Asia, has suffered from nervous necrosis virus (NNV) infection, resulting in massive mortality of Asian seabass larvae and enormous economic losses. Identification and characterization of disease resistance genes is important. Previous transcriptome analysis of Asians seabass epithelial cells after NNV infection revealed a highly inducible gene, receptor-transporting protein 3 ( rtp3 ), indicating it could play an important role in Asian seabass – NNV interaction. To characterize this gene, we determined its expression pattern and subcellular localization. The rtp3 was highly induced in most examined tissues and organs of Asian seabass after NNV infection, and protein Rtp3 was localized in cytoplasm. Further association study in multiple families revealed that a microsatellite marker, (GT)ntt(GT)n, in the 3′ UTR of rtp3 was significantly associated with VNN disease resistance in Asian seabass. Our results imply that rtp3 may be a novel disease resistance gene in Asian seabass. This data could improve our understanding of molecular interaction between Asian seabass and NNV, and has the potential to be applied in marker-assisted selection for disease resistance breeding in Asian seabass.

Hypoxia is one of the major challenges in aquaculture industry. Breeding of fish tolerant to hypoxia is an important task in genetic improvement of aquaculture species. The Asian seabass, Lates calcarifer, is an important foodfish species. We identified and characterized the hypoxia-inducible factor inhibitor (HIF1αn) gene in the Asian seabass. The full-length cDNA sequence of the HIF1αn was 3425 bp, with an ORF of 1065 bp, encoding 354 amino acids. The genomic sequence of the gene was 8667 bp in length, and contained eight exons and seven introns. Phylogenetic analysis of the gene in fish and tetrapods revealed that the HIF1αn in the Asian seabass was closely related to that of tilapia (Oreochromis niloticus). The HIF1αn was highly up-regulated in the gill, spleen and heart after 3.5-hours hypoxia treatment. We identified three SNPs in the third and fourth introns of the HIF1αn gene. The SNP (i.e. SNP 9332241 (C/T)) in the fourth intron was significantly (P 0.01) associated with hypoxia tolerance. This SNP might be useful in selecting Asian seabass for improved tolerance to hypoxia. Our data also provide useful information for further detailed analysis of the function of the HIF1αn gene in hypoxia tolerance.

Suv39h1 is a member of the Suv39h sub-family of histone lysine methyltransferases. The gene plays an important role in heterochromatin organization, chromosome segregation, and mitotic progression in humans and model organisms. However, its functions in food fish species remain largely unknown. We identified the full-length cDNA of the Suv39h1 in the Asian seabass, analysed its genomic structure and identified polymorphisms. The gene was expressed in all tissues examined, but highly expressed in the gill, spleen and kidney. After challenging Asian seabass juveniles with iridovirus, the expressions of the gene altered, suggesting that the gene is related to the infection by iridovirus. Three SNPs and one indel were detected in the 3’ UTR of the gene. The indel was significantly associated with the resistance against iridovirus. In silico-mapping revealed that the gene was mapped to chromosome 6, where a significant QTL for the resistance to iridovirus was located previously in the Asian seabass. Altogether, the data suggest that the gene Suv39h1 play an important role in iridovirus resistance in the Asian seabass. This study also laid the foundation for further detailed analysis of the function of the Suv39h1 gene in resistance to iridovirus. The indel may supply a useful DNA marker for selecting Asian seabass fingerlings resistant to iridovirus.

Oil palm is the world’s leading source of vegetable oil and fat. Dura, Pisifera and Tenera are three forms of oil palm. The genome sequence of Pisifera is available whereas the Dura form has not been sequenced yet. We sequenced the genome of one elite Dura palm, and re-sequenced 17 palm genomes. The assemble genome sequence of the elite Dura tree contained 10,971 scaffolds and was 1.701 Gb in length, covering 94.49% of the oil palm genome. 36,105 genes were predicted. Re-sequencing of 17 additional palm trees identified 18.1 million SNPs. We found high genetic variation among palms from different geographical regions, but lower variation among Southeast Asian Dura and Pisifera palms. We mapped 10,000 SNPs on the linkage map of oil palm. In addition, high linkage disequilibrium (LD) was detected in the oil palms used in breeding populations of Southeast Asia, suggesting that LD mapping is likely to be practical in this important oil crop. Our data provide a valuable resource for accelerating genetic improvement and studying the mechanism underlying phenotypic variations of important oil palm traits.

Diseases are one of the major challenges in tilapia aquaculture. Identification of DNA markers associated with disease resistance may facilitate the acceleration of the selection for disease resistance. Gadd45a (growth arrest and DNA damage 45 A), a stress-inducible gene in humans and mice, has not been studied in fish. We characterized the two prologues of Gadd45a genes in hybrid tilapia. Gadd45a1 and Gadd45a2 shared an identical gene structure and showed an amino acid sequence identity of 73.8%. Their expressions were detected in all 10 tissues examined, with the kidney and gill having high transcriptional expressions. The expression levels of Gadd45a1 were significantly lower than those of Gadd45a2 in all examined tissues. After a challenge with a bacterial pathogen Streptococcus agalactiae, the expressions of the two genes were up-regulated significantly in the spleen, kidney, liver and intestine. These findings suggest that the two Gadd45a genes play an important role in the resistance to S. agalactiae in tilapia. We identified 10 SNPs in the two genes. The SNP markers in the two Gadd45a genes could be used to examine whether they are associated with resistance to S. agalactiae.

Polymorphisms in the FTO gene are associated with obesity and body mass index in humans and livestock. Little information of whether FTO plays an important role in aquaculture fish species is available. We cloned and characterized the FTO gene in an economically important food fish species: Asian seabass (Lates calcarifer). The full-length cDNA of the gene is 3679 bp, containing an ORF of 1935 bp encoding 644 amino acids, a 216 bp 5′ UTR and a 1538 bp 3′ UTR. The gene consisted of nine exons and eight introns and was 117,679 bp in length. Phylogenetic analysis revealed that the gene in Asian seabass was closely related to those of Japanese flounder and Nile tilapia. Analysis of its expressions using qRT-PCR showed that it was expressed ubiquitously, but was higher in the liver, stomach and intestine. Comparative analysis of the genomic sequences of part of intron 1 of the gene among 10 unrelated individuals identified two SNPs. Analysis of associations between SNPs and traits (i.e. growth, oil content, omega-3 and -6 contents) in an F2 family demonstrated that the two SNPs were significantly associated with growth, oil content, omega-3 content and omega-3/-6 ratio. Altogether, our data suggest that the gene or/and its linked genes play an important role in growth and fatty acid synthesis, and that the SNPs associated with traits may be used as markers for selecting quicker growth and higher omega-3/-6 ratio at the fingerling stage. NRF (Natl Research Foundation, S’pore)

Resistance to pathogens is important for the sustainability and profitability of food fish production. In immune-related genes, the lipopolysaccharide-binding protein (LBP) gene is an important mediator of the inflammatory reaction. We analyzed the cDNA and genomic structure of the LBP gene in tilapia. The full-length cDNA (1901 bp) of the gene contained a 1416 bp open reading frame, encoding 471 amino acid residues. Its genomic sequence was 5577 bp, comprising 15 exons and 14 introns. Under normal conditions, the gene was constitutively expressed in all examined tissues. The highest expression was detected in intestine and kidney. We examined the responses of the gene to challenges with two bacterial pathogens Streptcoccus agalactiae and Aeromonas hydrophila. The gene was significantly upregulated in kidney and spleen post-infection with S. agalactiae and A. hydrophila, respectively. However, the expression profiles of the gene after the challenge with the two pathogens were different. Furthermore, we identified three SNPs in the gene. There were significant associations (p &lt, 0.05) of two of the three SNPs with the resistance to A. hydrophila, but not with the resistance to S. agalactiae or growth performance. These results suggest that the LBP gene is involved in the acute-phase immunologic response to the bacterial infections, and the responses to the two bacterial pathogens are different. The two SNPs associated with the resistance to A. hydrophila may be useful in the selection of tilapia resistant to A. hydrophila.

Jatropha curcas is a new promising bioenergy crop due to the high oil content in its seeds that can be converted into biodiesel. Seed size, a major determinant of Jatropha oil yield, is a target trait for Jatropha breeding. Due to the vital roles of phytohormone auxin in controlling seed and fruit development, we screened key genes in auxin pathway including ARF and IAA families and downstream effectors to identify candidate genes controlling seed size in Jatropha. As a result, JcARF19 was mapped in the major quantitative trait locus (QTL) region and significantly associated with seed length. By using expression QTL (eQTL) analysis to link variants with functional candidate genes, we provided evidences that seed traits were affected by the interaction of JcARF19 and JcIAA9. ARF19 and IAA9, involved in auxin signal transduction, were conserved in higher plants. These data including the single-nucleotide polymorphisms (SNPs) in the two genes could lead to utilization of the genes by integrating favored alleles into elite varieties through marker-assisted selection.

Mast cell proteases play an important role in the regulation of the immune response. We identified the cDNA of the mast cell protease 8 (MCP-8) gene and analyzed its genomic structure in tilapia. The ORF of the MCP-8 was 768 bp, encoding 255 amino acids. Quantitative real-time PCR revealed that the MCP-8 gene was expressed predominantly in spleen, moderately in liver, blood, brain, gill, intestine, skin, and weakly expressed in kidney, muscle and eye. After a challenge with Streptococcus agalactiae, the gene was induced significantly (p0.05) in intestine, kidney, spleen and liver. Furthermore, we identified five single nucleotide polymorphisms (SNPs) in the MCP-8 gene and found that three SNPs were significantly associated (p0.05) with resistance against S. agalactiae. However, we found no association between four SNPs and growth traits (p0.05). These results suggest that the MCP-8 gene play an important role in the resistance to S. agalactiae in tilapia. The SNP markers in the MCP-8 gene associated with the resistance to the bacterial pathogen may facilitate selection of tilapia resistant to the bacterial disease.

Sexual dimorphism is an interesting biological phenomenon. Previous studies showed that DNA methylation might play a role in sexual dimorphism. However, the overall picture of the genome-wide methylation landscape in sexually dimorphic species remains unclear. We analyzed the DNA methylation landscape and transcriptome in hybrid tilapia (Oreochromis spp.) using whole genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq). We found 4,757 sexually dimorphic differentially methylated regions (DMRs), with significant clusters of DMRs located on chromosomal regions associated with sex determination. CpG methylation in promoter regions was negatively correlated with the gene expression level. MAPK/ERK pathway was upregulated in male tilapia. We also inferred active cis-regulatory regions (ACRs) in skeletal muscle tissues from WGBS datasets, revealing sexually dimorphic cis-regulatory regions. These results suggest that DNA methylation contribute to sex-specific phenotypes and serve as resources for further investigation to analyze the functions of these regions and their contributions towards sexual dimorphisms.

Identification of differentially expressed genes (DEGs) and regulated pathways in response to stressors using a whole-genome approach is critical to understanding the mechanisms underlying stress responses. We challenged Asian seabass with lipopolysaccharide (LPS), Vibrio harveyi, high salinity and fasting, and sequenced six cDNA libraries of intestine samples using Roche 454 RNA-seq. Over 1 million reads (average size: 516 bp) were obtained. The de novo assembly obtained 83 911 unisequences with an average length of 747 bp. In total, 62.3% of the unisequences were annotated. We observed overall similar expression profiles among different challenges, while a number of DEGs and regulated pathways were identified under specific challenges. More than 1000 DEGs and over 200 regulated pathways for each stressor were identified. Thirty-seven genes were differentially expressed in response to all challenges. Our data suggest that there is a global coordination and fine-tuning of gene regulation during different challenges. In addition, we detected dramatic immune responses in intestines under different stressors. This study is the first step towards the comprehensive understanding of the mechanisms underlying stress responses and supplies significant transcriptome resources for studying biological questions in non-model fish species.

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R 2 > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.

Understanding the genomic signatures of population differentiation is fundamental to obtain a comprehensive view of the evolutionary process of organisms. Grass carp is one of the most important fish species in the world due to its significant value in aquaculture and world-wide vegetation biocontrol. However, little is known about the contemporary population structure and also the genetic basis of adaptation to a wide range of latitudinal environments. Using 43310 SNPs generated by genotyping by sequencing in 197 grass carps from nine populations, we examined the patterns of population differentiation, historical introduction and evidence of local selection. The overall genetic differentiation across all native populations was unexpectedly low. Nevertheless, these native populations were clearly differentiated into three genetic clusters, corresponding to the Yangtze River, the Pearl River and the Heilongjiang River System, respectively. Populations in Malaysia, India and Nepal, with the earliest introduction records, most likely have an origin from the Pearl River System. Using conceptually different approaches, 451 loci were detected under potential local selection, among which 84 were annotated to have a gene feature. 19.0% of the genes under putative selection were involved in immune responses, while 42.9% of the annotated loci showed a signature of latitudinal variation. This study provides valuable information for application of genomic tools in addressing questions concerning population differentiation that was influenced by both neutral and adaptive forces, as well as human activities.

Zinc finger AN1-type domain 3 (ZFAND3) is essential for spermatogenesis in mice. However, its function in teleosts remains unclear. In this study, we characterized the ZFAND3 gene (termed as OsZFAND3) in an important food fish, tilapia. The OsZFAND3 cDNA sequence is 1,050 bp in length, containing an ORF of 615 bp, which encodes a putative peptide of 204 amino acid residues. Quantitative real-time PCR revealed that the OsZFAND3 transcripts were exclusively expressed in the testis and ovary. In situ hybridization showed that the high expression of OsZFAND3 transcripts was predominantly localized in the spermatocyte and spermatid. These results suggest that OsZFAND3 is involved in male germ cell maturation. Three single nucleotide polymorphisms (SNPs) were detected in the introns of OsZFAND3. The OsZFAND3 gene was mapped in the sex-determining locus on linkage group 1 (LG1). The three SNPs in the OsZFAND3 gene were strictly associated with sex phenotype, suggesting that the OsZFAND3 gene is tightly linked to the sex-determining locus. Our study provides new insights into the functions of the OsZFAND3 gene in tilapia and a foundation for further detailed analysis of the OsZFAND3 gene in sex determination and differentiation.

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

MicroRNAs (miRNAs) are small noncoding RNAs that play crucial regulatory roles by targeting mRNAs for silencing. To identify miRNAs in Jatropha curcas L, a bioenergy crop, cDNA clones from two small RNA libraries of leaves and seeds were sequenced and analyzed using bioinformatic tools. Fifty-two putative miRNAs were found from the two libraries, among them six were identical to known miRNAs and 46 were novel. Differential expression patterns of 15 miRNAs in root, stem, leave, fruit and seed were detected using quantitative real-time PCR. Ten miRNAs were highly expressed in fruit or seed, implying that they may be involved in seed development or fatty acids synthesis in seed. Moreover, 28 targets of the isolated miRNAs were predicted from a jatropha cDNA library database. The miRNA target genes were predicted to encode a broad range of proteins. Sixteen targets had clear BLASTX hits to the Uniprot database and were associated with genes belonging to the three major gene ontology categories of biological process, cellular component, and molecular function. Four targets were identified for JcumiR004. By silencing JcumiR004 primary miRNA, expressions of the four target genes were up-regulated and oil composition were modulated significantly, indicating diverse functions of JcumiR004.

Asian seabass is an important food fish in Southeast Asia. Viral nervous necrosis (VNN) disease, triggered by nervous necrosis virus (NNV) infection, has caused mass mortality of Asian seabass larvae, resulting in enormous economic losses in the Asian seabass industry. In order to better understand the complex molecular interaction between Asian seabass and NNV, we investigated the transcriptome profiles of Asian seabass epithelial cells, which play an essential role in immune regulation, after NNV infection. Using the next generation sequencing (NGS) technology, we sequenced mRNA from eight samples (6, 12, 24, 48 h post-inoculation) of mock and NNV-infected Asian seabass epithelial cell line, respectively. Clean reads were de novo assembled into a transcriptome consisting of 89026 transcripts with a N50 of 2617 bp. Furthermore, 251 differentially expressed genes (DEGs) in response to NNV infection were identified. Top DEGs include protein asteroid homolog 1-like (ASTE1), receptor-transporting protein 3 (RTP3), heat shock proteins 30 (HSP30) and 70 (HSP70), Viperin, interferon regulatory factor 3 (IRF3) and other genes related to innate immunity. Our data suggest that abundant and diverse genes corresponding to NNV infection. The results of this study could also offer vital information not only for identification of novel genes involved in Asian seabass-NNV interaction, but also for our understanding of the molecular mechanism of Asian seabass' response to viral infection. In addition, 24807 simple sequence repeats (SSRs) were detected in the assembled transcriptome, providing valuable resources for studying genetic variations and accelerating quantitative trait loci (QTL) mapping for disease resistance in Asian seabass in the future.

Prolactin (encoded by PRL) is a multifunctional hormone involved in osmoregulation, reproduction, growth, development, immunomodulation, endocrine and metabolic regulation. We cloned the full-length cDNA of Asian seabass PRL, searched for polymorphism in the DNA sequence, and conducted association analyses. Twelve SNPs and one 4-bp deletion were identified in PRL. The SNP c.264+127C>G was used for linkage mapping, and this gene was mapped to linkage group 11. The c.264+980_983delTTGT, c.264+127C>G, c.264+138T>G, c.264+269T>C and c.330C>G polymorphisms were genotyped in 521 individuals with growth trait records. Association analyses between single markers and growth traits revealed that the c.264+269T>C SNP was significantly associated with body weight (BW), total length (TL), standard length (SL) and Fulton's condition factor (KTL and KSL), while the other four were not. Analysis of haplotypes showed that there were 10 haplotypes and 22 haplotype combinations in the population. The differences of BW, TL, KTL and KSL among different haplotype combinations were significant.

Liver-expressed antimicrobial peptide 2 (LEAP-2) is a novel antimicrobial peptide (AMP) recently found in vertebrates, and exhibits distinct amino acid sequence, secondary structure and expression pattern from other peptides. In this study, the LEAP-2 gene and its full-length cDNA were cloned from grass carp. Grass carp LEAP-2 gene consists of two introns and three exons. The translated product contains 92 amino acids, including a 26 amino acids signal peptide and a mature peptide of 41 amino acids. Grass carp LEAP-2 gene was expressed in a wide range of tissues except blood, with the highest level of transcripts found in liver. Upon induction by Aeromonas hydrophila, its expression was significantly up-regulated in liver, gill, skin, muscle, spleen, blood, head kidney, heart and intestine, but down-regulated in trunk kidney and brain. The transcript level was high in embryos at the 16-cell stage but declined gradually afterwards, suggesting that LEAP-2 transcripts in early embryos might be maternal. Mature peptides obtained by in vitro expression displayed selective antimicrobial activities. These results together further our understanding of the physiological function of LEAP-2 in vertebrates.

Summary Quantitative trait loci (QTL) affecting growth traits have previously been mapped in linkage groups (LG) 2, 3 and 23 of Barramundi (Lates carcalifer), but these QTL have not been verified in different genetic backgrounds and environments. Here, we report the identification and verification of QTL for growth traits on LG2, 3, 10 and 23 in F1 families constructed using brooders from the Singapore Marine Aquaculture Center (MAC) and from wild stocks collected in Thailand (THAI). The previously detected QTL for body weight and length linked to marker Lca371 on LG2 were confirmed in both the MAC and THAI families, whereas other QTL previously mapped to LG3 and 23 were only detected in one of the two families. QTL for body weight and length were identified in the MAC family, but not in the THAI family, in a region where the insulin-like growth factor 2 (IGF2) and tyrosine hydroxylase 1 (TH1) genes are located on LG10. Significant epistatic interactions were identified between markers Lca287 on LG2 and IGF2 on LG10 for growth trait QTL in the MAC family, but not in the THAI family. Effects of the IGF2, TH1 and parvalbumin 1 candidate genes were family-specific. Our results indicate that some but not all QTL are family-specific in Barramundi.

Barramundi (Lates calcarifer) is an important farmed marine food fish species. Its compact genome (∼700 Mb) is among the smallest genomes of food fish species. We established a first-generation genetic linkage map of Barramundi with a mapping panel containing three parents (two males and one female) and 93 progeny. A total of 240 microsatellite markers were mapped into 24 linkage groups. Among these markers, 10 were located in ESTs and known genes. The total lengths of the female and male maps were 873.8 and 414.5 cM with an average marker spacing of 6.20 and 4.70 cM, respectively. Comparing the flanking sequences of the 240 Barramundi microsatellites with the assembled whole-genome sequences of Tetraodon nigrovidiris revealed 55 homologous sequences located in 19 of the 21 chromosomes of T. nigrovidiris. The map will not only enable the mapping of quantitative trait loci, but also provide new resources for understanding the evolution of fish genomes.

Tilapia is the common name for a group of cichlid fishes. Identification of DNA markers significantly associated with important traits in candidate genes may speed up genetic improvement. L-Amino acid oxidase (LAO) plays a crucial role in the innate immune defences of animals. Previously, whether LAO variants were associated with economic traits had not been studied in fish. We characterized the cDNA sequence of the LAO gene of hybrid tilapia (Oreochromis spp.). Its ORF was 1536 bp, encoding a flavoenzyme of 511 amino acids. This gene consisted of seven exons and six introns. Its expression was detected in the intestine, blood, kidney, skin, liver. It was highly expressed in the intestine. After a challenge with a bacterial pathogen, Streptococcus agalactiae, its expression was up-regulated significantly in the liver, intestine and spleen (P

Tilapia is a group of cultured teleost fishes whose production is threatened by some diseases. Identification of DNA markers associated with disease resistance in candidate genes may facilitate to accelerate the selection of disease resistance. The gene encoding a duodenase, which can trigger immune response, has not been studied in fish. We characterized the cDNA of duodenase-1 gene of hybrid tilapia. Its ORF is 759 bp, encoding a serine protease of 252 amino acids. This gene consisted of five exons and four introns. Its expression was detected in all 10 tissues examined, and it was highly expressed in the intestine and kidney. After a challenge with the bacterial pathogen, Streptococcus agalactiae, its expression was up-regulated significantly in the intestine, liver and spleen. We identified seven SNPs in the gene and found that four of them were significantly associated with the resistance to S. agalactiae (P < 0.05). The CGTCC haplotype, CAGTC/CGGTC and CGTCC/CGTCC diplotype were significantly associated with the resistance to S. agalactiae (P = 0.00, 0.04 and < 0.0001, respectively). In addition, one SNP was associated significantly with growth traits (P < 0.05). These findings suggest that the duodenase-1 gene plays an important role in the resistance to S. agalactiae in tilapia. The SNP markers in the duodenase-1 gene associated with resistance to the bacterial pathogen, may facilitate the selection of tilapia resistant to the bacterial disease.

We determined the complete mtDNA nucleotide sequence of Lates calcarifer using the shotgun sequencing method. The mitochondrial DNA (mtDNA) was 16,535 base pairs (bp) in length, and contained 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and one major noncoding control region (CR). The CR was unusually short at only 768 bp. A striking feature of the mitochondrial genome was the high G+C content (46.1%), which is among the highest in fish. The gene order was identical to that of a typical vertebrate. Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes of 30 fish species representing 14 suborders clearly showed Lates calcarifer was located in the cluster of fish species from the order Perciformes, supporting the traditional systematic classification. We characterized single-nucleotide polymorphisms (SNPs) in the CR by sequencing the complete CR of 25 individuals obtained from Australia and Singapore. A total of 68 SNPs were detected. Eighteen SNPs were fixed with alternative nucleotides in Australian and Singapore seabass, and these SNPs could be used for differentiating fish from the two countries.

The fresh water pufferfish Tetraodon nigroviridis is a model organism for studying evolution of genome and gene functions, but its mitochondrial genome (mtDNA) sequence is still not available. We determined the complete nucleotide sequence of its mtDNA using shotgun sequencing. The T. nigroviridis mtDNA was 16,462 bp, and contained 13 protein coding genes, 22 tRNAs, 2 rRNAs and a major non-coding region. The gene order was identical to the common type of vertebrate mtDNA, whereas the G+C content in the sense strand was 46.9%, much higher than most other fish species. One hundred and three SNPs were detected in the control region of the mtDNA of 35 individuals, a majority of SNPs were detected in the 5′ end of the control region. A phylogenetic study including 21 fish species was performed on concatenated amino acid sequences of 12 protein coding genes, and revealed that the T. nigroviridis was clustered with Fugu rubripes into a group. The complete mtDNA sequence and SNPs in its control region will be use...

Growth is an important trait in animal breeding. However, the genetic effects underpinning fish growth variability are still poorly understood. QTL mapping and analysis of candidate genes are effective methods to address this issue. We conducted a genome-wide QTL analysis for growth in tilapia. A total of 10, 7 and 8 significant QTLs were identified for body weight, total length and standard length at 140 dph, respectively. The majority of these QTLs were sex-specific. One major QTL for growth traits was identified in the sex-determining locus in LG1, explaining 71.7%, 67.2% and 64.9% of the phenotypic variation (PV) of body weight, total length and standard length, respectively. In addition, a candidate gene GHR2 in a QTL was significantly associated with body weight, explaining 13.1% of PV. Real-time qPCR revealed that different genotypes at the GHR2 locus influenced the IGF-1 expression level. The markers located in the major QTL for growth traits could be used in marker-assisted selection of tilapia. The associations between GHR2 variants and growth traits suggest that the GHR2 gene should be an important gene that explains the difference in growth among tilapia species.

Studies on the zebrafish model have contributed to our understanding of several important developmental processes, especially those that can be easily studied in the embryo. However, our knowledge on late events such as gonad differentiation in the zebrafish is still limited. Here we provide an analysis on the gene sets expressed in the adult zebrafish testis and ovary in an attempt to identify genes with potential role in (zebra)fish gonad development and function. We produced 10 533 expressed sequence tags (ESTs) from zebrafish testis or ovary and downloaded an additional 23 642 gonad-derived sequences from the zebrafish EST database. We clustered these sequences together with over 13 000 kidney-derived zebrafish ESTs to study partial transcriptomes for these three organs. We searched for genes with gonad-specific expression by screening macroarrays containing at least 2600 unique cDNA inserts with testis-, ovary- and kidney-derived cDNA probes. Clones hybridizing to only one of the two gonad probes were selected, and subsequently screened with computational tools to identify 72 genes with potentially testis-specific and 97 genes with potentially ovary-specific expression, respectively. PCR-amplification confirmed gonad-specificity for 21 of the 45 clones tested (all without known function). Our study, which involves over 47 000 EST sequences and specialized cDNA arrays, is the first analysis of adult organ transcriptomes of zebrafish at such a scale. The study of genes expressed in adult zebrafish testis and ovary will provide useful information on regulation of gene expression in teleost gonads and might also contribute to our understanding of the development and differentiation of reproductive organs in vertebrates.

Leukocyte cell-derived chemotaxin-2 (LECT2) is an important protein of the innate immune system for the defense against bacterial infection. We cloned and characterized the LECT2 gene from Asian seabass (Lates calcarifer). Its complete cDNA consisted of an open reading frame of 459 bp encoding a protein of 152 amino acids. The genomic DNA sequence of this gene consists of four exons and three introns. Quantitative real-time PCR revealed that the LECT2 gene was expressed predominantly in liver while its expression was moderate in spleen and heart, and weak in other tissues. The LECT2 transcript was up-regulated in the kidney, spleen and liver in response to a challenge with a pathogenic bacterium Vibrio harveyi. In addition, we identified three single nucleotide polymorphisms (SNPs) in the LECT2 gene, and found significant associations between these polymorphisms and resistance to the big belly disease. These results suggest that the LECT2 gene play an important role in resistance to bacterial pathogens in fish. The SNP markers in the gene associated with the resistance to bacterial pathogens may facilitate selecting Asian seabass resistant to bacterial diseases.

Summary Parvalbumins are extremely abundant in fish muscle and play an important role in muscle relaxation. In this study, two parvalbumin genes (PVALB1 and PVALB2) were cloned from Asian seabass (Lates calcarifer). The cDNAs for PVALB1 and PVALB2 were 840 and 667 bp respectively. Both genes consisted of five exons and four introns, encoded 109 amino acids, and were of beta lineage. Using real-time polymerase chain reaction, expression of PVALB1 was detected in all 10 tissues tested, with expression in brain, kidney, muscle and small intestine being 15- to 322-fold higher than in the other tissues. Expression of PVALB2 was detected only in muscle, brain and intestine, and was up to 10-fold lower than PVALB1 expression. A (CT)17 microsatellite in the 3¢-untranslated region of PVALB1 and three single nucleotide polymorphisms (SNPs) in the third intron of PVALB2 were identified. The microsatellite in PVALB1 was significantly associated with body weight and body length at 90 days post-hatch (P < 0.01), whereas the SNPs in PVALVB2 were not associated with these traits.

Background Aquaculture is the quickest growing sector in agriculture. However, QTL for important traits have been only identified in a few aquaculture species. We conducted QTL mapping for growth traits in an Asian seabass F2 family with 359 individuals using 123 microsatellites and 22 SNPs, and performed association mapping in four populations with 881 individuals. Results Twelve and nine significant QTL, as well as 14 and 10 suggestive QTL were detected for growth traits at six and nine months post hatch, respectively. These QTL explained 0.9-12.0% of the phenotypic variance. For body weight, two QTL intervals at two stages were overlapped while the others were mapped onto different positions. The IFABP-a gene located in a significant QTL interval for growth on LG5 was cloned and characterized. A SNP in exon 3 of the gene was significantly associated with growth traits in different populations. Conclusions The results of QTL mapping for growth traits suggest that growth at different stages was controlled by some common QTL and some different QTL. Positional candidate genes and association mapping suggest that the IFABP-a is a strong candidate gene for growth. Our data supply a basis for fine mapping QTL, marker-assisted selection and further detailed analysis of the functions of the IFABP-a gene in fish growth.

Three F2 families have been established for QTL mapping by crossing genetically divergent Wild boar, Meishan and Pietrain pigs. More than 900 F2 animals are typed for a large number of marker loci and are recorded for over 100 quantitative traits. In a first analysis QTLs on chromosome 6 have been mapped in the region of the ryanodine receptor gene. Also several traits were significantly associated with different GH haplotyes.

Calreticulin (CRT) is a Ca2+-binding molecular chaperone in the endoplasmic reticulum. We cloned and characterized the CRT gene in an important marine food fish species Asian seabass (Lates calcarifer). The full-length DNA of the CRT gene was 2194 bp, including a complete open reading frame encoding 420 amino acid residues, a 113 bp 5′-untranslated region and an 818 bp 3′-untranslated region. The CRT gene contained nine exons and eight introns covering a total of 2772 bp genomic DNA from the start to stop codon. Ten single nucleotide polymorphisms (SNPs) were detected in introns and an exon in six individuals collected from five different locations. The CRT gene was assigned to linkage group 4 of the linkage map of Asian seabass. Quantitative real-time PCR revealed that the CRT gene was highly expressed in liver at the age of 1, 3 and 7 months under normal conditions, whereas its expression in liver reduced sharply after 0.5 to 2 h cold challenge at 16°C, and then increased slowly. A preliminary association analysis showed a significant (P < 0.001) association between the SNP6 in the CRT gene and the mortality after cold challenge at 16°C. Our results suggest that the CRT gene is associated with cold tolerance of Asian seabass and further investigation will be necessary to illustrate the underlying mechanisms.

Background Body weight and length are economically important traits in foodfish species influenced by quantitative trait loci (QTL) and environmental factors. It is usually difficult to dissect the genetic and environmental effects. Asian seabass (Lates calcarifer) is an important marine foodfish species with a compact genome (~700 Mb). The recent construction of a first generation linkage map of Asian seabass with 240 microsatellites provides a good opportunity to determine the number and position of QTL, and the magnitude of QTL effects with a genome scan. Results We conducted a genome scan for QTL affecting body weight, standard length and condition factors in an F1 family containing 380 full-sib individuals from a breeding stock by using 97 microsatellites evenly covering 24 chromosomes. Interval mapping and multiple QTL model mapping detected five significant and 27 suggestive QTL on ten linkage groups (LGs). Among the five significant QTL detected, three (qBW2-a, qTL2-a and qSL2-a) controlling body weight, total and standard length respectively, were mapped on the same region near Lca287 on LG2, and explained 28.8, 58.9 and 59.7% of the phenotypic variance. The other two QTL affecting body weight, qBW2-b and qBW3, were located on LG2 and 3, and accounted for 6.4 and 8.8% of the phenotypic variance. Suggestive QTL associated with condition factors are located on six different LGs. Conclusion This study presents the first example of QTL detection for growth-related traits in an F1 family of a marine foodfish species. The results presented here will enable further fine-mapping of these QTL for marker-assisted selection of the Asian seabass, eventually identifying individual genes responsible for growth-related traits.

Fgf21 is a newly discovered fibroblast growth factor. It is typically induced by fasting and plays important roles in the regulation of glucose and lipid metabolisms and energy balance in mammals, whereas potential functions of this gene in teleosts are still unknown. We identified the Fgf21 gene and studied its functions in Asian seabass (Lates calcarifer). The cDNA of the Fgf21 encoded a protein with 206 amino acids. Analysis of DNA and amino acid sequences of Fgf21 genes revealed that the sequences and structure of the Fgf21 genes were highly conserved in vertebrates. Real-time PCR revealed that Fgf21 was exclusively expressed in the intestine and kidney, which was different from the expression profiles of mammals. Fgf21 was down-regulated under fasting, whereas it was significantly increased during the LPS challenge. Exogenous recombinant FGF21 significantly suppressed the appetite of Asian seabass. Our data suggest that Fgf21 plays a role in energy regulation and acute phase response in Asian seabass, and may have different functions in fish and mammals. In addition, we identified one SNP in Fgf21. By using this SNP, the gene was mapped on the linkage group 23, where a suggestive QTL for growth was mapped previously. Association mapping identified significant associations between Fgf21 genotypes at the SNP and growth traits. These results not only provide important information of the functions of Fgf21, but also suggest that the SNP in this gene can be used as a marker in selecting fast-growing individuals of Asian seabass.

For detecting QTL in the whole swine genome, 1068 pigs from three F2 populations constructed by crossing European Wild boar and Pietrain (W x P), Meishan and Pietrain (M x P), and Wild Boar and Meishan (W x M) were genotyped for genetic markers evenly spaced at approximately 20 cM intervals. AQTL analysis was performed using a least-squares method. Here the results of the QTL analysis on the porcine chromosome 7 are presented. QTL for carcass composition (e.g. head weight, carcass length, backfat depth, abdominal fat and bacon meat) were mapped in the chromosomal region CYPA/CYPD-TNFB-S0102 in M x P and W x M, but not in W x P. The QTL explained 5.3%-27.2% of the F2 phenotypic variance in the two F2 populations. Most traits affected by the mapped QTL were related to carcass fatness. The mode of gene action of QTL was additive. Surprisingly, in contrast to the parental phenotype, the QTL alleles from fatty Meishan were associated with thinner backfat than Pietrain and Wild Boar alleles, suggesting that the genome of the fatty Meishan pig contains genes which can reduce fat content of carcass substantially.

Lysozymes are important proteins of the innate immune system for the defense against bacterial infection. We cloned and analyzed chicken-type (c-type) and goose-type (g-type) lysozymes from Asian seabass (Lates calcarifer). The deduced amino acid sequence of the c-type lysozyme contained 144 residues and possessed typical structure residues, conserved catalytic residues (Glu(50) and Asp(67)) and a "GSTDYGIFQINS" motif. The deduced g-type lysozyme contained 187 residues and possessed a goose egg white lysozyme (GEWL) domain containing three conserved catalytic residues (Glu(71), Asp(84), Asp(95)) essential for catalytic activity. Real time quantitative PCR (qRT-PCR) revealed that the two lysozyme genes were constitutively expressed in all the examined tissues. The c-type lysozyme was most abundant in liver, while the g-type lysozyme was predominantly expressed in intestine and weakly expressed in muscle. The c-type and g-type transcripts were up-regulated in the kidney, spleen and liver in response to a challenge with Vibrio harveyi. The up-regulation of the c-type lysozyme was much stronger than that of the g-type lysozyme in kidney and spleen. The recombinant proteins of the c-type and g-type lysozymes showed lytic activities against the bacterial pathogens Vibrio harveyi and Photobacterium damselae in a dosage-dependent manner. We identified single nucleotide polymorphisms (SNPs) in the two lysozyme genes. There were significant associations of these polymorphisms with resistance to the big belly disease. These results suggest that the c- and g-type genes play an important role in resistance to bacterial pathogens in fish. The SNP markers in the two genes associated with the resistance to bacterial pathogens may facilitate the selection of Asian seabass resistant to bacterial diseases.

The goldfish (Carassius auratus auratus; Cyprinidae) is not only an important ornamental fish species, but also a useful model for biological studies. The sequence of goldfish genes present in the public database was searched for short tandem repeats, and 11 polymorphic microsatellites were detected within eight important genes. Two microsatellites were located in coding regions of the c-myc and GAP-43 genes, respectively, whereas the others were mainly located in 5′ and 3′ untranscribed regions of other genes, such as gonadotropin and activin. The average allele number of these microsatellites was 5.5 per locus with a range between 3 and 9. These microsatellites will be useful for ecological and population genetic studies of this species, as well as for the ornamental fish industry.

Iron is one of the most important minor elements in the shells of bivalves. This study was designed to investigate the involvement of ferritin, the principal protein for iron storage, in shell growth. A novel ferritin subunit (Fth1) cDNA from the freshwater pearl mussel (Hyriopsis cumingii) was isolated and characterized. The complete cDNA contained 822 bp, with an open reading frame (ORF) of 525 bp, a 153 bp 5' untranslated region (UTR) and a 144 bp 3' UTR. The complete genomic DNA was 4125 bp, containing four exons and three introns. The ORF encoded a protein of 174 amino acids without a signal sequence. The deduced ferritin contained a highly conserved motif for the ferroxidase center comprising seven residues of a typical vertebrate heavy-chain ferritin. It contained one conserved iron associated residue (Try27) and iron-binding region signature 1 residues. The mRNA contained a 27 bp iron-responsive element with a typical stem-loop structure in the 5'-UTR position. Copy number variants (CNVs) of Fth1 in two populations (PY and JH) were detected using quantitative real-time PCR. Associations between CNVs and growth were also analyzed. The results showed that the copy number of the ferritin gene of in the diploid genome ranged from two to 12 in PY, and from two to six in JH. The copy number variation in PY was higher than that in JH. In terms of shell length, mussels with four copies of the ferritin gene grew faster than those with three copies (P

BackgroundMicroRNAs (miRNAs) play an important role in the regulation of many fundamental biological processes. So far miRNAs have been only identified in a few fish species, although there are over 30,000 fish species living under different environmental conditions on the earth. Here, we described an approach to identify conserved miRNAs and characterized their expression patterns in different tissues for the first time in a food fish species Asian seabass (Lates calcarifer).Methodology/principal findingsBy combining a bioinformatics analysis with an approach of homolog-based PCR amplification and sequencing, 63 novel miRNAs belonging to 29 conserved miRNA families were identified. Of which, 59 miRNAs were conserved across 10-86 species (E value ≤ 10⁻⁴) and 4 miRNAs were conserved only in fish species. qRT-PCR analysis showed that miR-29, miR-103, miR-125 and several let-7 family members were strongly and ubiquitously expressed in all tissues tested. Interestingly, miR-1, miR-21, miR-183, miR-184 and miR-192 showed highly conserved tissue-specific expression patterns. Exposure of the Asian seabass to lipopolysaccharide (LPS) resulted in up-regulation of over 50% of the identified miRNAs in spleen suggesting the importance of the miRNAs in acute inflammatory immune responses.Conclusions/significanceThe approach used in this study is highly effective for identification of conserved miRNAs. The identification of 63 miRNAs and determination of the spatial expression patterns of these miRNAs are valuable resources for further studies on post-transcriptional gene regulation in Asian seabass and other fish species. Further identification of the target genes of these miRNAs would shed new light on their regulatory roles of microRNAs in fish.

Background Fish diseases caused by pathogens are limiting their production and trade, affecting the economy generated by aquaculture. Innate immunity system is the first line of host defense in opposing pathogenic organisms or any other foreign material. For identification of immune-related genes in Asian seabass Lates calcarifer, an important marine foodfish species, we injected bacterial lipopolysaccharide (LPS), a commonly used elicitor of innate immune responses to eight individuals at the age of 35 days post-hatch and applied the suppression subtractive hybridization (SSH) technique to selectively amplify spleen cDNA of differentially expressed genes. Results Sequencing and bioinformatic analysis of 3351 ESTs from two SSH libraries yielded 1692 unique transcripts. Of which, 618 transcripts were unknown/novel genes and the remaining 1074 were similar to 743 known genes and 105 unannotated mRNA sequences available in public databases. A total of 161 transcripts were classified to the category "response to stimulus" and 115 to "immune system process". We identified 25 significantly up-regulated genes (including 2 unknown transcripts) and 4 down-regulated genes associated with immune-related processes upon challenge with LPS. Quantitative real-time PCR confirmed the differential expression of these genes after LPS challenge. Conclusions The present study identified 1692 unique transcripts upon LPS challenge for the first time in Asian seabass by using SSH, sequencing and bioinformatic analysis. Some of the identified transcripts are vertebrate homologues and others are hitherto unreported putative defence proteins. The obtained immune-related genes may allow for a better understanding of immunity in Asian seabass, carrying out detailed functional analysis of these genes and developing strategies for efficient immune protection against infections in Asian seabass.