Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish.

BACKGROUND: Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens. RESULTS: We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication. CONCLUSIONS: Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution.

Genomic Basis of Striking Fin Shapes and Colors in the Fighting Fish.

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.

Construction of high-resolution recombination maps in Asian seabass.

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.

Genome-wide identification of markers for selecting higher oil content in oil palm.

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 < 0.01) higher than that of trees without any beneficial QTL genotypes (average O/B of 28.24%). QTL combinations showed that the higher the number of QTL with beneficial genotypes, the higher the resulting average O/B in the breeding population. CONCLUSIONS: A linkage map with 1480 DNA markers was constructed and used to identify QTL for oil content traits. Pyramiding the identified QTL with beneficial genotypes associated with oil content traits using DNA markers has the potential to accelerate genetic improvement for oil yield in the breeding population of oil palm.

Genetic variability, local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass.

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.

RNA-Seq revealed the impairment of immune defence of tilapia against the infection of Streptococcus agalactiae with simulated climate warming.

Global warming is one of the causes of disease outbreaks in fishes. Understanding its mechanisms is critical in aquaculture and fisheries. We used tilapia to study the effects of a high temperature on the infection of a bacterial pathogen Streptococcus agalactiae using RNA-Seq. We found that the dissolved oxygen level in water at 32 °C is lower than at 22 °C, and tilapia infected with the pathogen died more rapidly at 32 °C. The gene expression profiles showed significant differences in fish raised under different conditions. We identified 126 and 576 differentially expressed genes (DEGs) at 4 and 24 h post infection at 22 °C, respectively, whereas at 32 °C, the data were 312 and 1670, respectively. Almost all responding pathways at 22 °C were involved in the immune responses, whereas at 32 °C, the enriched pathways were not only involved in immune responses but also involved in oxygen and energy metabolisms. We identified significant signals of immunosuppression of immune responses at 32 °C. In addition, many of the enriched transcription factors and DEGs under positive selection were involved in immune responses, oxygen and/or energy metabolisms. Our results suggest that global warming could reduce the oxygen level in water and impair the defence of tilapia against bacterial infection.

The intestinal microbiome of fish under starvation.

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.

The MCP-8 gene and its possible association with resistance to Streptococcus agalactiae in tilapia.

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 (p < 0.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 (p < 0.05) with resistance against S. agalactiae. However, we found no association between four SNPs and growth traits (p > 0.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.

Characterization of the LECT2 gene and its associations with resistance to the big belly disease in Asian seabass.

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.

Wittig/B─H insertion reaction: A unique access to trisubstituted Z-alkenes.

The Wittig reaction, which is one of the most effective methods for synthesizing alkenes from carbonyl compounds, generally gives thermodynamically stable E-alkenes, and synthesis of trisubstituted Z-alkenes from ketones presents notable challenges. Here, we report what we refer to as Wittig/B─H insertion reactions, which innovatively combine a Wittig reaction with carbene insertion into a B─H bond and constitute a promising method for the synthesis of thermodynamically unstable trisubstituted Z-boryl alkenes. Combined with the easy transformations of boryl group, this methodology provides efficient access to a variety of previously unavailable trisubstituted Z-alkenes and thus provides a platform for discovery of pharmaceuticals. The unique Z-selectivity of the reaction is determined by the maximum overlap of the orbitals between the B─H bond of the borane adduct and the alkylidene carbene intermediate in the transition state.

Saliva is more sensitive than nasopharyngeal or nasal swabs for diagnosis of asymptomatic and mild COVID-19 infection.

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.

Mapping and Validating QTL for Fatty Acid Compositions and Growth Traits in Asian Seabass.

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.

QTL Mapping for Resistance to Iridovirus in Asian Seabass Using Genotyping-by-Sequencing.

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.

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass.

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.

Genome-wide methylation analysis identified sexually dimorphic methylated regions in hybrid tilapia.

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.

Mapping QTL for Resistance Against Viral Nervous Necrosis Disease in Asian Seabass.

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.

Signatures of selection in tilapia revealed by whole genome resequencing.

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.

Construction of a high-density linkage map and fine mapping of QTL for growth in Asian seabass.

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.

Molecular characterization and mapping of Fgf21 gene in a foodfish species asian seabass.

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.