Endometriosis is characterized by ectopic endometrial-like epithelium and stroma, of which molecular characteristics remain to be fully elucidated. We sequenced 107 ovarian endometriotic and 82 normal uterine endometrial epithelium samples isolated by laser microdissection. In both endometriotic and normal epithelium samples, numerous somatic mutations were identified within genes frequently mutated in endometriosis-associated ovarian cancers. KRAS is frequently mutated in endometriotic epithelium, with a higher mutant allele frequency (MAF) accompanied by arm-level allelic imbalances. Analyses of MAF, combined with multiregional sequencing, illuminated spatiotemporal evolution of the endometriosis and uterine endometrium genomes. We sequenced 109 single endometrial glands and found that each gland carried distinct cancer-associated mutations, demonstrating the heterogeneity of the genomic architecture of endometrial epithelium. Remarkable increases in MAF of mutations in cancer-associated genes in endometriotic epithelium suggest retrograde flow of endometrial cells already harboring cancer-associated mutations, with selective advantages at ectopic sites, leading to the development of endometriosis.

Protein knockdown using the auxin-inducible degron (AID) technology is useful to study protein function in living cells because it induces rapid depletion, which makes it possible to observe an immediate phenotype. However, the current AID system has two major drawbacks: leaky degradation and the requirement for a high dose of auxin. These negative features make it difficult to control precisely the expression level of a protein of interest in living cells and to apply this method to mice. Here, we overcome these problems by taking advantage of a bump-and-hole approach to establish the AID version 2 (AID2) system. AID2, which employs an OsTIR1(F74G) mutant and a ligand, 5-Ph-IAA, shows no detectable leaky degradation, requires a 670-times lower ligand concentration, and achieves even quicker degradation than the conventional AID. We demonstrate successful generation of human cell mutants for genes that were previously difficult to deal with, and show that AID2 achieves rapid target depletion not only in yeast and mammalian cells, but also in mice.

ObjectiveGout, caused by hyperuricaemia, is a multifactorial disease. Although genome-wide association studies (GWASs) of gout have been reported, they included self-reported gout cases in which clinical information was insufficient. Therefore, the relationship between genetic variation and clinical subtypes of gout remains unclear. Here, we first performed a GWAS of clinically defined gout cases only.MethodsA GWAS was conducted with 945 patients with clinically defined gout and 1213 controls in a Japanese male population, followed by replication study of 1048 clinically defined cases and 1334 controls.ResultsFive gout susceptibility loci were identified at the genome-wide significance level (p<5.0×10−8), which contained well-known urate transporter genes (ABCG2 and SLC2A9) and additional genes: rs1260326 (p=1.9×10−12; OR=1.36) of GCKR (a gene for glucose and lipid metabolism), rs2188380 (p=1.6×10−23; OR=1.75) of MYL2-CUX2 (genes associated with cholesterol and diabetes mellitus) and rs4073582 (p=6.4×10−9; OR=1.66) of CNIH-2 (a gene for regulation of glutamate signalling). The latter two are identified as novel gout loci. Furthermore, among the identified single-nucleotide polymorphisms (SNPs), we demonstrated that the SNPs of ABCG2 and SLC2A9 were differentially associated with types of gout and clinical parameters underlying specific subtypes (renal underexcretion type and renal overload type). The effect of the risk allele of each SNP on clinical parameters showed significant linear relationships with the ratio of the case–control ORs for two distinct types of gout (r=0.96 [p=4.8×10−4] for urate clearance and r=0.96 [p=5.0×10−4] for urinary urate excretion).ConclusionsOur findings provide clues to better understand the pathogenesis of gout and will be useful for development of companion diagnostics.

Traditional approaches for gene mapping from candidate gene studies to positional cloning strategies have been applied for Mendelian disorders. Since 2005, next-generation sequencing (NGS) technologies are improving as rapid, high-throughput and cost-effective approaches to fulfill medical sciences and research demands. Using NGS, the underlying causative genes are directly distinguished via a systematic filtering, in which the identified gene variants are checked for novelty and functionality. During the past 2 years, the role of more than 100 genes has been distinguished in rare Mendelian disorders by means of whole-exome sequencing (WES). Combination of WES with traditional approaches, consistent with linkage analysis, has had the greatest impact on those disorders following autosomal mode of inheritance; in more than 60 identified genes, the causal variants have been transmitted at homozygous or compound heterozygous state. Recent literatures focusing on identified new causal genes in Mendelian disorders using WES are reviewed in the present survey.

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.

BackgroundThe human leukocyte antigen (HLA) region, the 3.8-Mb segment of the human genome at 6p21, has been associated with more than 100 different diseases, mostly autoimmune diseases. Due to the complex nature of HLA genes, there are difficulties in elucidating complete HLA gene sequences especially HLA gene haplotype structures by the conventional sequencing method. We propose a novel, accurate, and cost-effective method for generating phase-defined complete sequencing of HLA genes by using indexed multiplex next generation sequencing.ResultsA total of 33 HLA homozygous samples, 11 HLA heterozygous samples, and 3 parents-child families were subjected to phase-defined HLA gene sequencing. We applied long-range PCR to amplify six HLA genes (HLA-A, -C, -B, DRB1, -DQB1, and –DPB1) followed by transposase-based library construction and multiplex sequencing with the MiSeq sequencer. Paired-end reads (2 × 250 bp) derived from the sequencer were aligned to the six HLA gene segments of UCSC hg19 allowing at most 80 bases mismatch. For HLA homozygous samples, the six amplicons of an individual were pooled and simultaneously sequenced and mapped as an individual-tagging method. The paired-end reads were aligned to corresponding genes of UCSC hg19 and unambiguous, continuous sequences were obtained. For HLA heterozygous samples, each amplicon was separately sequenced and mapped as a gene-tagging method. After alignments, we detected informative paired-end reads harboring SNVs on both forward and reverse reads that are used to separate two chromosomes and to generate two phase-defined sequences in an individual. Consequently, we were able to determine the phase-defined HLA gene sequences from promoter to 3′-UTR and assign up to 8-digit HLA allele numbers, regardless of whether the alleles are rare or novel. Parent–child trio-based sequencing validated our sequencing and phasing methods.ConclusionsOur protocol generated phased-defined sequences of the entire HLA genes, resulting in high resolution HLA typing and new allele detection.

Gout is a common disease which results from hyperuricemia. We have reported that the dysfunction of urate exporter ABCG2 is the major cause of renal overload (ROL) hyperuricemia, but its involvement in renal underexcretion (RUE) hyperuricemia, the most prevalent subtype, is not clearly explained so far. In this study, the association analysis with 644 hyperuricemia patients and 1,623 controls in male Japanese revealed that ABCG2 dysfunction significantly increased the risk of RUE hyperuricemia as well as overall and ROL hyperuricemia, according to the severity of impairment. ABCG2 dysfunction caused renal urate underexcretion and induced hyperuricemia even if the renal urate overload was not remarkable. These results show that ABCG2 plays physiologically important roles in both renal and extra-renal urate excretion mechanisms. Our findings indicate the importance of ABCG2 as a promising therapeutic and screening target of hyperuricemia and gout.

ObjectiveA genome-wide association study (GWAS) of gout and its subtypes was performed to identify novel gout loci, including those that are subtype-specific.MethodsPutative causal association signals from a GWAS of 945 clinically defined gout cases and 1213 controls from Japanese males were replicated with 1396 cases and 1268 controls using a custom chip of 1961 single nucleotide polymorphisms (SNPs). We also first conducted GWASs of gout subtypes. Replication with Caucasian and New Zealand Polynesian samples was done to further validate the loci identified in this study.ResultsIn addition to the five loci we reported previously, further susceptibility loci were identified at a genome-wide significance level (p<5.0×10−8): urate transporter genes (SLC22A12 and SLC17A1) and HIST1H2BF-HIST1H4E for all gout cases, and NIPAL1 and FAM35A for the renal underexcretion gout subtype. While NIPAL1 encodes a magnesium transporter, functional analysis did not detect urate transport via NIPAL1, suggesting an indirect association with urate handling. Localisation analysis in the human kidney revealed expression of NIPAL1 and FAM35A mainly in the distal tubules, which suggests the involvement of the distal nephron in urate handling in humans. Clinically ascertained male patients with gout and controls of Caucasian and Polynesian ancestries were also genotyped, and FAM35A was associated with gout in all cases. A meta-analysis of the three populations revealed FAM35A to be associated with gout at a genome-wide level of significance (pmeta=3.58×10−8).ConclusionsOur findings including novel gout risk loci provide further understanding of the molecular pathogenesis of gout and lead to a novel concept for the therapeutic target of gout/hyperuricaemia.