DeepNeo: a webserver for predicting immunogenic neoantigens.

Non-self epitopes, whether originated from foreign substances or somatic mutations, trigger immune responses when presented by major histocompatibility complex (MHC) molecules and recognized by T cells. Identification of immunogenically active neoepitopes has significant implications in cancer and virus medicine. However, current methods are mostly limited to predicting physical binding of mutant peptides and MHCs. We previously developed a deep-learning based model, DeepNeo, to identify immunogenic neoepitopes by capturing the structural properties of peptide-MHC pairs with T cell reactivity. Here, we upgraded our DeepNeo model with up-to-date training data. The upgraded model (DeepNeo-v2) was improved in evaluation metrics and showed prediction score distribution that better fits known neoantigen behavior. The immunogenic neoantigen prediction can be conducted at https://deepneo.net.

Altered tumor signature and T-cell profile after chemotherapy reveal new therapeutic opportunities in high-grade serous ovarian carcinoma.

Chemotherapy combined with debulking surgery is the standard treatment protocol for high-grade serous ovarian carcinoma (HGSOC). Nonetheless, a significant number of patients encounter relapse due to the development of chemotherapy resistance. To better understand and address this resistance, we conducted a comprehensive study investigating the transcriptional alterations at the single-cell resolution in tissue samples from patients with HGSOC, using single-cell RNA sequencing and T-cell receptor sequencing techniques. Our analyses unveiled notable changes in the tumor signatures after chemotherapy, including those associated with epithelial-mesenchymal transition and cell cycle arrest. Within the immune compartment, we observed alterations in the T-cell profiles, characterized by naïve or pre-exhausted populations following chemotherapy. This phenotypic change was further supported by the examination of adjoining T-cell receptor clonotypes in paired longitudinal samples. These findings underscore the profound impact of chemotherapy on reshaping the tumor landscape and the immune microenvironment. This knowledge may provide clues for the development of future therapeutic strategies to combat treatment resistance in HGSOC.

Short tandem repeat expansions in cortical layer-specific genes implicate in phenotypic severity and adaptability of autism spectrum disorder.

AIM: Short tandem repeats (STRs) are repetitive DNA sequences and highly mutable in various human disorders. While the involvement of STRs in various genetic disorders has been extensively studied, their role in autism spectrum disorder (ASD) remains largely unexplored. In this study, we aimed to investigate genetic association of STR expansions with ASD using whole genome sequencing (WGS) and identify risk loci associated with ASD phenotypes. METHODS: We analyzed WGS data of 634 ASD families and performed genome-wide evaluation for 12,929 STR loci. We found rare STR expansions that exceeded normal repeat lengths in autism cases compared to unaffected controls. By integrating single cell RNA and ATAC sequencing datasets of human postmortem brains, we prioritized STR loci in genes specifically expressed in cortical development stages. A deep learning method was used to predict functionality of ASD-associated STR loci. RESULTS: In ASD cases, rare STR expansions predominantly occurred in early cortical layer-specific genes involved in neurodevelopment, highlighting the cellular specificity of STR-associated genes in ASD risk. Leveraging deep learning prediction models, we demonstrated that these STR expansions disrupted the regulatory activity of enhancers and promoters, suggesting a potential mechanism through which they contribute to ASD pathogenesis. We found that individuals with ASD-associated STR expansions exhibited more severe ASD phenotypes and diminished adaptability compared to non-carriers. CONCLUSION: Short tandem repeat expansions in cortical layer-specific genes are associated with ASD and could potentially be a risk genetic factor for ASD. Our study is the first to show evidence of STR expansion associated with ASD in an under-investigated population.

Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder.

Three-dimensional chromatin interactions regulate gene expressions. The significance of de novo mutations (DNMs) in chromatin interactions remains poorly understood for autism spectrum disorder (ASD). We generated 813 whole-genome sequences from 242 Korean simplex families to detect DNMs, and identified target genes which were putatively affected by non-coding DNMs in chromatin interactions. Non-coding DNMs in chromatin interactions were significantly involved in transcriptional dysregulations related to ASD risk. Correspondingly, target genes showed spatiotemporal expressions relevant to ASD in developing brains and enrichment in biological pathways implicated in ASD, such as histone modification. Regarding clinical features of ASD, non-coding DNMs in chromatin interactions particularly contributed to low intelligence quotient levels in ASD probands. We further validated our findings using two replication cohorts, Simons Simplex Collection (SSC) and MSSNG, and showed the consistent enrichment of non-coding DNM-disrupted chromatin interactions in ASD probands. Generating human induced pluripotent stem cells in two ASD families, we were able to demonstrate that non-coding DNMs in chromatin interactions alter the expression of target genes at the stage of early neural development. Taken together, our findings indicate that non-coding DNMs in ASD probands lead to early neurodevelopmental disruption implicated in ASD risk via chromatin interactions.

The emerging genetic diversity of hereditary spastic paraplegia in Korean patients.

Hereditary Spastic Paraplegias (HSP) are a group of rare inherited neurological disorders characterized by progressive loss of corticospinal motor-tract function. Numerous patients with HSP remain undiagnosed despite screening for known genetic causes of HSP. Therefore, identification of novel genetic variations related to HSP is needed. In this study, we identified 88 genetic variants in 54 genes from whole-exome data of 82 clinically well-defined Korean HSP families. Fifty-six percent were known HSP genes, and 44% were composed of putative candidate HSP genes involved in the HSPome and originally reported neuron-related genes, not previously diagnosed in HSP patients. Their inheritance modes were 39, de novo; 33, autosomal dominant; and 10, autosomal recessive. Notably, ALDH18A1 showed the second highest frequency. Fourteen known HSP genes were firstly reported in Koreans, with some of their variants being predictive of HSP-causing protein malfunction. SPAST and REEP1 mutants with unknown function induced neurite abnormality. Further, 54 HSP-related genes were closely linked to the HSP progression-related network. Additionally, the genetic spectrum and variation of known HSP genes differed across ethnic groups. These results expand the genetic spectrum for HSP and may contribute to the accurate diagnosis and treatment for rare HSP.

Immune signatures correlate with L1 retrotransposition in gastrointestinal cancers.

Long interspersed nuclear element-1 (LINE-1 or L1) retrotransposons are normally suppressed in somatic tissues mainly due to DNA methylation and antiviral defense. However, the mechanism to suppress L1s may be disrupted in cancers, thus allowing L1s to act as insertional mutagens and cause genomic rearrangement and instability. Whereas the frequency of somatic L1 insertions varies greatly among individual tumors, much remains to be learned about underlying genetic, cellular, or environmental factors. Here, we report multiple correlates of L1 activity in stomach, colorectal, and esophageal tumors through an integrative analysis of cancer whole-genome and matched RNA-sequencing profiles. Clinical indicators of tumor progression, such as tumor grade and patient age, showed positive association. A potential L1 expression suppressor, TP53, was mutated in tumors with frequent L1 insertions. We characterized the effects of somatic L1 insertions on mRNA splicing and expression, and demonstrated an increased risk of gene disruption in retrotransposition-prone cancers. In particular, we found that a cancer-specific L1 insertion in an exon of MOV10, a key L1 suppressor, caused exon skipping and decreased expression of the affected allele due to nonsense-mediated decay in a tumor with a high L1 insertion load. Importantly, tumors with high immune activity, for example, those associated with Epstein-Barr virus infection or microsatellite instability, tended to carry a low number of L1 insertions in genomes with high expression levels of L1 suppressors such as APOBEC3s and SAMHD1 Our results indicate that cancer immunity may contribute to genome stability by suppressing L1 retrotransposition in gastrointestinal cancers.

Convolutional neural network model to predict causal risk factors that share complex regulatory features.

Major progress in disease genetics has been made through genome-wide association studies (GWASs). One of the key tasks for post-GWAS analyses is to identify causal noncoding variants with regulatory function. Here, on the basis of >2000 functional features, we developed a convolutional neural network framework for combinatorial, nonlinear modeling of complex patterns shared by risk variants scattered among multiple associated loci. When applied for major psychiatric disorders and autoimmune diseases, neural and immune features, respectively, exhibited high explanatory power while reflecting the pathophysiology of the relevant disease. The predicted causal variants were concentrated in active regulatory regions of relevant cell types and tended to be in physical contact with transcription factors while residing in evolutionarily conserved regions and resulting in expression changes of genes related to the given disease. We demonstrate some examples of novel candidate causal variants and associated genes. Our method is expected to contribute to the identification and functional interpretation of potential causal noncoding variants in post-GWAS analyses.

Serine metabolism in the brain regulates starvation-induced sleep suppression in Drosophila melanogaster.

Sleep and metabolism are physiologically and behaviorally intertwined; however, the molecular basis for their interaction remains poorly understood. Here, we identified a serine metabolic pathway as a key mediator for starvation-induced sleep suppression. Transcriptome analyses revealed that enzymes involved in serine biosynthesis were induced upon starvation in Drosophila melanogaster brains. Genetic mutants of astray (aay), a fly homolog of the rate-limiting phosphoserine phosphatase in serine biosynthesis, displayed reduced starvation-induced sleep suppression. In contrast, a hypomorphic mutation in a serine/threonine-metabolizing enzyme, serine/threonine dehydratase (stdh), exaggerated starvation-induced sleep suppression. Analyses of double mutants indicated that aay and stdh act on the same genetic pathway to titrate serine levels in the head as well as to adjust starvation-induced sleep behaviors. RNA interference-mediated depletion of aay expression in neurons, using cholinergic Gal4 drivers, phenocopied aay mutants, while a nicotinic acetylcholine receptor antagonist selectively rescued the exaggerated starvation-induced sleep suppression in stdh mutants. Taken together, these data demonstrate that neural serine metabolism controls sleep during starvation, possibly via cholinergic signaling. We propose that animals have evolved a sleep-regulatory mechanism that reprograms amino acid metabolism for adaptive sleep behaviors in response to metabolic needs.

Selection on the regulation of sympathetic nervous activity in humans and chimpanzees.

Adrenergic α2C receptor (ADRA2C) is an inhibitory modulator of the sympathetic nervous system. Knockout mice for this gene show physiological and behavioural alterations that are associated with the fight-or-flight response. There is evidence of positive selection on the regulation of this gene during chicken domestication. Here, we find that the neuronal expression of ADRA2C is lower in human and chimpanzee than in other primates. On the basis of three-dimensional chromatin structure, we identified a cis-regulatory region whose DNA sequences have been significantly accelerated in human and chimpanzee. Active histone modification marks this region in rhesus macaque but not in human and chimpanzee; instead, repressive marks are enriched in various human brain samples. This region contains two neuron-restrictive silencer factor (NRSF) binding motifs, each of which harbours a polymorphism. Our genotyping and analysis of population genome data indicate that at both polymorphic sites, the derived allele has reached fixation in humans and chimpanzees but not in bonobos, whereas only the ancestral allele is present among macaques. Our CRISPR/Cas9 genome editing and reporter assays show that both derived nucleotides repress ADRA2C, most likely by increasing NRSF binding. In addition, we detected signatures of recent positive selection for lower neuronal ADRA2C expression in humans. Our findings indicate that there has been selective pressure for enhanced sympathetic nervous activity in the evolution of humans and chimpanzees.

Functional annotation of noncoding causal variants in autoimmune diseases.

Interpretation of noncoding disease variants, which comprise the vast majority of Genome-wide association studies (GWAS) hits, remains a momentous challenge due to haplotype structure and our limited understanding of the mechanisms and physiological contexts of noncoding elements. GWAS have identified loci underlying human diseases, but assigning the causal nucleotide changes still remain a controversial issue. Here we addressed these issues through the combination of high-density genotyping and epigenomic data using a random forest model to discover the noncoding causal variants. Focusing on autoimmune diseases, we triaged putative causal variants for atopic dermatitis and inflammatory bowel diseases. Making a filtering pipeline, we found three interesting single nucleotide polymorphisms (rs1800630, rs1799964 and rs4796793) in the upstream site of TNF and STAT3 genes, two frequent genes shared in some autoimmune diseases, and show how those variants affect on TNF and STAT3 expression levels. All data and source codes related to this manuscript are available at https://github.com/jieunjung511/Autoimmune-research.

Network perturbation by recurrent regulatory variants in cancer.

Cancer driving genes have been identified as recurrently affected by variants that alter protein-coding sequences. However, a majority of cancer variants arise in noncoding regions, and some of them are thought to play a critical role through transcriptional perturbation. Here we identified putative transcriptional driver genes based on combinatorial variant recurrence in cis-regulatory regions. The identified genes showed high connectivity in the cancer type-specific transcription regulatory network, with high outdegree and many downstream genes, highlighting their causative role during tumorigenesis. In the protein interactome, the identified transcriptional drivers were not as highly connected as coding driver genes but appeared to form a network module centered on the coding drivers. The coding and regulatory variants associated via these interactions between the coding and transcriptional drivers showed exclusive and complementary occurrence patterns across tumor samples. Transcriptional cancer drivers may act through an extensive perturbation of the regulatory network and by altering protein network modules through interactions with coding driver genes.

Intersection of genetics and epigenetics in monozygotic twin genomes.

As a final function of various epigenetic mechanisms, chromatin regulation is a transcription control process that especially demonstrates active interaction with genetic elements. Thus, chromatin structure has become a principal focus in recent genomics researches that strive to characterize regulatory functions of DNA variants related to diseases or other traits. Although researchers have been focusing on DNA methylation when studying monozygotic (MZ) twins, a great model in epigenetics research, interactions between genetics and epigenetics in chromatin level are expected to be an imperative research trend in the future. In this review, we discuss how the genome, epigenome, and transcriptome of MZ twins can be studied in an integrative manner from this perspective.

Global transcription network incorporating distal regulator binding reveals selective cooperation of cancer drivers and risk genes.

Global network modeling of distal regulatory interactions is essential in understanding the overall architecture of gene expression programs. Here, we developed a Bayesian probabilistic model and computational method for global causal network construction with breast cancer as a model. Whereas physical regulator binding was well supported by gene expression causality in general, distal elements in intragenic regions or loci distant from the target gene exhibited particularly strong functional effects. Modeling the action of long-range enhancers was critical in recovering true biological interactions with increased coverage and specificity overall and unraveling regulatory complexity underlying tumor subclasses and drug responses in particular. Transcriptional cancer drivers and risk genes were discovered based on the network analysis of somatic and genetic cancer-related DNA variants. Notably, we observed that the risk genes were functionally downstream of the cancer drivers and were selectively susceptible to network perturbation by tumorigenic changes in their upstream drivers. Furthermore, cancer risk alleles tended to increase the susceptibility of the transcription of their associated genes. These findings suggest that transcriptional cancer drivers selectively induce a combinatorial misregulation of downstream risk genes, and that genetic risk factors, mostly residing in distal regulatory regions, increase transcriptional susceptibility to upstream cancer-driving somatic changes.

HBx induces hypomethylation of distal intragenic CpG islands required for active expression of developmental regulators.

Epigenetic alterations caused by viral oncoproteins are strong initiation factors for cancer development, but their mechanisms are largely unknown. To identify the epigenetic effects of viral hepatitis B virus X (HBx) that lead to hepatocellular carcinoma (HCC), we profiled the DNA methylomes of normal and HBx transgenic mouse liver. Intriguingly, severe hypomethylation of intragenic CpG islands (CGIs) was observed in HBx liver before the full development of HCC. Normally, these CGIs were highly methylated (mCGIs) by the DNMT3L complex and marked with epigenetic signatures associated with active expression, such as H3K36me3. Hypomethylation of mCGI was caused by the downregulation of Dnmt3L and Dnmt3a due to HBx bound to their promoters, along with HDAC1. These events lead to the downregulation of many developmental regulators that could facilitate tumorigenesis. Here we provide an intriguing epigenetic regulation mediated by mCGI that is required for cell differentiation and describe a previously unidentified epigenetic role for HBx in promoting HCC development.

Predicting the recurrence of noncoding regulatory mutations in cancer.

BACKGROUND: One of the greatest challenges in cancer genomics is to distinguish driver mutations from passenger mutations. Whereas recurrence is a hallmark of driver mutations, it is difficult to observe recurring noncoding mutations owing to a limited amount of whole-genome sequenced samples. Hence, it is required to develop a method to predict potentially recurrent mutations. RESULTS: In this work, we developed a random forest classifier that predicts regulatory mutations that may recur based on the features of the mutations repeatedly appearing in a given cohort. With breast cancer as a model, we profiled 35 quantitative features describing genetic and epigenetic signals at the mutation site, transcription factors whose binding motif was disrupted by the mutation, and genes targeted by long-range chromatin interactions. A true set of mutations for machine learning was generated by interrogating publicly available pan-cancer genomes based on our statistical model of mutation recurrence. The performance of our random forest classifier was evaluated by cross validations. The variable importance of each feature in the classification of mutations was investigated. Our statistical recurrence model for the random forest classifier showed an area under the curve (AUC) of ~0.78 in predicting recurrent mutations. Chromatin accessibility at the mutation sites, the distance from the mutations to known cancer risk loci, and the role of the target genes in the regulatory or protein interaction network were among the most important variables. CONCLUSIONS: Our methods enable to characterize recurrent regulatory mutations using a limited number of whole-genome samples, and based on the characterization, to predict potential driver mutations whose recurrence is not found in the given samples but likely to be observed with additional samples.

Genome-wide reorganization of histone H2AX toward particular fragile sites on cell activation.

γH2AX formation by phosphorylation of the histone variant H2AX is the key process in the repair of DNA lesions including those arising at fragile sites under replication stress. Here we demonstrate that H2AX is dynamically reorganized to preoccupy γH2AX hotspots on increased replication stress by activated cell proliferation and that H2AX is enriched in aphidicolin-induced replisome stalling sites in cycling cells. Interestingly, H2AX enrichment was particularly found in genomic regions that replicate in early S phase. High transcription activity, a hallmark of early replicating fragile sites, was a determinant of H2AX localization. Subtelomeric H2AX enrichment was also attributable to early replication and high gene density. In contrast, late replicating and infrequently transcribed regions, including common fragile sites and heterochromatin, lacked H2AX enrichment. In particular, heterochromatin was inaccessible to H2AX incorporation, maybe partly explaining the cause of mutation accumulation in cancer heterochromatin. Meanwhile, H2AX in actively dividing cells was intimately colocalized with INO80. INO80 silencing reduced H2AX levels, particularly at the INO80-enriched sites. Our findings suggest that active DNA replication is accompanied with the specific localization of H2AX and INO80 for efficient damage repair or replication-fork stabilization in actively transcribed regions.

Regulation of the boundaries of accessible chromatin.

Regulatory regions maintain nucleosome-depleted, open chromatin status but simultaneously require the presence of nucleosomes for specific histone modifications. It remains unclear how these can be achieved for proper regulatory function. Here we demonstrate that nucleosomes positioned within accessible chromatin regions near the boundaries provide platforms for histone modifications while preventing the occlusion of regulatory elements. These boundary nucleosomes were particularly enriched for active or poised regulatory marks in human, such as histone acetylations, H3K4 methylations, H3K9me3, H3K79me2, and H4K20me1. Additionally, we found that based on a genome-wide profiling of ~100 recombinant yeast strains, the location of open chromatin borders tends to vary mostly within 150 bp upon genetic perturbation whereas this positional variation increases in proportion to the sequence preferences of the underlying DNA for nucleosome formation. More than 40% of the local boundary shifts were associated with genetic variation in cis- or trans-acting factors. A sizeable fraction of the identified genetic factors was also associated with nearby gene expression, which was correlated with the distance between the transcription start site (tss) and the boundary that faces the tss. Taken together, the variation in the width of accessible chromatin regions may arise in conjunction with the modulation of the boundary nucleosomes by post-translational modifications or by chromatin regulators and in association with the activity of nearby gene transcription.

Genetic landscape of open chromatin in yeast.

Chromatin regulation underlies a variety of DNA metabolism processes, including transcription, recombination, repair, and replication. To perform a quantitative genetic analysis of chromatin accessibility, we obtained open chromatin profiles across 96 genetically different yeast strains by FAIRE (formaldehyde-assisted isolation of regulatory elements) assay followed by sequencing. While 5∼10% of open chromatin region (OCRs) were significantly affected by variations in their underlying DNA sequences, subtelomeric areas as well as gene-rich and gene-poor regions displayed high levels of sequence-independent variation. We performed quantitative trait loci (QTL) mapping using the FAIRE signal for each OCR as a quantitative trait. While individual OCRs were associated with a handful of specific genetic markers, gene expression levels were associated with many regulatory loci. We found multi-target trans-loci responsible for a very large number of OCRs, which seemed to reflect the widespread influence of certain chromatin regulators. Such regulatory hotspots were enriched for known regulatory functions, such as recombinational DNA repair, telomere replication, and general transcription control. The OCRs associated with these multi-target trans-loci coincided with recombination hotspots, telomeres, and gene-rich regions according to the function of the associated regulators. Our findings provide a global quantitative picture of the genetic architecture of chromatin regulation.

Genome-wide methylation profiling of ADPKD identified epigenetically regulated genes associated with renal cyst development.

Autosomal dominant polycystic kidney disease (ADPKD) is a common human genetic disease characterized by the formation of multiple fluid-filled cysts in bilateral kidneys. Although mutations in polycystic kidney disease 1 (PKD1) are predominantly responsible for ADPKD, the focal and sporadic property of individual cystogenesis suggests another molecular mechanism such as epigenetic alterations. To determine the epigenomic alterations in ADPKD and their functional relevance, ADPKD and non-ADPKD individuals were analyzed by unbiased methylation profiling genome-wide and compared with their expression data. Intriguingly, PKD1 and other genes related to ion transport and cell adhesion were hypermethylated in gene-body regions, and their expressions were downregulated in ADPKD, implicating epigenetic silencing as the key mechanism underlying cystogenesis. Especially, in patients with ADPKD, PKD1 was hypermethylated in gene-body region and it was associated with recruitment of methyl-CpG-binding domain 2 proteins. Moreover, treatment with DNA methylation inhibitors retarded cyst formation of Madin-Darby Canine Kidney cells, accompanied with the upregulation of Pkd1 expression. These results are consistent with previous studies that knock-down of PKD1 was sufficient for cystogenesis. Therefore, our results reveal a critical role for hypermethylation of PKD1 and cystogenesis-related regulatory genes in cyst development, suggesting epigenetic therapy as a potential treatment for ADPKD.

Genome-wide profiles of H2AX and γ-H2AX differentiate endogenous and exogenous DNA damage hotspots in human cells.

Phosphorylation of the histone variant H2AX forms γ-H2AX that marks DNA double-strand break (DSB). Here, we generated the sequencing-based maps of H2AX and γ-H2AX positioning in resting and proliferating cells before and after ionizing irradiation. Genome-wide locations of possible endogenous and exogenous DSBs were identified based on γ-H2AX distribution in dividing cancer cells without irradiation and that in resting cells upon irradiation, respectively. γ-H2AX-enriched regions of endogenous origin in replicating cells included sub-telomeres and active transcription start sites, apparently reflecting replication- and transcription-mediated stress during rapid cell division. Surprisingly, H2AX itself, prior to phosphorylation, was specifically located at these endogenous hotspots. This phenomenon was only observed in dividing cancer cells but not in resting cells. Endogenous H2AX was concentrated on the transcription start site of actively transcribed genes but was irrelevant to pausing of RNA polymerase II (pol II), which precisely coincided with γ-H2AX of endogenous origin. γ-H2AX enrichment upon irradiation also coincided with actively transcribed regions, but unlike endogenous γ-H2AX, it extended into the gene body and was not specifically concentrated on the pausing site of pol II. Sub-telomeres were less responsive to external DNA damage than to endogenous stress. Our findings provide insight into DNA repair programs of cancer and may have implications for cancer therapy.