Tcf7l2 in hepatocytes regulates de novo lipogenesis in diet-induced non-alcoholic fatty liver disease in mice.

AIMS/HYPOTHESIS: Non-alcoholic fatty liver disease (NAFLD) associated with type 2 diabetes may more easily progress towards severe forms of non-alcoholic steatohepatitis (NASH) and cirrhosis. Although the Wnt effector transcription factor 7-like 2 (TCF7L2) is closely associated with type 2 diabetes risk, the role of TCF7L2 in NAFLD development remains unclear. Here, we investigated how changes in TCF7L2 expression in the liver affects hepatic lipid metabolism based on the major risk factors of NAFLD development. METHODS: Tcf7l2 was selectively ablated in the liver of C57BL/6N mice by inducing the albumin (Alb) promoter to recombine Tcf7l2 alleles floxed at exon 5 (liver-specific Tcf7l2-knockout [KO] mice: Alb-Cre;Tcf7l2f/f). Alb-Cre;Tcf7l2f/f and their wild-type (Tcf7l2f/f) littermates were fed a high-fat diet (HFD) or a high-carbohydrate diet (HCD) for 22 weeks to reproduce NAFLD/NASH. Mice were refed a standard chow diet or an HCD to stimulate de novo lipogenesis (DNL) or fed an HFD to provide exogenous fatty acids. We analysed glucose and insulin sensitivity, metabolic respiration, mRNA expression profiles, hepatic triglyceride (TG), hepatic DNL, selected hepatic metabolites, selected plasma metabolites and liver histology. RESULTS: Alb-Cre;Tcf7l2f/f essentially exhibited increased lipogenic genes, but there were no changes in hepatic lipid content in mice fed a normal chow diet. However, following 22 weeks of diet-induced NAFLD/NASH conditions, liver steatosis was exacerbated owing to preferential metabolism of carbohydrate over fat. Indeed, hepatic Tcf7l2 deficiency enhanced liver lipid content in a manner that was dependent on the duration and amount of exposure to carbohydrates, owing to cell-autonomous increases in hepatic DNL. Mechanistically, TCF7L2 regulated the transcriptional activity of Mlxipl (also known as ChREBP) by modulating O-GlcNAcylation and protein content of carbohydrate response element binding protein (ChREBP), and targeted Srebf1 (also called SREBP1) via miRNA (miR)-33-5p in hepatocytes. Eventually, restoring TCF7L2 expression at the physiological level in the liver of Alb-Cre;Tcf7l2f/f mice alleviated liver steatosis without altering body composition under both acute and chronic HCD conditions. CONCLUSIONS/INTERPRETATION: In mice, loss of hepatic Tcf7l2 contributes to liver steatosis by inducing preferential metabolism of carbohydrates via DNL activation. Therefore, TCF7L2 could be a promising regulator of the NAFLD associated with high-carbohydrate diets and diabetes since TCF7L2 deficiency may lead to development of NAFLD by promoting utilisation of excess glucose pools through activating DNL. DATA AVAILABILITY: RNA-sequencing data have been deposited into the NCBI GEO under the accession number GSE162449 ( www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE162449 ).

Epigenetic regulation of DHRS2 by SUV420H2 inhibits cell apoptosis in renal cell carcinoma.

Renal cell carcinoma (RCC), also known as kidney cancer, is a common malignant tumor of the urinary system. While surgical treatment is essential, novel therapeutic targets and corresponding drugs for RCC are still needed due to the high relapse rate and low five-year survival rate. In this study, we found that SUV420H2 is overexpressed in renal cancers and that high SUV420H2 expression is associated with a poor prognosis, as evidenced by RCC RNA-seq results derived from the TCGA. SUV420H2 knockdown using siRNA led to growth suppression and cell apoptosis in the A498 cell line. Furthermore, we identified DHRS2 as a direct target of SUV420H2 in the apoptosis process through a ChIP assay with a histone 4 lysine 20 (H4K20) trimethylation antibody. Rescue experiments showed that cotreatment with siSUV420H2 and siDHRS2 attenuated cell growth suppression induced by SUV420H2 knockdown only. Additionally, treatment with the SUV420H2 inhibitor A-196 induced cell apoptosis via upregulation of DHRS2. Taken together, our findings suggest that SUV420H2 may be a potential therapeutic target for the treatment of renal cancer.

Characterization of terminal-ileal and colonic Crohn's disease in treatment-naïve paediatric patients based on transcriptomic profile using logistic regression.

BACKGROUND: Inflammatory bowel disease (IBD) is a chronic and idiopathic inflammatory disorder of the gastrointestinal tract and comprises ulcerative colitis (UC) and Crohn's disease (CD). Crohn's disease can affect any part of the gastrointestinal tract, but mainly the terminal ileum and colon. In the present study, we aimed to characterize terminal-ileal CD (ICD) and colonic CD (CCD) at the molecular level, which might enable a more optimized approach for the clinical care and scientific research of CD. METHODS: We analyzed differentially expressed genes in samples from 23 treatment-naïve paediatric patients with CD and 25 non-IBD controls, and compared the data with previously published RNA-Seq data using multi-statistical tests and confidence intervals. We implemented functional profiling and proposed statistical methods for feature selection using a logistic regression model to identify genes that are highly associated in ICD or CCD. We also validated our final candidate genes in independent paediatric and adult cohorts. RESULTS: We identified 550 genes specifically expressed in patients with CD compared with those in healthy controls (p < 0.05). Among these DEGs, 240 from patients with CCD were mainly involved in mitochondrial dysfunction, whereas 310 from patients with ICD were enriched in the ileum functions such as digestion, absorption, and metabolism. To choose the most effective gene set, we selected the most powerful genes (p-value ≤ 0.05, accuracy ≥ 0.8, and AUC ≥ 0.8) using logistic regression. Consequently, 33 genes were identified as useful for discriminating CD location; the accuracy and AUC were 0.86 and 0.83, respectively. We then validated the 33 genes with data from another independent paediatric cohort (accuracy = 0.93, AUC = 0.92) and adult cohort (accuracy = 0.88, AUC = 0.72). CONCLUSIONS: In summary, we identified DEGs that are specifically expressed in CCD and ICD compared with those in healthy controls and patients with UC. Based on the feature selection analysis, 33 genes were identified as useful for discriminating CCD and ICD with high accuracy and AUC, for not only paediatric patients but also independent cohorts. We propose that our approach and the final gene set are useful for the molecular classification of patients with CD, and it could be beneficial in treatments based on disease location.

RNA-seq based transcriptome analysis of EHMT2 functions in breast cancer.

For breast cancer treatment, hormone therapy is effective for hormone receptor-positive breast cancer but not for TNBC (triple-negative breast cancer). Thus, many researchers have attempted to identify more effective therapeutic candidates for all subtypes of breast cancer. In this study, we established an RNA-seq analytical pipeline to analyze the subtype-specific functions of EHMT2 in the MB231 and MCF7 cell lines. After EHMT2 knockdown, we identified subtype-specific DEGs (differentially expressed genes) and overlapping DEGs. Through GO (Gene Ontology) analysis, GSEA (gene set enrichment analysis), and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis using the DEGs, we identified the subtype-specific functions of EHMT2 in the MB231 and MCF7 cell lines. Therefore, herein, we suggest that EHMT2 is an attractive therapeutic target for the treatment of all types of breast cancer.

Characterization of sexual size dimorphism and sex-biased genes expression profile in the olive flounder.

Sexual size dimorphism (SSD) is a widespread phenomenon in fish species, including in the olive flounder. Although it is well established that female olive flounders acquire more bone mass than males, the underlying mechanism and timing of this SSD remains controversial. Here, the gene expression profiles of adult male and female olive flounder fish were explored to better understand the SSD mechanisms. Using RNA sequencing, a total of 4784 sex-biased differentially expressed genes (DEGs) in the fin with asymptotic growth after maturity were identified, among which growth-related factors were found. Gene ontology and pathway enrichment studies were performed to predict potential SSD-related genes and their functions. According to functional analysis, negative regulation of cell proliferation was significantly enriched in males, and anabolism related genes were highly expressed in females. In addition, pathway analysis using the Kyoto Encyclopedia of Genes and Genomes database revealed that five sexual dimorphism-related candidate genes (bambia, smurf1, dvl2, cul1a, and dvl3) were enriched in osteogenesis-contributing pathways. These results suggest that these five candidate genes may be relevant for skeletal development in olive flounders. Altogether, this study adds new knowledge for a better understanding of SSD-related growth traits in olive flounder, which can be used for enhancing aquaculture productivity with reduced production costs.

Stabilization of E2-EPF UCP protein is implicated in hepatitis B virus-associated hepatocellular carcinoma progression.

Hepatitis B virus (HBV) X protein (HBx) is associated with hepatocarcinogenesis. E2-EPF ubiquitin carrier protein (UCP) catalyzes ubiquitination of itself and von Hippel-Lindau protein (pVHL) for degradation and associates with tumor growth and metastasis. However, it remains unknown whether HBx modulates the enzyme activity of UCP and thereby influences hepatocarcinogenesis. Here, we show that UCP is highly expressed in liver tissues of HBx-transgenic mice, but not non-transgenic mice. UCP was more frequently expressed in HBV-positive liver cancers than in HBV-negative liver cancers. HBx binds to UCP specifically and serotype independently, and forms a ternary complex with UCP and pVHL. HBx inhibits self-ubiquitination of UCP, but enhances UCP-mediated pVHL ubiquitination, resulting in stabilization of hypoxia-inducible factor-1α and -2α. HBx and UCP stabilize each other by mutually inhibiting their ubiquitination. HBx promotes cellular proliferation and metastasis via UCP. Our findings suggest that UCP plays a key role in HBV-related hepatocarcinogenesis.

Bioinformatic analysis of regulation of natural antisense transcripts by transposable elements in human mRNA.

Non-coding RNA is no longer considered to be "junk" DNA, based on evidence uncovered in recent decades. In particular, the important role played by natural antisense transcripts (NATs) in regulating the expression of genes is receiving increasing attention. However, the regulatory mechanisms of NATs remain incompletely understood. It is well-known that the insertion of transposable elements (TEs) can affect gene transcription. Using a bioinformatics approach, we identified NATs using human mRNA sequences from the UCSC Genome Browser Database. Our in silico analysis identified 1079 NATs and 700 sense-antisense gene pairs. We identified 179 NATs that showed evidence of having been affected by TEs during cellular gene expression. These findings may provide an understanding of the complex regulation mechanisms of NATs. If our understanding of NATs as modulators of gene expression is further enhanced, we can develop ways to control gene expression.

Generation of expandable human pluripotent stem cell-derived hepatocyte-like liver organoids.

BACKGROUND & AIMS: The development of hepatic models capable of long-term expansion with competent liver functionality is technically challenging in a personalized setting. Stem cell-based organoid technologies can provide an alternative source of patient-derived primary hepatocytes. However, self-renewing and functionally competent human pluripotent stem cell (PSC)-derived hepatic organoids have not been developed. METHODS: We developed a novel method to efficiently and reproducibly generate functionally mature human hepatic organoids derived from PSCs, including human embryonic stem cells and induced PSCs. The maturity of the organoids was validated by a detailed transcriptome analysis and functional performance assays. The organoids were applied to screening platforms for the prediction of toxicity and the evaluation of drugs that target hepatic steatosis through real-time monitoring of cellular bioenergetics and high-content analyses. RESULTS: Our organoids were morphologically indistinguishable from adult liver tissue-derived epithelial organoids and exhibited self-renewal. With further maturation, their molecular features approximated those of liver tissue, although these features were lacking in 2D differentiated hepatocytes. Our organoids preserved mature liver properties, including serum protein production, drug metabolism and detoxifying functions, active mitochondrial bioenergetics, and regenerative and inflammatory responses. The organoids exhibited significant toxic responses to clinically relevant concentrations of drugs that had been withdrawn from the market due to hepatotoxicity and recapitulated human disease phenotypes such as hepatic steatosis. CONCLUSIONS: Our organoids exhibit self-renewal (expandable and further able to differentiate) while maintaining their mature hepatic characteristics over long-term culture. These organoids may provide a versatile and valuable platform for physiologically and pathologically relevant hepatic models in the context of personalized medicine. LAY SUMMARY: A functionally mature, human cell-based liver model exhibiting human responses in toxicity prediction and drug evaluation is urgently needed for pre-clinical drug development. Here, we develop a novel human pluripotent stem cell-derived hepatocyte-like liver organoid that is critically advanced in terms of its generation method, functional performance, and application technologies. Our organoids can contribute to the better understanding of liver development and regeneration, and provide insights for metabolic studies and disease modeling, as well as toxicity assessments and drug screening for personalized medicine.

The Ancient Phosphatidylinositol 3-Kinase Signaling System Is a Master Regulator of Energy and Carbon Metabolism in Algae.

Algae undergo a complete metabolic transformation under stress by arresting cell growth, inducing autophagy and hyper-accumulating biofuel precursors such as triacylglycerols and starch. However, the regulatory mechanisms behind this stress-induced transformation are still unclear. Here, we use biochemical, mutational, and "omics" approaches to demonstrate that PI3K signaling mediates the homeostasis of energy molecules and influences carbon metabolism in algae. In Chlamydomonas reinhardtii, the inhibition and knockdown (KD) of algal class III PI3K led to significantly decreased cell growth, altered cell morphology, and higher lipid and starch contents. Lipid profiling of wild-type and PI3K KD lines showed significantly reduced membrane lipid breakdown under nitrogen starvation (-N) in the KD. RNA-seq and network analyses showed that under -N conditions, the KD line carried out lipogenesis rather than lipid hydrolysis by initiating de novo fatty acid biosynthesis, which was supported by tricarboxylic acid cycle down-regulation and via acetyl-CoA synthesis from glycolysis. Remarkably, autophagic responses did not have primacy over inositide signaling in algae, unlike in mammals and vascular plants. The mutant displayed a fundamental shift in intracellular energy flux, analogous to that in tumor cells. The high free fatty acid levels and reduced mitochondrial ATP generation led to decreased cell viability. These results indicate that the PI3K signal transduction pathway is the metabolic gatekeeper restraining biofuel yields, thus maintaining fitness and viability under stress in algae. This study demonstrates the existence of homeostasis between starch and lipid synthesis controlled by lipid signaling in algae and expands our understanding of such processes, with biotechnological and evolutionary implications.

Developing scalable cultivation systems of hepatic spheroids for drug metabolism via genomic and functional analyses.

Spheroids, a widely used three-dimensional (3D) culture model, are standard in hepatocyte culture as they preserve long-term hepatocyte functionality and enhance survivability. In this study, we investigated the effects of three operation modes in 3D culture - static, orbital shaking, and under vertical bidirectional flow using spheroid forming units (SFUs) on hepatic differentiation and drug metabolism to propose the best for mass production of functionally enhanced spheroids. Spheroids in SFUs exhibited increased hepatic gene expression, albumin secretion, and cytochrome P450 3A4 (CYP3A4) activity during the differentiation period (12 days). SFUs advantages include facilitated mass production and a relatively earlier peak of CYP3A4 activity. However, CYP3A4 activity was not well maintained under dimethyl sulfoxide (DMSO)-free conditions (13-18 days), dramatically reducing drug metabolism capability. Continued shear stimulation without differentiation stimuli in assay conditions markedly attenuated CYP3A4 activity, which was less severe in static conditions. In this condition, SFU spheroids exhibited dedifferentiation characteristics, such as increased proliferation and Notch signaling genes. We found that the dedifferentiation could be overcome by using the serum-free medium formulation. Therefore, we suggest that SFUs represent the best option for the mass production of functionally improved spheroids and so the serum-free conditions should be maintained during drug metabolism analysis.

EHMT2 is a metastasis regulator in breast cancer.

Various modes of epigenetic regulation of breast cancer proliferation and metastasis have been investigated, but epigenetic mechanisms involved in breast cancer metastasis remain elusive. Thus, in this study, EHMT2 (a histone methyltransferase) was determined to be significantly overexpressed in breast cancer tissues and in Oncomine data. In addition, knockdown of EHMT2 reduced cell migration/invasion and regulated the expression of EMT-related markers (E-cadherin, Claudin 1, and Vimentin). Furthermore, treatment with BIX-01294, a specific inhibitor of EHMT2, affected migration/invasion in MDA-MB-231 cells. Therefore, our findings demonstrate functions of EHMT2 in breast cancer metastasis and suggest that targeting EHMT2 may be an effective therapeutic strategy for preventing breast cancer metastasis.

A liver-specific gene expression panel predicts the differentiation status of in vitro hepatocyte models.

Alternative cell sources, such as three-dimensional organoids and induced pluripotent stem cell-derived cells, might provide a potentially effective approach for both drug development applications and clinical transplantation. For example, the development of cell sources for liver cell-based therapy has been increasingly needed, and liver transplantation is performed for the treatment for patients with severe end-stage liver disease. Differentiated liver cells and three-dimensional organoids are expected to provide new cell sources for tissue models and revolutionary clinical therapies. However, conventional experimental methods confirming the expression levels of liver-specific lineage markers cannot provide complete information regarding the differentiation status or degree of similarity between liver and differentiated cell sources. Therefore, in this study, to overcome several issues associated with the assessment of differentiated liver cells and organoids, we developed a liver-specific gene expression panel (LiGEP) algorithm that presents the degree of liver similarity as a "percentage." We demonstrated that the percentage calculated using the LiGEP algorithm was correlated with the developmental stages of in vivo liver tissues in mice, suggesting that LiGEP can correctly predict developmental stages. Moreover, three-dimensional cultured HepaRG cells and human pluripotent stem cell-derived hepatocyte-like cells showed liver similarity scores of 59.14% and 32%, respectively, although general liver-specific markers were detected. CONCLUSION: Our study describes a quantitative and predictive model for differentiated samples, particularly liver-specific cells or organoids; and this model can be further expanded to various tissue-specific organoids; our LiGEP can provide useful information and insights regarding the differentiation status of in vitro liver models. (Hepatology 2017;66:1662-1674).

Comparative in silico profiling of epigenetic modifiers in human tissues.

The technology of tissue differentiation from human pluripotent stem cells has attracted attention as a useful resource for regenerative medicine, disease modeling and drug development. Recent studies have suggested various key factors and specific culture methods to improve the successful tissue differentiation and efficient generation of human induced pluripotent stem cells. Among these methods, epigenetic regulation and epigenetic signatures are regarded as an important hurdle to overcome during reprogramming and differentiation. Thus, in this study, we developed an in silico epigenetic panel and performed a comparative analysis of epigenetic modifiers in the RNA-seq results of 32 human tissues. We demonstrated that an in silico epigenetic panel can identify epigenetic modifiers in order to overcome epigenetic barriers to tissue-specific differentiation.

Combinatory RNA-Sequencing Analyses Reveal a Dual Mode of Gene Regulation by ADAR1 in Gastric Cancer.

BACKGROUND: Adenosine deaminase acting on RNA 1 (ADAR1) is known to mediate deamination of adenosine-to-inosine through binding to double-stranded RNA, the phenomenon known as RNA editing. Currently, the function of ADAR1 in gastric cancer is unclear. AIMS: This study was aimed at investigating RNA editing-dependent and editing-independent functions of ADAR1 in gastric cancer, especially focusing on its influence on editing of 3' untranslated regions (UTRs) and subsequent changes in expression of messenger RNAs (mRNAs) as well as microRNAs (miRNAs). METHODS: RNA-sequencing and small RNA-sequencing were performed on AGS and MKN-45 cells with a stable ADAR1 knockdown. Changed frequencies of editing and mRNA and miRNA expression were then identified by bioinformatic analyses. Targets of RNA editing were further validated in patients' samples. RESULTS: In the Alu region of both gastric cell lines, editing was most commonly of the A-to-I type in 3'-UTR or intron. mRNA and protein levels of PHACTR4 increased in ADAR1 knockdown cells, because of the loss of seed sequences in 3'-UTR of PHACTR4 mRNA that are required for miRNA-196a-3p binding. Immunohistochemical analyses of tumor and paired normal samples from 16 gastric cancer patients showed that ADAR1 expression was higher in tumors than in normal tissues and inversely correlated with PHACTR4 staining. On the other hand, decreased miRNA-148a-3p expression in ADAR1 knockdown cells led to increased mRNA and protein expression of NFYA, demonstrating ADAR1's editing-independent function. CONCLUSIONS: ADAR1 regulates post-transcriptional gene expression in gastric cancer through both RNA editing-dependent and editing-independent mechanisms.

Profiling of the Major Phenolic Compounds and Their Biosynthesis Genes in Sophora flavescens Aiton.

Sophorae Radix (Sophora flavescens Aiton) has long been used in traditional medicine in East Asia due to the various biological activities of its secondary metabolites. Endogenous contents of phenolic compounds (phenolic acid, flavonol, and isoflavone) and the main bioactive compounds of Sophorae Radix were analyzed based on the qualitative HPLC analysis and evaluated in different organs and at different developmental stages. In total, 11 compounds were detected, and the composition of the roots and aerial parts (leaves, stems, and flowers) was significantly different. trans-Cinnamic acid and p-coumaric acid were observed only in the aerial parts. Large amounts of rutin and maackiain were detected in the roots. Four phenolic acid compounds (benzoic acid, caffeic acid, ferulic acid, and chlorogenic acid) and four flavonol compounds (kaempferol, catechin hydrate, epicatechin, and rutin) were higher in aerial parts than in roots. To identify putative genes involved in phenolic compounds biosynthesis, a total of 41 transcripts were investigated. Expression patterns of these selected genes, as well as the multiple isoforms for the genes, varied by organ and developmental stage, implying that they are involved in the biosynthesis of various phenolic compounds both spatially and temporally.

Vibrio injenensis sp. nov., isolated from human clinical specimens.

Vibrio species are well known as motile, mostly oxidase-positive, facultative anaerobic Gram-negative bacteria. They are abundant in aquatic environments and are a common cause of human infections including diarrhea, soft tissue diseases, and bacteremia. Here, two Gram-negative bacteria, designated M12-1144T and M12-1181, were isolated from human clinical specimens and identified using a polyphasic taxonomic approach. Phylogenetic study based on 16S rRNA gene sequence analysis revealed that the isolates belong to the genus Vibrio, and are closely related to Vibrio metschnikovii KCTC 32284T (98.3%) and Vibrio cincinnatiensis KCTC 2733T (97.8%). The major fatty acids were summed feature 3 (C16:1 ω7c/C16:1 ω6c, 38.0%), C16:0 (23.0%), and summed feature 8 (C18:1 ω7c or C18:1 ω6c, 19.3%) and major polar lipids were diphosphatidylglycerol, phosphatidylglycerol, and phosphatidylethanolamine. The G + C content of the genomic DNA was determined to be 44.1 mol%. DNA-DNA relatedness between the two newly isolated strains and V. metschnikovii KCTC 32284T and V. cincinnatiensis KCTC 2733T was between 42.6 to 47.5%. The similarities of genome-to-genome distance between M12-1144T and related species ranged from 18.4-54.8%. Based on these results, a new species of the genus Vibrio, Vibrio injenensis is proposed. The type strain is M12-1144 T(=KCTC 32233T =JCM 30011T).

Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway.

The mammalian Y chromosome has unique characteristics compared with the autosomes or X chromosomes. Here we report the finished sequence of the chimpanzee Y chromosome (PTRY), including 271 kb of the Y-specific pseudoautosomal region 1 and 12.7 Mb of the male-specific region of the Y chromosome. Greater sequence divergence between the human Y chromosome (HSAY) and PTRY (1.78%) than between their respective whole genomes (1.23%) confirmed the accelerated evolutionary rate of the Y chromosome. Each of the 19 PTRY protein-coding genes analyzed had at least one nonsynonymous substitution, and 11 genes had higher nonsynonymous substitution rates than synonymous ones, suggesting relaxation of selective constraint, positive selection or both. We also identified lineage-specific changes, including deletion of a 200-kb fragment from the pericentromeric region of HSAY, expansion of young Alu families in HSAY and accumulation of young L1 elements and long terminal repeat retrotransposons in PTRY. Reconstruction of the common ancestral Y chromosome reflects the dynamic changes in our genomes in the 5-6 million years since speciation.

Functional innovations of three chronological mesohexaploid Brassica rapa genomes.

BACKGROUND: The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa. RESULTS: Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were 'stress response' and 'development' both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization. CONCLUSION: We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Genome-wide analysis of DNA methylation patterns in horse.

BACKGROUND: DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known. RESULTS: We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density. CONCLUSIONS: We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.