Search

The nervous necrosis virus (NNV) causes the viral nervous necrosis (VNN) disease in aquatic animals and has been a major threat in aquaculture. Thus, it is essential for the development of a prevention method to minimize economic losses caused by NNV such as the identification of NNV resistance genes and application of these genes in molecular breeding to increase disease resistance. gab3 is an important NNV resistance gene in Asian seabass. However, the mechanism of gab3 in NNV resistance has not been elucidated. In this study, knockdown of gab3 in NNV-infected Asian seabass cells resulted in a significant decrease in viral RNA and virus titers. Knockout of gab3 in zebrafish led to an increased survival rate and resistant time after NNV infection. Cellular localization of the GAB3 and NNV by immunofluorescence staining showed that the GAB3 was translocated from the nucleus to the cytoplasm, and finally reached the cell membrane of SB cells after 48 h post NNV infection. Our study suggests that gab3 plays an important role in NNV replication and silencing gab3 can inhibit virus replication.

Understanding of the relationships between genotypes and phenotypes is a central problem in biology. Although teleosts have colorful phenotypes, not much is known about their underlying mechanisms. Our previous study showed that golden skin color in Mozambique tilapia was mapped in the major locus containing the Pmel gene, and an insertion in 3' UTR of Pmel17 was fully correlated with the golden color. However, the molecular mechanism of how Pmel17 determines the golden skin color is unknown. In this study, knockout of Pmel17 with CRISPR/Cas9 in blackish tilapias resulted in golden coloration, and rescue of Pmel17 in golden tilapias recovered the wild-type blackish color, indicating that Pmel17 is the gene determining the golden and blackish color. Functional analysis in vitro showed that the insertion in the 3' UTR of Pmel17 reduced the transcripts of Pmel17. Our data supplies more evidence to support that Pmel17 is the gene for blackish and golden colors, and highlights that the insertion in the 3' UTR of Pmel17 is the causative mutation for the golden coloration.

Chia (Salvia hispanica L.), originated in central and southern Mexico and Guatemala, is an emerging industry crop due to its high content of omega‐3 fatty acids and dietary fiber in its seeds. The seeds also have a high concentration of proteins and essential amino acids, and are becoming a promising source of bioactive peptides. Polymorphic DNA markers are essential tools to analyse genetic diversity and to accelerate genetic improvement. However, in chia, polymorphic and codominant DNA markers are still lacking. In this study, fourteen polymorphic microsatellites were identified from DNA sequences and were characterized. The average allele number was 4.8 while the expected and observed heterozygosity was 0.24 and 0.34, respectively. The average probability of identity (PI) was 0.50 while the combined PI was 9 × 10–6. These first 14 microsatellites in chia are useful in genetic analysis and traceability. These 14 polymorphic microsatellites were used in analysing genetic diversity and population relationships in six cultivars originating in Mexico, Australia and Bolivia. Results showed that allelic diversity and gene diversity were low and ranged from 2.79 to 3.64 and 0.27 to 0.38, respectively. The Mexico black cultivar showed the highest allelic (3.64) and gene diversity (0.38). The six cultivars were closely related with high identity (≥ 0.893). Taken together, these chia cultivars contain low genetic variation. Therefore, to initiate a breeding program for improving traits, it is essential to use seeds from multiple cultivars to enlarge genetic variation in the founder population.

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.

High-density genetic maps are essential for mapping QTL, improving genome assembly, comparative genomics, and studying sex chromosome evolution. The northern snakehead (Channa argus) is an economically important foodfish species with significant sexual dimorphism, where the males grow much faster and bigger than the females. However, to date, the sex determination pattern is still not clear, limiting identification of sex chromosomes, even sex determination genes and development of monosex populations that are valuable for both sex evolution of vertebrates and aquaculture practices. Here, a sex-averaged map and two sex-specific genetic maps were constructed with 2974, 2323, and 2338 SNPs, respectively. Little difference was observed in the pattern of sex-specific recombination between female- and male-specific genetic maps. Genome scan identified a major locus for sex determination at LG16. Females and males are, respectively, homogametic and heterogametic, suggesting an XY sex determination system for this species. By resequencing genomes, InDels in the sex-associated QTL region were discovered and used for developing sex-specific PCR assays for fast sexing of snakehead. These high-density genetic maps provide useful resources for future genomic studies in snakehead and its related species. The PCR assays for sexing are of importance in developing all male populations for aquaculture. This research is supported by the internal fund of the Institute of Fisheries, Hangzhou Academy of Agricultural Science.

The oil palm, originating from Africa, is the most productive oil crop species. Palm oil is an important source of edible oil. Its current global plantation area is over 23 million ha. The theoretical oil yield potential of the oil palm is 18.2 tons/ha/year. However, current average oil yield is only 3.8 tons/ha/year. In the past 100 years, conventional breeding and improvement of field management played important roles in increasing oil yield. However, conventional breeding for trait improvement was limited by its very long (10–20 years) phenotypic selection cycle, although it improved oil yield by ~10–20% per generation. Molecular breeding using novel molecular technologies will accelerate genetic improvement and may reduce the need to deforest and to use arable land for expanding oil palm plantations, which in turn makes palm oil more sustainable. Here, we comprehensively synthesize information from relevant literature of the technologies, achievements, and challenges of molecular approaches, including tissue culture, haploid breeding, mutation breeding, marker-assisted selection (MAS), genomic selection (GS), and genome editing (GE). We propose the characteristics of ideal palms and suggest a road map to breed ideal palms for sustainable palm oil.

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)

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), leads to mass mortality in mariculture. However, phenotypic selection for resistance against VNN is very difficult. To facilitate marker-assisted selection (MAS) for resistance against VNN and understanding of the genetic architecture underlying the resistance against this disease, we mapped quantitative trait loci (QTL) for resistance against VNN in Asian seabass. We challenged fingerlings at 37 days post-hatching (dph), from a single back-cross family, with NNV at a concentration of 9 × 10(6) TCID50/ml for 2 h. Daily mortalities were recorded and collected. A panel of 330 mortalities and 190 surviving fingerlings was genotyped using 149 microsatellites with 145 successfully mapped markers covering 24 linkage groups (LGs). Analysis of QTL for both resistance against VNN and survival time was conducted using interval mapping. Five significant QTL located in four LGs and eight suggestive QTL in seven LGs were identified for resistance. Another five significant QTL in three LGs and five suggestive QTL in three LGs were detected for survival time. One significant QTL, spanning 3 cM in LG20, was identified for both resistance and survival time. These QTL explained 2.2-4.1% of the phenotypic variance for resistance and 2.2-3.3% of the phenotypic variance for survival time, respectively. Our results suggest that VNN resistance in Asian seabass is controlled by many loci with small effects. Our data provide information for fine mapping of QTL and identification of candidate genes for a better understanding of the mechanism of disease resistance.

Leaf area is an important parameter in oil palm breeding as it is positively correlated with oil yield. However, measurement of leaf area is tedious and also destructive. In the present study, a breeding population with 145 palms derived from a cross between Deli Dura and Avros Pisifera was used to construct a high-density linkage map and identify quantitative trait loci (QTL) for leaf area in oil palm. Using genotyping by sequencing, a linkage map containing 2413 SNPs was constructed. The total length of the linkage map was 1161.89 cM, with an average marker spacing of 0.48 cM. Based on the continuous phenotyping of leaf area from 2010 to 2015, two suggestive QTL for leaf area were mapped on chromosomes (Chr) 3 and 9. The allelic effects of the QTL associated with leaf area in the mapping population were estimated by linear regression using ordinary least squares method. The QTL on Chr 9 explained 13.3% of phenotypic variation for leaf area. A candidate gene, ARC5, within the narrow interval of QTL on Chr 9 was identified. The gene was significantly higher expressed in leaf than root and fruit of oil palm. This high-quality and SNP-based map supplies a base to fine map QTL for agronomic traits in oil palm, and the markers closely linked to the stable QTL may be used in marker-assisted selection in oil palm breeding.

Oil palm (Elaeis guineensis Jacq.) is the leading oil-producing crops and the most important edible oil resource worldwide. DNA markers and genetic linkage maps are essential resources for marker-assisted selection to accelerate genetic improvement. We conducted RAD-seq on an Illumina NextSeq500 to discover genome-wide SNPs, and used the SNPs to construct a linkage map for an oil palm (Tenera) population derived from a cross between a Deli Dura and an AVROS Pisifera. The RAD-seq produced 1,076 million single-end reads across the breeding population containing 155 trees. Mining this dataset detected 510,251 loci. After filtering out loci with low accuracy and more than 20% missing data, 11,394 SNPs were retained. Using these SNPs, in combination with 188 anchor SNPs and 123 microsatellites, we constructed a linkage map containing 10,023 markers covering 16 chromosomes. The map length is 2,938.2 cM with an average marker space of 0.29 cM. The large number of SNPs will supply ample choices of DNA markers in analysing the genetic diversity, population structure and evolution of oil palm. This high-density linkage map will contribute to mapping quantitative trait loci (QTL) for important traits, thus accelerating oil palm genetic improvement.

Single nucleotide polymorphisms (SNPs) are the most preferred DNA markers for genetic studies. We employed bioinformatics tools to screen genome-wide gene-associated SNPs from more than 487 million raw sequencing reads of RNA-seq data available for Asian seabass (Lates calcarifer). 48,048 putative SNPs in genes were identified. We improved a SNP genotyping method and used it for genotyping a subset of potential SNPs. Over 82 % of the putative SNPs were validated in tested populations. The improved SNP genotyping method was cost-effective and easy to conduct. The genome-wide SNPs discovered in genes supply an important resource for genetic studies of Asian seabass and comparative genomic studies. The improved genotyping method is suitable for small to middle-scale genotyping of SNPs.

Oil palm is the most productive oil crop in the world and composes 36% of the world production. However, the molecular mechanisms of hybrids vigor (or heterosis) between Dura, Pisifera and their hybrid progeny Tenera has not yet been well understood. Here we compared the temporal and spatial compositions of lipids and transcriptomes for two oil yielding organs mesocarp and endosperm from Dura, Pisifera and Tenera. Multiple lipid biosynthesis pathways are highly enriched in all non-additive expression pattern in endosperm, while cytokinine biosynthesis and cell cycle pathways are highly enriched both in endosperm and mesocarp. Compared with parental palms, the high oil content in Tenera was associated with much higher transcript levels of EgWRI1, homolog of Arabidopsis thaliana WRINKLED1. Among 338 identified genes in lipid synthesis, 207 (61%) has been identified to contain the WRI1 specific binding AW motif. We further functionally identified EgWRI1-1, one of three EgWRI1 orthologs, by genetic complementation of the Arabidopsis wri1 mutant. Ectopic expression of EgWRI1-1 in plant produced dramatically increased seed mass and oil content, with oil profile changed. Our findings provide an explanation for EgWRI1 as an important gene contributing hybrid vigor in lipid biosynthesis in oil palm.

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.

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.

An (AC)n-microsatellite-enriched library for Larimichthys crocea was constructed in this study. Primers for fifty simple sequence repeat (SSR) loci were synthesized and genotyped on 30 L. crocea individuals from Guanjingyang wild population (WP) in Fujian Province and 38 individuals from Ningbo cultured population (CP) in Zhejiang Province. Only 21 loci were successfully amplified and polymorphic in two populations. In WP, the observed heterozygosity (HO) ranged from 0.233 to 0.900 and the expected heterozygosity (HE) ranged from 0.326 to 0.893, with an average of 7.8 alleles/locus, the polymorphism information content (PIC) ranged from 0.283 to 0.866 (mean 0.731). In CP, the HO ranged from 0.189 to 0.892 and the HE ranged from 0.333 to 0.800, with an average of 4.4 alleles/locus. The probability test showed significant departures from HWE in 9 and 2 of the 21 loci in WP and in CP, respectively. Deficiency of heterozygotes at four loci showed the presence of null alleles (P

A high-density genetic map is essential for comparative genomic studies and fine mapping of QTL and can also facilitate genome sequence assembly. Here, a high density genetic map of Asian seabass was constructed with 3321 SNPs generated by sequencing 144 individuals in a F2 family. The length of the map was 1577.67 cM with an average marker interval of 0.52 cM. A high level of genomic synteny among Asian seabass, European seabass, Nile tilapia and stickleback was detected. Using this map, one genome-wide significant and five suggestive QTL for growth traits were detected in six linkage groups (i.e. LG4, LG5, LG11, LG13, LG14 and LG15). These QTL explained 10.5–16.0% of phenotypic variance. A candidate gene, ACOX1 within the significant QTL on LG5 was identified. The gene was differentially expressed between fast- and slow-growing Asian seabass. The high-density SNP-based map provides an important tool for fine mapping QTL in molecular breeding and comparative genome analysis.

High genetic diversity is thought to characterize successful invasive species, as the potential to adapt to new environments is enhanced and inbreeding is reduced. The red swamp crayfish, Procambarus clarkii, native to northeastern Mexico and south-central USA was introduced to Nanjing, China from Japan in 1929. Little is known about the genetic diversity and population structure of this species in China. We examined the genetic diversity and population structure of six P.clarkii populations using nine polymorphic microsatellites. Among the six populations, Nanjing population showed the highest allele number, allele richness and gene diversity, which is consistent with records indicating Nanjing may be the first site of introduction. In all six populations, significant heterozygote deficit was observed, suggesting founder effects and non-random mating. Analysis of bottleneck under infinite allele model, stepwise mutation model and two-phased model of mutation revealed evidence of a recent bottleneck in all these populations. Pairwise genetic distance analysis, AMOVA and assignment tests demonstrated high genetic differentiation between populations. Pairwise genetic distance did not fit the pairwise geographic distance, suggesting that human mediated dispersal have played a role in the population expansion and genetic differentiation.

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.

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.

Tilapia is the common name for nearly a hundred species of cichlid fish which originated in Africa (El-Sayed 2006). It is one of the most important food fish genera in the world (El-Sayed 2006). A few red tilapia variants have been commonly seen in fish markets. While the genetic background of these varieties is not well documented, their derivation is generally attributed to the crossbreeding between mutant reddish-orange Mozambique tilapia (Oreochromis mossambicus) with other tilapia species like Nile tilapia (O. niloticus) and blue tilapia (O. aureus). Florida red tilapia is thought to be the result of blue tilapia mating with Mozambique tilapia (El-Sayed 2006). In Nile and Mozambique tilapia, microsatellites have been isolated (Kocher et al. 1998; Yue and Orban 2002) and applied in genetic studies (e.g. Kocher et al. 1998), whereas in Florida red tilapia, no microsatellites are available. Preliminary tests of cross-species amplification of 30 microsatellites isolated from Nile and Mozambique tilapia revealed that around 30% of tested microsatellites could not amplify products in red tilapia (Feng Liu and Gen Hua Yue, unpublished data). Therefore, to facilitate a better understanding of the genetic variations in tilapia species, we have isolated 19 polymorphic microsatellites from Florida red tilapia, and characterized them in Nile and Mozambique tilapia. These 19 microsatellites were all polymorphic in Florida red tilapia with an average allele number of 4.90 per locus, and average expected and observed heterozygosities of 0.62 and 0.74 respectively. Fourteen out of the 19 microsatellites amplified polymorphic products in Mozambique tilapia and all 19 amplified in Nile tilapia, respectively. Fourteen out of 19 and 13 out of 15 microsatellites were in HardyWeinberg equilibrium (HWE) in Nile and Mozambique

The brown mud mussel Glauconome rugosa is a delicious mussel food species. Due to over-harvesting for food in the wild, some populations may be endangered. Seventeen di- and tetranucleotide microsatellites were isolated and characterized using 24 wild caught individuals. The average allele number of the 17 microsatellites was 8/locus with a range of 5–14, and the expected heterozygosity averaged at 0.80, ranging from 0.65 to 0.91. Eleven out of 17 microsatellites agreed with the Hardy–Weinberg equilibrium (HWE) and showed independent segregation. Extreme size difference (over 50 bp) between alleles was seen at five of the 17 loci. These microsatellites will facilitate the studying of the genetic diversity and population structure of wild populations of G. rugosa.

During genotyping of 38 microsatellites for QTL (quantitative trait loci) mapping in three F2 swine populations, five mutant alleles were detected in a total of 66,436 parent-offspring transfers of microsatellite alleles, which gives an overall mutation rate of 7.52 x 10(-5) per locus per generation. No significant (P > 0.05) association between mutation rates and other factors (i.e., GC contents in the flanking regions, heterozygosity, and repeat number) was revealed. Detailed sequencing showed that four out of five mutant alleles were caused by insertions of one to five repeats, respectively. The other mutant allele was produced by either an insertion of three repeats or a change of 30 base pairs (a deletion of 16 CT repeats and an insertion of one CA repeat). An insertion of one base pair in the flanking region of a microsatellite was also detected. Together, these data indicate that expansions are more common than contractions among microsatellites and that the mutation processes are very complicated, do not fit with the strict stepwise mutation model and may vary from locus to locus.

The Asian swamp eel Monopterus albus is distributed in rice fields, rivers and ponds. Due to large scale use of pesticides in rice fields and destruction of natural habitats, the wild populations are decreasing very quickly. To facilitate studies on genetic diversity, we isolated 16 polymorphic novel microsatellites and characterized them using 48 individuals from two natural populations from two locations. The average allele number of the 16 microsatellites was 10.1/locus, ranging from 2 to 20. The expected heterozygosity was 0.79 whereas the observed heterozygosity was 0.65. In both populations, all markers conformed to HWE and showed no significant linkage between pairwise markers. The genetic differentiation between the two populations was significant (F ST = 0.24, P

The golden pompano (Trachinotus blochii) is a marine fish species in tropical regions. No information about genetic variation and population structure of wild populations is available. A first set of 12 polymorphic microsatellites isolated from this species were characterized. The number of alleles ranged from 5 to 14 with an average of 8.0 alleles per locus. All 12 markers conformed to HWE and were in linkage equilibrium. These 12 microsatellite markers supply the first set of co-dominant DNA markers for studying population genetics of the species T. blochii.

Twelve novel polymorphic microsatellites were isolated from the endangered Hawksbill turtle (Eretmochelys imbricate). Eight of 12 markers were used to study genetic diversity of two sea turtle species: E. imbricate and green sea turtle (C. mydas). In E. imbricate, the average allele number of the eight microsatellites was 6.25/locus with a range of 3–13. The average expected and observed heterozygosity was 0.66 and 0.63 respectively. In C. mydas, the average allele number of the eight markers was 11.63/locus. The observed heterozyosity (0.68) was lower than the expected heterozyosity (0.79). Most of 12 microsatellites amplified specific and polymorphic PCR products in other six turtle species. Hence, the developed microsatellites would facilitate studies on genetic diversity and population structure of E. imbricate and other marine turtle species.

Microsatellites are popular molecular markers in genetic and evolutionary studies. Their mutational dynamics have been extensively studied in humans and fruit flies, but few data were available in fish. By genotyping 55 individuals of a F1 pedigree, we investigated the mutation rates and patterns of 49 microsatellites in one of the most important fresh water fish species, the common carp (Cyprinus carpio L.). The overall mutation rate of the 49 loci was 5.56 x 10(-4)/locus/generation (95% confidence interval 1.52 x 10(-4) and 1.63 x 10(-3)). The change of allele size was between +2 to -5 repeat units, assuming that the mutation allele arose from the parental allele most similar in size to the mutant.

Asian arowana (Scleropages formosus) is a highly endangered fish species listed in Appendix 1 of CITES since 1980. Fourteen novel polymorphic microsatellites were isolated from a CA-enriched partial genomic DNA library, and were used for studying genetic diversity of 41 red arowanas from the wild population. Surprisingly the average allele number of the 14 polymorphic microsatellites was as high as 12.4/locus. The average observed heterozygosity was 0.78 ranging from 0.51 to 0.95, and the gene diversity was quite high (0.78). All these data suggest that high level of genetic diversity existed in the red Asian arowana population.

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.

Background High density linkage maps are essential for comparative analysis of synteny, fine mapping of quantitative trait loci (QTL), searching for candidate genes and facilitating genome sequence assembly. However, in most foodfish species, marker density is still low. We previously reported a first generation linkage map with 240 DNA markers and its application to preliminarily map QTL for growth traits in Asian seabass (Lates calcarifer). Here, we report a high-resolution linkage map with 790 microsatellites and SNPs, comparative analysis of synteny, fine-mapping of QTL and the identification of potential candidate genes for growth traits. Results A second generation linkage map of Asian seabass was developed with 790 microsatellite and SNP markers. The map spanned a genetic length of 2411.5 cM, with an average intermarker distance of 3.4 cM or 1.1 Mb. This high density map allowed for comparison of the map with Tetraodon nigroviridis genome, which revealed 16 synteny regions between the two species. Moreover, by employing this map we refined QTL to regions of 1.4 and 0.2 cM (or 400 and 50 kb) in linkage groups 2 and 3 in a population containing 380 progeny; potential candidate genes for growth traits in QTL regions were further identified using comparative genome analysis, whose effects on growth traits were investigated. Interestingly, a QTL cluster at Lca371 underlying growth traits of Asian seabass showed similarity to the cathepsin D gene of human, which is related to cancer and Alzheimer's disease. Conclusions We constructed a high resolution linkage map, carried out comparative mapping, refined the positions of QTL, identified candidate genes for growth traits and analyzed their effects on growth. Our study developed a framework that will be indispensable for further identification of genes and analysis of molecular variation within the refined QTL to enhance understanding of the molecular basis of growth and speed up genetic improvement of growth performance, and it also provides critical resource for future genome sequence assembly and comparative genomics studies on the evolution of fish genomes.

Four microsatellites were identified by screening the DNA sequences of Asian seabass (Lates calcarifer) deposited to GenBank. Two markers each are located in the growth hormone gene (GH) and in the insulin-like growth factor II gene ( IGF-2), respectively. The markers were characterized by genotyping 34 Asian seabass individuals. All 4 microsatellites showed polymorphism: the number of alleles per locus ranged from 2 to 11 (average, 5.0), while the expected heterozygosity ranged from 0.51 to 0.85 (average, 0.63) at the 4 loci. Cross-priming with all 4 primer pairs was tested in species belonging to 5 different genera, but no bands were amplified. These microsatellites are the first genomic DNA markers characterized in L. calcarifer; thus they may be valuable for research and aquaculture production of this species.

Background Grass carp (Ctenopharyngodon idella) belongs to the family Cyprinidae which includes more than 2000 fish species. It is one of the most important freshwater food fish species in world aquaculture. A linkage map is an essential framework for mapping traits of interest and is often the first step towards understanding genome evolution. The aim of this study is to construct a first generation genetic map of grass carp using microsatellites and SNPs to generate a new resource for mapping QTL for economically important traits and to conduct a comparative mapping analysis to shed new insights into the evolution of fish genomes. Results We constructed a first generation linkage map of grass carp with a mapping panel containing two F1 families including 192 progenies. Sixteen SNPs in genes and 263 microsatellite markers were mapped to twenty-four linkage groups (LGs). The number of LGs was corresponding to the haploid chromosome number of grass carp. The sex-specific map was 1149.4 and 888.8 cM long in females and males respectively whereas the sex-averaged map spanned 1176.1 cM. The average resolution of the map was 4.2 cM/locus. BLAST searches of sequences of mapped markers of grass carp against the whole genome sequence of zebrafish revealed substantial macrosynteny relationship and extensive colinearity of markers between grass carp and zebrafish. Conclusions The linkage map of grass carp presented here is the first linkage map of a food fish species based on co-dominant markers in the family Cyprinidae. This map provides a valuable resource for mapping phenotypic variations and serves as a reference to approach comparative genomics and understand the evolution of fish genomes and could be complementary to grass carp genome sequencing project.

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.

Background Microsatellites in cDNA are useful as molecular markers because they represent transcribed genes and can be used as anchor markers for linkage and comparative mapping, as well as for studying genome evolution. Microsatellites in cDNA can be detected in existing ESTs by data mining. However, in most fish species, no ESTs are available or the number of ESTs is limited, although fishes represent half of the vertebrates on the earth. We developed a simple and efficient method for isolation of microsatellites from cDNA in fish. Results The method included normalization of 150 ng cDNA using 0.5 U duplex-specific nuclease (DSN) at 65°C for 30 min, enrichment of microsatellites using biotinylated oligonucleotides and magnetic field, and directional cloning of cDNA into a vector. We tested this method to enrich CA- and GA-microsatellites from cDNA of Asian seabass, and demonstrated that enrichment of microsatellites from normalized cDNA could increased the efficiency of microsatellite isolation over 30 times as compared to direct sequencing of clones from cDNA libraries. One hundred and thirty-nine (36.2%) out of 384 clones from normalized cDNA contained microsatellites. Unique microsatellite sequences accounted for 23.6% (91/384) of sequenced clones. Sixty microsatellites isolated from cDNA were characterized, and 41 were polymorphic. The average allele number of the 41 microsatellites was 4.85 ± 0.54, while the expected heterozygosity was 0.56 ± 0.03. All the isolated microsatellites inherited in a Mendelian pattern. Conclusion Normalization of cDNA substantially increased the efficiency of enrichment of microsatellites from cDNA. The described method for isolation of microsatellites from cDNA has the potential to be applied to a wide range of fish species. The microsatellites isolated from cDNA could be useful for linkage and comparative mapping, as well as for studying genome evolution.