Clonal Expansion and Diversification of Cancer-Associated Mutations in Endometriosis and Normal Endometrium

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.

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.

P2434Digital zoom decreases radiation exposure dose up to 30% in percutaneous coronary intervention

Background

Interventional cardiology is gaining greater popularity worldwide with each passing year. Reduction of exposure dose is a very imminent and an important issue in cardiology procedure. Although a newer radiation reduction technique, device and procedure are very valuable and expected, we should consider about therapy technique, radiation technique, devices, and the way to protection. Digital zoom digitally enlarges images in real time by up to 2.5-fold at lower doses than those used with traditional field of view changes. In our phantom examination the average dose reduction of digital zoom was 27%.

Methods and results

This study is designated as single-center, retrospective, not-randomized, observation study. 2101 eligible cases were collected. We assigned the cases of PCI without the use of Digital zoom to the Conventional group and those involving the use of Digital zoom to the Digital zoom group. There were 806 patients in the Conventional group and 1195 in the Digital zoom group. Because we had begun using Digital zoom from January 2015 onwards, all patients in the Conventional group had undergone PCI from January 2013 to December 2014 and all patients in the Digital zoom group had undergone PCI from January 2015 to December 2016. In addition, we calculated the RAK/minute and DAP/minute for an accurate assessment. To minimize the difference of characteristics between two groups, propensity score including all baseline variables was performed. Furthermore, Predictors of radiation exposure were investigated using multivariable least square methods. Inter group differences were observed in DAP, RAK, DAP/min, and RAK/min (Digital zoom group vs conventional group: DAP, 16000 cGy cm2 [from 1st quartile to 3rd quartile; 10300–24400] vs 20700 [13400–29500], p<0.001; DAP/min, 557 cGy cm2/min [392–737] vs 782 [571–1010], p<0.01; RAK, 1590 Gy [990–2410] vs 1850 [1220–2720], p<0.01; RAK/min, 54.7 Gy/min [38.5–73.2] vs 71.2 [51.5–93.0], p<0.01). Even after propensity score matching, intergroup differences in DAP (810 cases), DAP/min (811 cases), RAK (746 cases), and RAK/min (744 cases) persisted. Furthermore, the least squares method showed that Digital zoom is an important predictor of DAP (β=0.17, p<0.01) and RAK (β=0.12, p<0.01).

Conclusion

Digital zoom is an old and cost-free technique, but one of most powerful reduction of exposure method. Propensity score adjustment and least square methods show that digital zoom is one of independent effective method.

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⁻⁵; 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 (P_Conditional < 8.3 × 10⁻⁴). 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⁻⁸, OR = 3.48, 95% CI = 2.23–5.43). MS risk at HLA-DQβ1 Phe9 was independent of HLA-DRB1*15:01 (P_Conditional = 1.5 × 10⁻⁵, OR = 2.91, 95% CI = 1.79–4.72), while HLA-DRB1*15:01 was just significant when conditioned on HLA-DQβ1 Phe9 (P_Conditional = 0.037). Regarding a case-control analysis for NMOSD, HLA-DQA1*05:03 had a significant association with NMOSD (P = 1.5 × 10⁻⁴, 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.

Electronic supplementary material

The online version of this article (10.1186/s12974-019-1551-z) contains supplementary material, which is available to authorized users.

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.

High Order Formation and Evolution of Hornerin in Primates

Genomic duplication or loss can accelerate evolution because the number of repeats could affect molecular pathways and phenotypes. We have previously reported that the repeated region of filaggrin (FLG), a crucial component of the outer layers of mammalian skin, had high levels of nucleotide diversity with species-specific divergence and expansion and that it evolved under the birth-and-death model. We focused on hornerin (HRNR), a member of the same gene family that harbor similar tandem repeats as FLG, and examined the formation process of repeated regions and the evolutional model that best fit the HRNR repeated region in the crab-eating macaque (Macaca fascicularis), orangutan (Pongo abelii), gorilla (Gorilla gorilla), and chimpanzee (Pan troglodytes) and compared them with the human (Homo sapiens) sequence. Paar et al. (2011) and Takaishi et al. (2005) have different theories as to the formation of the repeated region of HRNR; both groups share the longest repeat length of 1,404 bp (quartic or longest unit), but they differed in the process. We identified the formation described by Paar et al. {[(“39 bp (primary) × 9” × 2 (secondary)) × 2 (tertiary)] × 5 (quartic)} to be conserved in all species except the crab-eating macaque. We detected high nucleotide diversities between the longest repeats, which fits the birth-and-death model. We concluded that the high order repeat formation of HRNR was conserved in primates except the crab-eating macaque. As previously identified in FLG, the longest repeats have high levels of nucleotide diversity, which could contribute to phenotypic differences between closely related species.

Novel MXD4-NUTM1 fusion transcript identified in primary ovarian undifferentiated small round cell sarcoma

Primary ovarian sarcomas are extremely rare tumors, and their genomic and transcriptomic alterations remain to be elucidated. We performed whole exome sequencing of primary tumor and matched normal blood samples derived from one patient with ovarian undifferentiated small round cell sarcoma. We identified 8 nonsynonymous somatic mutations, and all mutations were missense or nonsense changes. Next, we performed RNA sequencing of the tumor sample and identified two in-frame fusion transcripts: MXD4-NUTM1 and ARL6-POT1. Most NUTM1 exons were retained in the MXD4-NUTM1 fusion transcript, and we confirmed an increase in NUTM1 mRNA and protein expression in tumor tissue. Further genomic and transcriptomic analyses might lead to the development of new therapeutic strategies based on the molecular characteristics of ovarian undifferentiated small round cell sarcoma.

Comprehensive microbiome analysis of tonsillar crypts in IgA nephropathy

Background

Immunoglobulin A nephropathy (IgAN) is the most prevalent primary chronic glomerular disease, in which the mucosal immune response elicited particularly in the tonsils or intestine has been estimated to be involved in the development of the disease. To explore the relationship between IgAN and bacterial flora in the tonsils, we conducted a comprehensive microbiome analysis.

Methods

We enrolled 48 IgAN patients, 21 recurrent tonsillitis (RT) patients without urine abnormalities and 30 children with tonsillar hyperplasia (TH) who had undergone tonsillectomy previously. Genomic DNA from tonsillar crypts of each patient was extracted, and V4 regions of the 16S ribosomal RNA gene were amplified and analysed using a high-throughput multiplexed sequencing approach. Differences in genus composition among the three study groups were statistically analysed by permutational multivariate analysis of variance and visualized by principal component analysis (PCA).

Results

Substantial diversity in bacterial composition was detected in each sample. Prevotella spp., Fusobacterium spp., Sphingomonas spp. and Treponema spp. were predominant in IgAN patients. The percentage of abundance of Prevotella spp., Haemophilus spp., Porphyromonas spp. and Treponema spp. in IgAN patients was significantly different from that in TH patients. However, there was no significant difference in the percentage of abundance of any bacterial genus between IgAN and RT patients. PCA did not distinguish IgAN from RT, although it discriminated TH. No significant differences in microbiome composition among the groups of IgAN patients according to clinicopathological parameters were observed.

Conclusions

Similar patterns of bacteria are present in tonsillar crypts of both IgAN and RT patients, suggesting that the host response to these bacteria might be important in the development of IgAN.

Identification of a novel variant of the RET proto-oncogene in a novel family with Hirschsprung's disease

PURPOSE

Hirschsprung's disease (HSCR) is a congenital disorder of the enteric nervous system characterized by the absence of ganglion cells in the Auerbach's and Meissner's plexuses. Although about 7% of cases are hereditary, the causal mutations have not been completely characterized. We encountered a novel family with inherited HSCR and screened them for causal mutations.

METHODS

A Japanese family of five female patients and six unaffected individuals was subjected to a whole-exome analysis with a next-generation sequencer.

RESULTS

After exome sequencing and the annotation of mutations, we identified co-segregated mutations with sequential filtering steps via a standard protocol. Eight mutations were identified: two on chromosome 10 and six on chromosome 11. We used pathogenicity prediction tools such as Genomic Evolutionary Rate Profiling, SIFT, and PolyPhen2 to predict the impact of mutations on the protein activity. S922Y, a novel mutation of RET, was identified as a likely causal mutation. In addition, a mutation of rs2435357T, known as enhancer of RET located in intron 1 of RET, was detected in this family.

CONCLUSION

The coexistence of RET mutations in both the exon (S922Y) and intron1 (rs2435357T) indicated a risk of HSCR in this family.

Multiple common and rare variants of ABCG2 cause gout

Objective

Previous studies have suggested an association between gout susceptibility and common dysfunctional variants in ATP-binding cassette transporter subfamily G member 2/breast cancer resistance protein (ABCG2/BCRP), including rs72552713 (Q126X) and rs2231142 (Q141K). However, the association of rare ABCG2 variants with gout is unknown. Therefore, we investigated the effects of rare ABCG2 variants on gout susceptibility in this study.

Methods

We sequenced the exons of ABCG2 in 480 patients with gout and 480 healthy controls (Japanese males). We also performed functional analyses of non-synonymous variants of ABCG2 and analysed the correlation between urate transport function and scores from the protein prediction algorithms (Sorting Intolerant from Tolerant (SIFT) and Polymorphism Phenotyping v2 (PolyPhen-2)). Stratified association analyses and multivariate logistic regression analysis were performed to evaluate the effects of rare and common ABCG2 variants on gout susceptibility.

Results

We identified 3 common and 19 rare non-synonymous variants of ABCG2. SIFT scores were significantly correlated with the urate transport function, although some ABCG2 variants showed inconsistent scores. When the effects of common variants were removed by stratified association analysis, the rare variants of ABCG2 were associated with a significantly increased risk of gout (OR=3.2, p=6.4×10⁻³). Multivariate logistic regression analysis revealed that the size effect of these rare ABCG2 variants (OR=2.7, p=3.0×10⁻³) was similar to that of the common variants, Q126X (OR=3.4, p=3.2×10⁻⁶) and Q141K (OR=2.3, p=2.7×10⁻¹⁶).

Conclusions

This study revealed that multiple common and rare variants of ABCG2 are independently associated with gout. These results could support both the ‘Common Disease, Common Variant’ and ‘Common Disease, Multiple Rare Variant’ hypotheses for the association between ABCG2 and gout susceptibility.

AMBRA1 is involved in T cell receptor-mediated metabolic reprogramming through an ATG7-independent pathway

Metabolic reprogramming contributes to dynamic alteration of cell functions and characteristics. In T cells, TCR-mediated signaling evokes metabolic reprogramming and autophagy. AMBRA1 is known to serve in the facilitation of autophagy and quality control of mitochondria, but the role of AMBRA1 in T cell metabolic alteration is unknown. Here, we show that AMBRA1, but not ATG7, plays a role in TCR-mediated control of glycolytic factors and mitochondrial mass, while both AMBRA1 and ATG7 are required for autolysosome formation. Our results suggested that AMBRA1 is a core factor that controls both autophagy and metabolic regulation.

Comprehensive genetic exploration of selective tooth agenesis of mandibular incisors by exome sequencing

Tooth agenesis is described as the absence of one or more teeth. It is caused by a failure in tooth development and is one of the most common human developmental anomalies. We herein report genomic analyses of selective mandibular incisor agenesis (SMIA) using exome sequencing. Two Japanese families with SMIA were subjected to exome sequencing, and family with sequence similarity 65 member A (FAM65), nuclear factor of activated T-cells 3 (NFATC3) and cadherin-related 23 gene (CDH23) were detected. In the follow-up study, 51 Japanese and 32 Korean sporadic patients with SMIA were subjected to exome analyses, and 18 reported variants in PAX9, AXIN2, EDA, EDAR, WNT10A, BMP2 and GREM2 and 27 variants of FAM65, NFATC3 and CDH23 were found in 38 patients. Our comprehensive genetic study of SMIA will pave the way for a full understanding of the genetic etiology of SMIA and provide targets for treatment.

Structure and evolution of the filaggrin gene repeated region in primates

BACKGROUND

The evolutionary dynamics of repeat sequences is quite complex, with some duplicates never having differentiated from each other. Two models can explain the complex evolutionary process for repeated genes-concerted and birth-and-death, of which the latter is driven by duplications maintained by selection. Copy number variations caused by random duplications and losses in repeat regions may modulate molecular pathways and therefore affect phenotypic characteristics in a population, resulting in individuals that are able to adapt to new environments. In this study, we investigated the filaggrin gene (FLG), which codes for filaggrin-an important component of the outer layers of mammalian skin-and contains tandem repeats that exhibit copy number variation between and within species. To examine which model best fits the evolutionary pathway for the complete tandem repeats within a single exon of FLG, we determined the repeat sequences in crab-eating macaque (Macaca fascicularis), orangutan (Pongo abelii), gorilla (Gorilla gorilla), and chimpanzee (Pan troglodytes) and compared these with the sequence in human (Homo sapiens).

RESULTS

In this study we compared concerted and birth-and-death evolution models, commonly used for gene copies. We found that there is high nucleotide diversity between filaggrin repeat regions, which fits the birth-and-death model. Phylogenetic analyses also suggested that independent duplication events created the repeat sequences in crab-eating macaques and orangutans, while different duplication and loss events created the repeats in gorillas, chimpanzees, and humans. Comparison of the repeat sequences detected purifying selection within species and lineage-specific duplications across species. We also found variation in the length of the repeated region within species such as chimpanzee and crab-eating macaque.

CONCLUSIONS

We conclude that the copy number variation in the repeat sequences of FLG between primates may be a consequence of species-specific divergence and expansion.

GWAS of clinically defined gout and subtypes identifies multiple susceptibility loci that include urate transporter genes

Objective

A genome-wide association study (GWAS) of gout and its subtypes was performed to identify novel gout loci, including those that are subtype-specific.

Methods

Putative 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.

Results

In addition to the five loci we reported previously, further susceptibility loci were identified at a genome-wide significance level (p<5.0×10⁻⁸): 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 (p_meta=3.58×10⁻⁸).

Conclusions

Our 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.

Germline Variants of Prostate Cancer in Japanese Families

Prostate cancer (PC) is the second most common cancer in men. Family history is the major risk factor for PC. Only two susceptibility genes were identified in PC, BRCA2 and HOXB13. A comprehensive search of germline variants for patients with PC has not been reported in Japanese families. In this study, we conducted exome sequencing followed by Sanger sequencing to explore responsible germline variants in 140 Japanese patients with PC from 66 families. In addition to known susceptibility genes, BRCA2 and HOXB13, we identified TRRAP variants in a mutually exclusive manner in seven large PC families (three or four patients per family). We also found shared variants of BRCA2, HOXB13, and TRRAP from 59 additional small PC families (two patients per family). We identified two deleterious HOXB13 variants (F127C and G132E). Further exploration of the shared variants in rest of the families revealed deleterious variants of the so-called cancer genes (ATP1A1, BRIP1, FANCA, FGFR3, FLT3, HOXD11, MUTYH, PDGFRA, SMARCA4, and TCF3). The germline variant profile provides a new insight to clarify the genetic etiology and heterogeneity of PC among Japanese men.

A partial nuclear genome of the Jomons who lived 3000 years ago in Fukushima, Japan

The Jomon period of the Japanese Archipelago, characterized by cord-marked ‘jomon' potteries, has yielded abundant human skeletal remains. However, the genetic origins of the Jomon people and their relationships with modern populations have not been clarified. We determined a total of 115 million base pair nuclear genome sequences from two Jomon individuals (male and female each) from the Sanganji Shell Mound (dated 3000 years before present) with the Jomon-characteristic mitochondrial DNA haplogroup N9b, and compared these nuclear genome sequences with those of worldwide populations. We found that the Jomon population lineage is best considered to have diverged before diversification of present-day East Eurasian populations, with no evidence of gene flow events between the Jomon and other continental populations. This suggests that the Sanganji Jomon people descended from an early phase of population dispersals in East Asia. We also estimated that the modern mainland Japanese inherited <20% of Jomon peoples' genomes. Our findings, based on the first analysis of Jomon nuclear genome sequence data, firmly demonstrate that the modern mainland Japanese resulted from genetic admixture of the indigenous Jomon people and later migrants.

Shared Genetic Risk Factors of Intracranial, Abdominal, and Thoracic Aneurysms

BACKGROUND

Intracranial aneurysms (IAs), abdominal aortic aneurysms (AAAs), and thoracic aortic aneurysms (TAAs) all have a familial predisposition. Given that aneurysm types are known to co-occur, we hypothesized that there may be shared genetic risk factors for IAs, AAAs, and TAAs.

METHODS AND RESULTS

We performed a mega-analysis of 1000 Genomes Project-imputed genome-wide association study (GWAS) data of 4 previously published aneurysm cohorts: 2 IA cohorts (in total 1516 cases, 4305 controls), 1 AAA cohort (818 cases, 3004 controls), and 1 TAA cohort (760 cases, 2212 controls), and observed associations of 4 known IA, AAA, and/or TAA risk loci (9p21, 18q11, 15q21, and 2q33) with consistent effect directions in all 4 cohorts. We calculated polygenic scores based on IA-, AAA-, and TAA-associated SNPs and tested these scores for association to case-control status in the other aneurysm cohorts; this revealed no shared polygenic effects. Similarly, linkage disequilibrium-score regression analyses did not show significant correlations between any pair of aneurysm subtypes. Last, we evaluated the evidence for 14 previously published aneurysm risk single-nucleotide polymorphisms through collaboration in extended aneurysm cohorts, with a total of 6548 cases and 16 843 controls (IA) and 4391 cases and 37 904 controls (AAA), and found nominally significant associations for IA risk locus 18q11 near RBBP8 to AAA (odds ratio [OR]=1.11; P=4.1×10(-5)) and for TAA risk locus 15q21 near FBN1 to AAA (OR=1.07; P=1.1×10(-3)).

CONCLUSIONS

Although there was no evidence for polygenic overlap between IAs, AAAs, and TAAs, we found nominally significant effects of two established risk loci for IAs and TAAs in AAAs. These two loci will require further replication.

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.

A large number of variants associated with human complex diseases have been discovered by genome-wide association studies (GWASs). These discoveries have been anticipated to be translated into the definitive understanding of disease pathogeneses; however, functional characterization of the disease-associated SNPs remains a formidable challenge. Here we explored regulatory mechanism of a variant on chromosome 9p21 associated with endometriosis, a common gynecological disorder. By scrutinizing linkage disequilibrium structure and DNase I hypersensitive sites across the risk locus, we prioritized rs17761446 as a candidate causal variant. The results of our “allele-specific” functional genomic approaches sheds light on regulatory mechanisms underlying 9p21 endometriosis risk locus, in which preferential bindings of TCF7L2 and its coactivator EP300 to the protective G allele of rs17761446 lead to stronger chromatin interaction with the promoter of ANRIL, which in turn activate transcription of the non-coding RNA. Motivated by the fact that TCF7L2 was a key transcription factor of Wnt signaling pathway, we postulated that the induction of Wnt signaling activated expression levels of ANRIL and cell cycle inhibitors, CDKN2A/2B. Functional genomics on common disease will unlock functional aspect of genotype-phenotype correlations in the post-GWAS stage.

Comprehensive genetic testing approach for major inherited kidney diseases, using next-generation sequencing with a custom panel

BACKGROUND

Gene identification of hereditary kidney diseases by DNA sequencing is important for precise diagnosis, treatment, and genetic consultations. However, the conventional Sanger sequencing is now practically powerless in the face of ever increasing numbers of reported causative genes of various hereditary diseases. The advent of next-generation sequencing technology has enabled large-scale, genome-wide, simultaneous sequence analyses of multiple candidate genes.

METHODS

We designed and verified a comprehensive diagnosis panel for approximately 100 major inherited kidney diseases, including 127 known genes. The panel was named Simple, sPEedy and Efficient Diagnosis of Inherited KIdney Diseases (SPEEDI-KID). We applied the panel to 73 individuals, clinically diagnosed with an inherited kidney disease, from 56 families.

RESULTS

The panel efficiently covered the candidate genes and allowed a prompt and accurate genetic diagnosis. Moreover, 18 unreported mutations suspected as the disease causes were detected. All these mutations were validated by Sanger sequencing, with 100 % concordance.

CONCLUSION

In conclusion, we developed a powerful diagnostic method, focusing on inherited kidney diseases, using a custom panel, SPEEDI-KID, allowing a fast, easy, and comprehensive diagnosis regardless of the disease type.

Novel kinase fusion transcripts found in endometrial cancer

Recent advances in RNA-sequencing technology have enabled the discovery of gene fusion transcripts in the transcriptome of cancer cells. However, it remains difficult to differentiate the therapeutically targetable fusions from passenger events. We have analyzed RNA-sequencing data and DNA copy number data from 25 endometrial cancer cell lines to identify potential therapeutically targetable fusion transcripts, and have identified 124 high-confidence fusion transcripts, of which 69% are associated with gene amplifications. As targetable fusion candidates, we focused on three in-frame kinase fusion transcripts that retain a kinase domain (CPQ-PRKDC, CAPZA2-MET, and VGLL4-PRKG1). We detected only CPQ-PRKDC fusion transcript in three of 122 primary endometrial cancer tissues. Cell proliferation of the fusion-positive cell line was inhibited by knocking down the expression of wild-type PRKDC but not by blocking the CPQ-PRKDC fusion transcript expression. Quantitative real-time RT-PCR demonstrated that the expression of the CPQ-PRKDC fusion transcript was significantly lower than that of wild-type PRKDC, corresponding to a low transcript allele fraction of this fusion, based on RNA-sequencing read counts. In endometrial cancers, the CPQ-PRKDC fusion transcript may be a passenger aberration related to gene amplification. Our findings suggest that transcript allele fraction is a useful predictor to find bona-fide therapeutic-targetable fusion transcripts.