Sex-inducing effects toward planarians widely present among parasitic flatworms.

Various parasitic flatworms infect vertebrates for sexual reproduction, often causing devastating diseases in their hosts. Consequently, flatworms are of great socioeconomic and biomedical importance. Although the cessation of parasitic flatworm sexual reproduction is a major target of anti-parasitic drug design, little is known regarding bioactive compounds controlling flatworm sexual maturation. Using the planarian Dugesia ryukyuensis, we observed that sex-inducing substances found in planarians are also widespread in parasitic flatworms, such as monogeneans and flukes (but not in tapeworms). Reverse-phase HPLC analysis revealed the sex-inducing substance(s) eluting around the tryptophan retention time in the fluke Calicophoron calicophorum, consistent with previous studies on the planarian Bipalium nobile, suggesting that the substance(s) is likely conserved among flatworms. Moreover, six of the 18 ovary-inducing substances identified via transcriptome and metabolome analyses are involved in purine metabolism. Our findings provide a basis for understanding and modifying the life cycles of various parasitic flatworms.

Establishment of the immunological self in juvenile Patiria pectinifera post-metamorphosis.

Ontogeny of the immune system is a fundamental immunology issue. One indicator of immune system maturation is the establishment of the immunological self, which describes the ability of the immune system to distinguish allogeneic individuals (allorecognition ability). However, the timing of immune system maturation during invertebrate ontogeny is poorly understood. In the sea star Patiria pectinifera, cells that have dissociated from the embryos and larvae are able to reconstruct larvae. This reconstruction phenomenon is possible because of a lack of allorecognition capability in the larval immune system, which facilitates the formation of an allogeneic chimera. In this study, we revealed that the adult immune cells of P. pectinifera (coelomocytes) have allorecognition ability. Based on a hypothesis that allorecognition ability is acquired before and after metamorphosis, we conducted detailed morphological observations and survival time analysis of metamorphosis-induced chimeric larvae. The results showed that all allogeneic chimeras died within approximately two weeks to one month of reaching the juvenile stage. In these chimeras, the majority of the epidermal cell layer was lost and the mesenchymal region expanded, but cell death appeared enhanced in the digestive tract. These results indicate that the immunological self of P. pectinifera is established post-metamorphosis during the juvenile stage. This is the first study to identify the timing of immune system maturation during echinodermal ontogenesis. As well as establishing P. pectinifera as an excellent model for studies on self- and non-self-recognition, this study enhances our understanding of the ontogeny of the immune system in invertebrates.

Evaluation of short-term epigenetic age fluctuation.

Considerable effort has been spent on lowering and maintaining the epigenetic age. However, the extent to which epigenetic age fluctuates under normal conditions is poorly understood. Therefore, we analyzed methylation data from monocytes and peripheral blood mononuclear cells collected from two Japanese men. The ranges of the Pan-tissue, Skin and blood, and DNAm PhenoAge epigenetic age during 3 months were ≥ 5.62, ≥ 3.04, and ≥ 8.23 years, and the maximum daily changes were 5.21, 3.20, and 6.53 years, respectively. These fluctuations were not suppressed by correcting for cell-type composition. Although the underlying biological mechanism remains unclear, there was a nonnegligible degree of age fluctuation which should inform personalized clinical applications.

Longitudinal DNA methylation dynamics as a practical indicator in clinical epigenetics.

BACKGROUND: One of the fundamental assumptions of DNA methylation in clinical epigenetics is that DNA methylation status can change over time with or without interplay with environmental and clinical conditions. However, little is known about how DNA methylation status changes over time under ordinary environmental and clinical conditions. In this study, we revisited the high frequency longitudinal DNA methylation data of two Japanese males (24 time-points within three months) and characterized the longitudinal dynamics. RESULTS: The results showed that the majority of CpGs on Illumina HumanMethylation450 BeadChip probe set were longitudinally stable over the time period of three months. Focusing on dynamic and stable CpGs extracted from datasets, dynamic CpGs were more likely to be reported as epigenome-wide association study (EWAS) markers of various traits, especially those of immune- and inflammatory-related traits; meanwhile, the stable CpGs were enriched in metabolism-related genes and were less likely to be EWAS markers, indicating that the stable CpGs are stable both in the short-term within individuals and under various environmental and clinical conditions. CONCLUSIONS: This study indicates that CpGs with different stabilities are involved in different functions and traits, and thus, they are potential indicators that can be applied for clinical epigenetic studies to outline underlying mechanisms.

D-Tryptophan enhances the reproductive organ-specific expression of the amino acid transporter homolog Dr-SLC38A9 involved in the sexual induction of planarian Dugesia ryukyuensis.

BACKGROUND: Many animals switch between asexual and sexual reproduction in nature. We previously established a system for the sexual induction of planarian Dugesia ryukyuensis by feeding asexual planarians with minced sexual planarians. We identified DL-tryptophan (Trp) as one of the sex-inducing substances. DL-Trp can induce ovarian development, the first and essential step of sexual induction. D-Trp must act as a principal bioactive compound in terms of ovarian development, because the ovary-inducing activity of D-Trp was 500 times more potent than that of L-Trp. However, how Trp controls sexual induction is still unknown. RESULTS: In this study, qRT-PCR analyses suggested that the putative amino acid transporter gene Dr-SLC38A9 is highly expressed in sexual worms, especially in the yolk glands. In situ hybridization analyses showed that Dr-SLC38A9 is expressed in the ovarian primordia of asexual worms and in the mature ovaries, testes, and yolk glands of sexual worms. In addition, Dr-SLC38A9 RNA interference during sexual induction resulted in the suppression of the development of reproductive organs. These results suggest that Dr-SLC38A9 is involved in the development of these organs. Moreover, we demonstrated that the reproductive organ-specific expression of Dr-SLC38A9 is enhanced by the addition of D-Trp. CONCLUSION: We propose that D-Trp activates the expression of Dr-SLC38A9 to promote sexual induction in the planarian D. ryukyuensis.

pmar1/phb homeobox genes and the evolution of the double-negative gate for endomesoderm specification in echinoderms.

In several model animals, the earliest phases of embryogenesis are regulated by lineage-specific genes, such as Drosophila bicoid Sea urchin (echinoid) embryogenesis is initiated by zygotic expression of pmar1, a paired-class homeobox gene that has been considered to be present only in the lineage of modern urchins (euechinoids). In euechinoids, Pmar1 promotes endomesoderm specification by repressing the hairy and enhancer of split C (hesC) gene. Here, we have identified the basal echinoid (cidaroid) pmar1 gene, which also promotes endomesoderm specification but not by repressing hesC A further search for related genes demonstrated that other echinoderms have pmar1-related genes named phb Functional analyses of starfish Phb proteins indicated that, similar to cidaroid Pmar1, they promote activation of endomesoderm regulatory gene orthologs via an unknown repressor that is not HesC. Based on these results, we propose that Pmar1 may have recapitulated the regulatory function of Phb during the early diversification of echinoids and that the additional repressor HesC was placed under the control of Pmar1 in the euechinoid lineage. This case provides an exceptional model for understanding how early developmental processes diverge.

Author Correction: Transcriptomic analysis reveals differences in the regulation of amino acid metabolism in asexual and sexual planarians.

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

Transcriptomic analysis reveals differences in the regulation of amino acid metabolism in asexual and sexual planarians.

Many flatworms can alternate between asexual and sexual reproduction. This is a powerful reproductive strategy enabling them to benefit from the features of the two reproductive modes, namely, rapid multiplication and genetic shuffling. The two reproductive modes are enabled by the presence of pluripotent adult stem cells (neoblasts), by generating any type of tissue in the asexual mode, and producing and maintaining germ cells in the sexual mode. In the current study, RNA sequencing (RNA-seq) was used to compare the transcriptomes of two phenotypes of the planarian Dugesia ryukyuensis: an asexual OH strain and an experimentally sexualized OH strain. Pathway enrichment analysis revealed striking differences in amino acid metabolism in the two worm types. Further, the analysis identified serotonin as a new bioactive substance that induced the planarian ovary de novo in a postembryonic manner. These findings suggest that different metabolic states and physiological conditions evoked by sex-inducing substances likely modulate stem cell behavior, depending on their different function in the asexual and sexual reproductive modes. The combination of RNA-seq and a feeding assay in D. ryukyuensis is a powerful tool for studying the alternation of reproductive modes, disentangling the relationship between gene expression and chemical signaling molecules.

iMETHYL: an integrative database of human DNA methylation, gene expression, and genomic variation.

We launched an integrative multi-omics database, iMETHYL (http://imethyl.iwate-megabank.org). iMETHYL provides whole-DNA methylation (~24 million autosomal CpG sites), whole-genome (~9 million single-nucleotide variants), and whole-transcriptome (>14 000 genes) data for CD4+ T-lymphocytes, monocytes, and neutrophils collected from approximately 100 subjects. These data were obtained from whole-genome bisulfite sequencing, whole-genome sequencing, and whole-transcriptome sequencing, making iMETHYL a comprehensive database.

Genome-wide identification of inter-individually variable DNA methylation sites improves the efficacy of epigenetic association studies.

Epigenome-wide association studies, which searches for blood-based DNA methylation signatures associated with environmental exposures and/or disease susceptibilities, is a promising approach to a better understanding of the molecular aetiology of common diseases. To carry out large-scale epigenome-wide association studies while avoiding false negative detection, an efficient strategy to determine target CpG sites for microarray-based or sequencing-based DNA methylation profiling is essentially needed. Here, we propose and validate a hypothesis that a strategy focusing on CpG sites with high DNA methylation level variability may attain an improved efficacy. Through whole-genome bisulfite sequencing of purified blood cells collected from > 100 apparently healthy subjects, we identified ~2.0 million inter-individually variable CpG sites as potential targets. The efficacy of our strategy was estimated to be 3.7-fold higher than that of the most frequently used strategy. Our catalogue of inter-individually variable CpG sites will accelerate the discovery of clinically relevant DNA methylation biomarkers in future epigenome-wide association studies.

Characterization of TRPA channels in the starfish Patiria pectinifera: involvement of thermally activated TRPA1 in thermotaxis in marine planktonic larvae.

The vast majority of marine invertebrates spend their larval period as pelagic plankton and are exposed to various environmental cues. Here we investigated the thermotaxis behaviors of the bipinnaria larvae of the starfish, Patiria pectinifera, in association with TRPA ion channels that serve as thermal receptors in various animal species. Using a newly developed thermotaxis assay system, we observed that P. pectinifera larvae displayed positive thermotaxis toward high temperatures, including toward temperatures high enough to cause death. In parallel, we identified two TRPA genes, termed PpTRPA1 and PpTRPA basal, from this species. We examined the phylogenetic position, spatial expression, and channel properties of each PpTRPA. Our results revealed the following: (1) The two genes diverged early in animal evolution; (2) PpTRPA1 and PpTRPA basal are expressed in the ciliary band and posterior digestive tract of the larval body, respectively; and (3) PpTRPA1 is activated by heat stimulation as well as by known TRPA1 agonists. Moreover, knockdown and rescue experiments demonstrated that PpTRPA1 is involved in positive thermotaxis in P. pectinifera larvae. This is the first report to reveal that TRPA1 channels regulate the behavioral response of a marine invertebrate to temperature changes during its planktonic larval period.

Genetic Predisposition to Ischemic Stroke: A Polygenic Risk Score.

BACKGROUND AND PURPOSE: The prediction of genetic predispositions to ischemic stroke (IS) may allow the identification of individuals at elevated risk and thereby prevent IS in clinical practice. Previously developed weighted multilocus genetic risk scores showed limited predictive ability for IS. Here, we investigated the predictive ability of a newer method, polygenic risk score (polyGRS), based on the idea that a few strong signals, as well as several weaker signals, can be collectively informative to determine IS risk. METHODS: We genotyped 13 214 Japanese individuals with IS and 26 470 controls (derivation samples) and generated both multilocus genetic risk scores and polyGRS, using the same derivation data set. The predictive abilities of each scoring system were then assessed using 2 independent sets of Japanese samples (KyushuU and JPJM data sets). RESULTS: In both validation data sets, polyGRS was shown to be significantly associated with IS, but weighted multilocus genetic risk scores was not. Comparing the highest with the lowest polyGRS quintile, the odds ratios for IS were 1.75 (95% confidence interval, 1.33-2.31) and 1.99 (95% confidence interval, 1.19-3.33) in the KyushuU and JPJM samples, respectively. Using the KyushuU samples, the addition of polyGRS to a nongenetic risk model resulted in a significant improvement of the predictive ability (net reclassification improvement=0.151; P<0.001). CONCLUSIONS: The polyGRS was shown to be superior to weighted multilocus genetic risk scores as an IS prediction model. Thus, together with the nongenetic risk factors, polyGRS will provide valuable information for individual risk assessment and management of modifiable risk factors.

Intraindividual dynamics of transcriptome and genome-wide stability of DNA methylation.

Cytosine methylation at CpG dinucleotides is an epigenetic mechanism that affects the gene expression profiles responsible for the functional differences in various cells and tissues. Although gene expression patterns are dynamically altered in response to various stimuli, the intraindividual dynamics of DNA methylation in human cells are yet to be fully understood. Here, we investigated the extent to which DNA methylation contributes to the dynamics of gene expression by collecting 24 blood samples from two individuals over a period of 3 months. Transcriptome and methylome association analyses revealed that only ~2% of dynamic changes in gene expression could be explained by the intraindividual variation of DNA methylation levels in peripheral blood mononuclear cells and purified monocytes. These results showed that DNA methylation levels remain stable for at least several months, suggesting that disease-associated DNA methylation markers are useful for estimating the risk of disease manifestation.

Two macrophage migration inhibitory factors regulate starfish larval immune cell chemotaxis.

Immune cell recruitment is critical step in the inflammatory response and associated diseases. However, the underlying regulatory mechanisms are poorly understood in invertebrates. Mesenchyme cells of the starfish larvae, which allowed Metchnikoff to complete his landmark experiments, are important model for analysis of immune cell migration. The present study investigated the role of macrophage migration inhibitory factor (MIF)--an evolutionarily conserved cytokine that is functionally similar to chemokines--in the larvae of the starfish Patiria (Asterina) pectinifera, which were found to possess two orthologs, ApMIF1 and ApMIF2. ApMIF1 and ApMIF2 clustered with mammalian MIF and its homolog D-dopachrome tautomerase (DDT), respectively, in the phylogenetic analysis. In contrast to the functional similarity between mammalian MIF and DDT, ApMIF1 knockdown resulted in the excessive recruitment of mesenchyme cells in vivo, whereas ApMIF2 deficiency inhibited the recruitment of these cells to foreign bodies. Mesenchyme cells migrated along a gradient of recombinant ApMIF2 in vitro, whereas recombinant ApMIF1 completely blocked ApMIF2-induced directed migration. Moreover, the expression patterns of ApMIF1 and ApMIF2 messenger RNA in bacteria-challenged mesenchyme cells were consistent with in vivo observations of cell behaviors. These results indicate that ApMIF1 and ApMIF2 act as chemotactic inhibitory and stimulatory factors, respectively, and coordinately regulate mesenchyme cell recruitment during the immune response in starfish larvae. This is the first report describing opposing functions for MIF- and DDT-like molecules. Our findings provide novel insight into the mechanisms underlying immune regulation in invertebrates.

Adjustment of Cell-Type Composition Minimizes Systematic Bias in Blood DNA Methylation Profiles Derived by DNA Collection Protocols.

Differences in DNA collection protocols may be a potential confounder in epigenome-wide association studies (EWAS) using a large number of blood specimens from multiple biobanks and/or cohorts. Here we show that pre-analytical procedures involved in DNA collection can induce systematic bias in the DNA methylation profiles of blood cells that can be adjusted by cell-type composition variables. In Experiment 1, whole blood from 16 volunteers was collected to examine the effect of a 24 h storage period at 4°C on DNA methylation profiles as measured using the Infinium HumanMethylation450 BeadChip array. Our statistical analysis showed that the P-value distribution of more than 450,000 CpG sites was similar to the theoretical distribution (in quantile-quantile plot, λ = 1.03) when comparing two control replicates, which was remarkably deviated from the theoretical distribution (λ = 1.50) when comparing control and storage conditions. We then considered cell-type composition as a possible cause of the observed bias in DNA methylation profiles and found that the bias associated with the cold storage condition was largely decreased (λ adjusted = 1.14) by taking into account a cell-type composition variable. As such, we compared four respective sample collection protocols used in large-scale Japanese biobanks or cohorts as well as two control replicates. Systematic biases in DNA methylation profiles were observed between control and three of four protocols without adjustment of cell-type composition (λ = 1.12-1.45) and no remarkable biases were seen after adjusting for cell-type composition in all four protocols (λ adjusted = 1.00-1.17). These results revealed important implications for comparing DNA methylation profiles between blood specimens from different sources and may lead to discovery of disease-associated DNA methylation markers and the development of DNA methylation profile-based predictive risk models.

Reduction of systematic bias in transcriptome data from human peripheral blood mononuclear cells for transportation and biobanking.

Transportation of samples is essential for large-scale biobank projects. However, RNA degradation during pre-analytical operations prior to transportation can cause systematic bias in transcriptome data, which may prevent subsequent biomarker identification. Therefore, to collect high-quality biobank samples for expression analysis, specimens must be transported under stable conditions. In this study, we examined the effectiveness of RNA-stabilizing reagents to prevent RNA degradation during pre-analytical operations with an emphasis on RNA from peripheral blood mononuclear cells (PBMCs) to establish a protocol for reducing systematic bias. To this end, we obtained PBMCs from 11 healthy volunteers and analyzed the purity, yield, and integrity of extracted RNA after performing pre-analytical operations for freezing PBMCs at -80°C. We randomly chose 7 samples from 11 samples individually, and systematic bias in expression levels was examined by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), RNA sequencing (RNA-Seq) experiments and data analysis. Our data demonstrated that omission of stabilizing reagents significantly lowered RNA integrity, suggesting substantial degradation of RNA molecules due to pre-analytical freezing. qRT-PCR experiments for 19 selected transcripts revealed systematic bias in the expression levels of five transcripts. RNA-Seq for 25,223 transcripts also suggested that about 40% of transcripts were systematically biased. These results indicated that appropriate reduction in systematic bias is essential in protocols for collection of RNA from PBMCs for large-scale biobank projects. Among the seven commercially available stabilizing reagents examined in this study, qRT-PCR and RNA-Seq experiments consistently suggested that RNALock, RNA/DNA Stabilization Reagent for Blood and Bone Marrow, and 1-Thioglycerol/Homogenization solution could reduce systematic bias. On the basis of the results of this study, we established a protocol to reduce systematic bias in the expression levels of RNA transcripts isolated from PBMCs. We believe that these data provide a novel methodology for collection of high-quality RNA from PBMCs for biobank researchers.

Starfish ApDOCK protein essentially functions in larval defense system operated by mesenchyme cells.

In larvae of the starfish, Asterina pectinifera, mesenchyme cells operate in the defense system through various behaviors. We have investigated mesenchyme cell dynamics during the immune response by identifying ApDOCK, a new member of the DOCK180 superfamily protein. In 4-day-old bipinnaria larvae processed for morpholino oligonucleotide-mediated knockdown of ApDOCK, injection of inorganic foreign substances revealed that (1) mesenchyme cells fail to undergo either directed migration toward a large oil-droplet or persistent spreading on the oil-droplet after contact; (2) neither uptake of micro-beads nor cell-to-cell fusion on the large oil-droplet differed from that of mesenchyme cells from control larvae. Similar behaviors were also recorded in experiments where bacteria were injected. Under culture conditions, the expression level of ApDOCK mRNA was significantly associated with the immunological behavior of mesenchyme cells. Apparently, the mesenchyme cells from ApDOCK loss-of-function larvae exhibited insufficient lamellipodium formation via lack of fibrous form of actin organization at the leading edge. These results suggest that the migratory congregation and persistence of encapsulation of larval mesenchyme cells are intracellularly regulated by ApDOCK protein, and this regulation is associated with organization of cytoskeletal actin.

Characterization of a scavenger receptor cysteine-rich-domain-containing protein of the starfish, Asterina pectinifera: ApSRCR1 acts as an opsonin in the larval and adult innate immune systems.

Proteins containing a scavenger receptor cysteine-rich (SRCR) domain (SRCR proteins) play an important role in the innate immune system of various metazoan animals. In the starfish Asterina pectinifera, mesenchyme cells and coelomocytes govern the two distinct innate immune systems of the larvae and adults, respectively. Here we identify a cDNA encoding a protein containing nine SRCR domains termed ApSRCR1, and present characterization of the molecular structure, expression, subcellular localization and function of ApSRCR1 protein during ontogenesis of this animal. ApSRCR1 protein is a membrane-type protein with a predicted molecular mass of approximately 120 kDa. During ontogenesis, ApSRCR1 protein is de novo synthesized and localizes to cytoplasmic vesicles in both mesenchyme cells and coelomocytes without translation of maternal mRNA; however, the net production and modification by N-glycosylation of ApSRCR1 protein differs in each cell type. In both types of cell, functional inhibition of ApSRCR1 protein leads to incompetent bacterial clearance and failure of aggregate formation. However, this inhibitory effect is weaker in the mesenchyme cells than in the coelomocytes. In the bacteria-sensitized adult, ApSRCR1 protein is up-regulated and digested to enable its secretion into the coelomic fluid. This secreted form of ApSRCR1 protein can apparently bind to bacteria. Overall, we show that ApSRCR1 protein is finely regulated for expression not only during development but also in a sensitive innate immunological situation, and thereupon acts as an opsonin for bacteria to different extents in the larvae and adults of A. pectinifera.

Defense system by mesenchyme cells in bipinnaria larvae of the starfish, Asterina pectinifera.

Here we characterize starfish larval mesenchyme cells, in terms of not only their phagocytic behavior, but also their structural and functional properties as a defense system. Our study reveals the following: (1) most mesenchyme cells construct a dynamic network structure beneath the body wall; (2) mesenchyme cells phagocytically respond to almost all foreign materials and form syncytial aggregates to conceal relatively large amounts and large sizes of foreign material; (3) the morphologies of the syncytial aggregates differ from one another depending on the species and the surface configuration of the cellular foreign material; (4) no mesenchyme cells respond to live mesenchyme cells even though they phagocytose chemically fixed cells; (5) mesenchyme cells phagocytose both cellular constituents effluxed from the ectodermal cells and foreign materials taken into the blastocoel through the body wall. Together, these results suggest that mesenchyme cells are equipped with a spectrum of abilities to engage in a defense system in starfish larva.