DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus).

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10,000 years, offers an interesting ecological model to address such effects. Using bisulphite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic-limnetic axis), both at single CpG sites and in 1000 bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

Genetic association of gemcitabine/carboplatin-induced leukopenia and neutropenia in non-small cell lung cancer patients using whole-exome sequencing.

OBJECTIVES: Gemcitabine/carboplatin treatment is known to cause severe adverse drug reactions which can lead to the need for reduction or cessation of chemotherapy. It would be beneficial to identify patients at risk of severe hematological toxicity in advance before treatment start. This study aims to identify genetic markers for gemcitabine/carboplatin-induced leukopenia and neutropenia in non-small cell lung cancer patients. MATERIAL AND METHODS: Whole-exome sequencing was performed on 215 patients. Association analysis was performed on single-nucleotide variants (SNVs) and genes, and the validation was based on an independent genome-wide association study (GWAS). Based on the association and validation analyses the genetic variants were then selected for and used in weighted genetic risk score (wGRS) prediction models for leukopenia and neutropenia. RESULTS: Association analysis identified 50 and 111 SNVs, and 12 and 20 genes, for leukopenia and neutropenia, respectively. Of these SNVS 20 and 19 were partially validated for leukopenia and neutropenia, respectively. The genes SVIL (p = 2.48E-06) and EFCAB2 (p = 4.63E-06) were significantly associated with leukopenia contain the partially validated SNVs rs3740003, rs10160013, rs1547169, rs10927386 and rs10927387. The wGRS prediction models showed significantly different risk scores for high and low toxicity patients. CONCLUSION: We have identified and partially validated genetic biomarkers in SNVs and genes correlated to gemcitabine/carboplatin-induced leukopenia and neutropenia and created wGRS models for predicting the risk of chemotherapy-induced hematological toxicity. These results provide a strong foundation for further studies of chemotherapy-induced toxicity.

Genes and variants in hematopoiesis-related pathways are associated with gemcitabine/carboplatin-induced thrombocytopenia.

Chemotherapy-induced myelosuppression, including thrombocytopenia, is a recurrent problem during cancer treatments that may require dose alterations or cessations that could affect the antitumor effect of the treatment. To identify genetic markers associated with treatment-induced thrombocytopenia, we whole-exome sequenced 215 non-small cell lung cancer patients homogeneously treated with gemcitabine/carboplatin. The decrease in platelets (defined as nadir/baseline) was used to assess treatment-induced thrombocytopenia. Association between germline genetic variants and thrombocytopenia was analyzed at single-nucleotide variant (SNV) (based on the optimal false discovery rate, the severity of predicted consequence, and effect), gene, and pathway levels. These analyses identified 130 SNVs/INDELs and 25 genes associated with thrombocytopenia (P-value < 0.002). Twenty-three SNVs were validated in an independent genome-wide association study (GWAS). The top associations include rs34491125 in JMJD1C (P-value = 9.07 × 10-5), the validated variants rs10491684 in DOCK8 (P-value = 1.95 × 10-4), rs6118 in SERPINA5 (P-value = 5.83 × 10-4), and rs5877 in SERPINC1 (P-value = 1.07 × 10-3), and the genes CAPZA2 (P-value = 4.03 × 10-4) and SERPINC1 (P-value = 1.55 × 10-3). The SNVs in the top-scoring pathway "Factors involved in megakaryocyte development and platelet production" (P-value = 3.34 × 10-4) were used to construct weighted genetic risk score (wGRS) and logistic regression models that predict thrombocytopenia. The wGRS predict which patients are at high or low toxicity risk levels, for CTCAE (odds ratio (OR) = 22.35, P-value = 1.55 × 10-8), and decrease (OR = 66.82, P-value = 5.92 × 10-9). The logistic regression models predict CTCAE grades 3-4 (receiver operator characteristics (ROC) area under the curve (AUC) = 0.79), and large decrease (ROC AUC = 0.86). We identified and validated genetic variations within hematopoiesis-related pathways that provide a solid foundation for future studies using genetic markers for predicting chemotherapy-induced thrombocytopenia and personalizing treatments.

Comprehensive RNA sequencing of healthy human endometrium at two time points of the menstrual cycle.

Endometrial receptivity is crucial for implantation and establishment of a normal pregnancy. The shift from proliferative to receptive endometrium is still far from being understood. In this paper, we comprehensively present the transcriptome of the human endometrium by comparing endometrial biopsies from proliferative phase with consecutive biopsies 7-9 days after ovulation. The results show a clear difference in expression between the two time points using both total and small RNA sequencing. A total of 3,297 messenger RNAs (mRNAs), 516 long noncoding RNAs (lncRNAs), and 102 small noncoding RNAs were identified as statistically differentially expressed between the two time points. We show a thorough description of the change in mRNA between the two time points and display lncRNAs, small nucleolar RNAs, and small nuclear RNAs not previously reported in the healthy human endometrium. In conclusion, this paper reports in detail the shift in RNA expression from the proliferative to receptive endometrium.

Proteogenomics produces comprehensive and highly accurate protein-coding gene annotation in a complete genome assembly of Malassezia sympodialis.

Complete and accurate genome assembly and annotation is a crucial foundation for comparative and functional genomics. Despite this, few complete eukaryotic genomes are available, and genome annotation remains a major challenge. Here, we present a complete genome assembly of the skin commensal yeast Malassezia sympodialis and demonstrate how proteogenomics can substantially improve gene annotation. Through long-read DNA sequencing, we obtained a gap-free genome assembly for M. sympodialis (ATCC 42132), comprising eight nuclear and one mitochondrial chromosome. We also sequenced and assembled four M. sympodialis clinical isolates, and showed their value for understanding Malassezia reproduction by confirming four alternative allele combinations at the two mating-type loci. Importantly, we demonstrated how proteomics data could be readily integrated with transcriptomics data in standard annotation tools. This increased the number of annotated protein-coding genes by 14% (from 3612 to 4113), compared to using transcriptomics evidence alone. Manual curation further increased the number of protein-coding genes by 9% (to 4493). All of these genes have RNA-seq evidence and 87% were confirmed by proteomics. The M. sympodialis genome assembly and annotation presented here is at a quality yet achieved only for a few eukaryotic organisms, and constitutes an important reference for future host-microbe interaction studies.

Analysis of stranded information using an automated procedure for strand specific RNA sequencing.

BACKGROUND: Strand specific RNA sequencing is rapidly replacing conventional cDNA sequencing as an approach for assessing information about the transcriptome. Alongside improved laboratory protocols the development of bioinformatical tools is steadily progressing. In the current procedure the Illumina TruSeq library preparation kit is used, along with additional reagents, to make stranded libraries in an automated fashion which are then sequenced on Illumina HiSeq 2000. By the use of freely available bioinformatical tools we show, through quality metrics, that the protocol is robust and reproducible. We further highlight the practicality of strand specific libraries by comparing expression of strand specific libraries to non-stranded libraries, by looking at known antisense transcription of pseudogenes and by identifying novel transcription. Furthermore, two ribosomal depletion kits, RiboMinus and RiboZero, are compared and two sequence aligners, Tophat2 and STAR, are also compared. RESULTS: The, non-stranded, Illumina TruSeq kit can be adapted to generate strand specific libraries and can be used to access detailed information on the transcriptome. The RiboZero kit is very effective in removing ribosomal RNA from total RNA and the STAR aligner produces high mapping yield in a short time. Strand specific data gives more detailed and correct results than does non-stranded data as we show when estimating expression values and in assembling transcripts. Even well annotated genomes need improvements and corrections which can be achieved using strand specific data. CONCLUSIONS: Researchers in the field should strive to use strand specific data; it allows for more confidence in the data analysis and is less likely to lead to false conclusions. If faced with analysing non-stranded data, researchers should be well aware of the caveats of that approach.

The age and genomic integrity of neurons after cortical stroke in humans.

It has been unclear whether ischemic stroke induces neurogenesis or neuronal DNA rearrangements in the human neocortex. Using immunohistochemistry; transcriptome, genome and ploidy analyses; and determination of nuclear bomb test-derived (14)C concentration in neuronal DNA, we found neither to be the case. A large proportion of cortical neurons displayed DNA fragmentation and DNA repair a short time after stroke, whereas neurons at chronic stages after stroke showed DNA integrity, demonstrating the relevance of an intact genome for survival.

Sequencing degraded RNA addressed by 3' tag counting.

RNA sequencing has become widely used in gene expression profiling experiments. Prior to any RNA sequencing experiment the quality of the RNA must be measured to assess whether or not it can be used for further downstream analysis. The RNA integrity number (RIN) is a scale used to measure the quality of RNA that runs from 1 (completely degraded) to 10 (intact). Ideally, samples with high RIN (> 8) are used in RNA sequencing experiments. RNA, however, is a fragile molecule which is susceptible to degradation and obtaining high quality RNA is often hard, or even impossible when extracting RNA from certain clinical tissues. Thus, occasionally, working with low quality RNA is the only option the researcher has. Here we investigate the effects of RIN on RNA sequencing and suggest a computational method to handle data from samples with low quality RNA which also enables reanalysis of published datasets. Using RNA from a human cell line we generated and sequenced samples with varying RINs and illustrate what effect the RIN has on the basic procedure of RNA sequencing; both quality aspects and differential expression. We show that the RIN has systematic effects on gene coverage, false positives in differential expression and the quantification of duplicate reads. We introduce 3' tag counting (3TC) as a computational approach to reliably estimate differential expression for samples with low RIN. We show that using the 3TC method in differential expression analysis significantly reduces false positives when comparing samples with different RIN, while retaining reasonable sensitivity.

The ontogeny of sleep-wake cycles in zebrafish: a comparison to humans.

Zebrafish (Danio rerio) are used extensively in sleep research; both to further understanding of sleep in general and also as a model of human sleep. To date, sleep studies have been performed in larval and adult zebrafish but no efforts have been made to document the ontogeny of zebrafish sleep-wake cycles. Because sleep differs across phylogeny and ontogeny it is important to validate the use of zebrafish in elucidating the neural substrates of sleep. Here we describe the development of sleep and wake across the zebrafish lifespan and how it compares to humans. We find power-law distributions to best fit wake bout data but demonstrate that exponential distributions, previously used to describe sleep bout distributions, fail to adequately account for the data in either species. Regardless, the data reveal remarkable similarities in the ontogeny of sleep cycles in zebrafish and humans. Moreover, as seen in other organisms, zebrafish sleep levels are highest early in ontogeny and sleep and wake bouts gradually consolidate to form the adult sleep pattern. Finally, sleep percentage, bout duration, bout number, and sleep fragmentation are shown to allow for meaningful comparisons between zebrafish and human sleep.

Sleep-wake dynamics under extended light and extended dark conditions in adult zebrafish.

We characterize the effects of sleep deprivation on sleep-wake behavior, neurogenesis and stress in adult zebrafish, and describe light-induced changes in gene expression. Sleep deprivation was performed using two stimuli: mild electroshock and light. Comparisons were made between five groups of fish: naïve; electroshock sleep-deprived and yoked-control; fish exposed to constant light (increasing wakefulness); and fish exposed to constant darkness (increasing sleep). Behavioral parameters assessed were sleep percentage, number of sleep-wake transitions, and sleep and wake bout length. Using microarray technology, light-dark modulation of gene expression was examined. In parallel with gene expression, neurogenesis was measured and stress following sleep deprivation was assessed behaviorally and physiologically. Our results indicate that sleep duration is most effectively altered by varying exposure to ambient light. Further, while the sleep-wake dynamics are comparable to those observed in mammals, zebrafish may exhibit weaker sleep homeostasis and sleep pressure than do mammals; and sleep deprivation does not significantly alter their stress responses. Finally, modulation of gene expression by light and dark was observed. Genes upregulated during the dark period are broadly related to growth, morphogenesis, energy balance, and lipid synthesis. Genes upregulated during light are broadly related to synaptic plasticity and cell proliferation.

Effects of modafinil on sleep-wake cycles in larval zebrafish.

We describe, for the first time, the effects of the wakefulness-promoting drug modafinil on sleep and wakefulness in larval zebrafish. Modafinil is currently used to treat narcolepsy, hypersomnia, and shift-work disorder by increasing wakefulness. Tolerance and dependence are limited with modafinil use, differentiating it from common stimulants; however, the neural mechanisms of action of modafinil are still unknown. Zebrafish, a low-cost, prolific, and genetically tractable animal model, have recently become a key model in sleep research. Zebrafish express circadian rhythms, sleep homeostasis, and sleep pressure, and, in addition, respond to common hypnotics and stimulants in a manner similar to mammals. Therefore, in the current experiment we characterize the effects of modafinil on sleep-wake cycles in larval zebrafish as a first step to gaining further insight into the neural mechanisms underlying the effects of modafinil. We show that modafinil modulates sleep-wake activity in larval zebrafish in a manner consistent with what would be predicted from mammalian data. Modafinil increases wakefulness by lengthening wake-bouts, an effect that likely restricted to the night (lights-off). These results validate the use of zebrafish as an animal model for the study of sleep and provide a means for dissecting the neural mechanisms of modafinil, and, more broadly, sleep disorders.