Thousands of Qatari genomes inform human migration history and improve imputation of Arab haplotypes.

Arab populations are largely understudied, notably their genetic structure and history. Here we present an in-depth analysis of 6,218 whole genomes from Qatar, revealing extensive diversity as well as genetic ancestries representing the main founding Arab genealogical lineages of Qahtanite (Peninsular Arabs) and Adnanite (General Arabs and West Eurasian Arabs). We find that Peninsular Arabs are the closest relatives of ancient hunter-gatherers and Neolithic farmers from the Levant, and that founder Arab populations experienced multiple splitting events 12-20 kya, consistent with the aridification of Arabia and farming in the Levant, giving rise to settler and nomadic communities. In terms of recent genetic flow, we show that these ancestries contributed significantly to European, South Asian as well as South American populations, likely as a result of Islamic expansion over the past 1400 years. Notably, we characterize a large cohort of men with the ChrY J1a2b haplogroup (n = 1,491), identifying 29 unique sub-haplogroups. Finally, we leverage genotype novelty to build a reference panel of 12,432 haplotypes, demonstrating improved genotype imputation for both rare and common alleles in Arabs and the wider Middle East.

c-Rel employs multiple mechanisms to promote the thymic development and peripheral function of regulatory T cells in mice.

The NF-κB transcription factor c-Rel is a critical regulator of Treg ontogeny, controlling multiple points of the stepwise developmental pathway. Here, we found that the thymic Treg defect in c-Rel-deficient (cRel-/- ) mice is quantitative, not qualitative, based on analyses of TCR repertoire and TCR signaling strength. However, these parameters are altered in the thymic Treg-precursor population, which is also markedly diminished in cRel-/- mice. Moreover, c-Rel governs the transcriptional programme of both thymic and peripheral Tregs, controlling a core of genes involved with immune signaling, and separately in the periphery, cell cycle progression. Last, the immune suppressive function of peripheral cRel-/- tTregs is diminished in a lymphopenic model of T cell proliferation and is associated with decreased stability of Foxp3 expression. Collectively, we show that c-Rel is a transcriptional regulator that controls multiple aspects of Treg development, differentiation, and function via distinct mechanisms.

Reprogramming roadmap reveals route to human induced trophoblast stem cells.

The reprogramming of human somatic cells to primed or naive induced pluripotent stem cells recapitulates the stages of early embryonic development1-6. The molecular mechanism that underpins these reprogramming processes remains largely unexplored, which impedes our understanding and limits rational improvements to reprogramming protocols. Here, to address these issues, we reconstruct molecular reprogramming trajectories of human dermal fibroblasts using single-cell transcriptomics. This revealed that reprogramming into primed and naive pluripotency follows diverging and distinct trajectories. Moreover, genome-wide analyses of accessible chromatin showed key changes in the regulatory elements of core pluripotency genes, and orchestrated global changes in chromatin accessibility over time. Integrated analysis of these datasets revealed a role for transcription factors associated with the trophectoderm lineage, and the existence of a subpopulation of cells that enter a trophectoderm-like state during reprogramming. Furthermore, this trophectoderm-like state could be captured, which enabled the derivation of induced trophoblast stem cells. Induced trophoblast stem cells are molecularly and functionally similar to trophoblast stem cells derived from human blastocysts or first-trimester placentas7. Our results provide a high-resolution roadmap for the transcription-factor-mediated reprogramming of human somatic cells, indicate a role for the trophectoderm-lineage-specific regulatory program during this process, and facilitate the direct reprogramming of somatic cells into induced trophoblast stem cells.

RNA-seq Reveals Dysregulation of Novel Melanocyte Genes upon Oxidative Stress: Implications in Vitiligo Pathogenesis.

Oxidative stress is known to induce melanocyte death, but the underlying mechanisms are incompletely understood. To identify oxidative stress-induced global gene expression changes in melanocytes, we treated PIG1 melanocytes with H2O2 in a dose- and time-dependent manner and performed RNA-seq. This approach allowed us to capture the events occurring early as well as late phase after treatment with H2O2. Our bioinformatics analysis identified differentially expressed genes involved in various biological processes of melanocytes which are known to contribute to the vitiligo development, such as apoptosis, autophagy, cell cycle regulation, cell adhesion, immune and inflammatory responses, melanocyte pluripotency, and developmental signaling such as WNT and NOTCH pathways. We uncovered several novel genes that are not previously described to be involved in melanocytic response to stress nor in vitiligo pathogenesis. Quantitative PCR and western blot analysis of selected proteins, performed on PIG1 and primary human epidermal melanocytes, confirmed the RNA-seq data. Interestingly, we discovered an aberrant regulation of several transcription factors that are involved in diabetes, neurological, and psychiatric diseases, all of which are comorbid conditions in patients with vitiligo. Our results may lead to a better understanding of the molecular mechanisms underlying vitiligo pathogenesis and help developing new drug targets for effective treatment.

Interplay between Follistatin, Activin A, and BMP4 Signaling Regulates Postnatal Thymic Epithelial Progenitor Cell Differentiation during Aging.

A key feature of immune functional impairment with age is the progressive involution of thymic tissue responsible for naive T cell production. In this study, we identify two major phases of thymic epithelial cell (TEC) loss during aging: a block in mature TEC differentiation from the pool of immature precursors, occurring at the onset of puberty, followed by impaired bipotent TEC progenitor differentiation and depletion of Sca-1lo cTEC and mTEC lineage-specific precursors. We reveal that an increase in follistatin production by aging TECs contributes to their own demise. TEC loss occurs primarily through the antagonism of activin A signaling, which we show is required for TEC maturation and acts in dissonance to BMP4, which promotes the maintenance of TEC progenitors. These results support a model in which an imbalance of activin A and BMP4 signaling underpins the degeneration of postnatal TEC maintenance during aging, and its reversal enables the transient replenishment of mature TECs.

Gestational diabetes and maternal obesity are associated with epigenome-wide methylation changes in children.

Offspring of women with gestational diabetes mellitus (GDM) are at increased risk of developing metabolic disease, potentially mediated by epigenetic mechanisms. We recruited 608 GDM and 626 control offspring from the Danish National Birth Cohort, aged between 9 and 16 years. DNA methylation profiles were measured in peripheral blood of 93 GDM offspring and 95 controls using the Illumina HumanMethylation450 BeadChip. Pyrosequencing was performed for validation/replication of putative GDM-associated, differentially methylated CpGs in additional 905 offspring (462 GDM, 444 control offspring). We identified 76 differentially methylated CpGs in GDM offspring compared with controls in the discovery cohort (FDR, P < 0.05). Adjusting for offspring BMI did not affect the association between methylation levels and GDM status for any of the 76 CpGs. Most of these epigenetic changes were due to confounding by maternal prepregnancy BMI; however, 13 methylation changes were independently associated with maternal GDM. Three prepregnancy BMI-associated CpGs (cg00992687 and cg09452568 of ESM1 and cg14328641 of MS4A3) were validated in the replication cohort, while cg09109411 (PDE6A) was found to be associated with GDM status. The identified methylation changes may reflect developmental programming of organ disease mechanisms and/or may serve as disease biomarkers.

Age-Related Decline in Primary CD8+ T Cell Responses Is Associated with the Development of Senescence in Virtual Memory CD8+ T Cells.

Age-associated decreases in primary CD8+ T cell responses occur, in part, due to direct effects on naive CD8+ T cells to reduce intrinsic functionality, but the precise nature of this defect remains undefined. Aging also causes accumulation of antigen-naive but semi-differentiated "virtual memory" (TVM) cells, but their contribution to age-related functional decline is unclear. Here, we show that TVM cells are poorly proliferative in aged mice and humans, despite being highly proliferative in young individuals, while conventional naive T cells (TN cells) retain proliferative capacity in both aged mice and humans. Adoptive transfer experiments in mice illustrated that naive CD8 T cells can acquire a proliferative defect imposed by the aged environment but age-related proliferative dysfunction could not be rescued by a young environment. Molecular analyses demonstrate that aged TVM cells exhibit a profile consistent with senescence, marking an observation of senescence in an antigenically naive T cell population.

The Use of Targeted Next Generation Sequencing to Explore Candidate Regulators of TGF-ß1's Impact on Kidney Cells.

Aims/Hypothesis: Transforming growth factor-beta (TGF-ß1) plays an important regulatory role in the progression of chronic kidney failure. Further, damage to kidney glomerular mesangial cells is central to the progression of diabetic nephropathy. The aim of this study was to explore the genetic associations between mRNA, microRNA, and epigenetics in mesangial cells in response to TGF-ß1. Methods: The regulatory effects of TGF-ß1 on mesangial cells were investigated at different molecular levels by treating mesangial cells with TGF-ß1 for 3 days followed by genome-wide miRNA, RNA, DNA methylation, and H3K27me3 expression profiling using next generation sequencing (NGS). Results: Our results provide the first comprehensive, computationally integrated report of RNA-Seq, miRNA-Seq, and epigenomic analyses across all genetic variations, confirming the occurrence of DNA methylation and H3K27me3 in response to TGF-ß1. Our findings show that the expression of KLF7 and Gja4 are involved in TGF-ß1 regulated DNA methylation. Our data also provide evidence of the association between epigenetic changes and the expression of genes closely related to TGF-ß1 regulation. Conclusion: This study has advanced our current knowledge of mechanisms that contribute to the expression of TGF-ß1-regulated genes involved in the pathogenesis of kidney disease. The molecular underpinnings of TGF-ß1 stimulation of kidney cells was determined, thereby providing a robust platform for further target exploration.

Targeting of Polycomb Repressive Complex 2 to RNA by Short Repeats of Consecutive Guanines.

Polycomb repressive complex 2 (PRC2) is a histone methyltransferase that trimethylates H3K27, a mark of repressed chromatin. Mammalian PRC2 binds RNA promiscuously, with thousands of target transcripts in vivo. But what does PRC2 recognize in these RNAs? Here we show that purified human PRC2 recognizes G > C,U ≫ A in single-stranded RNA and has a high affinity for folded guanine quadruplex (G4) structures but little binding to duplex RNAs. Importantly, G-tract motifs are significantly enriched among PRC2-binding transcripts in vivo. DNA sequences coding for PRC2-binding RNA motifs are enriched at PRC2-binding sites on chromatin and H3K27me3-modified nucleosomes. Collectively, the abundance of PRC2-binding RNA motifs rationalizes the promiscuous RNA binding of PRC2, and their enrichment at Polycomb target genes provides a means for RNA-mediated regulation.

Genome-Wide Measures of Peripheral Blood Dna Methylation and Prostate Cancer Risk in a Prospective Nested Case-Control Study.

BACKGROUND: Global measures of peripheral blood DNA methylation have been associated with risk of some malignancies, including breast, bladder, and gastric cancer. Here, we examined genome-wide measures of peripheral blood DNA methylation in prostate cancer and its non-aggressive and aggressive disease forms. METHODS: We used a matched, case-control study of 687 incident prostate cancer samples, nested within a larger prospective cohort study. DNA methylation was measured in pre-diagnostic, peripheral blood samples using the Illumina Infinium HM450K BeadChip. Genome-wide measures of DNA methylation were computed as the median M-value of all CpG sites and according to CpG site location and regulatory function. We used conditional logistic regression to test for associations between genome-wide measures of DNA methylation and risk of prostate cancer and its subtypes, and by time between blood draw and diagnosis. RESULTS: We observed no associations between the genome-wide measure of DNA methylation based on all CpG sites and risk of prostate cancer or aggressive disease. Risk of non-aggressive disease was associated with higher methylation of CpG islands (OR = 0.80; 95%CI = 0.68-0.94), promoter regions (OR = 0.79; 95%CI = 0.66-0.93), and high density CpG regions (OR = 0.80; 95%CI = 0.68-0.94). Additionally, higher methylation of all CpGs (OR = 0.66; 95%CI = 0.48-0.89), CpG shores (OR = 0.62; 95%CI = 0.45-0.84), and regulatory regions (OR = 0.68; 95% CI = 0.51-0.91) was associated with a reduced risk of overall prostate cancer within 5 years of blood draw but not thereafter. CONCLUSIONS: A reduced risk of overall prostate cancer within 5 years of blood draw and non-aggressive prostate cancer was associated with higher genome-wide methylation of peripheral blood DNA. While these data have no immediate clinical utility, with further work they may provide insight into the early events of prostate carcinogenesis. Prostate 77:471-478, 2017. © 2017 Wiley Periodicals, Inc.

The NF-κB transcription factor RelA is required for the tolerogenic function of Foxp3(+) regulatory T cells.

The properties of CD4(+) regulatory T cell (Treg) subsets are dictated by distinct patterns of gene expression determined by FOXP3 and different combinations of various transcription factors. Here we show the NF-κB transcription factor RelA is constitutively active in naïve and effector Tregs. The conditional inactivation of Rela in murine FOXP3(+) cells induces a rapid onset, multi-focal autoimmune disease that depends on RelA being expressed in conventional T cells. In addition to promoting Treg lineage stability, RelA determines the size of the effector Treg population, a function influenced by the presence or absence of RelA in conventional T cells. These findings showing that RelA controls Treg stability and promotes the competitive fitness of effector Tregs highlight the importance of RelA activity in peripheral Treg induced tolerance.

A urinary microRNA signature can predict the presence of bladder urothelial carcinoma in patients undergoing surveillance.

BACKGROUND: The objective of this study was to determine whether microRNA (miRNA) profiling of urine could identify the presence of urothelial carcinoma of the bladder (UCB) and to compare its performance characteristics to that of cystoscopy. METHODS: In the discovery cohort we screened 81 patients, which included 21 benign controls, 30 non-recurrers and 30 patients with active cancer (recurrers), using a panel of 12 miRNAs. Data analysis was performed using a machine learning approach of a Support Vector Machine classifier with a Student's t-test feature selection procedure. This was trained using a three-fold cross validation approach and performance was measured using the area under the receiver operator characteristic curve (AUC). The miRNA signature was validated in an independent cohort of a further 50 patients. RESULTS: The best predictor to distinguish patients with UCB from non-recurrers was achieved using a combination of six miRNAs (AUC=0.85). This validated in an independent cohort (AUC=0.74) and detected UCB with a high sensitivity (88%) and sufficient specificity (48%) with all significant cancers identified. The performance of the classifier was best in detecting clinically significant disease such as presence of T1 Stage disease (AUC=0.92) and high-volume disease (AUC=0.81). Cystoscopy rates in the validation cohort would have been reduced by 30%. CONCLUSIONS: Urinary profiling using this panel of miRNAs shows promise for detection of tumour recurrence in the surveillance of UCB. Such a panel may be useful in reducing the morbidity and costs associated with cystoscopic surveillance, and now merits prospective evaluation.

Tracking the origins and drivers of subclonal metastatic expansion in prostate cancer.

Tumour heterogeneity in primary prostate cancer is a well-established phenomenon. However, how the subclonal diversity of tumours changes during metastasis and progression to lethality is poorly understood. Here we reveal the precise direction of metastatic spread across four lethal prostate cancer patients using whole-genome and ultra-deep targeted sequencing of longitudinally collected primary and metastatic tumours. We find one case of metastatic spread to the surgical bed causing local recurrence, and another case of cross-metastatic site seeding combining with dynamic remoulding of subclonal mixtures in response to therapy. By ultra-deep sequencing end-stage blood, we detect both metastatic and primary tumour clones, even years after removal of the prostate. Analysis of mutations associated with metastasis reveals an enrichment of TP53 mutations, and additional sequencing of metastases from 19 patients demonstrates that acquisition of TP53 mutations is linked with the expansion of subclones with metastatic potential which we can detect in the blood.

Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array.

BACKGROUND: The Illumina HumanMethylation450 BeadChip (HM450K) measures the DNA methylation of 485,512 CpGs in the human genome. The technology relies on hybridization of genomic fragments to probes on the chip. However, certain genomic factors may compromise the ability to measure methylation using the array such as single nucleotide polymorphisms (SNPs), small insertions and deletions (INDELs), repetitive DNA, and regions with reduced genomic complexity. Currently, there is no clear method or pipeline for determining which of the probes on the HM450K bead array should be retained for subsequent analysis in light of these issues. RESULTS: We comprehensively assessed the effects of SNPs, INDELs, repeats and bisulfite induced reduced genomic complexity by comparing HM450K bead array results with whole genome bisulfite sequencing. We determined which CpG probes provided accurate or noisy signals. From this, we derived a set of high-quality probes that provide unadulterated measurements of DNA methylation. CONCLUSIONS: Our method significantly reduces the risk of false discoveries when using the HM450K bead array, while maximising the power of the array to detect methylation status genome-wide. Additionally, we demonstrate the utility of our method through extraction of biologically relevant epigenetic changes in prostate cancer.

Rigorous assessment of gene set enrichment tests.

MOTIVATION: Several statistical tests are available to detect the enrichment of differential expression in gene sets. Such tests were originally proposed for analyzing gene sets associated with biological processes. The objective evaluation of tests on real measurements has not been possible as it is difficult to decide a priori, which processes will be affected in given experiments. RESULTS: We present a first large study to rigorously assess and compare the performance of gene set enrichment tests on real expression measurements. Gene sets are defined based on the targets of given regulators such as transcription factors (TFs) and microRNAs (miRNAs). In contrast to processes, TFs and miRNAs are amenable to direct perturbations, e.g. regulator over-expression or deletion. We assess the ability of 14 different statistical tests to predict the perturbations from expression measurements in Escherichia coli, Saccharomyces cerevisiae and human. We also analyze how performance depends on the quality and comprehensiveness of the regulator targets via a permutation approach. We find that ANOVA and Wilcoxons test consistently perform better than for instance Kolmogorov-Smirnov and hypergeometric tests. For scenarios where the optimal test is not known, we suggest to combine all evaluated tests into an unweighted consensus, which also performs well in our assessment. Our results provide a guide for the selection of existing tests as well as a basis for the development and assessment of novel tests.

MIRTFnet: analysis of miRNA regulated transcription factors.

BACKGROUND: Several expression datasets of miRNA transfection experiments are available to analyze the regulatory mechanisms downstream of miRNA effects. The miRNA induced regulatory effects can be propagated via transcription factors (TFs). We propose the method MIRTFnet to identify miRNA controlled TFs as active regulators if their downstream target genes are differentially expressed. METHODOLOGY/PRINCIPAL FINDINGS: MIRTFnet enables the determination of active transcription factors (TFs) and is sensitive enough to exploit the small expression changes induced by the activity of miRNAs. For this purpose, different statistical tests were evaluated and compared. Based on the identified TFs, databases, computational predictions and the literature we construct regulatory models downstream of miRNA actions. Transfecting miRNAs are connected to active regulators via a network of miRNA-TF, miRNA-kinase-TF as well as TF-TF relationships. Based on 43 transfection experiments involving 17 cancer relevant miRNAs we show that MIRTFnet detects active regulators reliably. CONCLUSIONS/SIGNIFICANCE: The consensus of the individual regulatory models shows that the examined miRNAs induce activity changes in a common core of transcription factors involved in cancer related processes such as proliferation or apoptosis.

miRSel: automated extraction of associations between microRNAs and genes from the biomedical literature.

BACKGROUND: MicroRNAs have been discovered as important regulators of gene expression. To identify the target genes of microRNAs, several databases and prediction algorithms have been developed. Only few experimentally confirmed microRNA targets are available in databases. Many of the microRNA targets stored in databases were derived from large-scale experiments that are considered not very reliable. We propose to use text mining of publication abstracts for extracting microRNA-gene associations including microRNA-target relations to complement current repositories. RESULTS: The microRNA-gene association database miRSel combines text-mining results with existing databases and computational predictions. Text mining enables the reliable extraction of microRNA, gene and protein occurrences as well as their relationships from texts. Thereby, we increased the number of human, mouse and rat miRNA-gene associations by at least three-fold as compared to e.g. TarBase, a resource for miRNA-gene associations. CONCLUSIONS: Our database miRSel offers the currently largest collection of literature derived miRNA-gene associations. Comprehensive collections of miRNA-gene associations are important for the development of miRNA target prediction tools and the analysis of regulatory networks. miRSel is updated daily and can be queried using a web-based interface via microRNA identifiers, gene and protein names, PubMed queries as well as gene ontology (GO) terms. miRSel is freely available online at http://services.bio.ifi.lmu.de/mirsel.