Ecohydrologic model with satellite-based data for predicting streamflow in ungauged basins.

Information on water availability in basins can be crucial for making decisions for effective water resource management in basins. As the operation of hydrometric stations in Korea is mainly focused on flood season and large rivers, most basins have lack or no observed data. Consequently, this complicates water resource planning and management. Remote sensing data is emerging as a powerful alternative to hydrological information in ungauged basins. This study investigated the applicability of Satellite-Remote Sensed Data (SRSD) as a source for model calibration in Prediction in Ungauged Basins (PUB) through modeling. Remote sensed leaf area index (LAI), actual evapotranspiration, and soil moisture data were used. Each SRSD was used alone to calibrate a hydrologic model to predict the daily streamflow for 28 basins in Korea. A vegetation module was added to the existing hydrologic model to use LAI. Among the SRSDs tested, the model calibrated with LAI had the most robust performance, predicting streamflow with acceptable accuracy compared to the traditional calibration based on streamflow. In particular, since the model account for vegetation actively interacting with evapotranspiration and soil moisture in the season of low flow, the LAI-calibrated model showed an advantage in improving the flow prediction performance. Although further research is required to utilize evapotranspiration and soil moisture data, the overall results of the LAI-based calibration were promising for predicting streamflow in ungauged basins where observations are scarce or absent, given that the satellite-derived LAI data were used alone without any preprocessing such as a bias correction. However, the prediction performance of the LAI-calibrated model was found to have a statistically significant relationship with local conditions. Therefore, by evaluating and improving the potential of SRSD in different region and climatic conditions, it is expected that the application of the SRSD-only calibration method can be extended to various ungauged basins.

A copula model to identify the risk of river water temperature stress for meteorological drought.

Drought is a natural phenomenon that can occur in all climatic zones, and is persistent and regionally widespread. Extreme drought caused by climate change can have serious consequences for freshwater ecosystems, which can have significant social and economic impacts. In this study, the effect of meteorological drought on river water temperature was analyzed probabilistically in order to identify the risk of river water temperature stress experienced by the aquatic ecosystem when a meteorological drought occurs. Meteorological drought is divided into a situation in which moisture is insufficiently supplied from the atmosphere and a situation in which the atmosphere requires excessive moisture from the earth's surface. Using the copula theory, a joint probabilistic model between the river water temperature and each meteorological drought caused by two causes is proposed. In order to consider the propagation time from meteorological drought to river water temperature, the optimal time-scale meteorological drought index is adopted through correlation analysis between the meteorological drought index calculated at various time-scales and the river water temperature. The optimal copula function of the drought index and river water temperature is determined using AIC analysis. Using the proposed model, a risk map is drawn for the river water temperature stress experienced by the aquatic ecosystem under the user-defined meteorological drought severity. The risk map identifies the stream sections where the river water temperature is relatively more sensitive to meteorological drought. The identified stream sections appear differently depending on the cause of the meteorological drought, the region, and the season.

Chromosome-level genome assembly of the shuttles hoppfish, Periophthalmus modestus.

BACKGROUND: The shuttles hoppfish (mudskipper), Periophthalmus modestus, is one of the mudskippers, which are the largest group of amphibious teleost fishes, which are uniquely adapted to live on mudflats. Because mudskippers can survive on land for extended periods by breathing through their skin and through the lining of the mouth and throat, they were evaluated as a model for the evolutionary sea-land transition of Devonian protoamphibians, ancestors of all present tetrapods. RESULTS: A total of 39.6, 80.2, 52.9, and 33.3 Gb of Illumina, Pacific Biosciences, 10X linked, and Hi-C data, respectively, was assembled into 1,419 scaffolds with an N50 length of 33 Mb and BUSCO score of 96.6%. The assembly covered 117% of the estimated genome size (729 Mb) and included 23 pseudo-chromosomes anchored by a Hi-C contact map, which corresponded to the top 23 longest scaffolds above 20 Mb and close to the estimated one. Of the genome, 43.8% were various repetitive elements such as DNAs, tandem repeats, long interspersed nuclear elements, and simple repeats. Ab initio and homology-based gene prediction identified 30,505 genes, of which 94% had homology to the 14 Actinopterygii transcriptomes and 89% and 85% to Pfam familes and InterPro domains, respectively. Comparative genomics with 15 Actinopterygii species identified 59,448 gene families of which 12% were only in P. modestus. CONCLUSIONS: We present the high quality of the first genome assembly and gene annotation of the shuttles hoppfish. It will provide a valuable resource for further studies on sea-land transition, bimodal respiration, nitrogen excretion, osmoregulation, thermoregulation, vision, and mechanoreception.

Copula-based Joint Drought Index using SPI and EDDI and its application to climate change.

The drought index, which mainly focuses on the moisture supply side of the atmosphere, which has been mainly used in the field of drought monitoring, has limitations that cannot reflect drought caused by changes in various climate variables such as an increase in surface air temperature due to global warming. To overcome these limitations, various evaporation demand-based drought indices have been proposed, focusing on the aspect of atmospheric moisture demand. However, drought indices that consider only precipitation or the demand for atmospheric evaporation are difficult to comprehensively interpret drought caused by various climatic factors. The novelty of this study is to propose a new drought index to simultaneously monitor droughts occurring in terms of atmospheric moisture supply and demand. The proposed Copula-based Joint Drought Index (CJDI) combines the Standardized Precipitation Index and the Evaporative Demand Drought Index using copula. Since CJDI reflects the correlation between the two drought indices, it is shown that CJDI can better monitor Korea's past droughts than other drought indices. It is found that quantification of past drought using CJDI can be used to objectively recognize the level of drought currently in progress by combining with drought severity-duration-frequency curves derived from partial duration series. As a result of analyzing the future drought pattern in Korea, it was revealed that the drought would be alleviated by about 11% in the case of SPI and SPEI, but the drought would intensify by about 89% in the case of EDDI. In the case of CJDI, it is projected that the drought is likely to intensify to about 17%. From the perspective of better reproducing past droughts and projecting a more convincing future drought than other drought indices, CJDI is expected to be fully utilized as a drought index to monitor droughts and establish climate change adaptation policies.

First draft genome for the sand-hopper Trinorchestia longiramus.

Crustacean amphipods are important trophic links between primary producers and higher consumers. Although most amphipods occur in or around aquatic environments, the family Talitridae is the only family found in terrestrial and semi-terrestrial habitats. The sand-hopper Trinorchestia longiramus is a talitrid species often found in the sandy beaches of South Korea. In this study, we present the first draft genome assembly and annotation of this species. We generated ~380.3 Gb of sequencing data assembled in a 0.89 Gb draft genome. Annotation analysis estimated 26,080 protein-coding genes, with 89.9% genome completeness. Comparison with other amphipods showed that T. longiramus has 327 unique orthologous gene clusters, many of which are expanded gene families responsible for cellular transport of toxic substances, homeostatic processes, and ionic and osmotic stress tolerance. This first talitrid genome will be useful for further understanding the mechanisms of adaptation in terrestrial environments, the effects of heavy metal toxicity, as well as for studies of comparative genomic variation across amphipods.

MYC deregulates TET1 and TET2 expression to control global DNA (hydroxy)methylation and gene expression to maintain a neoplastic phenotype in T-ALL.

BACKGROUND: While aberrant DNA methylation is a characteristic feature of tumor cells, our knowledge of how these DNA methylation patterns are established and maintained is limited. DNA methyltransferases and ten-eleven translocation methylcytosine dioxygenases (TETs) function has been found altered in a variety of cancer types. RESULTS: Here, we report that in T cell acute lymphoblastic leukemia (T-ALL) the MYC oncogene controls the expression of TET1 and TET2 to maintain 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) patterns, which is associated with tumor cell-specific gene expression. We found that cellular senescence and tumor regression upon MYC inactivation in T-ALL was associated with genome-wide changes in 5mC and 5hmC patterns. Correlating with the changes in DNA (hydroxy)methylation, we found that T-ALL overexpress TET1, while suppressing TET2 in a MYC-dependent fashion. Consequently, MYC inactivation led to an inverse expression pattern, decreasing TET1, while increasing TET2 levels. Knockdown of TET1 or ectopic expression of TET2 in T-ALL was associated with genome-wide changes in 5mC and 5hmC enrichment and decreased cell proliferation, suggesting a tumor promoting function of TET1, and a tumor suppressing role for TET2. Among the genes and pathways controlled by TET1, we found ribosomal biogenesis and translational control of protein synthesis highly enriched. CONCLUSIONS: Our finding that MYC directly deregulates the expression of TET1 and TET2 in T-ALL provides novel evidence that MYC controls DNA (hydroxy)methylation in a genome-wide fashion. It reveals a coordinated interplay between the components of the DNA (de)methylating machinery that contribute to MYC-driven tumor maintenance, highlighting the potential of specific TET enzymes for therapeutic strategies.

Metabolic Reprogramming by MYCN Confers Dependence on the Serine-Glycine-One-Carbon Biosynthetic Pathway.

MYCN amplification drives the development of neuronal cancers in children and adults. Given the challenge in therapeutically targeting MYCN directly, we searched for MYCN-activated metabolic pathways as potential drug targets. Here we report that neuroblastoma cells with MYCN amplification show increased transcriptional activation of the serine-glycine-one-carbon (SGOC) biosynthetic pathway and an increased dependence on this pathway for supplying glucose-derived carbon for serine and glycine synthesis. Small molecule inhibitors that block this metabolic pathway exhibit selective cytotoxicity to MYCN-amplified cell lines and xenografts by inducing metabolic stress and autophagy. Transcriptional activation of the SGOC pathway in MYCN-amplified cells requires both MYCN and ATF4, which form a positive feedback loop, with MYCN activation of ATF4 mRNA expression and ATF4 stabilization of MYCN protein by antagonizing FBXW7-mediated MYCN ubiquitination. Collectively, these findings suggest a coupled relationship between metabolic reprogramming and increased sensitivity to metabolic stress, which could be exploited as a strategy for selective cancer therapy. SIGNIFICANCE: This study identifies a MYCN-dependent metabolic vulnerability and suggests a coupled relationship between metabolic reprogramming and increased sensitivity to metabolic stress, which could be exploited for cancer therapy.See related commentary by Rodriguez Garcia and Arsenian-Henriksson, p. 3818.

Discovery and Validation of a Novel Neutrophil Activation Marker Associated with Obesity.

Obesity is accompanied by low-grade systemic inflammation that etiologically contributes to obesity-induced cardiovascular disease (CVD). Growing evidence supports that neutrophil, the most abundant type of leukocytes in human, is most likely to be the target peripheral leukocyte subtype initiating the inflammatory cascade in obesity. However, few studies have systematically assessed the genome wide changes in neutrophils associated with obesity. In this study, a hypothesis-free OMIC approach (i.e. the discovery phase) and a target approach (i.e. the validation phase) were used to identify obesity related neutrophil activation markers and their roles on CVD risks. In the discovery phase, genome wide DNA methylation, RNA-sequencing and quantitative proteomics were obtained from purified neutrophils (12 obese vs. 12 lean). In the validation phase, gene expression levels of the promising genes from the OMIC platforms were measured in 81 obese cases vs. 83 lean controls, and the association between the expression levels and CVD risks were evaluated. Significant difference was found for one gene, alkaline phosphatase, liver/bone/kidney (ALPL), across 3 OMIC platforms. In the validation phase, the gene expression levels of ALPL in leukocytes were significantly higher in obese compared with lean subjects (p < 0.05). Within the obese population, we observed that ALPL expression level showed significantly positive association with CVD risk factors (p < 0.05) including systolic blood pressure, diastolic blood pressure, mean arterial pressure, carotid intima-media thickness and borderline significance with fasting insulin (p = 0.08). This study identified one novel marker ALPL of neutrophil activation in response to obesity and provided evidence that obesity induced change in ALPL expression was associated with CVD risk factors.

The complete mitochondrial genome of the sand-hopper Trinorchestia longiramus (Amphipoda: Talitridae).

The complete mitochondrial genome of sand-hopper Trinorchestia longiramus was analyzed in this study, which is the first for the genus within the family Talitridae. The mitogenome sequence is 15,401 bp in length containing two ribosomal RNA genes, 22 transfer RNA genes, 13 protein-coding genes, and a control region as found in most amphipods. The gene order showed that T. longiramus has a unique control region location compared to other amphipods. Phylogenetic analysis using the maximum likelihood method positioned T. longiramus within the monophyletic clades of the family Talitridae.

Hydroxyurea differentially modulates activator and repressors of γ-globin gene in erythroblasts of responsive and non-responsive patients with sickle cell disease in correlation with Index of Hydroxyurea Responsiveness.

Hydroxyurea (HU), the first of two drugs approved by the US Food and Drug Administration for treating patients with sickle cell disease (SCD), produces anti-sickling effect by re-activating fetal γ-globin gene to enhance production of fetal hemoglobin. However, approximately 30% of the patients do not respond to HU therapy. The molecular basis of non-responsiveness to HU is not clearly understood. To address this question, we examined HU-induced changes in the RNA and protein levels of transcription factors NF-Y, GATA-1, -2, BCL11A, TR4, MYB and NF-E4 that assemble the γ-globin promoter complex and regulate transcription of γ-globin gene. In erythroblasts cultured from peripheral blood CD34+ cells of patients with SCD, we found that HU-induced changes in the protein but not the RNA levels of activator GATA-2 and repressors GATA-1, BCL11A and TR4 correlated with HU-induced changes in fetal hemoglobin (HbF) levels in the peripheral blood of HU high and low responders. However, HU did not significantly induce changes in the protein or RNA levels of activators NF-Y and NF-E4. Based on HU-induced changes in the protein levels of GATA-2, -1 and BCL11A, we calculated an Index of Hydroxyurea Responsiveness (IndexHU-3). Compared to the HU-induced fold changes in the individual transcription factor protein levels, the numerical values of IndexHU-3 statistically correlated best with the HU-induced peripheral blood HbF levels of the patients. Thus, IndexHU-3 can serve as an appropriate indicator for inherent HU responsiveness of patients with SCD.

Promoter Methylation Modulates Indoleamine 2,3-Dioxygenase 1 Induction by Activated T Cells in Human Breast Cancers.

Triple-negative breast cancer (TNBC) cells are modulated in reaction to tumor-infiltrating lymphocytes. However, their specific responses to this immune pressure are unknown. In order to address this question, we first used mRNA sequencing to compare the immunophenotype of the TNBC cell line MDA-MB-231 and the luminal breast cancer cell line MCF7 after both were cocultured with activated human T cells. Despite similarities in the cytokine-induced immune signatures of the two cell lines, MDA-MD-231 cells were able to transcribe more IDO1 than MCF7 cells. The two cell lines had similar upstream JAK/STAT1 signaling and IDO1 mRNA stability. However, using a series of breast cancer cell lines, IFNγ stimulated IDO1 protein expression and enzymatic activity only in ER-, not ER+, cell lines. Treatment with 5-aza-deoxycytidine reversed the suppression of IDO1 expression in MCF7 cells, suggesting that DNA methylation was potentially involved in IDO1 induction. By analyzing several breast cancer datasets, we discovered subtype-specific mRNA and promoter methylation differences in IDO1, with TNBC/basal subtypes exhibiting lower methylation/higher expression and ER+/luminal subtypes exhibiting higher methylation/lower expression. We confirmed this trend of IDO1 methylation by bisulfite pyrosequencing breast cancer cell lines and an independent cohort of primary breast tumors. Taken together, these findings suggest that IDO1 promoter methylation regulates anti-immune responses in breast cancer subtypes and could be used as a predictive biomarker for IDO1 inhibitor-based immunotherapy. Cancer Immunol Res; 5(4); 330-44. ©2017 AACR.

Long non-coding RNAs transcribed by ERV-9 LTR retrotransposon act in cis to modulate long-range LTR enhancer function.

LTR retrotransposons are repetitive DNA elements comprising ∼10% of the human genome. However, LTR sequences are disproportionately present in human long, non-coding RNAs (lncRNAs). Whether and how the LTR lncRNAs serve biological functions are largely unknown. Here we show that in primary human erythroblasts, lncRNAs transcribed from the LTR retrotransposons of ERV-9 human endogenous retrovirus activated transcription of key erythroid genes and modulated ex vivo erythropoiesis. To dissect the functional mechanism of ERV-9 lncRNAs, we performed genome-wide RNA and ChIRP analyses before and after global knockdown or locus-specific deletion of ERV-9 lncRNAs in human erythroblasts carrying ∼4000 copies of the ERV-9 LTRs and in transgenic mouse erythroblasts carrying a single copy of the primate-specific ERV-9 LTR in the 100 kb human ß-globin gene locus. We found that ERV-9 lncRNAs acted in cis to stabilize assembly of the ERV-9 LTR enhancer complex and facilitate long-range LTR enhancer function in activating transcription of downstream, cis-linked globin genes. Our findings suggested that LTR lncRNAs transcribed from many of the 4000 copies of ERV-9 LTR retrotransposons acted by a similar cis mechanism to modulate LTR enhancer function in activating transcription of downstream genes critical to cellular processes including erythropoiesis.

CAME: identification of chromatin accessibility from nucleosome occupancy and methylome sequencing.

Motivation: Chromatin accessibility plays a key role in epigenetic regulation of gene activation and silencing. Open chromatin regions allow regulatory elements such as transcription factors and polymerases to bind for gene expression while closed chromatin regions prevent the activity of transcriptional machinery. Recently, Methyltransferase Accessibility Protocol for individual templates-Bisulfite Genome Sequencing (MAPit-BGS) and nucleosome occupancy and methylome sequencing (NOMe-seq) have been developed for simultaneously profiling chromatin accessibility and DNA methylation on single molecules. Therefore, there is a great demand in developing computational methods to identify chromatin accessibility from MAPit-BGS and NOMe-seq. Results: In this article, we present CAME (Chromatin Accessibility and Methylation), a seed-extension based approach that identifies chromatin accessibility from NOMe-seq. The efficiency and effectiveness of CAME were demonstrated through comparisons with other existing techniques on both simulated and real data, and the results show that our method not only can precisely identify chromatin accessibility but also outperforms other methods. Availability and Implementation: CAME is implemented in java and the program is freely available online at http://sourceforge.net/projects/came/. Contacts: jechoi@gru.edu or khryu@dblab.chungbuk.ac.kr. Supplementary information: Supplementary data are available at Bioinformatics online.

High fat diet and exercise lead to a disrupted and pathogenic DNA methylome in mouse liver.

High-fat diet consumption and sedentary lifestyle elevates risk for obesity, non-alcoholic fatty liver disease, and cancer. Exercise training conveys health benefits in populations with or without these chronic conditions. Diet and exercise regulate gene expression by mediating epigenetic mechanisms in many tissues; however, such effects are poorly documented in the liver, a central metabolic organ. To dissect the consequences of diet and exercise on the liver epigenome, we measured DNA methylation, using reduced representation bisulfite sequencing, and transcription, using RNA-seq, in mice maintained on a fast food diet with sedentary lifestyle or exercise, compared with control diet with and without exercise. Our analyses reveal that genome-wide differential DNA methylation and expression of gene clusters are induced by diet and/or exercise. A combination of fast food and exercise triggers extensive gene alterations, with enrichment of carbohydrate/lipid metabolic pathways and muscle developmental processes. Through evaluation of putative protective effects of exercise on diet-induced DNA methylation, we show that hypermethylation is effectively prevented, especially at promoters and enhancers, whereas hypomethylation is only partially attenuated. We assessed diet-induced DNA methylation changes associated with liver cancer-related epigenetic modifications and identified significant increases at liver-specific enhancers in fast food groups, suggesting partial loss of liver cell identity. Hypermethylation at a subset of gene promoters was associated with inhibition of tissue development and promotion of carcinogenic processes. Our study demonstrates extensive reprogramming of the epigenome by diet and exercise, emphasizing the functional relevance of epigenetic mechanisms as an interface between lifestyle modifications and phenotypic alterations.

Transcriptional Profiling Reveals a Common Metabolic Program in High-Risk Human Neuroblastoma and Mouse Neuroblastoma Sphere-Forming Cells.

High-risk neuroblastoma remains one of the deadliest childhood cancers. Identification of metabolic pathways that drive or maintain high-risk neuroblastoma may open new avenues of therapeutic interventions. Here, we report the isolation and propagation of neuroblastoma sphere-forming cells with self-renewal and differentiation potential from tumors of the TH-MYCN mouse, an animal model of high-risk neuroblastoma with MYCN amplification. Transcriptional profiling reveals that mouse neuroblastoma sphere-forming cells acquire a metabolic program characterized by transcriptional activation of the cholesterol and serine-glycine synthesis pathways, primarily as a result of increased expression of sterol regulatory element binding factors and Atf4, respectively. This metabolic reprogramming is recapitulated in high-risk human neuroblastomas and is prognostic for poor clinical outcome. Genetic and pharmacological inhibition of the metabolic program markedly decreases the growth and tumorigenicity of both mouse neuroblastoma sphere-forming cells and human neuroblastoma cell lines. These findings suggest a therapeutic strategy for targeting the metabolic program of high-risk neuroblastoma.

CD133+CD24lo defines a 5-Fluorouracil-resistant colon cancer stem cell-like phenotype.

The chemotherapeutic agent 5-Fluorouracil (5-FU) is the most commonly used drug for patients with advanced colon cancer. However, development of resistance to 5-FU is inevitable in almost all patients. The mechanism by which colon cancer develops 5-FU resistance is still unclear. One recently proposed theory is that cancer stem-like cells underlie colon cancer 5-FU resistance, but the phenotypes of 5-FU-resistant colon cancer stem cells are still controversial. We report here that 5-FU treatment selectively enriches a subset of CD133+ colon cancer cells in vitro. 5-FU chemotherapy also increases CD133+ tumor cells in human colon cancer patients. However, sorted CD133+ colon cancer cells exhibit no increased resistance to 5-FU, and CD133 levels exhibit no correlation with colon cancer patient survival or cancer recurrence. Genome-wide analysis of gene expression between sorted CD133+ colon cancer cells and 5-FU-selected colon cancer cells identifies 207 differentially expressed genes. CD24 is one of the genes whose expression level is lower in the CD133+ and 5-FU-resistant colon cancer cells as compared to CD133+ and 5-FU-sensitive colon cancer cells. Consequently, CD133+CD24lo cells exhibit decreased sensitivity to 5-FU. Therefore, we determine that CD133+CD24lo phenotype defines 5-FU-resistant human colon cancer stem cell-like cells.

OGS2: genome re-annotation of the jewel wasp Nasonia vitripennis.

BACKGROUND: Nasonia vitripennis is an emerging insect model system with haplodiploid genetics. It holds a key position within the insect phylogeny for comparative, evolutionary and behavioral genetic studies. The draft genomes for N. vitripennis and two sibling species were published in 2010, yet a considerable amount of transcriptiome data have since been produced thereby enabling improvements to the original (OGS1.2) annotated gene set. We describe and apply the EvidentialGene method used to produce an updated gene set (OGS2). We also carry out comparative analyses showcasing the usefulness of the revised annotated gene set. RESULTS: The revised annotation (OGS2) now consists of 24,388 genes with supporting evidence, compared to 18,850 for OGS1.2. Improvements include the nearly complete annotation of untranslated regions (UTR) for 97 % of the genes compared to 28 % of genes for OGS1.2. The fraction of RNA-Seq validated introns also grow from 85 to 98 % in this latest gene set. The EST and RNA-Seq expression data provide support for several non-protein coding loci and 7712 alternative transcripts for 4146 genes. Notably, we report 180 alternative transcripts for the gene lola. Nasonia now has among the most complete insect gene set; only 27 conserved single copy orthologs in arthropods are missing from OGS2. Its genome also contains 2.1-fold more duplicated genes and 1.4-fold more single copy genes than the Drosophila melanogaster genome. The Nasonia gene count is larger than those of other sequenced hymenopteran species, owing both to improvements in the genome annotation and to unique genes in the wasp lineage. We identify 1008 genes and 171 gene families that deviate significantly from other hymenopterans in their rates of protein evolution and duplication history, respectively. We also provide an analysis of alternative splicing that reveals that genes with no annotated isoforms are characterized by shorter transcripts, fewer introns, faster protein evolution and higher probabilities of duplication than genes having alternative transcripts. CONCLUSIONS: Genome-wide expression data greatly improves the annotation of the N. vitripennis genome, by increasing the gene count, reducing the number of missing genes and providing more comprehensive data on splicing and gene structure. The improved gene set identifies lineage-specific genomic features tied to Nasonia's biology, as well as numerous novel genes. OGS2 and its associated search tools are available at http://arthropods.eugenes.org/EvidentialGene/nasonia/ , www.hymenopteragenome.org/nasonia/ and waspAtlas: www.tinyURL.com/waspAtlas . The EvidentialGene pipeline is available at https://sourceforge.net/projects/evidentialgene/ .

Targeted next generation sequencing approach identifies eighteen new candidate genes in normosmic hypogonadotropic hypogonadism and Kallmann syndrome.

The genetic basis is unknown for ∼60% of normosmic hypogonadotropic hypogonadism (nHH)/Kallmann syndrome (KS). DNAs from (17 male and 31 female) nHH/KS patients were analyzed by targeted next generation sequencing (NGS) of 261 genes involved in hypothalamic, pituitary, and/or olfactory pathways, or suggested by chromosome rearrangements. Selected variants were subjected to Sanger DNA sequencing, the gold standard. The frequency of Sanger-confirmed variants was determined using the ExAC database. Variants were classified as likely pathogenic (frameshift, nonsense, and splice site) or predicted pathogenic (nonsynonymous missense). Two novel FGFR1 mutations were identified, as were 18 new candidate genes including: AMN1, CCKBR, CRY1, CXCR4, FGF13, GAP43, GLI3, JAG1, NOS1, MASTL, NOTCH1, NRP2, PALM2, PDE3A, PLEKHA5, RD3, and TRAPPC9, and TSPAN11. Digenic and trigenic variants were found in 8/48 (16.7%) and 1/48 (2.1%) patients, respectively. NGS with confirmation by Sanger sequencing resulted in the identification of new causative FGFR1 gene mutations and suggested 18 new candidate genes in nHH/KS.

Phenotypic alteration of CD8+ T cells in chronic lymphocytic leukemia is associated with epigenetic reprogramming.

Immunosuppression is a prevalent clinical feature in chronic lymphocytic leukemia (CLL) patients, with many patients demonstrating increased susceptibility to infections as well as increased failure of an antitumor immune response. However, much is currently not understood regarding the precise mechanisms that attribute to this immunosuppressive phenotype in CLL. To provide further clarity to this particular phenomenon, we analyzed the T-cell profile of CLL patient samples within a large cohort and observed that patients with an inverted CD4/CD8 ratio had a shorter time to first treatment as well as overall survival. These observations coincided with higher expression of the immune checkpoint receptor PD-1 in CLL patient CD8+ T cells when compared to age-matched healthy donors. Interestingly, we discovered that increased PD-1 expression in CD8+ T cells corresponds with decreased DNA methylation levels in a distal upstream locus of the PD-1 gene PDCD1. Further analysis using luciferase reporter assays suggests that the identified PDCD1 distal upstream region acts as an enhancer for PDCD1 transcription and this region becomes demethylated during activation of naïve CD8+ T cells by anti-CD3/anti-CD28 antibodies and IL2. Finally, we conducted a genome-wide DNA methylation analysis comparing CD8+ T cells from CLL patients against healthy donors and identified additional differentially methylated genes with known immune regulatory functions including CCR6 and KLRG1. Taken together, our findings reveal the occurrence of epigenetic reprogramming taking place within CLL patient CD8+ T cells and highlight the potential mechanism of how immunosuppression is accomplished in CLL.