Clonal origin and genomic diversity in Lynch syndrome-associated endometrial cancer with multiple synchronous tumors: Identification of the pathogenicity of MLH1 p.L582H.

Lynch syndrome-associated endometrial cancer patients often present multiple synchronous tumors and this assessment can affect treatment strategies. We present a case of a 27-year-old woman with tumors in the uterine corpus, cervix, and ovaries who was diagnosed with endometrial cancer and exhibited cervical invasion and ovarian metastasis. Her family history suggested Lynch syndrome, and genetic testing identified a variant of uncertain significance, MLH1 p.L582H. We conducted immunohistochemical staining, microsatellite instability analysis, and Sanger sequencing for Lynch syndrome-associated cancers in three generations of the family and identified consistent MLH1 loss. Whole-exome sequencing for the corpus, cervical, and ovarian tumors of the proband identified a copy-neutral loss of heterozygosity (LOH) occurring at the MLH1 position in all tumors. This indicated that the germline variant and the copy-neutral LOH led to biallelic loss of MLH1 and was the cause of cancer initiation. All tumors shared a portion of somatic mutations with high mutant allele frequencies, suggesting a common clonal origin. There were no mutations shared only between the cervix and ovary samples. The profiles of mutant allele frequencies shared between the corpus and cervix or ovary indicated that two different subclones originating from the corpus independently metastasized to the cervix or ovary. Additionally, all tumors presented unique mutations in endometrial cancer-associated genes such as ARID1A and PIK3CA. In conclusion, we demonstrated clonal origin and genomic diversity in a Lynch syndrome-associated endometrial cancer, suggesting the importance of evaluating multiple sites in Lynch syndrome patients with synchronous tumors.

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.

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.

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.

Biased expression of mutant alleles in cancer-related genes in esophageal squamous cell carcinoma.

BACKGROUND: Recent progress of large-scale international studies has provided comprehensive catalogs of somatic mutations in cancers. Additionally, it has become evident that allelic imbalance in the abundance of somatic mutations between DNA and RNA were pervasive in various types of cancer. However, the allelic imbalance of the abundance of somatic mutations in esophageal squamous cell carcinoma (ESCC) has not been fully analyzed. METHODS: We performed exome sequencing for 25 Japanese patients with ESCC to detect a comprehensive catalog of somatic mutations in ESCC. Additionally, we performed mRNA sequencing to evaluate the allelic imbalance of the identified somatic mutations at the transcriptional level by comparing the mutant allele frequencies between RNA and DNA. RESULTS: The exome sequencing showed that TP53 and ZNF750 were significantly mutated genes. The expression levels of TP53 and ZNF750 were different depending on the mutation status. In almost all the tumors with missense mutations in TP53 and ZNF750, the mutant allele frequencies were higher in the RNA sequencing than those in the exome sequencing, indicating that the mutant alleles were preferentially expressed. By examining the allelic imbalances for all the identified missense mutations, we demonstrated that genes showing preferential expressions of the mutant alleles were involved in the pathways including cell cycle, cell death, and chromatin modification. CONCLUSIONS: The results of this study suggest that the allelic imbalance of the abundance of somatic mutations plays important roles in the initiation and progression of ESCC by modulating cancer-related biological pathways.

Biological significance of KRAS mutant allele expression in ovarian endometriosis.

KRAS is the most frequently mutated in ovarian endometriosis. However, it is unclear whether the KRAS mutant allele's mRNA is expressed and plays a biological role in ovarian endometriosis. Here, we performed mutation-specific RNA in situ hybridization to evaluate mutant allele expression of KRAS p.G12V, the most frequently detected mutation in ovarian endometriosis in our previous study, in formalin-fixed paraffin-embedded tissue (FFPE) samples of ovarian endometriosis, cancer cell lines, and ovarian cancers. First, we verified that mutant or wild-type allele of KRAS were expressed in all 5 cancer cell lines and 9 ovarian cancer cases corresponding to the mutation status. Next, we applied this assay to 26 ovarian endometriosis cases, and observed mutant allele expression of KRAS p.G12V in 10 cases. Mutant or wild-type allele of KRAS were expressed in line with mutation status in 12 available endometriosis cases for which KRAS gene sequence was determined. Comparison of clinical features between ovarian endometriosis with KRAS p.G12V mutant allele expression and with KRAS wild-type showed that KRAS p.G12V mutant allele expression was significantly associated with inflammation in ovarian endometriosis. Finally, we assessed the spatial distribution of KRAS mutant allele expression in 5 endometriosis cases by performing multiregional sampling. Intratumor heterogeneity of KRAS mutant allele expression was observed in two endometriosis cases, whereas the spatial distribution of KRAS p.G12V mutation signals were diffuse and homogenous in ovarian cancer. In conclusion, evaluation of oncogene mutant expression will be useful for clarifying the biological significance of oncogene mutations in benign tumors.

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.

Subtype-specific gout susceptibility loci and enrichment of selection pressure on ABCG2 and ALDH2 identified by subtype genome-wide meta-analyses of clinically defined gout patients.

OBJECTIVES: Genome-wide meta-analyses of clinically defined gout were performed to identify subtype-specific susceptibility loci. Evaluation using selection pressure analysis with these loci was also conducted to investigate genetic risks characteristic of the Japanese population over the last 2000-3000 years. METHODS: Two genome-wide association studies (GWASs) of 3053 clinically defined gout cases and 4554 controls from Japanese males were performed using the Japonica Array and Illumina Array platforms. About 7.2 million single-nucleotide polymorphisms were meta-analysed after imputation. Patients were then divided into four clinical subtypes (the renal underexcretion type, renal overload type, combined type and normal type), and meta-analyses were conducted in the same manner. Selection pressure analyses using singleton density score were also performed on each subtype. RESULTS: In addition to the eight loci we reported previously, two novel loci, PIBF1 and ACSM2B, were identified at a genome-wide significance level (p<5.0×10-8) from a GWAS meta-analysis of all gout patients, and other two novel intergenic loci, CD2-PTGFRN and SLC28A3-NTRK2, from normal type gout patients. Subtype-dependent patterns of Manhattan plots were observed with subtype GWASs of gout patients, indicating that these subtype-specific loci suggest differences in pathophysiology along patients' gout subtypes. Selection pressure analysis revealed significant enrichment of selection pressure on ABCG2 in addition to ALDH2 loci for all subtypes except for normal type gout. CONCLUSIONS: Our findings on subtype GWAS meta-analyses and selection pressure analysis of gout will assist elucidation of the subtype-dependent molecular targets and evolutionary involvement among genotype, phenotype and subtype-specific tailor-made medicine/prevention of gout and hyperuricaemia.

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.

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.

Different mutation profiles between epithelium and stroma in endometriosis and normal endometrium.

STUDY QUESTION: Are there common mutation profiles between epithelial and stromal cells in ovarian endometriotic tissue and the normal endometrium? SUMMARY ANSWER: Our study revealed no common mutations between epithelial and stromal cells in ovarian endometriotic tissue and the normal endometrium. WHAT IS KNOWN ALREADY: Epithelial cells in both ovarian endometriotic tissue and the normal endometrium harbor somatic mutations in cancer-associated genes such as phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and KRAS proto-oncogene, GTPase (KRAS). STUDY DESIGN, SIZE, DURATION: We performed a retrospective study to identify the mutation profiles of stromal cells in endometriotic tissue and the normal endometrium. We collected 11 endometriotic stroma samples and 10 normal endometrial stroma samples between 2013 and 2017 at a tertiary care center. PARTICIPANTS/MATERIALS, SETTING, METHODS: The laser microdissection method was used to obtain stromal cells in ovarian endometriotic and normal endometrial tissues from patients with ovarian endometriosis and/or other non-invasive gynecological diseases. Target gene sequencing was performed to assess and compare the mutation profiles of stromal cells with those of epithelial cells obtained in our previous study. For target gene sequencing, 76 genes were selected based on previous genomic analyses for ovarian endometriosis, normal endometrium, endometriosis-related ovarian cancer and endometrial cancer. MAIN RESULTS AND THE ROLE OF CHANCE: Stromal samples in ovarian endometrioma and normal endometrium harbor somatic mutations (18 mutations in 11 endometriosis samples and 16 mutations in 10 normal endometrial samples) but did not share any mutations with paired epithelial samples. The mutant allele frequency of stromal samples was significantly lower than that of epithelial samples in ovarian endometrioma (P = 6.0 × 10-11) and normal endometrium (P = 1.4 × 10-7). LIMITATIONS, REASONS FOR CAUTION: The number of genes evaluated in the mutational analysis was limited. Additionally, the functional roles of somatic mutations in stromal cells remain unclear. WIDER IMPLICATIONS OF THE FINDINGS: Different mutation profiles between paired epithelial and stromal cells in both ovarian endometrioma and normal endometrium suggest that origins of epithelial and stromal cells would be independent of each other in both normal endometrium and ovarian endometrioma; however, the theory of epithelial-mesenchymal transition is proposed in ovarian endometrioma. STUDY FUNDING/COMPETING INTEREST(S): This work was supported in part by the Japan Society for the Promotion of Science KAKENHI grant number JP15H02373 (Grant-in-Aid for Scientific Research A for I.I.), JP16H06267 (Grant-in-Aid for Young Scientists A for K.Y.), JP17K08688 (Grant-in-Aid for Scientific Research C for H.N.) and JP16H06279 (Grant-in-Aid for Scientific Research on Innovative Areas-Platforms for Advanced Technologies and Research Resources for H.N. and K.Y). There are no conflicts of interest to declare. TRIAL REGISTRATION NUMBER: Not applicable.

Allelic Imbalance in Regulation of ANRIL through Chromatin Interaction at 9p21 Endometriosis Risk Locus.

Genome-wide association studies (GWASs) have discovered numerous single nucleotide polymorphisms (SNPs) associated with human complex disorders. However, functional characterization of the disease-associated SNPs remains a formidable challenge. Here we explored regulatory mechanism of a SNP on chromosome 9p21 associated with endometriosis by leveraging "allele-specific" functional genomic approaches. By re-sequencing 1.29 Mb of 9p21 region and scrutinizing DNase-seq data from the ENCODE project, we prioritized rs17761446 as a candidate functional variant that was in perfect linkage disequilibrium with the original GWAS SNP (rs10965235) and located on DNase I hypersensitive site. Chromosome conformation capture followed by high-throughput sequencing revealed that the protective G allele of rs17761446 exerted stronger chromatin interaction with ANRIL promoter. We demonstrated that the protective allele exhibited preferential binding affinities to TCF7L2 and EP300 by bioinformatics and chromatin immunoprecipitation (ChIP) analyses. ChIP assays for histone H3 lysine 27 acetylation and RNA polymerase II reinforced the enhancer activity of the SNP site. The allele specific expression analysis for eutopic endometrial tissues and endometrial carcinoma cell lines showed that rs17761446 was a cis-regulatory variant where G allele was associated with increased ANRIL expression. Our work illuminates the allelic imbalances in a series of transcriptional regulation from factor binding to gene expression mediated by chromatin interaction underlie the molecular mechanism of 9p21 endometriosis risk locus. Functional genomics on common disease will unlock functional aspect of genotype-phenotype correlations in the post-GWAS stage.