New genes helped acorn barnacles adapt to a sessile lifestyle.

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

Rapid discovery of a new antifoulant: From in silico studies targeting barnacle chitin synthase to efficacy against barnacle settlement.

Due to the adverse environmental impacts of toxic heavy metal-based antifoulants, the screening of environmentally friendly antifoulants has become important for the development of marine antifouling technology. Compared with the traditional lengthy and costly screening method, computer-aided drug design (CADD) offers a promising and efficient solution that can accelerate the screening process of green antifoulants. In this study, we selected barnacle chitin synthase (CHS, an important enzyme for barnacle settlement and development) as the target protein for docking screening. Three CHS genes were identified in the barnacle Amphibalanus amphitrite, and their encoded proteins were found to share a conserved glycosyltransferase domain. Molecular docking of 31,561 marine natural products with AaCHSs revealed that zoanthamine alkaloids had the best binding affinity (-11.8 to -12.6 kcal/mol) to AaCHSs. Considering that the low abundance of zoanthamine alkaloids in marine organisms would limit their application as antifoulants, a marine fungal-derived natural product, mycoepoxydiene (MED), which has a similar chemical structure to zoanthamine alkaloids and the potential for large-scale production by fermentation, was selected and validated for stable binding to AaCHS2L2 using molecular docking and molecular dynamics simulations. Finally, the efficacy of MED in inhibiting cyprid settlement of A. amphitrite was confirmed by a bioassay that demonstrated an EC50 of 1.97 µg/mL, suggesting its potential as an antifoulant candidate. Our research confirmed the reliability of using AaCHSs as antifouling targets and has provided insights for the efficient discovery of green antifoulants by CADD.

OmicsSuite: a customized and pipelined suite for analysis and visualization of multi-omics big data.

With the advancements in high-throughput sequencing technologies such as Illumina, PacBio, and 10X Genomics platforms, and gas/liquid chromatography-mass spectrometry, large volumes of biological data in multiple formats can now be obtained through multi-omics analysis. Bioinformatics is constantly evolving and seeking breakthroughs to solve multi-omics problems; however, it is challenging for most experimental biologists to analyse data using command-line interfaces, coding, and scripting. Based on experience with multi-omics, we have developed OmicsSuite, a desktop suite that comprehensively integrates statistics and multi-omics analysis and visualization. The suite has 175 sub-applications in 12 categories, including Sequence, Statistics, Algorithm, Genomics, Transcriptomics, Enrichment, Proteomics, Metabolomics, Clinical, Microorganism, Single Cell, and Table Operation. We created the user interface with Sequence View, Table View, and intelligent components based on JavaFX and the popular Shiny framework. The multi-omics analysis functions were developed based on BioJava and 300+ packages provided by the R CRAN and Bioconductor communities, and it encompasses over 3000 adjustable parameter interfaces. OmicsSuite can directly read multi-omics raw data in FastA, FastQ, Mutation Annotation Format, mzML, Matrix, and HDF5 formats, and the programs emphasize data transfer directions and pipeline analysis functions. OmicsSuite can produce pre-publication images and tables, allowing users to focus on biological aspects. OmicsSuite offers multi-omics step-by-step workflows that can be easily applied to horticultural plant breeding and molecular mechanism studies in plants. It enables researchers to freely explore the molecular information contained in multi-omics big data (Source: https://github.com/OmicsSuite/, Website: https://omicssuite.github.io, v1.3.9).

Metabolic and transcriptional responses demonstrating enhanced thermal tolerance in domesticated abalone.

Global warming threatens aquatic systems and organisms. Many studies have focused on the vulnerability and stress responses of aquaculture organisms to future thermal conditions. However, it may be of more practical significance to reveal their acclimation potential and mechanisms. In this study, the physiological, metabolic, and transcriptional responses to long-term temperature acclimation of northern and southern populations of Pacific abalone Haliotis discus hannai, a commercially important gastropod sensitive to environmental changes, were compared. This study conducted two common-garden experiments, including a thermostatic experiment in the lab and an aquaculture experiment on the farm. The abalone population cultured in warmer southern waters was tolerant of ongoing high temperatures, whereas the abalone population originally cultured in cooler northern waters exhibited vulnerability to high temperatures but could enhance its thermal tolerance through the process of natural selection in warmer southern waters. This difference was linked to divergence in the metabolic and transcriptional processes of the two populations. The tolerant population exhibited a greater capacity for carbohydrate and amino acid metabolism regulation and energy redistribution to cope with heat stress. This capacity may have been selected for, and accumulated, over many generations because the tolerant population originated from the intolerant population over two decades ago. This work provides insight into the vulnerability and acclimation potential of abalone to heat stress and discloses the molecular and metabolic traits underlying this phenomenon. Future research on the ability of abalone and other commercial shellfish species to acclimate to global warming should take this potential into account.

A computer-aided diagnosis system using white-light endoscopy for the prediction of conventional adenoma with high grade dysplasia.

OBJECTIVES: We developed a computer-aided diagnosis system called ECRCCAD using standard white-light endoscopy (WLE) for predicting conventional adenomas with high-grade dysplasia (HGD) to optimise the patients' management decisions during colonoscopy. METHODS: Pretraining model was used to fine-tune the model parameters by transfer learning. 2,397 images of HGD and 2,487 low-grade dysplasia (LGD) images were randomly assigned (8:1:1) to the training, optimising, and internal validation dataset. The prospective validation dataset is the frames accessed from colonoscope videoes. One independent rural hospital provided an external validation dataset. Histopathological diagnosis was used as the standard criterion. The capability of the ECRCCAD to distinguish HGD was assessed and compared with two expert endoscopists. RESULTS: The accuracy, sensitivity and specificity for diagnosis of HGD in the internal validation set were 90.5%, 93.2%, 87.9%, respectively. While 88.2%, 85.4%, 89.8%, respectively, for the external validation set. For the prospective validation set, ECRCCAD achieved an AUC of 93.5% in diagnosing HGD. The performance of ECRCCAD in diagnosing HGD was better than that of the expert endoscopist in the external validation set (88.2% vs. 71.5%, P < 0.0001). CONCLUSION: ECRCCAD had good diagnostic capability for HGD and enabled a more convenient and accurate diagnosis using WLE.

Genomic insights into the adaptation and evolution of the nautilus, an ancient but evolving "living fossil".

The nautilus, commonly known as a "living fossil," is endangered and may be at risk of extinction. The lack of genomic information hinders a thorough understanding of its biology and evolution, which can shed light on the conservation of this endangered species. Here, we report the first high-quality chromosome-level genome assembly of Nautilus pompilius. The assembled genome size comprised 785.15 Mb. Comparative genomic analyses indicated that transposable elements (TEs) and large-scale genome reorganizations may have driven lineage-specific evolution in the cephalopods. Remarkably, evolving conserved genes and recent TE insertion activities were identified in N. pompilius, and we speculate that these findings reflect the strong adaptability and long-term survival of the nautilus. We also identified gene families that are potentially responsible for specific adaptation and evolution events. Our study provides unprecedented insights into the specialized biology and evolution of N. pompilius, and the results serve as an important resource for future conservation genomics of the nautilus and closely related species.

Evaluation of Genomic Selection for Seven Economic Traits in Yellow Drum (Nibea albiflora).

Yellow drum (Nibea albiflora) is an important maricultural fish in China, and genetic improvement is necessary for this species. This research evaluated the application of genomic selection methods to predict the genetic values of seven economic traits for yellow drum. Using genome-wide single-nucleotide polymorphisms (SNPs), we estimated the genetic parameters for seven traits, including body length (BL), swimming bladder index (SBI), swimming bladder weight (SBW), body thickness (BT), body height (BH), body length/body height ratio (LHR), and gonad weight index (GWI). The heritability estimates ranged from 0.309 to 0.843. We evaluated the prediction performance of various statistical methods, and no one method provided the highest predictive ability for all traits. We then evaluated and compared the use of genome-wide association study (GWAS)-informative SNPs and random SNPs for prediction and found that GWAS-informative SNPs obviously increased. It only needed 5 and 100 informative SNPs for LHR and BT to achieve almost the same predictive abilities as using genome-wide SNPs, and for BL, SBI, SBW, BH, and GWI, about 1000 to 3000 informative SNPs were needed to achieve whole-genome level predictive abilities. It can be concluded from the test results that breeders can use fewer SNPs to save the breeding costs of genomic selection for some traits.

Molecular characterization of carbapenem-resistant Klebsiella pneumoniae isolates with focus on antimicrobial resistance.

BACKGROUND: The enhancing incidence of carbapenem-resistant Klebsiella pneumoniae (CRKP)-mediated infections in Mengchao Hepatobiliary Hospital of Fujian Medical University in 2017 is the motivation behind this investigation to study gene phenotypes and resistance-associated genes of emergence regarding the CRKP strains. In current study, seven inpatients are enrolled in the hospital with complete treatments. The carbapenem-resistant K. pneumoniae whole genome is sequenced using MiSeq short-read and Oxford Nanopore long-read sequencing technology. Prophages are identified to assess genetic diversity within CRKP genomes. RESULTS: The investigation encompassed eight CRKP strains that collected from the patients enrolled as well as the environment, which illustrate that blaKPC-2 is responsible for phenotypic resistance in six CRKP strains that K. pneumoniae sequence type (ST11) is informed. The plasmid with IncR, ColRNAI and pMLST type with IncF[F33:A-:B-] co-exist in all ST11 with KPC-2-producing CRKP strains. Along with carbapenemases, all K. pneumoniae strains harbor two or three extended spectrum ß-lactamase (ESBL)-producing genes. fosA gene is detected amongst all the CRKP strains. The single nucleotide polymorphisms (SNP) markers are indicated and validated among all CRKP strains, providing valuable clues for distinguishing carbapenem-resistant strains from conventional K. pneumoniae. CONCLUSIONS: ST11 is the main CRKP type, and blaKPC-2 is the dominant carbapenemase gene harbored by clinical CRKP isolates from current investigations. The SNP markers detected would be helpful for characterizing CRKP strain from general K. pneumoniae. The data provides insights into effective strategy developments for controlling CRKP and nosocomial infection reductions.

Chromosome assembly of Collichthys lucidus, a fish of Sciaenidae with a multiple sex chromosome system.

Collichthys lucidus (C. lucidus) is a commercially important marine fish species distributed in coastal regions of East Asia with the X1X1X2X2/X1X2Y multiple sex chromosome system. The karyotype for female C. lucidus is 2n = 48, while 2n = 47 for male ones. Therefore, C. lucidus is also an excellent model to investigate teleost sex-determination and sex chromosome evolution. We reported the first chromosome genome assembly of C. lucidus using Illumina short-read, PacBio long-read sequencing and Hi-C technology. An 877 Mb genome was obtained with a contig and scaffold N50 of 1.1 Mb and 35.9 Mb, respectively. More than 97% BUSCOs genes were identified in the C. lucidus genome and 28,602 genes were annotated. We identified potential sex-determination genes along chromosomes and found that the chromosome 1 might be involved in the formation of Y specific metacentric chromosome. The first C. lucidus chromosome-level reference genome lays a solid foundation for the following population genetics study, functional gene mapping of important economic traits, sex-determination and sex chromosome evolution studies for Sciaenidae and teleosts.

Long-read nanopore sequencing-based draft genome of a carbapenem-resistant Pseudomonas aeruginosa isolate.

OBJECTIVES: Pseudomonas aeruginosa is a common Gram-negative bacterium causing various serious infections, such as lower respiratory tract infection and urinary tract infection in catheterised patients. Here we report the draft genome sequence of a carbapenem-resistant P. aeruginosa (CRPA) isolate. METHODS: The genome of the CRPA isolate was sequenced using a combination of short, highly accurate Illumina reads and additional coverage in very long Oxford Nanopore reads. RESULTS: The resulting assembly was highly contiguous, containing a total of 6624003bp with a GC content of 66.21%. Annotation identified 6389 protein-coding genes. Mutations in the oprD and mexR genes conferred resistance to carbapenems in the CRPA isolate. CONCLUSION: The draft genome sequence of this CRPA isolate could provide a solid basis for further research on the resistance mechanisms and the development of drug therapy for drug resistance genes.

Near-complete genome assembly and annotation of the yellow drum (Nibea albiflora) provide insights into population and evolutionary characteristics of this species.

Yellow drum (Nibea albiflora) is an important fish species in capture fishery and aquaculture in East Asia. We herein report the first and near-complete genome assembly of an ultra-homologous gynogenic female yellow drum using Illumina short sequencing reads. In summary, a total of 154.2 Gb of raw reads were generated via whole-genome sequencing and were assembled to 565.3 Mb genome with a contig N50 size of 50.3 kb and scaffold N50 size of 2.2 Mb (BUSCO completeness of 97.7%), accounting for 97.3%-98.6% of the estimated genome size of this fish. We further identified 22,448 genes using combined methods of ab initio prediction, RNAseq annotation, and protein homology searching, of which 21,614 (96.3%) were functionally annotated in NCBI nr, trEMBL, SwissProt, and KOG databases. We also investigated the nucleotide diversity (around 1/390) of aquacultured individuals and found the genetic diversity of the aquacultured population decreased due to inbreeding. Evolutionary analyses illustrated significantly expanded and extracted gene families, such as myosin and sodium: neurotransmitter symporter (SNF), could help explain swimming motility of yellow drum. The presented genome will be an important resource for future studies on population genetics, conservation, understanding of evolutionary history and genetic breeding of the yellow drum and other Nibea species.

A high-density genetic linkage map and QTL mapping for growth and sex of yellow drum (Nibea albiflora).

A high-density genetic linkage map is essential for the studies of comparative genomics and gene mapping, and can facilitate assembly of reference genome. Herein, we constructed a high-density genetic linkage map with 8,094 SNPs selected from 113 sequenced fish of a F1 family. Ultimately, the consensus map spanned 3818.24 cM and covered nearly the whole genome (99.4%) with a resolution of 0.47 cM. 1,457 scaffolds spanning 435.15 Mb were anchored onto 24 linkage groups, accounting for 80.7% of the draft genome assembly of the yellow drum. Comparative genomic analyses with medaka and zebrafish genomes showed superb chromosome-scale synteny between yellow drum and medaka. QTL mapping and association analysis congruously revealed 22 QTLs for growth-related traits and 13 QTLs for sex dimorphism. Some important candidate genes such as PLA2G4A, BRINP3 and P2RY1 were identified from these growth-related QTL regions. A gene family including DMRT1, DMRT2 and DMRT3 was identified from these sex-related QTL regions on the linkage group LG9. We demonstrate that this linkage map can facilitate the ongoing marker-assisted selection and genomic and genetic studies for yellow drum.

Genomewide association study for economic traits in the large yellow croaker with different numbers of extreme phenotypes.

A traditional genomewide association study (GWAS) detects genotype-phenotype associations by the vast number of genotyped individuals. This method requires large-scale samples and considerable sequencing costs. Extreme phenotypic sampling proposes make GWAS more cost-efficient and are applied more widely. With extreme phenotypic sampling, we performed a GWAS for n-3 highly unsaturated fatty acids (HUFA) and eviscerated weight (EW) traits in the large yellowcroaker population. Of the 32,249 and 29,748 detected SNPs for the two traits, three candidate regions were found in each trait. Three candidate regions associated with HUFA were known near genes on chromosomes 4 and 11, and three candidate regions were on chromosome 6, and 15 for the EW trait. By combing through our GWAS results and the biological functional analysis of the genes, we suggest that the FABP, DGAT, ATP8B1, FAF2 and CERS2 genes, as well as the IGF2, BORA, CYP1A1, GRTP1 and HOX genes are promising candidate genes for n-3 HUFA and EW, respectively, in the large yellow croaker.Moreover, compared with the different numbers of the extreme phenotypic sampling, we conclude that 60% of the extreme phenotypic subsample can obtain a similar result as GWAS with whole phenotypes. Thus, extreme phenotypic sampling could save 40% of the cost for genotyping and DNA extraction without loss of the candidate regions and functional genes. Our study may provide a basis for further genomic breeding and a reference for others who want to perform GWAS with extreme phenotypes.

The identification of growth, immune related genes and marker discovery through transcriptome in the yellow drum (Nibea albiflora).

Yellow drum (Nibea albiflora) is a commercially important marine fish, which is widely distributed in the coastal waters of China, Japan and Korea. Wild yellow drum resources have dramatically declined due to overfishing and ocean pollution. Genetic data can contribute to biodiversity conservation and protection. And molecular markers can play important roles in genetic breeding and aid in germplasm preservation in fish. In this study, 11 tissues (brain, heart, liver, kidney, muscle, head kidney, skin, fin, spleen, gonad and air bladder) were collected for pooled RNA sequencing. The unigenes were assembled using Trinity and EvidentialGene, and were then aligned to nr, nt, Swiss-Prot GO, KEGG, and KOG for annotation. Molecular markers (e.g. simple sequence repeat, SSR and single nucleotide polymorphism, SNP) were detected using MIcroSAtellite identification tool (MISA) and Genome Analysis Tool Kit (GATK). All clean reads were assembled into 109,209 transcripts, and 31,183 unigenes were generated after pruning and classifying, ranging from 201 to 19,857 bp in length (1230 bp in average), and 26,728 (85.7%) assembled unigenes had significant hits in public databases. Total of 27 and 103 unigenes were respectively identified as involved in growth- and immune-related pathways in the N. albiflora transcriptome. In addition, we identified a considerable quantity of molecular markers, including 11,484 SSRs and 56,186 SNPs. The growth- and immune-relevant genes and the molecular markers identified here provided a meaningful reference gene set and laid a foundation for future genetic selection and breeding for this species.

Whole-genome single-nucleotide polymorphism (SNP) marker discovery and association analysis with the eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) content in Larimichthys crocea.

Whole-genome single-nucleotide polymorphism (SNP) markers are valuable genetic resources for the association and conservation studies. Genome-wide SNP development in many teleost species are still challenging because of the genome complexity and the cost of re-sequencing. Genotyping-By-Sequencing (GBS) provided an efficient reduced representative method to squeeze cost for SNP detection; however, most of recent GBS applications were reported on plant organisms. In this work, we used an EcoRI-NlaIII based GBS protocol to teleost large yellow croaker, an important commercial fish in China and East-Asia, and reported the first whole-genome SNP development for the species. 69,845 high quality SNP markers that evenly distributed along genome were detected in at least 80% of 500 individuals. Nearly 95% randomly selected genotypes were successfully validated by Sequenom MassARRAY assay. The association studies with the muscle eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) content discovered 39 significant SNP markers, contributing as high up to ∼63% genetic variance that explained by all markers. Functional genes that involved in fat digestion and absorption pathway were identified, such as APOB, CRAT and OSBPL10. Notably, PPT2 Gene, previously identified in the association study of the plasma n-3 and n-6 polyunsaturated fatty acid level in human, was re-discovered in large yellow croaker. Our study verified that EcoRI-NlaIII based GBS could produce quality SNP markers in a cost-efficient manner in teleost genome. The developed SNP markers and the EPA and DHA associated SNP loci provided invaluable resources for the population structure, conservation genetics and genomic selection of large yellow croaker and other fish organisms.

Functional marker detection and analysis on a comprehensive transcriptome of large yellow croaker by next generation sequencing.

Large yellow croaker (Larimichthys crocea) is an important economic fish in China and Eastern Asia. Because of the exhaustive fishing and overdense aquaculture, the wild population and the mariculture of the species are facing serious challenges on germplasm degeneration and susceptibility to infectious disease agents. However, a comprehensive transcriptome from multi-tissues of the species has not been reported and functional molecular markers have not yet been detected and analyzed. In this work, we applied RNA-seq with the Illumina Hiseq2000 platform for a multi-tissue sample of large yellow croaker and assembled the transcriptome into 88,103 transcripts. Of them, 52,782 transcripts have been successfully annotated by nt/nr, InterPro, GO and KEGG database. Comparing with public fish proteins, we have found that 34,576 protein coding transcripts are shared in large yellow croaker with zebrafish, medaka, pufferfish, and stickleback. For functional markers, we have discovered 1,276 polymorphic SSRs and 261, 000 SNPs. The functional impact analysis of SNPs showed that the majority (~75%) of small variants cause synonymous mutations in proteins, followed by variations in 3' UTR region. The functional enrichment analysis illuminated that transcripts involved in DNA bindings, enzyme activities, and signal pathways prominently exhibit less single-nucleotide variants but genes for the constituent of the muscular tissue, the cytoskeleton, and the immunity system contain more frequent SNP mutations, which may reflect the structural and functional selections of the translated proteins. This is the first work for the high-throughput detection and analysis of functional polymorphic SSR and SNP markers in a comprehensive transcriptome of large yellow croaker. Our study provides valuable transcript sequence and functional marker resources for the quantitative trait locus (QTL) identification and molecular selection of the species in the research community.