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.

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.

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.

Discovery of Gene Sources for Economic Traits in Hanwoo by Whole-genome Resequencing.

Hanwoo, a Korean native cattle (Bos taurus coreana), has great economic value due to high meat quality. Also, the breed has genetic variations that are associated with production traits such as health, disease resistance, reproduction, growth as well as carcass quality. In this study, next generation sequencing technologies and the availability of an appropriate reference genome were applied to discover a large amount of single nucleotide polymorphisms (SNPs) in ten Hanwoo bulls. Analysis of whole-genome resequencing generated a total of 26.5 Gb data, of which 594,716,859 and 592,990,750 reads covered 98.73% and 93.79% of the bovine reference genomes of UMD 3.1 and Btau 4.6.1, respectively. In total, 2,473,884 and 2,402,997 putative SNPs were discovered, of which 1,095,922 (44.3%) and 982,674 (40.9%) novel SNPs were discovered against UMD3.1 and Btau 4.6.1, respectively. Among the SNPs, the 46,301 (UMD 3.1) and 28,613 SNPs (Btau 4.6.1) that were identified as Hanwoo-specific SNPs were included in the functional genes that may be involved in the mechanisms of milk production, tenderness, juiciness, marbling of Hanwoo beef and yellow hair. Most of the Hanwoo-specific SNPs were identified in the promoter region, suggesting that the SNPs influence differential expression of the regulated genes relative to the relevant traits. In particular, the non-synonymous (ns) SNPs found in CORIN, which is a negative regulator of Agouti, might be a causal variant to determine yellow hair of Hanwoo. Our results will provide abundant genetic sources of variation to characterize Hanwoo genetics and for subsequent breeding.

Single nucleotide polymorphism analysis of Korean native chickens using next generation sequencing data.

There are five native chicken lines in Korea, which are mainly classified by plumage colors (black, white, red, yellow, gray). These five lines are very important genetic resources in the Korean poultry industry. Based on a next generation sequencing technology, whole genome sequence and reference assemblies were performed using Gallus_gallus_4.0 (NCBI) with whole genome sequences from these lines to identify common and novel single nucleotide polymorphisms (SNPs). We obtained 36,660,731,136 ± 1,257,159,120 bp of raw sequence and average 26.6-fold of 25-29 billion reference assembly sequences representing 97.288 % coverage. Also, 4,006,068 ± 97,534 SNPs were observed from 29 autosomes and the Z chromosome and, of these, 752,309 SNPs are the common SNPs across lines. Among the identified SNPs, the number of novel- and known-location assigned SNPs was 1,047,951 ± 14,956 and 2,948,648 ± 81,414, respectively. The number of unassigned known SNPs was 1,181 ± 150 and unassigned novel SNPs was 8,238 ± 1,019. Synonymous SNPs, non-synonymous SNPs, and SNPs having character changes were 26,266 ± 1,456, 11,467 ± 604, 8,180 ± 458, respectively. Overall, 443,048 ± 26,389 SNPs in each bird were identified by comparing with dbSNP in NCBI. The presently obtained genome sequence and SNP information in Korean native chickens have wide applications for further genome studies such as genetic diversity studies to detect causative mutations for economic and disease related traits.

Antimicrobial peptides in the centipede Scolopendra subspinipes mutilans.

The centipede Scolopendra subspinipes mutilans is an environmentally beneficial and medically important arthropod species. Although this species is increasingly applied as a reliable source of new antimicrobial peptides, the transcriptome of this species is a prerequisite for more rational selection of antimicrobial peptides. In this report, we isolated total RNA from the whole body of adult centipedes, S. subspinipes mutilans, that were nonimmunized and immunized against Escherichia coli, and we generated a total of 77,063 pooled contigs and singletons using high-throughput sequencing. To screen putative antimicrobial peptides, in silico analyses of the S. subspinipes mutilans transcriptome were performed based on the physicochemical evidence of length, charge, isoelectric point, and in vitro and in vivo aggregation scores together with the existence of continuous antimicrobial peptide stretches. Moreover, we excluded some transcripts that showed similarity with both previously known antimicrobial peptides and the human proteome, had a proteolytic cleavage site, and had downregulated expression compared with the nonimmunized sample. As a result, we selected 17 transcripts and tested their antimicrobial activity with a radial diffusion assay. Among them, ten synthetic peptides experimentally showed antimicrobial activity against microbes and no toxicity to mouse erythrocytes. Our results provide not only a useful set of antimicrobial peptide candidates and an efficient strategy for novel antimicrobial peptide development but also the transcriptome data of a big centipede as a valuable resource.

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.

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.

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.

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.

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.

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.

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.

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.

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.