Association Between Muscle Growth and Transcription of a Mutant MSTN Gene in Olive Flounder (Paralichthys olivaceus).

Myostatin (MSTN, also known as growth differentiation factor-8 (GDF-8)), a member of the transforming growth factor ß (TGF-ß) superfamily, functions as a negative regulator of skeletal muscle development and growth. However, it is also expressed in a wide range of tissues in fish and thus may have more diverse roles in this group than in mammals. In this study, we assessed the genome-wide transcriptional expression pattern associated with the CRISPR/Cas9-mutated MSTN gene in the olive flounder (Paralichthys olivaceus) in association with changes in cell proliferation and transportation processes. There were no differences in the hepatosomatic index, and the growth of male and female fish increased in the F1 progeny of the MSTN mutants. Furthermore, the histopathological analysis showed that myostatin editing resulted in a 41.24% increase in back muscle growth and 46.92% increase in belly muscle growth in male flounder compared with normal flounder, and a 16.01% increase in back muscle growth and 14.26% increase in belly muscle growth in female flounder compared with normal flounder. This study demonstrates that editing of the myostatin gene enhances muscle growth in olive flounder, with a notably more pronounced effect observed in males. Consequently, myostatin-edited male flounder could represent a valuable asset for the flounder aquaculture industry.

Genome-Wide Differential Methylation Profiles from Two Terpene-Rich Medicinal Plant Extracts Administered in Osteoarthritis Rats.

Extracts from the plants Phlomis umbrosa and Dipsacus asperoides-which are widely used in Korean and Chinese traditional medicine to treat osteoarthritis and other bone diseases-were used to treat experimental osteoarthritis (OA) rats. Genome-wide differential methylation regions (DMRs) of these medicinal-plant-treated rats were profiled as therapeutic evidence associated with traditional medicine, and they need to be investigated further using detailed molecular research to extrapolate traditional practices to modern medicine. In total, 49 protein-encoding genes whose expression is differentially regulated during disease progression and recovery have been discovered via systematic bioinformatic analysis and have been approved/proposed as druggable targets for various bone diseases by the US food and drug administration. Genes encoding proteins involved in the PI3K/AKT pathway were found to be enriched, likely as this pathway plays a crucial role during OA progression as well as during the recovery process after treatment with the aforementioned plant extracts. The four sub-networks of PI3K/AKT were highly regulated by these plant extracts. Overall, 29 genes were seen in level 2 (51-75%) DMRs and were correlated highly with OA pathogenesis. Here, we propose that these genes could serve as targets to study OA; moreover, the iridoid and triterpenoid phytochemicals obtained from these two plants may serve as potential therapeutic agents.

Draft genome of the Korean smelt Hypomesus nipponensis and its transcriptomic responses to heat stress in the liver and muscle.

Pond smelt (Hypomesus nipponensis) is a cold-freshwater fish species and a winter economic aquaculture resource in South Korea. Because of its high susceptibility to abnormal water temperature from global warming, a large number of smelt die in hot summers. Here, we present the first draft genome of H. nipponensis and transcriptomic changes in molecular mechanisms or intracellular responses under heat stress. We combined Illumina and PacBio sequencing technologies to generate the draft genome of H. nipponensis. Based on the reference genome, we conducted transcriptome analysis of liver and muscle tissues under normal (NT, 5°C) vs. warm (HT, 23°C) conditions to identify heat stress-induced genes and gene categories. We observed a total of 1987 contigs with N50 of 0.46 Mbp, with the largest contig (3.03 Mbp) in the assembled genome. A total of 20,644 protein-coding genes were predicted, and 19,224 genes were functionally annotated: 15,955 genes for Gene Ontology terms and 11,560 genes for KEGG Orthology. We conducted the lost and gained genes analysis compared with three species that: human, zebrafish, and salmon. In the lost genes analysis, we detected that smelt lost 4461 (22.16%), 2825 (10.62%), and 1499 (3.09%) genes compare with above three species, respectively. In the gained genes analysis, we observed that smelt gained 1133 (5.49%), 1670 (8.09%), and 229 (1.11%) genes compared with the above species, respectively. From transcriptome analysis, a total of 297 and 331 differentially expressed genes (DEGs) with a false discovery rate <0.05 were identified in the liver and muscle tissues, respectively. Gene enrichment analysis of DEGs indicates that upregulated genes were significantly enriched for lipid biosynthetic process (GO:0008610, P < 0.001) and regulation of apoptotic process (GO:0042981, P < 0.01), and genes were downregulated by immune responses such as myeloid cell differentiation (GO:0030099, P < 0.001) in the liver under heat stress. In muscle tissue, upregulated genes were enriched for hypoxia (GO:0001666, P < 0.05), transcription regulator activity (GO:0140110, P < 0.001), and calcium-release channel activity (GO:0015278, P < 0.01), and genes were downregulated for a nicotinamide nucleotide biosynthetic process (GO:0019359, P < 0.01). The results of KEGG pathway analysis were similar to that of gene enrichment analysis. The draft genome and transcriptomic of H. nipponensis will be a useful genetic resource for functional and evolutionary studies. Our findings will improve understanding of molecular mechanisms and heat responses and be useful for predicting survival of the smelt and its closely related species under global warming.

Machine learning, transcriptome, and genotyping chip analyses provide insights into SNP markers identifying flower color in Platycodon grandiflorus.

Bellflower is an edible ornamental gardening plant in Asia. For predicting the flower color in bellflower plants, a transcriptome-wide approach based on machine learning, transcriptome, and genotyping chip analyses was used to identify SNP markers. Six machine learning methods were deployed to explore the classification potential of the selected SNPs as features in two datasets, namely training (60 RNA-Seq samples) and validation (480 Fluidigm chip samples). SNP selection was performed in sequential order. Firstly, 96 SNPs were selected from the transcriptome-wide SNPs using the principal compound analysis (PCA). Then, 9 among 96 SNPs were later identified using the Random forest based feature selection method from the Fluidigm chip dataset. Among six machines, the random forest (RF) model produced higher classification performance than the other models. The 9 SNP marker candidates selected for classifying the flower color classification were verified using the genomic DNA PCR with Sanger sequencing. Our results suggest that this methodology could be used for future selection of breeding traits even though the plant accessions are highly heterogeneous.

A genome-wide investigation of the effect of farming and human-mediated introduction on the ubiquitous seaweed Undaria pinnatifida.

Human activity is an important driver of ecological and evolutionary change on our planet. In particular, domestication and biological introductions have important and long-lasting effects on species' genomic architecture and diversity. However, genome-wide analysis of independent domestication and introduction events within a single species has not previously been performed. The Pacific kelp Undaria pinnatifida provides such an opportunity because it has been cultivated in its native range in Northeast Asia but also introduced to four other continents in the past 50 years. Here we present the results of a genome-wide analysis of natural, cultivated and introduced populations of U. pinnatifida to elucidate human-driven evolutionary change. We demonstrate that these three categories of origin can be distinguished at the genome level, reflecting the combined influence of neutral (demography and migration) and non-neutral (selection) processes.

The transcriptome data from the leaves of four Papaver species captured at the plant's three developmental life cycles.

The plants in the Papaver genus are widely known as Poppies, which is used for ornamental and medicinal purposes, to utilize its plants derived alkaloids and attractive flowers. From this genus, we have sequenced the transcriptomes of four species's (Papaver rhoeas (two cultivar), Papaver nudicaule (five cultivar), Papaver fauriei, and Papaver somniferum) leaves at three developmental stages (i.e., leaf rosette (30 days), elongation and branching (60 days), and blossom and seed formations (90 days)), to elucidate the secondary metabolite biosynthesis gene expression profiles at respective plant stages.

Author Correction: Whole genome sequencing reveals the impact of recent artificial selection on red sea bream reared in fish farms.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Expression of Heat Shock Proteins in Thermally Challenged Pacific Abalone Haliotis discus hannai.

Summer mortality, caused by thermal conditions, is the biggest threat to abalone aquaculture production industries. Various measures have been taken to mitigate this issue by adjusting the environment; however, the cellular processes of Pacific abalone (Haliotis discus hannai) have been overlooked due to the paucity of genetic information. The draft genome of H. discus hannai has recently been reported, prompting exploration of the genes responsible for thermal regulation in Pacific abalone. In this study, 413 proteins were systematically annotated as members of the heat shock protein (HSP) super families, and among them 26 HSP genes from four Pacific abalone tissues (hemocytes, gill, mantle, and muscle) were differentially expressed under cold and heat stress conditions. The co-expression network revealed that HSP expression patterns were tissue-specific and similar to those of other shellfish inhabiting intertidal zones. Finally, representative HSPs were selected at random and their expression patterns were identified by RNA sequencing and validated by qRT-PCR to assess expression significance. The HSPs expressed in hemocytes were highly similar in both analyses, suggesting that hemocytes could be more reliable samples for validating thermal condition markers compared to other tissues.

Identification of transcriptome-wide, nut weight-associated SNPs in Castanea crenata.

Nut weight is one of the most important traits that can affect a chestnut grower's returns. Due to the long juvenile phase of chestnut trees, the selection of desired characteristics at early developmental stages represents a major challenge for chestnut breeding. In this study, we identified single nucleotide polymorphisms (SNPs) in transcriptomic regions, which were significantly associated with nut weight in chestnuts (Castanea crenata), using a genome-wide association study (GWAS). RNA-sequencing (RNA-seq) data were generated from large and small nut-bearing trees, using an Illumina HiSeq. 2000 system, and 3,271,142 SNPs were identified. A total of 21 putative SNPs were significantly associated with chestnut weight (false discovery rate [FDR] < 10-5), based on further analyses. We also applied five machine learning (ML) algorithms, support vector machine (SVM), C5.0, k-nearest neighbour (k-NN), partial least squares (PLS), and random forest (RF), using the 21 SNPs to predict the nut weights of a second population. The average accuracy of the ML algorithms for the prediction of chestnut weights was greater than 68%. Taken together, we suggest that these SNPs have the potential to be used during marker-assisted selection to facilitate the breeding of large chestnut-bearing varieties.

Whole genome sequencing reveals the impact of recent artificial selection on red sea bream reared in fish farms.

Red sea bream, a popular fish resource in Korea and Japan, is being bred in fish farms of the two countries. It is hypothesized that the genomes of red sea bream are influenced by decades of artificial selection. This study investigates the impact of artificial selection on genomes of red sea bream. Whole genome sequencing was conducted for 40 samples of red sea bream either from Ehime, Nagasaki and Tongyeong fish farms or from the wild. Population stratification based on whole genome data was investigated and the genomic regions of fish farm populations under selection were identified using XP-EHH and relative nucleotide diversity. Gene ontology analysis revealed that different functions were enriched in different fish farms. In conclusion, this study highlights the difference between independently cultured red sea bream populations by showing that influence of artificial selection acted upon completely different genes related to different functions including metabolic and developmental processes.

First draft genome sequence of the rock bream in the family Oplegnathidae.

The rock bream (Oplegnathus fasciatus) is one of the most economically valuable marine fish in East Asia, and due to various environmental factors, there is substantial revenue loss in the production sector. Therefore, knowledge of its genome is required to uncover the genetic factors and the solutions to these problems. In this study, we constructed the first draft genome of O. fasciatus as a reference for the family Oplegnathidae. The genome size is estimated to be 749 Mb, and it was assembled into 766 Mb by combining Illumina and PacBio sequences. A total of 24,053 transcripts (23,338 genes) are predicted, and among those transcripts, 23,362 (97%), are annotated with functional terms. Finally, the completeness of the genome assembly was assessed by CEGMA, which resulted in the complete mapping of 220 (88.7%) core genes in the genome. To the best of our knowledge, this is the first draft genome for the family Oplegnathidae.

Transcriptome Profiling of Two Ornamental and Medicinal Papaver Herbs.

The Papaver spp. (Papaver rhoeas (Corn poppy) and Papaver nudicaule (Iceland poppy)) genera are ornamental and medicinal plants that are used for the isolation of alkaloid drugs. In this study, we generated 700 Mb of transcriptome sequences with the PacBio platform. They were assembled into 120,926 contigs, and 1185 (82.2%) of the benchmarking universal single-copy orthologs (BUSCO) core genes were completely present in our assembled transcriptome. Furthermore, using 128 Gb of Illumina sequences, the transcript expression was assessed at three stages of Papaver plant development (30, 60, and 90 days), from which we identified 137 differentially expressed transcripts. Furthermore, three co-occurrence heat maps are generated from 51 different plant genomes along with the Papaver transcriptome, i.e., secondary metabolite biosynthesis, isoquinoline alkaloid biosynthesis (BIA) pathway, and cytochrome. Sixty-nine transcripts in the BIA pathway along with 22 different alkaloids (quantified with LC-QTOF-MS/MS) were mapped into the BIA KEGG map (map00950). Finally, we identified 39 full-length cytochrome transcripts and compared them with other genomes. Collectively, this transcriptome data, along with the expression and quantitative metabolite profiles, provides an initial recording of secondary metabolites and their expression related to Papaver plant development. Moreover, these profiles could help to further detail the functional characterization of the various secondary metabolite biosynthesis and Papaver plant development associated problems.

Analysis of the Draft Genome of the Red Seaweed Gracilariopsis chorda Provides Insights into Genome Size Evolution in Rhodophyta.

Red algae (Rhodophyta) underwent two phases of large-scale genome reduction during their early evolution. The red seaweeds did not attain genome sizes or gene inventories typical of other multicellular eukaryotes. We generated a high-quality 92.1 Mb draft genome assembly from the red seaweed Gracilariopsis chorda, including methylation and small (s)RNA data. We analyzed these and other Archaeplastida genomes to address three questions: 1) What is the role of repeats and transposable elements (TEs) in explaining Rhodophyta genome size variation, 2) what is the history of genome duplication and gene family expansion/reduction in these taxa, and 3) is there evidence for TE suppression in red algae? We find that the number of predicted genes in red algae is relatively small (4,803-13,125 genes), particularly when compared with land plants, with no evidence of polyploidization. Genome size variation is primarily explained by TE expansion with the red seaweeds having the largest genomes. Long terminal repeat elements and DNA repeats are the major contributors to genome size growth. About 8.3% of the G. chorda genome undergoes cytosine methylation among gene bodies, promoters, and TEs, and 71.5% of TEs contain methylated-DNA with 57% of these regions associated with sRNAs. These latter results suggest a role for TE-associated sRNAs in RNA-dependent DNA methylation to facilitate silencing. We postulate that the evolution of genome size in red algae is the result of the combined action of TE spread and the concomitant emergence of its epigenetic suppression, together with other important factors such as changes in population size.

Genome sequence of pacific abalone (Haliotis discus hannai): the first draft genome in family Haliotidae.

Background: Abalones are large marine snails in the family Haliotidae and the genus Haliotis belonging to the class Gastropoda of the phylum Mollusca. The family Haliotidae contains only one genus, Haliotis, and this single genus is known to contain several species of abalone. With 18 additional subspecies, the most comprehensive treatment of Haliotidae considers 56 species valid [ 1 ]. Abalone is an economically important fishery and aquaculture animal that is considered a highly prized seafood delicacy. The total global supply of abalone has increased 5-fold since the 1970s and farm production increased explosively from 50 mt to 103 464 mt in the past 40 years. Additionally, researchers have recently focused on abalone given their reported tumor suppression effect. However, despite the valuable features of this marine animal, no genomic information is available for the Haliotidae family and related research is still limited. To construct the H . discus hannai genome, a total of 580-G base pairs using Illumina and Pacbio platforms were generated with 322-fold coverage based on the 1.8-Gb estimated genome size of H . discus hannai using flow cytometry. The final genome assembly consisted of 1.86 Gb with 35 450 scaffolds (>2 kb). GC content level was 40.51%, and the N50 length of assembled scaffolds was 211 kb. We identified 29 449 genes using Evidence Modeler based on the gene information from ab initio prediction, protein homology with known genes, and transcriptome evidence of RNA-seq. Here we present the first Haliotidae genome, H . discus hannai , with sequencing data, assembly, and gene annotation information. This will be helpful for resolving the lack of genomic information in the Haliotidae family as well as providing more opportunities for understanding gastropod evolution.

Identification of a gene encoding a membrane-anchored toll-like receptor 5 (TLR5M) in Oplegnathus fasciatus that responds to flagellin challenge and activates NF-κB.

Toll-like receptor 5 (TLR5) recognizes bacterial flagellin and induces the downstream signaling through the myeloid differentiation primary response gene 88 (MyD88) protein to produce proinflammatory cytokines. In this study, we describe a TLR5 membrane form (OfTLR5M) and its adaptor protein MyD88 (OfMyD88) in rock bream, Oplegnathus fasciatus. Both Oftlr5m (6.7 kb) and Ofmyd88 (3.7 kb) genes displayed a quinquepartite structure with five exons and four introns. Protein structure of OfTLR5M revealed the conventional architecture of TLRs featured by an extracellular domain with 22 leucine rich repeats (LRR), a transmembrane domain and an endodomain with TIR motif. Primary OfTLR5M sequence shared a higher homology with teleost TLR5M. The evolutional analysis confirmed that TLR5 identified in the current study is a membrane receptor and the data further suggested the co-evolution of the membrane-anchored and soluble forms of TLR5 in teleosts. Inter-lineage comparison of gene structures in vertebrates indicated that the tlr5m gene has evolved with extensive rearrangement; whereas, the myd88 gene has maintained a stable structure throughout the evolution. Inspection of 5' flanking region of these genes disclosed the presence of several transcription factor binding sites including NF-κB. Quantitative real-time PCR (qPCR) detected Oftlr5m mRNA in eleven tissues with the highest abundance in liver. In vivo flagellin administration strongly induced the transcripts of both Oftlr5m and Ofmyd88 in gills and head kidney tissues suggesting their ligand-mediated upregulation. In a luciferase assay, HEK293T cells transiently transfected with Oftlr5m and Ofmyd88 demonstrated a higher NF-κB activity than the mock control, and the luciferase activity was intensified when cells were stimulated with flagellin. Collectively, our study represents the genomic, evolutional, expressional and functional insights into a receptor and adaptor molecules of teleost origin that are involved in flagellin sensing.