HIF1α-AS1 is a DNA:DNA:RNA triplex-forming lncRNA interacting with the HUSH complex.

DNA:DNA:RNA triplexes that are formed through Hoogsteen base-pairing of the RNA in the major groove of the DNA duplex have been observed in vitro, but the extent to which these interactions occur in cells and how they impact cellular functions remains elusive. Using a combination of bioinformatic techniques, RNA/DNA pulldown and biophysical studies, we set out to identify functionally important DNA:DNA:RNA triplex-forming long non-coding RNAs (lncRNA) in human endothelial cells. The lncRNA HIF1α-AS1 was retrieved as a top hit. Endogenous HIF1α-AS1 reduces the expression of numerous genes, including EPH Receptor A2 and Adrenomedullin through DNA:DNA:RNA triplex formation by acting as an adapter for the repressive human silencing hub complex (HUSH). Moreover, the oxygen-sensitive HIF1α-AS1 is down-regulated in pulmonary hypertension and loss-of-function approaches not only result in gene de-repression but also enhance angiogenic capacity. As exemplified here with HIF1α-AS1, DNA:DNA:RNA triplex formation is a functionally important mechanism of trans-acting gene expression control.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing.

BACKGROUND: Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 µg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits. RESULTS: We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits' fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection. CONCLUSIONS: The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Whole-exome sequencing in syndromic craniosynostosis increases diagnostic yield and identifies candidate genes in osteogenic signaling pathways.

Craniosynostosis (CS) is a common congenital anomaly defined by premature fusion of one or more cranial sutures. Syndromic CS involves additional organ anomalies or neurocognitive deficits and accounts for 25%-30% of the cases. In a recent population-based study by our group, 84% of the syndromic CS cases had a genetically verified diagnosis after targeted analyses. A number of different genetic causes were detected, confirming that syndromic CS is highly heterogeneous. In this study, we performed whole-exome sequencing of 10 children and parents from the same cohort where previous genetic results were negative. We detected pathogenic, or likely pathogenic, variants in four additional genes (NFIA, EXTL3, POLR2A, and FOXP2) associated with rare conditions. In two of these (POLR2A and FOXP2), CS has not previously been reported. We further detected a rare predicted damaging variant in SH3BP4, which has not previously been related to human disease. All findings were clustered in genes involved in the pathways of osteogenesis and suture patency. We conclude that whole-exome sequencing expands the list of genes associated with syndromic CS, and provides new candidate genes in osteogenic signaling pathways.

Multimodal human thymic profiling reveals trajectories and cellular milieu for T agonist selection.

To prevent autoimmunity, thymocytes expressing self-reactive T cell receptors (TCRs) are negatively selected, however, divergence into tolerogenic, agonist selected lineages represent an alternative fate. As thymocyte development, selection, and lineage choices are dependent on spatial context and cell-to-cell interactions, we have performed Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) and spatial transcriptomics on paediatric human thymu​​s. Thymocytes expressing markers of strong TCR signalling diverged from the conventional developmental trajectory prior to CD4+ or CD8+ lineage commitment, while markers of different agonist selected T cell populations (CD8αα(I), CD8αα(II), T(agonist), Treg(diff), and Treg) exhibited variable timing of induction. Expression profiles of chemokines and co-stimulatory molecules, together with spatial localisation, supported that dendritic cells, B cells, and stromal cells contribute to agonist selection, with different subsets influencing thymocytes at specific developmental stages within distinct spatial niches. Understanding factors influencing agonist T cells is needed to benefit from their immunoregulatory effects in clinical use.

Rapid SARS-CoV-2 variant monitoring using PCR confirmed by whole genome sequencing in a high-volume diagnostic laboratory.

OBJECTIVES: The emerging SARS-CoV-2 variants of concern (VoC), B.1.1.7, B.1.351 and P.1, with increased transmission and/or immune evasion, emphasise the need for broad and rapid variant monitoring. Our high-volume laboratory introduced a PCR variant assay (Variant PCR) in January 2021 based on the protocol by Vogels et al. STUDY DESIGN: To assess whether Variant PCR could be used for rapid B.1.1.7, B.1.351 and P.1 screening, all positive SARS-CoV-2 airway samples were prospectively tested in parallel using both the Variant PCR and whole genome sequencing (WGS). RESULTS: In total 1,642 SARS-CoV-2 positive samples from individual patients were tested within a time span of 4 weeks. For all samples with valid results from both Variant PCR and WGS, no VoC was missed by Variant PCR (totalling 399 VoC detected). Conversely, all of the samples identified as "other lineages" (i.e., "non-VoC lineages") by the Variant PCR, were confirmed by WGS. CONCLUSIONS: The Variant PCR based on the protocol by Vogels et al., is an effective method for rapid screening for VoC, applicable for most diagnostic laboratories within a pandemic setting. WGS is still required to confirm the identity of certain variants and for continuous surveillance of emerging VoC.

Isolation and genome-wide characterization of cellular DNA:RNA triplex structures.

RNA can directly bind to purine-rich DNA via Hoogsteen base pairing, forming a DNA:RNA triple helical structure that anchors the RNA to specific sequences and allows guiding of transcription regulators to distinct genomic loci. To unravel the prevalence of DNA:RNA triplexes in living cells, we have established a fast and cost-effective method that allows genome-wide mapping of DNA:RNA triplex interactions. In contrast to previous approaches applied for the identification of chromatin-associated RNAs, this method uses protein-free nucleic acids isolated from chromatin. High-throughput sequencing and computational analysis of DNA-associated RNA revealed a large set of RNAs which originate from non-coding and coding loci, including super-enhancers and repeat elements. Combined analysis of DNA-associated RNA and RNA-associated DNA identified genomic DNA:RNA triplex structures. The results suggest that triplex formation is a general mechanism of RNA-mediated target-site recognition, which has major impact on biological functions.

A histone acetylome-wide association study of Alzheimer's disease identifies disease-associated H3K27ac differences in the entorhinal cortex.

We quantified genome-wide patterns of lysine H3K27 acetylation (H3K27ac) in entorhinal cortex samples from Alzheimer's disease (AD) cases and matched controls using chromatin immunoprecipitation and highly parallel sequencing. We observed widespread acetylomic variation associated with AD neuropathology, identifying 4,162 differential peaks (false discovery rate < 0.05) between AD cases and controls. Differentially acetylated peaks were enriched in disease-related biological pathways and included regions annotated to genes involved in the progression of amyloid-ß and tau pathology (for example, APP, PSEN1, PSEN2, and MAPT), as well as regions containing variants associated with sporadic late-onset AD. Partitioned heritability analysis highlighted a highly significant enrichment of AD risk variants in entorhinal cortex H3K27ac peak regions. AD-associated variable H3K27ac was associated with transcriptional variation at proximal genes including CR1, GPR22, KMO, PIM3, PSEN1, and RGCC. In addition to identifying molecular pathways associated with AD neuropathology, we present a framework for genome-wide studies of histone modifications in complex disease.

Native Chromatin Immunoprecipitation-Sequencing (ChIP-Seq) from Low Cell Numbers.

ChIP-seq is the current method of choice for genome-wide protein location analysis. Here, we present a native (non-cross-linked) ChIP procedure suitable for histone proteins, coupled with an efficient library preparation technique for subsequent next-generation sequencing. The method enables ChIP-seq starting with 50,000 or more cells.

Deregulation of an imprinted gene network in prostate cancer.

Multiple epigenetic alterations contribute to prostate cancer progression by deregulating gene expression. Epigenetic mechanisms, especially differential DNA methylation at imprinting control regions (termed DMRs), normally ensure the exclusive expression of imprinted genes from one specific parental allele. We therefore wondered to which extent imprinted genes become deregulated in prostate cancer and, if so, whether deregulation is due to altered DNA methylation at DMRs. Therefore, we selected presumptive deregulated imprinted genes from a previously conducted in silico analysis and from the literature and analyzed their expression in prostate cancer tissues by qRT-PCR. We found significantly diminished expression of PLAGL1/ZAC1, MEG3, NDN, CDKN1C, IGF2, and H19, while LIT1 was significantly overexpressed. The PPP1R9A gene, which is imprinted in selected tissues only, was strongly overexpressed, but was expressed biallelically in benign and cancerous prostatic tissues. Expression of many of these genes was strongly correlated, suggesting co-regulation, as in an imprinted gene network (IGN) reported in mice. Deregulation of the network genes also correlated with EZH2 and HOXC6 overexpression. Pyrosequencing analysis of all relevant DMRs revealed generally stable DNA methylation between benign and cancerous prostatic tissues, but frequent hypo- and hyper-methylation was observed at the H19 DMR in both benign and cancerous tissues. Re-expression of the ZAC1 transcription factor induced H19, CDKN1C and IGF2, supporting its function as a nodal regulator of the IGN. Our results indicate that a group of imprinted genes are coordinately deregulated in prostate cancers, independently of DNA methylation changes.

Concomitant downregulation of the imprinted genes DLK1 and MEG3 at 14q32.2 by epigenetic mechanisms in urothelial carcinoma.

BACKGROUND: The two oppositely imprinted and expressed genes, DLK1 and MEG3, are located in the same gene cluster at 14q32. Previous studies in bladder cancer have suggested that tumor suppressor genes are located in this region, but these have not been identified. RESULTS: We observed that both DLK1 and MEG3 are frequently silenced in urothelial cancer tissues and cell lines. The concomitant downregulation of the two genes is difficult to explain by known mechanisms for inactivating imprinted genes, namely deletion of active alleles or epitype switching. Indeed, quantitative PCR revealed more frequent copy number gains than losses in the gene cluster that were, moreover, consistent within each sample, excluding gene losses as the cause of downregulation. Instead, we observed distinctive epigenetic alterations at the three regions controlling DLK1 and MEG3 expression, namely the DLK1 promoter; the intergenic (IG) and MEG3 differentially methylated regions (DMRs). Bisulfite sequencing and pyrosequencing revealed novel patterns of DNA methylation in tumor cells, which were distinct from that of either paternal allele. Furthermore, chromatin immunoprecipitation demonstrated loss of active and gain of repressive histone modifications at all regulatory sequences. CONCLUSIONS: Our data support the idea that the main cause of the prevalent downregulation of DLK1 and MEG3 in urothelial carcinoma is epigenetic silencing across the 14q32 imprinted gene cluster, resulting in the unusual concomitant inactivation of oppositely expressed and imprinted genes.

The SNP rs6441224 influences transcriptional activity and prognostically relevant hypermethylation of RARRES1 in prostate cancer.

Epigenetic aberrations are frequent in prostate cancer and could be useful for detection and prognostication. However, the underlying mechanisms and the sequence of these changes remain to be fully elucidated. The tumor suppressor gene RARRES1 (TIG1) is frequently hypermethylated in several cancers. Having noted changes in the expression of its paralogous neighbor gene LXN at 3q25.32, we used pyrosequencing to quantify DNA methylation at both genes and determine its relationship with clinicopathological parameters in 86 prostate cancer tissues from radical prostatectomies. Methylation at LXN and RARRES1 was highly correlated. Increasing methylation was associated with worse clinical features, including biochemical recurrence, and decreased expression of both genes. However, expression of three neighboring genes was unaffected. Intriguingly, RARRES1 methylation was influenced by the genotype of the rs6441224 single-nucleotide polymorphism (SNP) in its promoter. We found that this SNP is located within an ETS-family-response element and that the more strongly methylated allele confers lower activity in reporter assays. Concomitant methylation of RARRES1 and LXN in cancerous tissues was also detected in prostate cancer cell lines and was shown to be associated with repressive histone modifications and transcriptional downregulation. In conclusion, we found that genotype-associated hypermethylation of the ETS-family target gene RARRES1 influences methylation at its neighbor gene LXN and could be useful as a prognostic biomarker.

Specific changes in the expression of imprinted genes in prostate cancer--implications for cancer progression and epigenetic regulation.

Epigenetic dysregulation comprising DNA hypermethylation and hypomethylation, enhancer of zeste homologue 2 (EZH2) overexpression and altered patterns of histone modifications is associated with the progression of prostate cancer. DNA methylation, EZH2 and histone modifications also ensure the parental-specific monoallelic expression of at least 62 imprinted genes. Although it is therefore tempting to speculate that epigenetic dysregulation may extend to imprinted genes, expression changes in cancerous prostates are only well documented for insulin-like growth factor 2 (IGF2). A literature and database survey on imprinted genes in prostate cancer suggests that the expression of most imprinted genes remains unchanged despite global disturbances in epigenetic mechanisms. Instead, selective genetic and epigenetic changes appear to lead to the inactivation of a sub-network of imprinted genes, which might function in the prostate to limit cell growth induced via the PI3K/Akt pathway, modulate androgen responses and regulate differentiation. Whereas dysregulation of IGF2 may constitute an early change in prostate carcinogenesis, inactivation of this imprinted gene network is rather associated with cancer progression.

Selective changes of retroelement expression in human prostate cancer.

Retroelements constitute a large part of the human genome. These sequences are mostly silenced in normal cells, but genome-wide DNA hypomethylation in cancers might lead to their re-expression. Whether this re-expression really occurs in human cancers is largely unkown. We therefore investigated expression and DNA methylation of several classes of retroelements in human prostate cancer tissues and cell lines by quantitative reverse transcription-polymerase chain reaction and pyrosequencing, respectively. The most striking finding was strong and generalized increased expression of the HERV-K_22q11.23 provirus in cancers, including de novo expression of a spliced accessory Np9 transcript in some tumors. In parallel, DNA methylation in the long terminal repeat (LTR) decreased. Conversely, HERVK17 expression was significantly diminished in cancer tissues, but this decrease was unrelated to LTR methylation. Expression of both proviruses was restricted to androgen-responsive prostate cancer cell lines and LTRs sequences containing steroid hormone-responsive elements bound the androgen receptor and conferred androgen responsiveness to reporter constructs. Expression of LINE-1 5'-untranslated region (UTR) and 3'-UTR sequences in prostate cancers rather decreased, despite significant hypomethylation of the internal LINE-1 promoter. Increased expression of the young AluYa5 and AluYb8 families was restricted to individual tumors. Our findings demonstrate a surprising specificity of changes in expression and DNA methylation of retroelements in prostate cancer. In particular, LINE-1 hypomethylation does not lead to generalized overexpression, but specific human endogenous retrovirus-K proviruses display conspicuous changes in their expression hinting at significant functions during prostate carcinogenesis.

Insights into cancer mechanisms from genomic research on urological cancers.

Molecular mechanisms driving cancer development and progression are rarely unique to one cancer type. Rather, recent genomic studies of urological cancers suggest that common mechanisms recur with variations. Examples include alterations in hypoxia response regulation, epigenetic regulator proteins, and signal transduction pathways in renal, prostatic and urothelial carcinomas. Consideration of these variations alongside the common basic cancer mechanisms might be important for the successful development of targeted therapies.

Epigenetic inactivation of the placentally imprinted tumor suppressor gene TFPI2 in prostate carcinoma.

BACKGROUND: Imprinted genes are often arranged in clusters epigenetically controlled by differentially methylated regions (DMR) containing bivalent histone modifications. Both DNA hypermethylation and hypomethylation in cancer can therefore disturb imprinted gene expression. We have studied expression, DNA methylation and histone modifications of TFPI2, a presumed tumor suppressor, and that of other genes in the 7q21 imprinted gene cluster in prostate cancer. MATERIALS AND METHODS: TFPI, TFPI2, SGCE and PON2 expression were assessed by qRT-PCR in prostate cancer tissues and cell lines. DNA methylation and histone modifications were investigated by bisulphite sequencing and chromatin immunoprecipatation. RESULTS: TFPI2 was highly variably expressed in cancer tissues, in contrast to TFPI, and did not correlate to unchanged SGCE and significantly elevated PON2 expression. TFPI2 expression variations were unrelated to global DNA hypomethylation, but were associated with promoter methylation. PC3 cells with high expression retained normal methylation and bivalent histone modifications at DMR and promoter, whereas low-expressing LNCaP cells presented aberrant DNA methylation and more repressive histone modifications. CONCLUSION: Epigenetic disturbances in the 7q21 cluster affect imprinted genes in a non-coordinate manner suggesting an unstable epigenetic state prone to selection for specific expression changes.

Systemic reduction of soluble complement receptor II/CD21 during pregnancy to levels reminiscent of autoimmune disease.

Complement receptor type II/CD21 is the functional receptor for complement fragments such as C3d, iC3b and the Epstein Barr Virus. A soluble form of CD21 (sCD21) is shed from lymphocytes surface and is able to bind to its ligands found in the plasma. The amount of sCD21 in serum may modulate immunity as the plasma levels are correlated with autoimmune conditions, such as Systemic Lupus Erythematosus, Rheumatoid Arthritis and Sjoegren's Syndrome. Because of the fact that pregnancy may lead to remission of autoimmune diseases we determined the serum levels of sCD21 during pregnancy and postpartum. The serum sCD21 levels during pregnancy are significantly lower as compared to that of the healthy controls. There were no significant differences in sCD21 levels between the mother and the cord blood also immediately after parturition. Restoration of sCD21 levels to normal values takes between 6 weeks and 1 year after childbirth. Our study indicates that CD21-shedding is affected during pregnancy comparable to that of autoimmunity.