Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies. In order to provide the best genome coverage for the analysis of performance and production traits, a large number of relatively evenly distributed SNPs are needed. Gene-associated SNPs may fulfill these requirements of large numbers and genome wide distribution. In addition, gene-associated SNPs could themselves be causative SNPs for traits. The objective of this project was to identify large numbers of gene-associated SNPs using high-throughput next generation sequencing.

Yellowfin seabream (Acanthopagrus latus), a protandrous hermaphroditic fish, is a good model for studying the mechanism of sex reversal. However, limited knowledge is known about the genetic information related to reproduction and sex differentiation in this species. Here, we performed de novo transcriptome sequencing analysis of the testis, ovotestis, and ovary to identify sex-related genes in yellowfin seabream. The results assembled 71,765 unigenes in which 16,126 and 17,560 unigenes were differentially expressed in the ovotestis and ovary compared to the testis, respectively. The most differentially expressed gene (DEG)-enriched Kyoto Encyclopedia of Genes and Genomes and GO pathways were closely associated with the synthesis of sex steroid hormones. Functional analyses identified 55 important sex-related DEGs, including 32 testis-biased DEGs (dmrt1, amh, and sox9, etc.), 20 ovary-biased DEGs (cyp19a, foxl2, and wnt4, etc.), and 3 ovotestis-biased DEGs (lhb, dmrt2, and foxh1). Furthermore, the testis-specific expression of dmrt1 and the brain-pituitary-ovary axis expression of foxl2 were characterized, suggesting that they might play important roles in sex differentiation in yellowfin seabream. Our present work provided an important molecular basis for elucidating the mechanisms underlying sexual transition and reproductional regulation in yellowfin seabream.

Hepatic stellate cells (HSCs) play a crucial role in the progression of liver fibrosis, which can be considered as the specific therapeutic target of anti-fibrotic treatment. Targeted induction of HSCs to hepatocytes via delivery of clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (dCas9) system holds promise for hepatic fibrosis treatment. Our study here revealed that CRISPR/dCas9-VP64 system encapsulated in AML12 cell-derived exosomes could efficiently and successfully be delivered into the HSCs. In turn, the CRISPR/dCas9-VP64 system loaded in the exosomes can be efficiently released into the HSCs. As a proof-of-concept study, gRNA against hepatocyte nuclear factor 4α (HNF4α) together with the delivery of CRISPR/dCas9-VP64 system induced the HSCs to hepatocyte-like phenotype. In conclusion, our study here revealed that CRISPR/dCas9-VP64 system encapsulated in AML12 cell-derived exosomes could be functional in HSCs, emerging as a gene therapy strategy for hepatic fibrosis.

Duplicated genes prevail in vertebrates and are important in the acquisition of new genes and novelties. Whole genome duplication (WGD) is one of the sources of duplicated genes. It can provide raw materials for natural selection by increasing the flexibility and complexity of the genome. WGDs are the driving force for the evolution of vertebrates and contribute greatly to their species diversity, especially in fish species with complicated WGD patterns. Here, we constructed the DupScan database (https://dupscan.sysumeg.com/) by integrating 106 chromosomal-level genomes, which can analyze and visualize synteny at both the gene and genome scales, visualize the Ka, Ks, and 4DTV values, and browse genomes. DupScan was used to perform functional adaptation for the intricate WGD investigation based on synteny matching. DupScan supports the analysis of five WGD rounds (R): VGD2 (vertebrate genome duplication 2), Ars3R (Acipenser-ruthenus-specific 3R), Pss3R (Polyodon-spathula-specific 3R), Ts3R (teleost-specific duplication 3R), Ss4R (salmonid-specific 4R), and Cs4R (carp-specific 4R). DupScan serves as one-stop analysis platform for synteny and WGD research in which users can analyze and predict synteny and WGD patterns across 106 species of whole genome sequences. This further aided us in elucidating genome evolutionary patterns across over 60,000 vertebrate species with synteny and WGD events.

Invertebrate shrimp are one of the dominant benthic macrofaunae in the deep-sea environment. The microbiota of shrimp intestine can contribute to the adaptation of their host. The impact of surrounding sediment on intestinal microbiota has been observed in cultured shrimp species, but needs to be further investigated in deep-sea shrimp. The characterization of bacterial, archaeal, and fungal community structure and their interrelationships is also limited. In this study, wild-type deep-sea shrimp and the surrounding sediment were sampled. Shrimp individuals incubated in a sediment-absent environment were also used in this study. Microbial community structure of the shrimp intestine and sediment was investigated through amplicon sequencing targeting bacterial 16S rRNA genes, archaeal 16S rRNA genes, and fungal ITS genes. The results demonstrate distinct differences in community structure between shrimp intestine and the surrounding sediment and between surface and deep (5 mbsf) sediment. The composition of the intestinal microbiota in shrimp living without sediment was different from that of wild-type shrimp, indicating that the presence or absence of sediment could influence the shrimp intestinal microbiota. Carbohydrate metabolism, energy metabolism (carbon fixation, methane metabolism, nitrogen metabolism, and sulfur metabolism), amino acid metabolism, and xenobiotic biodegradation were the most commonly predicted microbial functionalities and they interacted closely with one another. Overall, this study provided comprehensive insights into bacterial, archaeal, and fungal community structure of deep-sea shrimp intestine as well as potential ecological interactions with the surrounding sediment. This study could update our understanding of the microbiota characteristics in shrimp and sediment in deep-sea ecosystems.

The ordering and orientation of genomic scaffolds to reconstruct chromosomes is an essential step during de novo genome assembly. Because this process utilizes various mapping techniques that each provides an independent line of evidence, a combination of multiple maps can improve the accuracy of the resulting chromosomal assemblies. We present ALLMAPS, a method capable of computing a scaffold ordering that maximizes colinearity across a collection of maps. ALLMAPS is robust against common mapping errors, and generates sequences that are maximally concordant with the input maps. ALLMAPS is a useful tool in building high-quality genome assemblies. ALLMAPS is available at: https://github.com/tanghaibao/jcvi/wiki/ALLMAPS .

Increasing utilization of stabilized iron sulfides (FeS) nanoparticles implies an elevated release of the materials into the environment. To understand potential impacts and underlying mechanisms of nanoparticle-induced stress, we used the transcriptome sequencing (RNA-seq) technique to characterize the transcriptomes from adult zebrafish exposed to 10 mg/L carboxymethyl cellulose (CMC) stabilized FeS nanoparticles for 96 h, demonstrating striking differences in the gene expression profiles in liver. The exposure caused significant expression alterations in genes related to immune and inflammatory responses, detoxification, oxidative stress and DNA damage/repair. The complement and coagulation cascades Kyoto encyclopedia of genes and genomes (KEGG) pathway was found significantly up-regulated under nanoparticle exposure. The quantitative real-time polymerase chain reaction using twelve genes confirmed the RNA-seq results. We identified several candidate genes commonly regulated in liver, which may serve as gene indicators when exposed to the nanoparticles. Hepatic inflammation was further confirmed by histological observation of pyknotic nuclei, and vacuole formation upon exposure. Tissue accumulation tests showed a 2.2 times higher iron concentration in the fish tissue upon exposure. This study provides preliminary mechanistic insights into potential toxic effects of organic matter stabilized FeS nanoparticles, which will improve our understanding of the genotoxicity caused by stabilized nanoparticles.

As the global market for fisheries and aquaculture products expands, mislabeling of these products has become a growing concern in the food safety arena. Molecular species identification techniques hold the potential for rapid, accurate assessment of proper labeling. Here we developed and evaluated DNA barcodes for use in differentiating United States domestic and imported catfish species. First, we sequenced 651 base-pair barcodes from the cytochrome oxidase I (COI) gene from individuals of 9 species (and an Ictalurid hybrid) of domestic and imported catfish in accordance with standard DNA barcoding protocols. These included domestic Ictalurid catfish, and representative imported species from the families of Clariidae and Pangasiidae. Alignment of individual sequences from within a given species revealed highly consistent barcodes (98% similarity on average). These alignments allowed the development and analyses of consensus barcode sequences for each species and comparison with limited sequences in public databases (GenBank and Barcode of Life Data Systems). Validation tests carried out in blinded studies and with commercially purchased catfish samples (both frozen and fresh) revealed the reliability of DNA barcoding for differentiating between these catfish species. The developed protocols and consensus barcodes are valuable resources as increasing market and governmental scrutiny is placed on catfish and other fisheries and aquaculture products labeling in the United States.

The lectin pathway of the complement system is characterized by two groups of soluble pattern recognition molecules, mannose-binding lectins (MBLs) and ficolins. These molecules recognize and bind carbohydrates in pathogens and activate complement leading to opsonization, leukocyte activation, and direct pathogen killing. While MBLs have been reported in many fish species, ficolins do not appear to be present in the teleost lineage, despite their importance in invertebrate and higher vertebrate innate immunity. A protein with a similar fibrinogen-like domain, microfibrillar-associated protein 4, MFAP4, is present in fish, albeit with no described immune function. We examined whether MFAP4 genes in fish may potentially act as pathogen receptors in the absence of ficolin. We isolated and characterized five MFAP4 genes from channel catfish. Linkage mapping and phylogenetic analysis indicated that at least three of the catfish MFAP4 genes are tightly clustered on a single chromosome, suggesting that they may have arisen through tandem duplication. Divergent, duplicated families of MFAP4 genes are also present in other teleost species. Expression analysis of the catfish MFAP4 transcripts revealed unique patterns of homeostatic expression among the genes in gill, spleen, skin, liver, and muscle. Expression of the five MFAP4 transcripts showed significant changes in expression as soon as 4h after infection with either Edwardsiella ictaluri or Flavobacterium columnare with modulation of expression continuing up to 7 d following pathogen exposure. Several different tissues and gene-specific patterns were captured and transcript expression changes of >30-fold were observed over the course of the bacterial challenges. Our results suggest a novel role for MFAP4 in teleost immune responses.

Upon the completion of whole genome sequencing, thorough genome annotation that associates genome sequences with biological meanings is essential. Genome annotation depends on the availability of transcript information as well as orthology information. In teleost fish, genome annotation is seriously hindered by genome duplication. Because of gene duplications, one cannot establish orthologies simply by homology comparisons. Rather intense phylogenetic analysis or structural analysis of orthologies is required for the identification of genes. To conduct phylogenetic analysis and orthology analysis, full-length transcripts are essential. Generation of large numbers of full-length transcripts using traditional transcript sequencing is very difficult and extremely costly.

Recent advances in next-generation sequencing technologies have drastically increased throughput and significantly reduced sequencing costs. However, the average read lengths in next-generation sequencing technologies are short as compared with that of traditional Sanger sequencing. The short sequence reads pose great challenges for de novo sequence assembly. As a pilot project for whole genome sequencing of the catfish genome, here we attempt to determine the proper sequence coverage, the proper software for assembly, and various parameters used for the assembly of a BAC physical map contig spanning approximately a million of base pairs.

Sampling challenges in deep-sea ecosystems lead to a lack of knowledge about the distribution of microbes in different submarine canyons. To study microbial diversity and community turnover under different ecological processes, we performed 16S/18S rRNA gene amplicon sequencing for sediment samples from a submarine canyon in the South China Sea. Bacteria, archaea, and eukaryotes made up 57.94% (62 phyla), 41.04% (12 phyla), and 1.02% (4 phyla) of the sequences, respectively. Thaumarchaeota, Planctomycetota, Proteobacteria, Nanoarchaeota, and Patescibacteria are the five most abundant phyla. Heterogeneous community composition was mainly observed in vertical profiles rather than horizontal geographic locations, and microbial diversity in the surface layer was much lower than that in deep layers. According to the null model tests, homogeneous selection dominated community assembly within each sediment layer, whereas heterogeneous selection and dispersal limitation dominated community assembly between distant layers. Different sedimentation processes of sediments, i.e., rapid deposition caused by turbidity currents or slow sedimentation, seem to be primarily responsible for these vertical variations. Finally, functional annotation through shotgun-metagenomic sequencing found that glycosyl transferases and glycoside hydrolases are the most abundant carbohydrate-active enzyme categories. The most likely expressed sulfur cycling pathways include assimilatory sulfate reduction, the link between inorganic and organic sulfur transformation, and organic sulfur transformation, while the potentially activated methane cycling pathways include aceticlastic methanogenesis and aerobic and anaerobic oxidation of methane. Overall, our study revealed high levels of microbial diversity and putative functions in canyon sediments and the important influence of sedimentary geology on microbial community turnover between vertical sediment layers. IMPORTANCE Deep-sea microbes have received growing attention due to their contribution to biogeochemical cycles and climate change. However, related research lags due to the difficulty of collecting samples. Based on our previous study, which revealed the formation of sediments under the dual action of turbidity currents and seafloor obstacles in a submarine canyon in the South China Sea, this interdisciplinary research provides new insights into how sedimentary geology influences microbial community assembly in sediments. We proposed some uncommon or new findings, including the following: (i) microbial diversity was much lower on the surface than in deeper layers (ii) archaea and bacteria dominated the surface and deep layers, respectively; (iii) sedimentary geology played key roles in vertical community turnover; and (iv) the microbes have great potential to catalyze sulfur, carbon, and methane cycling. This study may lead to extensive discussion of the assembly and function of deep-sea microbial communities in the context of geology.

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

In China, laparoscopic splenectomy and esophagogastric devascularization (LSED) are effective and safe tools that are used to treat esophageal-fundic variceal bleeding with portal hypertension (PHT) with minimal trauma; however, due to the increased difficulty of operation, their application in massive splenomegaly (MS) remains limited. This study sought to determine the efficacy and safety of LSED in treating MS patients with PHT.

In this study, we performed a comprehensive analysis of the transcriptome of one- and two-year-old male and female brains of Cynoglossus semilaevis by high-throughput Illumina sequencing. A total of 77,066 transcripts, corresponding to 21,475 unigenes, were obtained with a N50 value of 4349 bp. Of these unigenes, 33 genes were found to have significant differential expression and potentially associated with growth, from which 18 genes were down-regulated and 12 genes were up-regulated in two-year-old males, most of these genes had no significant differences in expression among one-year-old males and females and two-year-old females. A similar analysis was conducted to look for genes associated with reproduction; 25 genes were identified, among them, five genes were found to be down regulated and 20 genes up regulated in two-year-old males, again, most of the genes had no significant expression differences among the other three. The performance of up regulated genes in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was significantly different between two-year-old males and females. Males had a high gene expression in genetic information processing, while female's highly expressed genes were mainly enriched on organismal systems. Our work identified a set of sex-biased genes potentially associated with growth and reproduction that might be the candidate factors affecting sexual dimorphism of tongue sole, laying the foundation to understand the complex process of sex determination of this economic valuable species.

Teleost fish of the genus Danio are excellent models to study the genetic and cellular bases of pigment pattern variation in vertebrates. The two sister species Danio rerio and Danio aesculapii show divergent patterns of horizontal stripes and vertical bars that are partly caused by the divergence of the potassium channel gene kcnj13. Here, we show that kcnj13 is required only in melanophores for interactions with xanthophores and iridophores, which cause location-specific pigment cell shapes and thereby influence colour pattern and contrast in D. rerio. Cis-regulatory rather than protein coding changes underlie kcnj13 divergence between the two Danio species. Our results suggest that homotypic and heterotypic interactions between the pigment cells and their shapes diverged between species by quantitative changes in kcnj13 expression during pigment pattern diversification.

Yellow catfish (Pelteobagrus fulvidraco) is a pivotal freshwater aquaculture species in China. It shows sexual size dimorphism favoring male in growth. Whole transcriptome approach is required to get the overview of genetic toolkit for understanding the sex determination mechanism aiming at devising its monosex production. Beside gonads, the brain is also considered as a major organ for vertebrate reproduction. Transcriptomic analyses on the brain and of different developmental stages will provide the dynamic view necessary for better understanding its sex determination. In this regard, we have performed a de novo assembly of yellow catfish brain transcriptome by high throughput Illumina sequencing. A total number of 154,507 contigs were obtained with the lengths ranging from 201 to 27,822 bp and N50 of 2,101 bp, as well as 20,699 unigenes were identified. Of these unigenes, 13 and 54 unigenes were detected to be XY-specifically expressed genes (SEGs) for one and 2-year-old yellow catfish, while the corresponding numbers of XX-SEGs for those two stages were 19 and 13, respectively. Our work identifies a set of annotated genes that are candidate factors affecting sexual dimorphism as well as simple sequence repeat (SSR) and single nucleotide variation (SNV) in yellow catfish. To validate the expression patterns of the sex-related genes, we performed quantitative real-time PCR (qRT-PCR) indicating the reliability and accuracy of our analysis. The results in our study may enhance our understanding of yellow catfish sex determination and potentially help to improve the production of all-male yellow catfish for aquaculture.

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Forkhead box P3 (FOXP3), an X-linked tumor suppressor gene, plays an important role in breast cancer. However, the biological functions of FOXP3 in breast cancer apoptosis remain unclear. To investigate the underlying genes and networks regulated by FOXP3 in breast cancer, RNA sequencing was performed to compare FOXP3-overexpressing MDA-MB-231 cells and control MDA-MB-231 cells. Differentially expressed genes were identified, and functional enrichment analysis comparing the two groups was performed. The differentially expressed genes were mainly enriched in phagosomes, oxytocin, serotonergic synapses and the phospholipase D signaling pathway. Furthermore, gene set enrichment analysis revealed the enrichment of a gene signature associated with apoptosis in FOXP3-overexpressing MDA-MB-231 cells compared with wild-type cells. Further analysis showed that programmed cell death 4 (PDCD4), a key molecule involved in apoptosis, was overexpressed in FOXP3-MDA-MB-231 cells. Reverse transcription-quantitative PCR and western blotting showed that FOXP3 upregulated the expression of PDCD4 in breast cancer cells. Clinical sample analysis using a public database showed that the expression level of PDCD4 was associated with breast cancer clinical stages. Overall, the present study suggested that FOXP3 can promote the apoptosis of breast cancer cells by upregulating the expression of PDCD4, thus exerting a tumor suppressive function.

Glypican-3(GPC3) is a transmembrane protein which has been found to be frequently overexpressed on the surfaces of liver cancer (LC) cells, which contributes to both the growth and metastasis of LC cells. Recently, the expression of GPC3 has been reported to be inversely associated with glucose metabolism activity in LC patients, suggesting that GPC3 may play a role in the regulation of glucose metabolism in LC. However, the role of GPC3 in glucose metabolism reprogramming, as well as in LC cell growth and metastasis, is unknown. Here, we found that GPC3 significantly contributed to the reprogramming of glucose metabolism in LC cells. On the one hand, GPC3 enhanced the glycolysis of LC cells through upregulation of the glycolytic genes of Glut1, HK2, and LDH-A. On the other hand, GPC3 repressed mitochondrial respiration through downregulation of peroxisome proliferator-activated receptor-gamma coactivator 1-alpha (PGC-1α), which has been well known as a crucial regulator in mitochondrial biogenesis. Mechanistic investigations revealed that HIF-1α was involved in both GPC3-regulated upregulation of glycolytic genes of HK2, PKM2, and Glut1 and downregulation of mitochondrial biogenesis regulator PGC-1α in LC cells. Additionally, GPC3-regulated reprogramming of glucose metabolism played a critical role in the growth and metastasis of LC cells. Conclusion. Our findings demonstrate that GPC3 is a critical regulator of glucose metabolism reprogramming in LC cells, which provides a strong line of evidence for GPC3 as an important therapeutic target to normalize glucose metabolic aberrations responsible for LC progression.