Search

The palm family (Arecaceae), consisting of ∼ 2600 species, is the third most economically important family of plants. The African oil palm (Elaeis guineensis) is one of the most important palms. However, the genome sequences of palms available are still limited and fragmented. Here, we report a high-quality chromosome-level reference genome of an oil palm Dura. The genome of 1.7 Gb was assembled by integrating long reads with ∼ 150 × genome coverage. The assembled genome covered 94.5% of the estimated genome size, within which 91.6% were assigned into 16 pseudochromosomes and 73.7% were repetitive sequences. Relying on the conserved synteny with oil palm, the existing draft genome sequences of both date palm and coconut were further assembled into chromosomal level. Transposon burst, particularly long terminal repeat retrotransposons (LTRs) retrotransposons, following the last whole-genome duplication (WGD), likely explains genome size variation across palms. Convergent evolution of fruit colors tends to eliminate the roles of the virescens gene in controlling accumulation of anthocyanins in exocarp of ripe fruit of palms. Recent duplications of high tandemly repeated pathogenesis-related proteins (PRs) from the same tandem arrays played an important role in defense responses to Ganoderma. Whole genome re-sequencing of both ancestral African and introduced oil palms in Southeast Asia revealed that genes under putative selection were notably associated with stress responses, suggesting adaptation to stresses in the new habitat. The genomic resources and insights gained in this study could be exploited for accelerating genetic improvement and understanding the evolution of palms.

Sexual dimorphism is common across the animal kingdom. Knowledge of the mechanisms of sexual size dimorphism is limited although it is important in biology and aquaculture. Tilapia is the common name for ~ 100 species of cichlid fish. Some are important aquaculture species and males outgrow females. To gain novel insights into the mechanisms underlying sexual size dimorphism, we analyzed the differences of brain transcriptomes between males and females in Mozambique tilapia and studied the function of the pro-opiomelanocortin (Pomc) gene in tilapia and zebrafish. The transcriptome analysis identified 123, 55, and 2706 sex-biased genes at 5, 30, and 90 dph (days post-hatch), respectively, indicating sexual dimorphism of gene expressions in the brain. The expression of Pomc in the tilapia brain was a female-biased at 30, 90, and 120 dph. An analysis of the DNA sequence located upstream of the tilapia Pomc transcriptional start site identified two estrogenic response elements. In vitro luciferase assay of the two elements revealed that β-estradiol significantly enhanced the expression of luciferase activity, suggesting that the expression of Pomc is mediated by estrogen. We knocked out Pomc in zebrafish using Crispr/Cas-9. The Pomc-knockout zebrafish showed faster growth and higher sensitivity to feeding as compared to the wild-type fish. Taken together, our results indicate that Pomc contributes to sexual size dimorphism and suggest that the high estrogen level in females promotes the expression of Pomc and suppresses feeding in female tilapias, which leads to the slower growth of female tilapias. Nanyang Technological University Submitted/Accepted version This study is part of the PhD thesis of Zi Yi Wan, which was supported by Temasek Life Sciences Laboratory and Nanyang Technological University, Singapore.

We aimed to test the sensitivity of naso-oropharyngeal saliva and self-administered nasal (SN) swab compared to nasopharyngeal (NP) swab for COVID-19 testing in a large cohort of migrant workers in Singapore. We also tested the utility of next-generation sequencing (NGS) for diagnosis of COVID-19. Saliva, NP and SN swabs were collected from subjects who presented with acute respiratory infection, their asymptomatic roommates, and prior confirmed cases who were undergoing isolation at a community care facility in June 2020. All samples were tested using RT-PCR. SARS-CoV-2 amplicon-based NGS with phylogenetic analysis was done for 30 samples. We recruited 200 subjects, of which 91 and 46 were tested twice and thrice respectively. In total, 62.0%, 44.5%, and 37.7% of saliva, NP and SN samples were positive. Cycle threshold (Ct) values were lower during the earlier period of infection across all sample types. The percentage of test-positive saliva was higher than NP and SN swabs. We found a strong correlation between viral genome coverage by NGS and Ct values for SARS-CoV-2. Phylogenetic analyses revealed Clade O and lineage B.6 known to be circulating in Singapore. We found saliva to be a sensitive and viable sample for COVID-19 diagnosis.

The viral load of SARS-CoV-2 in clinical samples as measured by the primary diagnostic tool of RT-PCR is an imperfect readout for infection potential as most targeted assays designed for sensitivity, indiscriminately detect short and long RNA fragments, although infectivity is embodied only in the whole virus and its intact genome. Here, we used next-generation sequencing (NGS) to characterize 155 clinical samples and show sensitive and quantitative detection of viral RNA which confirmed subgenomic RNA in 57.6% of samples and provided a novel method to determine relative integrity of viral RNA in samples. The relative abundance of long fragments quantified as a viral fragmentation score was positively associated with viral load and inversely related to time from disease onset. An empirically determined score cut-off for presence of substantially fragmented RNA was able to identify 100% of samples collected after 8 days of illness with poor infection potential in line with current clinical understanding of infectiousness of SARS-CoV-2. The quantification of longer fragments in addition to existing short targets in an NGS or RT-PCR-based assay could provide a valuable readout of infection potential simultaneous to the detection of any fragments of SARS-CoV-2 RNA in test samples.

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)

We aimed to test the sensitivity of naso-oropharyngeal saliva and self-administered nasal (SN) swab compared to nasopharyngeal (NP) swab for COVID-19 testing in a large cohort of migrant workers in Singapore. We also tested the utility of next-generation sequencing (NGS) for diagnosis of COVID-19. Saliva, NP and SN swabs were collected from subjects who presented with acute respiratory infection, their asymptomatic roommates, and prior confirmed cases who were undergoing isolation at a community care facility in June 2020. All samples were tested using RT-PCR. SARS-CoV-2 amplicon-based NGS with phylogenetic analysis was done for 30 samples. We recruited 200 subjects, of which 91 and 46 were tested twice and thrice respectively. In total, 62.0%, 44.5%, and 37.7% of saliva, NP and SN samples were positive. Cycle threshold (Ct) values were lower during the earlier period of infection across all sample types. The percentage of test-positive saliva was higher than NP and SN swabs. We found a strong correlation between viral genome coverage by NGS and Ct values for SARS-CoV-2. Phylogenetic analyses revealed Clade O and lineage B.6 known to be circulating in Singapore. We found saliva to be a sensitive and viable sample for COVID-19 diagnosis.

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.

Growth is an economically important trait in aquaculture. To improve growth trait of the Asian seabass (Lates calcarifer) we have been carrying out, since 2004, a selective breeding program. This study focuses on growth traits in the F2 fish generation, comprised of offspring from 23 mass crosses from 383 F1 brooders. Using genotyping analysis for 10 microsatellites from both brood stock and progeny we have reconstructed the pedigree of each mass-cross. For F2 generation at 90 days post hatch (dph), we have recorded body weight (BW) for 12,117 individuals and total length (TL), standard length (SL) and condition factor (Ktl and Ksl) for 3530 individuals; and all five traits for 2136 individuals at 270 dph. At 90 dph the average BW was 46.88 ± 20.95 g. Combining pedigree information, recorded growth traits and Restricted Maximum Likelihood method, we have estimated that the narrow sense heritability (h2) in F1 fish for BW, TL, SL, Ktl and Ksl was at, 90 dph, 0.12 ± 0.03, 0.11 ± 0.03, 0.10 ± 0.03, 0.20 ± 0.04 and 0.11 ± 0.03, respectively and, at 270 dph, 0.34 ± 0.07, 0.32 ± 0.07, 0.30 ± 0.06, 0.13 ± 0.04 and 0.11 ± 0.04, respectively. At 90 dph the realised heritability for BW was 0.13. Comparing with F1 generation, the growth performance of F2 fish was increased by 14.4%. Heritability of growth traits will be useful for future genetic improvement programmes of the Asian seabass. NRF (Natl Research Foundation, S’pore) Published version

Background Oil palm (Elaeis guineensis, Jacq.) is the most important source of edible oil. The improvement of oil yield is currently slow in conventional breeding programs due to long generation intervals. Marker-assisted selection (MAS) has the potential to accelerate genetic improvement. To identify DNA markers associated with oil content traits for MAS, we performed quantitative trait loci (QTL) mapping using genotyping by sequencing (GBS) in a breeding population derived from a cross between Deli Dura and Ghana Pisifera, containing 153 F1 trees. Results We constructed a high-density linkage map containing 1357 SNPs and 123 microsatellites. The 16 linkage groups (LGs) spanned 1527 cM, with an average marker space of 1.03 cM. One significant and three suggestive QTL for oil to bunch (O/B) and oil to dry mesocarp (O/DM) were mapped on LG1, LG8, and LG10 in a F1 breeding population, respectively. These QTL explained 7.6–13.3% of phenotypic variance. DNA markers associated with oil content in these QTL were identified. Trees with beneficial genotypes at two QTL for O/B showed an average O/B of 30.97%, significantly (P

Asian seabass is an important food fish species. While improving growth, increasing the nutritional value is important, omega-3 fatty acids are indispensable to human health. Identifying and validating DNA markers associated with traits is the first step towards marker-assisted selection (MAS). We quantified 13 different fatty acids and three growth traits in 213 F2 Asian seabass from a family at the age 270 days post hatch, and screened QTL for these traits. The content of total fatty acids in 100 g flesh was 2.57 ± 0.80 g, while the proportions of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) were 16.96 ± 2.20% and 5.42 ± 0.90%, respectively. A linkage map with 2424 SNPs was constructed and used for QTL mapping. For fatty acid compositions, 14 significant QTL were identified on three linkage groups (LG5, LG11 and LG14), with phenotypic variance explained (PVE) from 12.8 to 24.6%. Thirty-nine suggestive QTL were detected on 16 LGs. Two significant QTL for EPA were identified on LG5 and LG14, with PVE of 15.2% and 15.1%, respectively. No significant QTL was identified for DHA. For growth traits, six significant and 13 suggestive QTL were identified on two and seven LGs, respectively. Only a few significant QTL for fatty acids overlapped with previously mapped QTL for these traits, suggesting that most QTL detected in a family are family-specific and could only be used in MAS in the family per se. To facilitate population-wide molecular breeding, more powerful methods (e.g. GWAS) should be used to identify SNPs for genomic selection.

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.

Asian seabass is an important marine foodfish species in Asia-Pacific region. Previous studies using microsatellites and mitochondrial DNA sequences showed that the populations in Australia are significantly different from those in Southeast Asia. However, within Australian and Southeast Asian populations, there are some conflicts regarding to population differentiations. It is essential for us to use more DNA markers to examine the population structure of Asian seabass. We collected Asian seabass samples from seven wild populations distributed across Indonesian Archipelago from Southeast Asia to Australia and Papua New Guinea. We measured the genetic diversity within each population and analyzed the genetic relationships among the wild populations using 96 SNPs. We detected significant population differentiations among the seven wild populations of Asian seabass. The genetic diversity in Southeast Asian populations was significantly higher than that in Australian and Papua New Guinean populations. There were significant genetic differentiations among the Australian and Papua New Guinean populations while the differentiation in Southeast Asian populations was minor. In addition, we identified footprints of directional selection at three loci using outlier tests. Our data supply useful information about the management of wild stocks of Asian seabass and the use of wild Asian seabass to start a breeding program.

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.

Identifying quantitative trait loci (QTL) for viral disease resistance is of particular importance in selective breeding programs of fish species. Genetic markers linked to QTL can be useful in marker-assisted selection (MAS) for elites resistant to specific pathogens. Here, we conducted a genome scan for QTL associated with Singapore grouper iridovirus (SGIV) resistance in an Asian seabass (Lates calcarifer) family, using a high-density linkage map generated with genotyping-by-sequencing. One genome-wide significant and three suggestive QTL were detected at LG21, LG6, LG13, and LG15, respectively. The phenotypic variation explained (PVE) by the four QTL ranged from 7.5 to 15.6%. The position of the most significant QTL at LG21 was located between 31.88 and 36.81 cM. The SNP marker (SNP130416) nearest to the peak of this QTL was significantly associated with SGIV resistance in an unrelated multifamily population. One candidate gene, MECOM, close to the peak of this QTL region, was predicted. Evidence of alternative splicing was observed for MECOM and one specific category of splicing variants was differentially expressed at 5 days post-SGIV infection. The QTL detected in this study are valuable resources and can be used in the selective breeding programs of Asian seabass with regard to resistance to SGIV.

Background: A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing. Results: A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males. Conclusions: These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass. NRF (Natl Research Foundation, S’pore) Published version

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.

Single nucleotide polymorphisms (SNPs) are important DNA markers for genetic studies. However, the number of SNPs is limited in important aquaculture species. We have identified 23,535 high quality SNPs in genes of Nile tilapia ( Oreochromis niloticus ) by analyzing the available RNA-seq data. Among them, 17,699 SNPs were mapped onto 22 linkage groups (LG) using in silico mapping. The number of mapped SNPs ranged from 332 on LG2 to 1148 on LG16_21 with an average number of 805 SNPs/LG. The mapped SNPs were located in 7146 genes. This SNP resource will be useful in monitoring the genetic diversity and mapping QTL for genetic improvement of tilapia species. We further studied the relationships among individuals from four populations (GIFT, O. niloticus , Red tilapia and O. mossambicus ) using 101 SNPs from LG23. Ninety-four SNPs (~ 93%) were polymorphic in the tested populations. Red tilapia (77%) had the highest number of informative markers than GIFT (69%), Nile tilapia (56%), and Mozambique tilapia (12%). Significant population structure was detected through phylogenic and population structure analysis. The Red tilapia individuals were closely related to Mozambique tilapia, but were split into several subclades, confirming that red tilapia is a hybrid admixture of different resources of tilapias.

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.

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.

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South-East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South-East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. FST and principal components analysis revealed striking divergence between South-East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South-East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South-East Asia during mid-Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the ‘genomic islands’ scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid-Pleistocene.

Natural selection and selective breeding for genetic improvement have left detectable signatures within the genome of a species. Identification of selection signatures is important in evolutionary biology and for detecting genes that facilitate to accelerate genetic improvement. However, selection signatures, including artificial selection and natural selection, have only been identified at the whole genome level in several genetically improved fish species. Tilapia is one of the most important genetically improved fish species in the world. Using next-generation sequencing, we sequenced the genomes of 47 tilapia individuals. We identified a total of 1.43 million high-quality SNPs and found that the LD block sizes ranged from 10–100 kb in tilapia. We detected over a hundred putative selective sweep regions in each line of tilapia. Most selection signatures were located in non-coding regions of the tilapia genome. The Wnt signaling, gonadotropin-releasing hormone receptor and integrin signaling pathways were under positive selection in all improved tilapia lines. Our study provides a genome-wide map of genetic variation and selection footprints in tilapia, which could be important for genetic studies and accelerating genetic improvement of tilapia.

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.

The polyacrylamide hydrogel (PAMG) has been used in cosmetology in China, Ukraine and Russia since 1990s. Because the monomer acrylamide(AM) used to produce PAMG has been implicated as a potential mutagen and reproductive toxicant[1,2], it is important to accurately determine the amount of residual AM monomer in the PAMG. In this study, a quick, practical and simple method to determine AM is presented with respect to the hydrogel. AM is analysed quantitatively by ODS-3 column with ultraviolet (UV) absorbance detector. AM is separated from interferential component with an aqueous solution of 0.9%NaCl (NS) adjusted at pH~3.7 using hydrochloric acid and then detected at a UV wavelength of 210 nm. The results show that ODS-3 is effective approach for quantifying AM concentrations in PAMG. This method has a lower detection limit of 0.003µg/ml and a linear response range of 0.003 and 0.9 µg/ml (depending on the range required for analysis). Precision studies give coefficients of variation of

Polyacrylamide (PAM) was usually atoxic, stable. Its hydrogel (PAMG) has been used in plastic and aesthetic surgery more than 10 years in P.R.China, Ukraine and Russia. But there were some complications with PAMG injected in patients. Considering the complicacy in vivo, it was necessary to study the PAMG’s stability. In this paper, H2O2 was added in PAMG in vitro and acrylamide (AM) content after oxidating was determined using HPLC method. Detection limit for AM can be achieved in the parts-per-billion rang. The AM content in supernatant at 0.5, 1,2 and 3 hr after oxidating of PAMG was 1.27e-8g/ml, 1.35e-8g/ml, 1.03e-7g/ml and 2.74e-7g/ml, respectively. The AM content in supernatant of PAMG was 7.70e-9g/ml. These results indicated that PAMG can be degraded to AM by hydroxyl radical. Under the same condition, the AM content was stable. So, AM can exist for a while when it was produce by degrading of PAMG. That can increase the toxicity of PAMG.

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.

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.

Background Starvation not only affects the nutritional and health status of the animals, but also the microbial composition in the host’s intestine. Next-generation sequencing provides a unique opportunity to explore gut microbial communities and their interactions with hosts. However, studies on gut microbiomes have been conducted predominantly in humans and land animals. Not much is known on gut microbiomes of aquatic animals and their changes under changing environmental conditions. To address this shortcoming, we determined the microbial gene catalogue, and investigated changes in the microbial composition and host-microbe interactions in the intestine of Asian seabass in response to starvation. Results We found 33 phyla, 66 classes, 130 orders and 278 families in the intestinal microbiome. Proteobacteria (48.8%), Firmicutes (15.3%) and Bacteroidetes (8.2%) were the three most abundant bacteria taxa. Comparative analyses of the microbiome revealed shifts in bacteria communities, with dramatic enrichment of Bacteroidetes, but significant depletion of Betaproteobacteria in starved intestines. In addition, significant differences in clusters of orthologous groups (COG) functional categories and orthologous groups were observed. Genes related to antibiotic activity in the microbiome were significantly enriched in response to starvation, and host genes related to the immune response were generally up-regulated. Conclusions This study provides the first insights into the fish intestinal microbiome and its changes under starvation. Further detailed study on interactions between intestinal microbiomes and hosts under dynamic conditions will shed new light on how the hosts and microbes respond to the changing environment.