Super enhancer loci of EGFR regulate EGFR variant 8 through enhancer RNA and strongly associate with survival in HNSCCs.

Epidermal growth factor receptor (EGFR) has been shown to be overexpressed in human cancers due to mutation, amplification, and epigenetic hyperactivity, which leads to deregulated transcriptional mechanism. Among the eight different EGFR isoforms, the mechanism of regulation of full-length variant 1 is well-known, no studies have examined the function & factors regulating the expression of variant 8. This study aimed to understand the function of EGFR super-enhancer loci and its associated transcription factors regulating the expression of EGFR variant 8. Our study shows that overexpression of variant 8 and its transcription was more prevalent than variant 1 in many cancers and positively correlated with the EGFR-AS1 expression in oral cancer and HNSCC. Notably, individuals overexpressing variant 8 showed shorter overall survival and had a greater connection with other clinical traits than patients with overexpression of variant 1. In this study, TCGA enhancer RNA profiling on the constituent enhancer (CE1 and CE2) region revealed that the multiple enhancer RNAs formed from CE2 by employing CE1 as a promoter. Our bioinformatic analysis further supports the enrichment of enhancer RNA specific chromatin marks H3K27ac, H3K4me1, POL2 and H2AZ on CE2. GeneHancer and 3D chromatin capture analysis showed clustered interactions between CE1, CE2 loci and this interaction may regulates expression of both EGFR-eRNA and variant 8. Moreover, increased expression of SNAI2 and its close relationship to EGFR-AS1 and variant 8 suggest that SNAI2 could regulates variant 8 overexpression by building a MegaTrans complex with both EGFR-eRNA and EGFR-AS1. Our findings show that EGFR variant 8 and its transcriptional regulation & chromatin modification by eRNAs may provide a rationale for targeting RNA splicing in combination with targeted EGFR therapies in cancer.

Noninvasive graft monitoring using donor-derived cell-free DNA in Japanese liver transplantation.

AIM: To evaluate the use of donor-derived cell-free DNA (dd-cfDNA) in diagnosing graft injuries in Japanese liver transplantation (LTx), including family-related living donors. METHODS: A total of 321 samples from 10 newly operated LTx recipients were collected to monitor the early dynamics of dd-cfDNA levels after LTx. Fifty-five samples from 55 recipients were collected during protocol biopsies (PB), whereas 36 samples from 27 recipients were collected during event biopsies, consisting of 11 biopsy-proven acute rejection (AR), 20 acute dysfunctions without rejection (ADWR), and 5 chronic rejections. The levels of dd-cfDNA were quantified using a next-generation sequencer based on single nucleotide polymorphisms. RESULTS: The dd-cfDNA levels were elevated significantly after LTx, followed by a rapid decline to the baseline in patients without graft injury within 30 days post-LTx. The dd-cfDNA levels were significantly higher in the 11 samples obtained during AR than those obtained during PB (p < 0.0001), which decreased promptly after treatment. The receiver operator characteristic curve analysis of diagnostic ability yielded areas under the curve of 0.975 and 0.897 for AR (rejection activity index [RAI] ≥3) versus PB and versus non-AR (ADWR + PB). The dd-cfDNA levels during AR were elevated earlier and correlated more strongly with the RAI (r = 0.740) than aspartate aminotransferase/alanine aminotransferase. The dd-cfDNA levels were neither associated with graft fibrosis based on histology nor the status of donor-specific antibodies in PB samples. CONCLUSIONS: Donor-derived cell-free DNA serves as a sensitive biomarker for detecting graft injuries in LTx. Further large-scale cohort studies are warranted to optimize its use in differentiating various post-LTx etiologies.

Spatial genomic diversity associated with APOBEC mutagenesis in squamous cell carcinoma arising from ovarian teratoma.

Although the gross and microscopic features of squamous cell carcinoma arising from ovarian mature cystic teratoma (MCT-SCC) vary from case to case, the spatial spreading of genomic alterations within the tumor remains unclear. To clarify the spatial genomic diversity in MCT-SCCs, we performed whole-exome sequencing by collecting 16 samples from histologically different parts of two MCT-SCCs. Both cases showed histological diversity within the tumors (case 1: nonkeratinizing and keratinizing SCC and case 2: nonkeratinizing SCC and anaplastic carcinoma) and had different somatic mutation profiles by histological findings. Mutation signature analysis revealed a significantly enriched apolipoprotein B mRNA editing enzyme catalytic subunit (APOBEC) signature at all sites. Intriguingly, the spread of genomic alterations within the tumor and the clonal evolution patterns from nonmalignant epithelium to cancer sites differed between cases. TP53 mutation and copy number alterations were widespread at all sites, including the nonmalignant epithelium, in case 1. Keratinizing and nonkeratinizing SCCs were differentiated by the occurrence of unique somatic mutations from a common ancestral clone. In contrast, the nonmalignant epithelium showed almost no somatic mutations in case 2. TP53 mutation and the copy number alteration similarities were observed only in nonkeratinizing SCC samples. Nonkeratinizing SCC and anaplastic carcinoma shared almost no somatic mutations, suggesting that each locally and independently arose in the MCT. We demonstrated that two MCT-SCCs with different histologic findings were highly heterogeneous tumors with clearly different clones associated with APOBEC-mediated mutagenesis, suggesting the importance of evaluating intratumor histological and genetic heterogeneity among multiple sites of MCT-SCC.

Movements of Ancient Human Endogenous Retroviruses Detected in SOX2-Expressing Cells.

Human endogenous retroviruses (HERVs) occupy approximately 8% of the human genome. HERVs, transcribed in early embryos, are epigenetically silenced in somatic cells, except under pathological conditions. HERV-K is thought to protect embryos from exogenous viral infection. However, uncontrolled HERV-K expression in somatic cells has been implicated in several diseases. Here, we show that SOX2, which plays a key role in maintaining the pluripotency of stem cells, is critical for HERV-K LTR5Hs. HERV-K undergoes retrotransposition within producer cells in the absence of Env expression. Furthermore, we identified new HERV-K integration sites in long-term culture of induced pluripotent stem cells that express SOX2. These results suggest that the strict dependence of HERV-K on SOX2 has allowed HERV-K to protect early embryos during evolution while limiting the potentially harmful effects of HERV-K retrotransposition on host genome integrity in these early embryos. IMPORTANCE Human endogenous retroviruses (HERVs) account for approximately 8% of the human genome; however, the physiological role of HERV-K remains unknown. This study found that HERV-K LTR5Hs and LTR5B were transactivated by SOX2, which is essential for maintaining and reestablishing pluripotency. HERV-K can undergo retrotransposition within producer cells without env expression, and new integration sites may affect cell proliferation. In induced pluripotent stem cells (iPSCs), genomic impairment due to HERV-K retrotransposition has been identified, but it is a rare event. Considering the retention of SOX2-responsive elements in the HERV-K long terminal repeat (LTR) for over 20 million years, we conclude that HERV-K may play important physiological roles in SOX2-expressing cells.

Genome-wide meta-analysis between renal overload type and renal underexcretion type of clinically defined gout in Japanese populations.

Despite progress in understanding of the genetic basis of gout, the precise factors affecting differences in gout susceptibility among different gout subtypes remain unclear. Using clinically diagnosed gout patients, we conducted a genome-wide meta-analysis of two distinct gout subtypes: the renal overload type and the renal underexcretion type. We provide genetic evidence at a genome-wide level of significance that supports a positive association between ABCG2 dysfunction and acquisition of the renal overload type.

APOBEC mediated mutagenesis drives genomic heterogeneity in endometriosis.

Endometriosis is a benign gynecologic condition, acting as a precursor of certain histological subtypes of ovarian cancers. The epithelial cells of endometriotic tissues and normal uterine endometrium accumulated somatic mutations in cancer-associated genes such as phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and Kirsten rat sarcoma (KRAS) proto-oncogene. To determine the genomic characteristic of endometriotic epithelial cells and normal uterine endometrium and to identify the predominant mutational process acting on them, we studied the somatic mutation profiles obtained from whole exome sequencing of 14 endometriotic epithelium and 11 normal uterine endometrium tissues and classified them into mutational signatures. We observed that single base substitutions 2/13 (SBS), attributed to Apolipoprotein B mRNA Editing Enzyme Catalytic Subunit (APOBEC) induced mutagenesis, were significant in endometriotic tissues, but not in the normal uterine endometrium. Additionally, the larger number and wider allele frequency distribution of APOBEC signature mutations, compared to cancer-associated driver mutations in endometriotic epithelium suggested APOBEC mutagenesis as an important source of mutational burden and heterogeneity in endometriosis. Further, the relative risk of enriched APOBEC signature mutations was higher in endometriosis patients who were carriers of APOBEC3A/3B germline deletion, a common polymorphism in East Asians which involves the complete loss of APOBEC3B coding region. Our results illustrate the significance of APOBEC induced mutagenesis in driving the genomic heterogeneity of endometriosis.

The landscape of genetic alterations of UVB-induced skin tumors in DNA repair-deficient mice.

Non-melanoma skin cancer (NMSC) is mainly caused by ultraviolet (UV)-induced somatic mutations and is characterized by UV signature modifications. Xeroderma pigmentosum group A (Xpa) knockout mice exhibit extreme UV-induced photo-skin carcinogenesis, along with a photosensitive phenotype. We performed whole-exome sequencing (WES) of squamous cell carcinoma (SCC) samples after repetitive ultraviolet B (UVB) exposure to investigate the differences in the landscape of somatic mutations between Xpa knockout and wild-type mice. Although the tumors that developed in mice harboured UV signature mutations in a similar set of cancer-related genes, the pattern of transcriptional strand asymmetry was largely different; UV signature mutations in Xpa knockout and wild-type mice preferentially occurred in transcribed and non-transcribed strands, respectively, reflecting a deficiency in transcription-coupled nucleotide excision repair in Xpa knockout mice. Serial time point analyses of WES for a tumor induced by only a single UVB exposure showed pathogenic mutations in Kras, Fat1, and Kmt2c, which may be driver genes for the initiation and promotion of SCC in Xpa knockout mice. Furthermore, the inhibitory effects on tumor production in Xpa knockout mice by the anti-inflammatory CXCL1 monoclonal antibody affected the pattern of somatic mutations, wherein the transcriptional strand asymmetry was attenuated and the activated signal transduction was shifted from the RAS/RAF/MAPK to the PIK3CA pathway.

AMBRA1 controls antigen-driven activation and proliferation of naive T cells.

AMBRA1 (activating molecule in Beclin1-regulated autophagy) is a member of the BECN1 (BECLIN1) protein complex, and it plays a role in autophagy, cell death, tumorigenesis and proliferation. We recently reported that on T-cell receptor (TCR) stimulation, AMBRA1 controlled both autophagy and the cell cycle with metabolic regulation. Accumulating evidence has shown that autophagy and metabolic control are pivotal for T-cell activation, clonal expansion and effector/memory cell fate decision. However, it is unknown whether AMBRA1 is involved in T-cell function under physiological conditions. We found that T cells in Ambra1-conditional knockout (cKO) mice induced an exacerbated graft versus host response when they were transplanted into allogeneic BALB/c mice. Furthermore, Ambra1-deficient T cells showed increased proliferation and cytotoxic capability toward specific antigens in response to in vivo stimulation using allogeneic spleen cells. This enhanced immune response mainly contributed to naive T-cell hyperactivity. The T-cell hyperactivity observed in this study was similar to those in some metabolic factor-deficient mice, but not those in other pro-autophagic factor-deficient mice. Under the static condition, however, naive T cells were reduced in Ambra1-cKO mice, the same as in pro-autophagic factor-deficient mice. Collectively, these results suggested that AMBRA1 was involved in regulating T cell-mediated immune responses through autophagy-dependent and -independent mechanisms.

Comprehensive discovery of CRISPR-targeted terminally redundant sequences in the human gut metagenome: Viruses, plasmids, and more.

Viruses are the most numerous biological entity, existing in all environments and infecting all cellular organisms. Compared with cellular life, the evolution and origin of viruses are poorly understood; viruses are enormously diverse, and most lack sequence similarity to cellular genes. To uncover viral sequences without relying on either reference viral sequences from databases or marker genes that characterize specific viral taxa, we developed an analysis pipeline for virus inference based on clustered regularly interspaced short palindromic repeats (CRISPR). CRISPR is a prokaryotic nucleic acid restriction system that stores the memory of previous exposure. Our protocol can infer CRISPR-targeted sequences, including viruses, plasmids, and previously uncharacterized elements, and predict their hosts using unassembled short-read metagenomic sequencing data. By analyzing human gut metagenomic data, we extracted 11,391 terminally redundant CRISPR-targeted sequences, which are likely complete circular genomes. The sequences included 2,154 tailed-phage genomes, together with 257 complete crAssphage genomes, 11 genomes larger than 200 kilobases, 766 genomes of Microviridae species, 56 genomes of Inoviridae species, and 95 previously uncharacterized circular small genomes that have no reliably predicted protein-coding gene. We predicted the host(s) of approximately 70% of the discovered genomes at the taxonomic level of phylum by linking protospacers to taxonomically assigned CRISPR direct repeats. These results demonstrate that our protocol is efficient for de novo inference of CRISPR-targeted sequences and their host prediction.

AMBRA1 p.Gln30Arg Mutation, Identified in a Cowden Syndrome Family, Exhibits Hyperproliferative Potential in hTERT-RPE1 Cells.

Cowden syndrome (CS) is a rare autosomal dominant disorder associated with multiple hamartomatous and neoplastic lesions in various organs. Most CS patients have been found to have germline mutations in the PTEN tumor suppressor. In the present study, we investigated the causative gene of CS in a family of PTEN (phosphatase and tensin homolog deleted on chromosome 10) -negative CS patients. Whole exome sequencing analysis revealed AMBRA1 (Autophagy and Beclin 1 Regulator 1) as a novel candidate gene harboring two germline variants: p.Gln30Arg (Q30R) and p.Arg1195Ser (R1195S). AMBRA1 is a key regulator of the autophagy signaling network and a tumor suppressor. To functionally validate the role of AMBRA1 in the clinical manifestations of CS, we generated AMBRA1 depletion and Q30R mutation in hTERT-RPE1 (humanTelomerase Reverse Transcriptase-immortalized Retinal Pigmented Epithelial cells) using the CRISPR-Cas9 gene editing system. We observed that both AMBRA1-depleted and mutant cells showed accumulation in the S phase, leading to hyperproliferation, which is a characteristic of hamartomatous lesions. Specifically, the AMBRA1 Q30R mutation disturbed the G1/S transition of cells, leading to continuous mitotic entry of mutant cells, irrespective of the extracellular condition. From our analysis of primary ciliogenesis in these cells, we speculated that the mitotic entry of AMBRA1 Q30R mutants could be due to non-functional primary cilia that lead to impaired processing of extracellular sensory signals. Additionally, we observed a situs inversus phenotype in ambra1-depleted zebrafish, a developmental abnormality resulting from dysregulated primary ciliogenesis. Taken together, we established that the AMBRA1 Q30R mutation that we observed in CS patients might play an important role in inducing the hyperproliferative potential of cells through regulating primary ciliogenesis.

Population genetics: past, present, and future.

We present selected topics of population genetics and molecular phylogeny. As several excellent review articles have been published and generally focus on European and American scientists, here, we emphasize contributions by Japanese researchers. Our review may also be seen as a belated 50-year celebration of Motoo Kimura's early seminal paper on the molecular clock, published in 1968.

HLA-B*39:01:01 is a novel risk factor for antithyroid drug-induced agranulocytosis in Japanese population.

Antithyroid drug (ATD) is a mainstay of Graves' disease (GD). About 0.1-0.5% of patients with GD treated with ATD exhibit ATD-induced agranulocytosis, which is characterized by severe reduction of circulating neutrophils. Immune-mediated responses have been proposed as a possible mechanism for the pathogenesis of ATD-induced agranulocytosis. Although it has been reported that the HLA class II allele (HLA-DRB1*08:03) was associated with ATD-induced agranulocytosis in multiple populations, the entire HLA region have not been explored in Japanese. Therefore, we performed HLA sequencing for 10 class I and 11 class II genes in 87 patients with ATD-induced agranulocytosis and 384 patients with GD who did not show ATD-induced agranulocytosis. By conducting case-control association studies at the HLA allele and haplotype levels, we replicated the association between HLA-DRB1*08:03:02 and ATD-induced agranulocytosis (P = 5.2 × 10-7, odds ratio = 2.80), and identified HLA-B*39:01:01 as an independent risk factor (P = 1.4 × 10-3, odds ratio = 3.35). To verify reproducibility of the novel association of HLA-B*39:01:01, we reanalyzed allele frequency data for HLA-B*39:01:01 from previous case-control association studies. The association of HLA-B*39:01:01 was significantly replicated in Chinese (P = 9.0 × 10-3), Taiwanese (P = 1.1 × 10-3), and European populations (P = 5.2 × 10-4). A meta-analysis combining results from the previous and current studies reinforced evidence of the association between HLA-B*39:01:01 and ATD-induced agranulocytosis (Pmeta = 1.2 × 10-9, pooled OR = 3.66, 95% CI; 2.41-5.57). The results of this study will provide a better understanding of the pathogenesis of ATD-induced agranulocytosis in the context of HLA-mediated hypersensitivity reactions.

Identification of ancient viruses from metagenomic data of the Jomon people.

Ancient DNA studies provide genomic information about the origins, population structures, and physical characteristics of ancient humans that cannot be solely examined by archeological studies. The DNAs extracted from ancient human bones, teeth, or tissues are often contaminated with coexisting bacterial and viral genomes that contain DNA from ancient microbes infecting those of ancient humans. Information on ancient viral genomes is useful in making inferences about the viral evolution. Here, we have utilized metagenomic sequencing data from the dental pulp of five Jomon individuals, who lived on the Japanese archipelago more than 3000 years ago; this is to detect ancient viral genomes. We conducted de novo assembly of the non-human reads where we have obtained 277,387 contigs that were longer than 1000 bp. These contigs were subjected to homology searches against a collection of modern viral genome sequences. We were able to detect eleven putative ancient viral genomes. Among them, we reconstructed the complete sequence of the Siphovirus contig89 (CT89) viral genome. The Jomon CT89-like sequence was determined to contain 59 open reading frames, among which five genes known to encode phage proteins were under strong purifying selection. The host of CT89 was predicted to be Schaalia meyeri, a bacterium residing in the human oral cavity. Finally, the CT89 phylogenetic tree showed two clusters, from both of which the Jomon sequence was separated. Our results suggest that metagenomic information from the dental pulp of the Jomon people is essential in retrieving ancient viral genomes used to examine their evolution.

Substantial anti-gout effect conferred by common and rare dysfunctional variants of URAT1/SLC22A12.

OBJECTIVES: Gout, caused by chronic elevation of serum uric acid levels, is the commonest form of inflammatory arthritis. The causative effect of common and rare variants of ATP-binding cassette transporter G2 (ABCG2/BCRP) on gout risk has been studied, but little attention has been paid to the effect of common (rs121907892, p.W258X) and rare variants of urate transporter 1 (URAT1/SLC22A12) on gout, despite dysfunctional variants of URAT1 having been identified as pathophysiological causes of renal hypouricaemia. METHODS: To address this important but overlooked issue, we investigated the effects of these URAT1 variants on gout susceptibility, using targeted exon sequencing on 480 clinically defined gout cases and 480 controls of Japanese males in combination with a series of functional analyses of newly identified URAT1 variants. RESULTS: Our results show that both common and rare dysfunctional variants of URAT1 markedly decrease the risk of gout (OR 0.0338, reciprocal OR 29.6, P = 7.66 × 10-8). Interestingly, we also found that the URAT1-related protective effect on gout eclipsed the ABCG2-related causative effect (OR 2.30-3.32). Our findings reveal only one dysfunctional variant of URAT1 to have a substantial anti-gout effect, even in the presence of causative variants of ABCG2, a 'gout gene'. CONCLUSION: Our findings provide a better understanding of gout/hyperuricaemia and its aetiology that is highly relevant to personalized health care. The substantial anti-gout effect of common and rare variants of URAT1 identified in the present study support the genetic concept of a 'Common Disease, Multiple Common and Rare Variant' model.

Proposing a molecular classification associated with hypercoagulation in ovarian clear cell carcinoma.

BACKGROUND: Although ovarian clear cell carcinoma (CCC) is associated with high incidence of thromboembolism, the clinicopathological and biological significance of hypercoagulable status in CCC remains unclear. MATERIALS AND METHODS: We retrospectively analyzed pretreatment D-dimer levels, thromboembolic status, and clinical outcome of 125 CCCs in the discovery set and 143 CCCs in two other independent validation sets. Next, we performed RNA sequencing of 93 CCCs and compared coagulation-related gene profiles with 2492 pan-cancer data. We investigated differences in molecular characteristics of CCC subclasses based on coagulation status. RESULTS: In the discovery dataset, D-dimer elevation above the normal range was significantly associated with shorter progression-free and overall survival, irrespective to thromboembolic status. Multivariate analysis identified D-dimer elevation and clinical stage as an independent prognostic factors. We confirmed the prognostic significance of D-dimer elevation in the validation sets. Tissue factor and IL6, which are considered key elements of cancer-induced hypercoagulation, were highly expressed in CCC than in other cancers regardless of D-dimer level. Higher activity of various oncogenic pathways was observed in CCC with compared to without D-dimer elevation. Moreover, hierarchical cluster analysis divided 57 CCCs with D-dimer elevation into immunologically hot and cold tumor subtypes. Hot tumors were characterized by enrichment of T-cell inflamed phenotype, inflammation, the epithelial-mesenchymal transition, and high serum levels of CRP, and cold tumors by enrichment of cell cycle and MYC pathways. CONCLUSIONS: CCC represents hypercoagulable disease and elevate D-dimer is a prognostic factor for decreased survival in CCC. D-dimer high CCC has distinct molecular characteristics into the inflammatory-driven pathway (hot tumor) and the immune-suppressive pathway (cold tumor). Treatment implication of our proposed molecular classification merits further investigation.

Clonal lineage from normal endometrium to ovarian clear cell carcinoma through ovarian endometriosis.

Clear cell carcinoma of the ovary is thought to arise from endometriosis. In addition, retrograde menstruation of shed endometrium is considered the origin of endometriosis. However, little evidence supports cellular continuity from uterine endometrium to clear cell carcinoma through endometriosis at the genomic level. Here, we performed multiregional whole-exome sequencing to clarify clonal relationships among uterine endometrium, ovarian endometriosis and ovarian clear cell carcinoma in a 56-year-old patient. Many somatic mutations including cancer-associated gene mutations in ARID1A, ATM, CDH4, NRAS and PIK3CA were shared among epithelium samples from uterine endometrium, endometriotic lesions distant from and adjacent to the carcinoma, and the carcinoma. The mutant allele frequencies of shared mutations increased from uterine endometrium to distant endometriosis, adjacent endometriosis, and carcinoma. Although a splice site mutation of ARID1A was shared among the four epithelium samples, a frameshift insertion in ARID1A was shared by adjacent endometriosis and carcinoma samples, suggesting that the biallelic mutations triggered malignant transformation. Somatic copy number alterations, including loss of heterozygosity events at PIK3CA and ATM, were identified only in adjacent endometriosis and carcinoma, suggesting that mutant allele-specific imbalance is another key factor driving malignant transformation. By reconstructing a clonal evolution tree based on the somatic mutations, we showed that the epithelium samples were derived from a single ancestral clone. Although the study was limited to a single patient, the results from this illustrative case could suggest the possibility that epithelial cells of ovarian endometriosis and clear cell carcinoma were descendants of uterine endometrial epithelium.

Germline mutations of multiple breast cancer-related genes are differentially associated with triple-negative breast cancers and prognostic factors.

Genetic testing for BRCA1/2 mutations has become the standard clinical practice. Recent findings suggest the clinical significance of multigene panel testing of BRCA1/2 and other cancer-related genes. However, the clinical features of patients with breast cancer with germline mutations identified using multigene panels remain unclear. In this study, DNA samples from 583 Chinese women with breast cancer were subjected to target sequencing for 54 cancer-related genes using a pre-capture pooling method followed by next-generation sequencing. We identified 79 pathogenic germline mutations in 21 cancer-related genes. Forty-five patients (7.7%) harbored BRCA1/2 mutations, and 38 patients (6.5%) carried pathogenic mutations in the remaining 19 genes. PALB2 was the most commonly (1.2%) mutated gene other than BRCA1/2. Most of the identified pathogenic mutations were novel, suggesting mutation screening by using multigene panel testing is important particularly for non-European populations. Mutations in BRCA1/2 and the other cancer-related genes were differentially associated with clinical features. BRCA1 mutation carriers were strongly associated with triple-negative breast cancer (TNBC), whereas BRCA2 mutation carriers were not. Tumors in BRCA1-mutation carriers had a high histological grade. Patients with BRCA2-mutated breast cancers were likely to develop E-cadherin-negative tumors with bone metastases. Furthermore, mutations in PALB2 were strongly associated with TNBC. We demonstrated the usefulness of multigene panel testing and observed that a substantial proportion of patients with breast cancer had hereditary risk factors. Identifying differential associations between mutation status and clinical features will advance our understanding regarding the pathologies of this heterogeneous disease.

Systematic identification and characterization of regulatory elements derived from human endogenous retroviruses.

Human endogenous retroviruses (HERVs) and other long terminal repeat (LTR)-type retrotransposons (HERV/LTRs) have regulatory elements that possibly influence the transcription of host genes. We systematically identified and characterized these regulatory elements based on publicly available datasets of ChIP-Seq of 97 transcription factors (TFs) provided by ENCODE and Roadmap Epigenomics projects. We determined transcription factor-binding sites (TFBSs) using the ChIP-Seq datasets and identified TFBSs observed on HERV/LTR sequences (HERV-TFBSs). Overall, 794,972 HERV-TFBSs were identified. Subsequently, we identified "HERV/LTR-shared regulatory element (HSRE)," defined as a TF-binding motif in HERV-TFBSs, shared within a substantial fraction of a HERV/LTR type. HSREs could be an indication that the regulatory elements of HERV/LTRs are present before their insertions. We identified 2,201 HSREs, comprising specific associations of 354 HERV/LTRs and 84 TFs. Clustering analysis showed that HERV/LTRs can be grouped according to the TF binding patterns; HERV/LTR groups bounded to pluripotent TFs (e.g., SOX2, POU5F1, and NANOG), embryonic endoderm/mesendoderm TFs (e.g., GATA4/6, SOX17, and FOXA1/2), hematopoietic TFs (e.g., SPI1 (PU1), GATA1/2, and TAL1), and CTCF were identified. Regulatory elements of HERV/LTRs tended to locate nearby and/or interact three-dimensionally with the genes involved in immune responses, indicating that the regulatory elements play an important role in controlling the immune regulatory network. Further, we demonstrated subgroup-specific TF binding within LTR7, LTR5B, and LTR5_Hs, indicating that gains or losses of the regulatory elements occurred during genomic invasions of the HERV/LTRs. Finally, we constructed dbHERV-REs, an interactive database of HERV/LTR regulatory elements (http://herv-tfbs.com/). This study provides fundamental information in understanding the impact of HERV/LTRs on host transcription, and offers insights into the transcriptional modulation systems of HERV/LTRs and ancestral HERVs.

Next-generation sequencing identifies contribution of both class I and II HLA genes on susceptibility of multiple sclerosis in Japanese.

BACKGROUND: The spectrum of classical and non-classical HLA genes related to the risk of multiple sclerosis (MS) and neuromyelitis optica spectrum disorder (NMOSD) in the Japanese population has not been studied in detail. We conducted a case-control analysis of classical and non-classical HLA genes. METHODS: We used next-generation sequencing (NGS)-based HLA genotyping methods for mapping risk for 45 MS patients, 31 NMOSD patients, and 429 healthy controls. We evaluated the association of the HLA variants with the risk of MS and NMOSD using logistic regression analysis and Fisher's exact test. RESULTS: We confirmed that HLA-DRB1*15:01 showed the strongest association with MS (P = 2.1 × 10-5; odds ratio [OR] = 3.44, 95% confidence interval [95% CI] = 1.95-6.07). Stepwise conditional analysis identified HLA-DRB1*04:05, HLA-B*39:01, and HLA-B*15:01 as being associated with independent MS susceptibility (PConditional < 8.3 × 10-4). With respect to amino acid polymorphisms in HLA genes, we found that phenylalanine at HLA-DQß1 position 9 had the strongest effect on MS susceptibility (P = 3.7 × 10-8, OR = 3.48, 95% CI = 2.23-5.43). MS risk at HLA-DQß1 Phe9 was independent of HLA-DRB1*15:01 (PConditional = 1.5 × 10-5, OR = 2.91, 95% CI = 1.79-4.72), while HLA-DRB1*15:01 was just significant when conditioned on HLA-DQß1 Phe9 (PConditional = 0.037). Regarding a case-control analysis for NMOSD, HLA-DQA1*05:03 had a significant association with NMOSD (P = 1.5 × 10-4, OR = 6.96, 95% CI = 2.55-19.0). CONCLUSIONS: We identified HLA variants associated with the risk of MS and NMOSD. Our study contributes to the understanding of the genetic architecture of MS and NMOSD in the Japanese population.

Genome-wide association study revealed novel loci which aggravate asymptomatic hyperuricaemia into gout.

OBJECTIVE: The first ever genome-wide association study (GWAS) of clinically defined gout cases and asymptomatic hyperuricaemia (AHUA) controls was performed to identify novel gout loci that aggravate AHUA into gout. METHODS: We carried out a GWAS of 945 clinically defined gout cases and 1003 AHUA controls followed by 2 replication studies. In total, 2860 gout cases and 3149 AHUA controls (all Japanese men) were analysed. We also compared the ORs for each locus in the present GWAS (gout vs AHUA) with those in the previous GWAS (gout vs normouricaemia). RESULTS: This new approach enabled us to identify two novel gout loci (rs7927466 of CNTN5 and rs9952962 of MIR302F) and one suggestive locus (rs12980365 of ZNF724) at the genome-wide significance level (p<5.0×10-8). The present study also identified the loci of ABCG2, ALDH2 and SLC2A9. One of them, rs671 of ALDH2, was identified as a gout locus by GWAS for the first time. Comparing ORs for each locus in the present versus the previous GWAS revealed three 'gout vs AHUA GWAS'-specific loci (CNTN5, MIR302F and ZNF724) to be clearly associated with mechanisms of gout development which distinctly differ from the known gout risk loci that basically elevate serum uric acid level. CONCLUSIONS: This meta-analysis is the first to reveal the loci associated with crystal-induced inflammation, the last step in gout development that aggravates AHUA into gout. Our findings should help to elucidate the molecular mechanisms of gout development and assist the prevention of gout attacks in high-risk AHUA individuals.